单自由度的振动响应及相图程序
各位好！我做了一个关于但自由度的振动响应及相图程序。它是在三种阻尼情况下运行的，但只出来阻尼为0.1的，效果也不好。另俩个没有显示出来，请大家帮忙看看原因。公式见附件。function xdot=free1(t,x,dump,zeta,alpha)
xdot=;
zeat=;alpha=;
tspan=linspace(0,40,400);
for i=1:3
=ode45('free1',tspan,',[],zeta(i),alpha(i));
subplot(2,1,1);
plot(t,x(:,1));hold on
subplot(2,1,2);
plot(x(:,1),x(:,2));
hold on
end
subplot(2,1,1)
xlabel('time(\tau)');
ylabel('displacement x(\tau)');
axis();
text(2.7 ,-1.3,'\zeta=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 portant');
axis([-2.0 2.0 -2.0 2.0]);
text();
text(0.8 ,-0.65,'1.0');text(0.8,0.1,'5.0');
页:
[1]