- 积分
- 0
- 注册时间
- 2009-3-26
- 仿真币
-
- 最后登录
- 1970-1-1
|
我做了一个两自由度振动模型,以其中一个参数作为变量,通过不同的值比较曲线,但是图像不能显示在一张图上,出现以下提示
??? Error using ==> plot
Vectors must be the same lengths.
程序如下:
function dy=matrixode(t,y)
dy=zeros(4,1)
dy(1)=y(3)
dy(2)=y(4)
dy(3)=-0.5*y(1)-0.5*(-1)*y(3)+y(4)-y(2)
dy(4)=1-0.5*y(1)-0.5*(-1)*y(3)
function dz=myfun(t,z)
dz=zeros(4,1)
dz(1)=z(3)
dz(2)=z(4)
dz(3)=-0.5*z(1)-0.5*(-10)*z(3)+z(4)-z(2)
dz(4)=1-0.5*z(1)-0.5*(-10)*z(3)
tspan=[0 1];
x_init=[0;0;-0.12;0];
[t,x]=ode45(@m1_equal,tspan,x_init);
plot(t,x(:,1))
tspan=[0 1];
y_init=[0;0;-0.12;0];
[t,y]=ode45(@m2_equal,tspan,y_init)
plot(t,y(:,1))
plot(t,x(:,1),'g-',t,y(:,1),'r-') |
|