如何求解这样的微分方程组

鸾凤阁

dE/dt=-k*E+g*E*D+x
  dD/dt=-[(1+Ip+E*E)*D-Ip+1]/y
  k=k0[1+sin(wt])

ko=3.3*10^7; g=6.0*10^7; y=10*10^-3;
x=1.0*10^-4; w=3.5*10^5; Ip=2*10^-3

E,D的初值‘你做主’（可任意设）
谢谢您的帮助！！！

[ 本帖最后由 ChaChing 于 2009-5-10 23:52 编辑 ]

liushuiwuxin

这样可以解除结果，但是不知道结果是不是真确的，是不是你要滴
function weiwen
M0=[1;2];
[t,M]=ode15s(@wei,[0:0.2:1],M0)
function F=wei(t,M)
E=M(1); D=M(2);
ko=3.3*10^7; g=6.0*10^7; y=10*10^-3;
x=1.0*10^-4; w=3.5*10^5; Ip=2*10^-3
k=k0*(1+sin(w*t));
DE=-k*E+g*E*D+x;
DD=-[(1+Ip+E*E)*D-Ip+1]/y;
F=[DE;DD]

[ 本帖最后由 ChaChing 于 2009-5-10 23:55 编辑 ]

鸾凤阁

谢谢你的帮助！！！
我用的是（引用你的函数名）
f=zeros(2,1)
E=M(1);
D=M(2);
f(1)=-k*E+g*E*D+x;
f(2)=-[(1+Ip+E*E)*D-Ip+1]/y;
这样就不需要f=[DE;DD]这样一步了！！！
非常给你写你的回复！！！
以后请多多指教！我是刚学不到五天！
你能告诉我一下spectrum一些基本语法吗？
不胜感激！！！