第6章综合实例(可作为第8次上机练习,或者第8次上机时测试).pptVIP

* * 第6章 综合实例 用 MATLAB 求解问题时，一般要经历建模和编程两个过程，只有在建模正确的前提下，方能得出正确的结果。 一、单自由度系统有阻尼自由振动 1.建立计算模型 由动力学可知，单自由度有阻尼自由振动的振动方程为： 无量刚化后有： 其中 , 上述方程的解为： 其中 x0 表示初始位置， v0 表示初始速度。 参数ωn ＝10， x0 ＝1， v0 ＝0，计算的终止时间 t=2。试求ξ 从 0.1 到 1运动方程的解，并画出波形。 2. MATLAB 编程 编写 M 文件 ex1.m %首先清空 MATLAB 的工作空间 clear; %给定初值 wn=10; tf=2; x0=1; v0=0; %计算不同的ξ 值所对应的振型 for j=1:10; eta(j)=0.1*j; wd(j)=wn*sqrt(1-eta(j)^2); %求振幅 A a=sqrt((wn*x0*eta(j)+v0)^2+(x0*wd(j))^2)/wd(j); 第6章 综合实例 %求相位角 phi=atan2(wd(j)*x0,v0+eta(j)*wn*x0); %设定自变量数组 t t=0:tf/1000:tf; %求过渡过程 x(j,:)=a*exp