jiangchunlei 发表于 2008-11-3 10:34

谁能帮我求解这个微分方程组不?

syms a r w j real
=dsolve('Dy1=r*(y5+y6)-a*(y2-y3)*i,Dy2=-r*y2-(a*y1-a*y4-2*w*y2)*i,Dy3=-r*y3-(a*y4-a*y1+2*w*y3)*i,Dy4=-2*r*y4-a*(y3-y2)*i,Dy5=r*(y4-y5)-j*(y7-y8)*i,Dy6=r*(y4-y6)+j*(y7-y8)*i,Dy7=r*y7-j*(y5-y6)*i,Dy8=r*y8+j*(y5-y6)*i','y1(0)=0,y2(0)=0,y3(0)=0,y4(0)=1,y5(0)=0,y6(0)=0,y7(0)=0,y8(0)=0')

sogooda 发表于 2008-11-3 16:29

回复 楼主 jiangchunlei 的帖子

看的好晕,贴个方程的图吧。
另外,你的方程如果很复杂的话建议数值求解。

jiangchunlei 发表于 2008-11-4 22:38

怎么做数值求解啊,你能帮我看一下吗
我的有八个要求的量就是y1,y2,y3,y4,y5,y6,y7,y8
r,a,w,j都是参数
我的qq是419271324不胜感谢哦

yufeng 发表于 2008-11-5 09:07

里面的系数a ij 是不是不知道?

jiangchunlei 发表于 2008-11-5 11:15

我的都是不大于1 的矩阵元i 是虚数单位a,w,r,j都是参数,我运行的时候加了初试条件就有带c1,c2。。。c8的结果,加了初试条件就提示??? Error using ==> dsolve
Error, (in dsolve/IC) The 'implicit' option is not available when giving Initial Conditions.
Error in ==> ded at 2
=Dsolve('Dy1=r*(y5+y6)-a*(y2-y3)*i,Dy2=-r*y2-(a*y1-a*y4-2*w*y2)*i,Dy3=-r*y3-(a*y4-a*y1+2*w*y3)*i,Dy4=-2*r*y4-a*(y3-y2)*i,Dy5=r*(y4-y5)-j*(y7-y8)*i,Dy6=r*(y4-y6)+j*(y7-y8)*i,D
>>
页: [1]
查看完整版本: 谁能帮我求解这个微分方程组不?