[新手求助]SOS，大家帮帮忙阿，谢谢啦

woai534 · 发表于 2006-5-12 15:36

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

function sys tstop=[0 8]; y0=[0;0;2]; ode45(@f,tstop,yo) plot(t,y) function dydt=f(t,y) a=[2*cos(t-y1); 3*sin(t-y1); 0]; b=[2 9 y2; 8 2 4;3 7 5]; dydt=a-b*y; y1是微分方程组中的一个量，这个y1直接这么放在程序中是会显示y1未定义的， y2也是同样的问题 我该怎么处理呢，我不想将这个矩阵乘出来？

[此贴子已经被作者于2006-5-12 18:31:03编辑过]

woai534 · 发表于 2006-5-12 18:18

难道我一定要把矩阵乘出来，写出各个方程的表达式不成？

woai534 · 发表于 2006-5-12 22:33

没人指导我吗？

[此贴子已经被作者于2006-5-12 23:09:22编辑过]

feichang61 · 发表于 2006-5-13 09:17

ode45(@f,tstop,yo)中yo应该是y0,y1应该是y(1),同理y2 不知你是不是这个意思

woai534 · 发表于 2006-5-13 13:20

是的，就是这个意思，我该怎么处理啊？求高人指点阿

happy · 发表于 2006-5-13 15:49

function sys tstop=[0 8]; y0=[0;0;2]; ode45(@f,tstop,y0) plot(t,y) function dydt=f(t,y) a=[2*cos(t-y(1)); 3*sin(t-y(1)); 0]; b=[2 9 y(2); 8 2 4;3 7 5]; dydt=a-b*y;

woai534 · 发表于 2006-5-13 16:02

难道我又敲错程序了？看来就是这样子的，那我那题其实我昨天就已经做出来了。谢谢HAPPY教授

happy · 发表于 2006-5-13 16:05

<DIV class=quote>以下是引用woai534在2006-5-13 16:02:49的发言： 
难道我又敲错程序了？看来就是这样子的，那我那题其实我昨天就已经做出来了。谢谢HAPPY教授</DIV>
 写程序的是自己要注意检查 有些错误只要小心点就能检查出来 问之间自己先检查清楚

woai534 · 发表于 2006-5-13 16:20

谢谢热心帮助我的人，谢谢

账号		自动登录	找回密码
密码			我要加入

[新手求助]SOS，大家帮帮忙阿，谢谢啦

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

回复：（woai534）[新手求助]SOS，大家帮帮忙阿，谢谢...

回复：（woai534）难道我又敲错程序了？看来就是这样子...