MATLAB第7章振动解析.ppt

下载文档 降价啦

40
0
约9.27千字
约 72页
2017-01-10 发布于湖北
举报
版权申诉
保障服务

MATLAB第7章振动解析.ppt

1、本文档共72页，可阅读全部内容。
2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。

MATLAB第7章振动解析

1 第七章动力学与振动 7.1 轨迹 7.2 单自由度系统 7.3 多自由度系统由初始条件建立执行文件menu71.m initcond=[2 0 0 1.5;1 0 0 2*pi;2 0 0 4]; tspan=linspace(0,5,1000); options=odeset(RelTol,1e-6,AbsTol,[1e-6 1e-6 1e-6 1e-6]); lintype=[-. -. -.]; for i=1:3 [t,x]=ode45(orbit,tspan,[initcond(i,:)],options); polar(x(:,3),x(:,1),lintype(2*(i-1)+1:2*i)); hold on end text(0.5,-1.2,椭圆轨迹); text(-1.2,1,圆轨迹); text(1.75,2,双曲线轨迹); 运行结果由初始条件（位移和速度均为1时）建立执行文件menu72.m zeta=[0.1 1.0 5.0]; tspan=linspace(0,40,400); lintype=[-b --r ---r]; for i=1:3 [t,x]=ode45(FreeOcillation,tspan,[1 1],[],zeta(i)); subplot(2,1,1); plot(t,x(:,1),lintype(2*(i-1)+1:2*i)); hold on subplot(2,1,2); plot(x(:,1),x(:,2),lintype(2*(i-1)+1:2*i)); hold on end subplot(2,1,1); xlabel(Time( \tau)); ylabel(Displacement x( \tau)); title(Displacement as a function of( \tau)); axis([0 40 -2.0 2.0]); text(2.7,-1.3,阻尼系数=0.1); text(3.6,-0.1,1.0); text(3.6,1.0,5.0); subplot(2,1,2); xlabel(Displacement); ylabel(Velocity); title(Phase portrait); axis([-2.0 2.0 -2.0 2.0]); text(0.7,-1.25,阻尼系数=0.1); text(0.8,-0.65,1.0); text(0.8,0.1,5.0); 运行结果 t = 0:0.001:0.6; x = sin(2*pi*50*t)+sin(2*pi*120*t); y = x + 2*randn(size(t)); plot(y(1:50)) title(Signal Corrupted with Zero-Mean Random Noise) xlabel(time (seconds)) 编制执行文件menu72f.m zeta=0.4;Omega=3.0;x0=50; tspan=linspace(0,30,6000); options=odeset(RelTol,1e-8,AbsTol,1e-8); lintype=[-b]; [t,x]=ode45(ForceOcillation,tspan,[0 0],options,zeta,Omega,x0); subplot(2,1,1); plot(t,x(:,1)); axis([0 30 -8 8]); hold on subplot(2,1,2); yy=x(:,1); N=2048;Nstart=3200;Fs=200; [f,Amplitude]=AmplitudeSpectrum(yy,Fs,Nstart,N); semilogy(f(1:40),2*Amplitude(1:40)); xlabel(Frequency); ylabel(Amplitude); title(Response spectrum of a linear system); hold on subplot(2,1,1); xlabel(Time( \tau)); ylabel(Displacement x( \tau)); title(Response of a linear system); hold on 运行结果频率响应阶跃响应脉冲响应 Bode函数: g = tf([1 0.1 7.5],[1 0.12 9 0 0]); bode(g) 编制m文件. k=[2 -1 0;-1 2 -1;