机械振动大作业可研报告.doc

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