哈工大 《机械振动基础》大作业.doc
文本预览下载声明
《机械振动基础》大作业
（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,
显示全部