当变量乘以常量时，dsolve出错（R2011a）

Question

我正在尝试使用以下代码：

ode1='D2y1=-1256.4*y1-5*Dy1+255.1*y2+182.781';   
ode2='D2y2=-151.5*y2-5*Dy2+255.1*y1-14.0459';   
CI='y1(0)=2,y2(0)=-2,Dy1(0)=0,Dy2(0)=0';

sol=dsolve(ode1,ode2,CI,'t');
sol.y1
sol.y2

和matlab返回错误：

  

???使用==＆gt;时出错mupadengine.mupadengine> mupadengine.feval at 144   MuPAD错误：错误：除以零[_invert];

     

在评估'stdlib :: normalNoExpand'

期间      

==＆gt;中的错误dsolve＆gt; mupadDsolve在215 T =   feval（symengine， 'symobj :: dsolve'，SYS中，x，ignoreConstraints）;

     

==＆gt;中的错误溶于96溶胶=   mupadDsolve（ignoreConstraints，varargin {1：narg}）;

     

==＆gt;中的错误maglevsol at 7 sol = dsolve（ode1，ode2，CI，'t'）;

但是，如果我运行这个

ode1='D2y1=-y1-5*Dy1+255.1*y2+182.781';   
ode2='D2y2=-y2-5*Dy2+255.1*y1-14.0459';   
CI='y1(0)=2,y2(0)=-2,Dy1(0)=0,Dy2(0)=0';

sol=dsolve(ode1,ode2,CI,'t');
sol.y1
sol.y2

我完全没有错误。 怎么了？我的版本是r2011a

Answer 1

由于IC的问题，似乎Matlab无法解决这个问题。 （在尝试数值求解器时给出的t = 0时导数的非数值）。我无法找到解决方法（Matlab 2015a实际上是悬而未决的，用于分析解决方案）。

但Maple能够解决它（Maple在ODE中非常好）。这是解决方案和图（对于所有那些硬编码的数值使用符号，然后在末尾使用subs也是一个好主意。不需要在整个地方都有魔术数字）。您可能希望将此报告给Mathworks

restart;
ode1 := diff(y1(t),t$2) = a*y1(t) + b*diff(y1(t),t)+ c*y2(t) + d;
ode2 := diff(y2(t),t$2)= e*y2(t) + b*diff(y2(t),t)+ c*y1(t) + f;
ic := y1(0) = 2, y2(0) = -2, D(y1)(0) = 0, D(y2)(0) = 0;
sol:=dsolve({ode1, ode2, ic}, {y1(t), y2(t)}):
sol:=subs( {a = -1256.4, b = -5, c= 182.781, e = -151.5, f = -14.0459,d=182.781},sol);


sol := {y1(t) = (-.1922320350+0.5088536230e-1*I)*
exp((-2.500000000+9.444367980*I)*t)+(1.096005949-0.7581365270e-1*I)*
exp((-2.500000000+36.14144316*I)*t)+(-.1922320331-0.5088536185e-1*I)*
exp((-2.500000000-9.444367980*I)*t)+(1.096005952+0.7581365750e-1*I)*
exp((-2.500000000-36.14144316*I)*t)+.1924521734, 

y2(t) = (-.8748433300+.2315780499*I)*exp((-2.500000000+9.444367980*I)*t)+
(-.2408287839+0.1665876887e-1*I)*exp((-2.500000000+36.14144316*I)*t)+
(-.8748433212-.2315780478*I)*exp((-2.500000000-9.444367980*I)*t)+
(-.2408287844-0.1665876992e1*I)*
exp((-2.500000000-36.14144316*I)*t)+.2313442207}

这是y1（t）和y2（t）的解决方案图2秒