※ 引述《qqazw5789tw (嘟嘟)》之銘言:
: 求解微分方程組
: p0'(t)=-0.01*e^(-0.001t)*p0(t)+0.021(t)
: p1'(t)=0.01*e^(-0.001t)*p0(t)-0.02*p1(t
: 初p0(0)=1,p1(0)=0
: 如何求出p0(t)和p1(t)
: 程式跑出來的結果是
: P[0](t) =
: 2000*exp(10*exp(-(1/1000)*t)-(1/50)*t)*(Int<0~t>((1/100000)*exp(-10*exp(-1/1000)
: t)+(1/50)*t)dt)+1/exp(10)
: P[1](t) =
: -2000*exp(10*exp(-(1/1000)*t)-(1/50)*t)*(Int<0~t>((1/100000)*exp(-10*exp(-1/1000
: t)+(1/50)*t)dt)+1-exp(10*exp(-(1/1000)*t)-(1/50)*t)
: 不知道用甚麼方法可以解這種問題
: 以及如何解這題
Clear[p0, p1]
eqn = NDSolve[{p0'[t] == -(1/100) E^(-(t/1000)) p0[t] + 2/100 p1[t],
p1'[t] == 1/100 E^(-(t/1000)) p0[t] - 2/100 p1[t], p0[0] == 1,
p1[0] == 0}, {p0[t], p1[t]}, t]
P0[t_] = eqn[[1, 1, 2]];
P1[t_] = eqn[[1, 2, 2]];
tMin = -1; tMax = 1;
Plot[{P0[t], P1[t]}, {t, tMin, tMax}]
ParametricPlot[{P0[t], P1[t]}, {t, tMin, tMax}]
--
By Mathematica 9.0
--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 140.109.112.96
※ 文章網址: http://www.ptt.cc/bbs/Math/M.1403706892.A.77A.html