请教一个很简单的问题
对于微分方程x''-x=cos(0.5*t)的解为x=-0.8*cos(0.5*t)。如果用龙哥库塔法解，却得到发散的结果，请问大家是怎么回事，我的程序有错么？function dy=wu(t,y)
dy=zeros(2,1);
dy(1)=y(2);
dy(2)=y(1)+cos(0.5*t)
=ode45('wu',,);
plot(t,x(:,1)); x=-0.8*cos(0.5*t)
只是它的一个特解吧？ x''-x=f(t):本身就是不稳定的，那种方法都一样。很有可能方程是x''++++++++++++++++x=cos(0.5*t)，符号被搞错了 本方程的通解为x=c1*exp(t)+c2*exp(-t) -0.8cos(0.5*t), c1, c2是任意常数
楼主的初始条件为x(0)=0.1, x'(0)=0，可得c1=c2=0.9
这样t=20并不会发散，另有原因
页:
[1]