求助：ode求解微分方程组
用于求解方程如下:function dx=ehpcs(t,x)
w=15*pi;
k=300;
J1=0.191;
J2=0.097;
dx=zeros(4,1);
T0=40;
dx(1)=x(2);
dx(2)=T0*sin(4*w*t)/J1-(k/J1)*(x(1)-x(3));
dx(3)=x(4);
dx(4)=-(k/J2)*(x(3)-x(1));
%ode求解程序如下：
tspan=;
x0=;
=ode45('ehpcs',tspan,x0);
plot(t,x(:,1),'r',t,x(:,3),'b-.')
grid
微分方程就是最简单的线性微分方程组，是实际应用的减振器的简化模型，但是求解出来的结果不是围绕x轴的波动数据，本人刚接触ode求解，自己没有检查出程序错误，希望有会的大神帮我看看是哪里出了错，不甚感激！
页:
[1]