《机械振动基础》大作业（2016 年春季学期）题目多自由度振动系统的固有频率和固有阵型姓名学号班级专业报告提交日期哈尔滨工业大学《机械振动基础》大作业1 报告要求1. 请根据课堂布置的2 道大作业题，任选其一，拒绝雷同和抄袭；2. 报告最好包含自己的心得、体会或意见、建议等；3. 报告统一用该模板撰写，字数不少于3000 字，上限不限；4. 正文格式： 小四号 字体，行距为 1.25 倍行距；5. 用 A4 纸单面打印 ；左侧装订， 1 枚钉；6. 课程报告需同时提交打印稿和电子文档予以存档，电子文档由班长收齐，统一发送至： liuyingxiang868@hit.edu.cn。7. 此页不得删除。评语：成绩（15 分）：教师签名：年月日《机械振动基础》大作业2 基于 MATLAB 的对多自由度振动系统的固有频率和固有阵型的分析一、利用 MATLAB 编程求七自由度系统的固有频率和固有振型模型如下图所示系统中各质量为（kg ）： m1=2;m2=4;m3=5;m4=6;m5=5;m6=4;m7=2; 各处弹性系数为（N/m ）k1=5;k2=5;k3=5;k4=5;k5=5;k6=5;k7=5; 二、实验程序如下：clear close m1=2; m2=4; m3=5; m4=6; m5=5; m6=4; m7=2; k1=5; k2=5; k3=5; k4=5; k5=5; k6=5; 《机械振动基础》大作业3 k7=5; m=[m1,0,0,0,0,0,0; 0,m2,0,0,0,0,0; 0,0,m3,0,0,0,0; 0,0,0,m4,0,0,0; 0,0,0,0,m5,0,0; 0,0,0,0,0,m6,0; 0,0,0,0,0,0,m7]; k=[k1+k2,-k2,0,0,0,0,0; -k2,k2+k3,-k3,0,0,0,0; 0,-k3,k3+k4,-k4,0,0,0; 0,0,-k4,k4+k5,-k5,0,0; 0,0,0,-k5,k5+k6,-k6,0; 0,0,0,0,-k6,k6+k7,-k7; 0,0,0,0,0,-k7,k7]; [V,D]=eig(k,m) for j=1:1:7 w(j)=sqrt(D(j,j)); for i=1:1:7 absV(i,j)=abs(V(i,j)); end end max=(absV); for j=1:1:7 for i=1:1:7 V(i,j)=V(i,j)/max(j); end 《机械振动基础》大作业4 end w figure x=1:7; for a=1:7 subplot(3,4,a),plot(x,V(x,a)); hold on; grid on; title(' 振型图 '); end 三、将程序在 MATLAB 上运行，得到结果如下图（ 1）多自由度系统的固有频率得到振型图如下《机械振动基础》大作业5 图（ 2） 多自由度系统的固有阵型图（ 2）中分别为本实验条件下，模拟出来小车的振动的阵型图。四、心得、体会1、对多自由度系统的认识和学习多自由度系统，是指必须通过两个以上的独立广义坐标才能够描述系统运动特性的系统， 或者说是自由度数目多于一个，但又不属于连续弹性体 （...
