hack3702
2014-09-10, 19:57
各位高手好,我需要用符号运算工具箱,求解一个四元的一阶线性非齐次微分方程组的稳态解析解,但是在运行之后一直显示busy,得不到有效的结果,请问是什么原因,该如何修改,谢谢
代码如下:
clear all;
syms R L Rs Ls Rr Lr Lm w1 w U0 t;
S=dsolve('(R+Rs)*ids+Ls*Dids+Lm*Didr-w*Ls*iqs-w*Lm*iqr=0', ...
'(R+Rs)*iqs+Ls*Diqs+Lm*Diqr+w*Ls*ids+w*Lm*idr=0', ...
'Rr*idr+Lr*Didr+Lm*Dids=2/3*U0*cos(w1*t)', ...
'Rr*iqr+Lr*Diqr+Lm*Diqs=0', ...
'ids(0)=0','iqs(0)=0','idr(0)=0','iqr(0)=0', ...
'Dids(0)=0','Diqs(0)=0','Didr(0)=0','Diqr(0)=0','t');
ids=S.ids,
iqs=S.iqs,
idr=S.idr,
iqr=S.iqr,
代码如下:
clear all;
syms R L Rs Ls Rr Lr Lm w1 w U0 t;
S=dsolve('(R+Rs)*ids+Ls*Dids+Lm*Didr-w*Ls*iqs-w*Lm*iqr=0', ...
'(R+Rs)*iqs+Ls*Diqs+Lm*Diqr+w*Ls*ids+w*Lm*idr=0', ...
'Rr*idr+Lr*Didr+Lm*Dids=2/3*U0*cos(w1*t)', ...
'Rr*iqr+Lr*Diqr+Lm*Diqs=0', ...
'ids(0)=0','iqs(0)=0','idr(0)=0','iqr(0)=0', ...
'Dids(0)=0','Diqs(0)=0','Didr(0)=0','Diqr(0)=0','t');
ids=S.ids,
iqs=S.iqs,
idr=S.idr,
iqr=S.iqr,