使用matlab绘制非线性受迫振动？

在百edit编辑器中输入度并保存:

clcclear

T=0:.1:100w=input('输入w:')if ~(0.8<=w&&w<=1.5)    error('不在范围内')endf=@回(t,x)[x(2)-(0.2*x(2)-x(1)+x(1)^3)-0.3*cos(w*t)][t,y]=ode45(f,T,[00])plot(t,y(:,1))text(2,max(y(:,1)/1.1),strcat('\omega=',(num2str(w))))

运行后在主窗口中输入 0.9

图像答:

你看看满意吗？

clearclc

x=0:0.01:5

w=10

y1=(1-exp(-x)).*(sin(w*x)/2+.5)

y2=1-exp(-x)

plot(x,y1)

hold on

plot(x,y2)

首先要告诉你的是自由振动能反应的是结构的固有特性，而动力特性这里应该是受迫振动体现的性质，所以根据振动理论，任意振型都可以用自由振动特征向量组表示出来，因此只要知道动力参数，就能求出结构的动力响应。

---