哈工大《机械振动基础》大作业教程.docVIP

《机械振动基础》大作业 （2015年春季学期） 题目 基于MATLAB求系统特性 姓 名 学 号 班 级 专 业 机械设计制造及其自动化 报告提交日期 哈尔滨工业大学 报告要求 请根据课堂布置的2道大作业题，任选其一，拒绝雷同和抄袭； 报告最好包含自己的心得、体会或意见、建议等； 报告统一用该模板撰写，字数不少于3000字，上限不限； 正文格式：小四号字体，行距为1.25倍行距； 用A4纸单面打印；左侧装订，1枚钉； 课程报告需同时提交打印稿和电子文档予以存档，电子文档由班长收齐，统一发送至：shanxiaobiao@hit.edu.cn。 此页不得删除。 评语： 成绩（15分）： 教师签名： 年 月 日 基于MATLAB求系统特性 题目简介： 已知9自由度无阻尼振动系统: 利用MATLAB编程（所有参数自己定，若有雷同，均为零分）求9系统的固有频率和固有振型,（报告中需给出MATLAB程序）； 针对本题的练习过程，写出自己对9自由度系统以及利用MATLAB求解多自由度矩阵的认识体会。 MATLAB程序图： m=[]; k1=[]; k=[]; c=[]; c1=[]; for i=1:9 a=input(输入质量矩阵m：); m(i,i)=a; end ; for j=1:9 b=input(输入刚度系数k：); k1(1,j)=b; end for l=1:8 k(l,l)=k1(l)+k1(l+1); k(9,9)=k1(9); k(l+1,l)=-k1(l+1); k(l,l+1)=-k1(l+1); k(9,8)=-k1(9); k(8,9)=-k1(9); end ; syms w; B=k-w^2*m %系统的特征矩阵B Y=det(B); %展开行列式 W=solve(Y); %求解wh lW=length(W); [V,D]=eig(k,m); for I=1:9 for J=1:9 V(J,I)=V(J,I)/V(5,I); end end V W 三 MATLAB结果输入输出： 3.1 程序输入内容： 输入质量矩阵m：1 输入质量矩阵m：2 输入质量矩阵m：3 输入质量矩阵m：4 输入质量矩阵m：5 输入质量矩阵m：6 输入质量矩阵m：7 输入质量矩阵m：8 输入质量矩阵m：9 输入刚度系数k：10 输入刚度系数k：11 输入刚度系数k：12 输入刚度系数k：13 输入刚度系数k：14 输入刚度系数k：15 输入刚度系数k：16 输入刚度系数k：17 输入刚度系数k：18 3.2 Matlab 输出界面截图： 3.3 输出结果： B = [ 21-w^2, -11, 0, 0, 0, 0, 0, 0, 0] [ -11, 23-2*w^2, -12, 0, 0, 0, 0, 0, 0] [ 0, -12, 25-3*w^2, -13, 0, 0, 0, 0, 0] [ 0, 0, -13, 27-4*w^2, -14, 0, 0, 0, 0] [ 0, 0, 0, -14, 29-5*w^2, -15, 0, 0, 0] [ 0, 0,