查看: 70|回复: 1|关注: 0

[未答复] ode45微分方程求解出错,请大佬指教

[复制链接]

新手

16 麦片

财富积分


050


4

主题

12

帖子

0

最佳答案
发表于 2019-12-2 16:15:11 | 显示全部楼层 |阅读模式
function o=rigid(t,x)
o=zeros(5,1);
o(1)=7.5+(10.2/104.6)*x(5)-((1.7+16)/15.2)*x(1);
o(2)=0+(16/15.2)*x(1)-(16/12.6)*x(2);
o(3)=2.7+(1.6/12.6)*x(2)-((5.9+9.2+3.6)/3.2)*x(3);
o(4)=0+(9.2/3.2)*x(3)-((2.6+6.6)/1194.0)*x(4);
o(5)=0+(3.6/6.2)*x(3)+(6.6/1194.0)*x(4)-(10.2/104.6)*x(5);

>> rigid(t, x)

ans =

-120.9368
  -24.7364
-671.2806
  338.1148
   57.8148

>>
[t,x]=ode45(@rigid,[0 10],[120 180 100 150 156]);
plot(t,x(:,1),'*',t,x(:,2),'.',t,x(:,3),'○',t,x(:,4),'△',t,x(:,5),'◇')
错误使用 plot
color/linetype 参数出错。

新手

16 麦片

财富积分


050


4

主题

12

帖子

0

最佳答案
 楼主| 发表于 2019-12-2 16:23:56 | 显示全部楼层
plot (t,x(:,1),'*',t,x(:,2),'.',t,x(:,3),'+',t,x(:,4),'-',t,x(:,5),'-.')

已经解决,谢谢
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

站长推荐上一条 /3 下一条

快速回复 返回顶部 返回列表