鸾凤阁 发表于 2009-5-8 12:56

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

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

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 发表于 2009-5-9 10:37

这样可以解除结果，但是不知道结果是不是真确的，是不是你要滴
function weiwen
M0=;
=ode15s(@wei,,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=

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

鸾凤阁 发表于 2009-5-10 11:52

回复 沙发 liushuiwuxin 的帖子

谢谢你的帮助！！！
我用的是（引用你的函数名）
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=这样一步了！！！
非常给你写你的回复！！！
以后请多多指教！我是刚学不到五天！
你能告诉我一下spectrum一些基本语法吗？
不胜感激！！！
页: [1]
查看完整版本: 如何求解这样的微分方程组