matlab解微分方程组，图像特别密集，是怎么回事？？

jiaoxiaolei2010

M文件如下：
%%CIM
function dx = funt_CIM(t,x)
dx = [x(1);x(2);x(3);x(4)]
w = 100
dx(1) = x(2)
dx(3) = x(4)
dx(2) = (-.3077e9+w^2)*x(1)+3.299*w*x(4)-.1477e10*x(3)^2
dx(4) = (-.3265e5-.1044e10*x(1)-.3580e10*x(3)^2+w^2)*x(3)-1.167*w*x(2)-.3580e10*x(3)^3
在命令窗口输入下面命令：
[t,x] = ode45('funt_CIM',[0,0.2],[0.001;0;0.001;0])
figure(1)
plot(t,x(:,1));grid on;
xlabel('time s');
ylabel('qf1 m');
figure(2)
plot(t,x(:,2));grid on;
xlabel('time s');
ylabel('qf1 dt m/s');
figure(3)
plot(t,x(:,3));grid on;
xlabel('time s');
ylabel('qf2 m');
figure(4)
plot(t,x(:,4));grid on;
xlabel('time s');
ylabel('qf2 dt m/s');
图像如下：

请问这是怎么回事啊，是刚性方程的问题吗？？？

blackbutter

应该是X 点太多，所以密集

jiaoxiaolei2010

点如果少的话，就不能反映真实变化趋势呀

lengyunfeng

jiaoxiaolei2010 发表于 2013-2-22 12:59
点如果少的话，就不能反映真实变化趋势呀

看你的图片，可能是函数对应的波状变幅的频率太大了，可以试试y－log(t)绘图形式来表达，或者其它方式增大绘图x坐标值之间的间距。

jiaoxiaolei2010

lengyunfeng 发表于 2013-2-22 22:47
看你的图片，可能是函数对应的波状变幅的频率太大了，可以试试y－log(t)绘图形式来表达，或者其它方式增 ...

嗯嗯，谢谢啦