在微分方程表示式輸入中,以大寫字母D來表示微分
dsolve呼叫格式:dsolve("equ1","equ2",.........."equN")
另外要注意:在微分方程表示式輸入中,以大寫字母D來表示微分,D,D2,.......Dn分別表示一階,二階和n階
2dx/dt+dy/dt-y=exp(-t)
dx/dt+x+y=0 其中初始條件:x(0)=1.5,y(0)=0
首先求解微分方程的通解:
s=dsolve("2*Dx+Dy-y=exp(-t)","Dx+x+y=0");%求解微分方程組的通解
>> s.x
%微分方程組變數x的通解
ans =
-C1*exp((1+2^(1/2))*t)-C2*exp(-(2^(1/2)-1)*t)+1/2*C1*exp((1+2^(1/2))*t)*2^(1/2)-1/2*C2*exp(-(2^(1/2)-1)*t)*2^(1/2)-1/2*exp(-t)
>> s.y
%微分方程組變數y的通解
C1*exp((1+2^(1/2))*t)+C2*exp(-(2^(1/2)-1)*t)
然後根據初始條件,求解微分方程組的特解:
>> s=dsolve("2*Dx+Dy-y=exp(-t)","dx+x+y=0","x(0)=1.5","y(0)=0");%微分方程組在給定初始條件下的特解
在微分方程表示式輸入中,以大寫字母D來表示微分
dsolve呼叫格式:dsolve("equ1","equ2",.........."equN")
另外要注意:在微分方程表示式輸入中,以大寫字母D來表示微分,D,D2,.......Dn分別表示一階,二階和n階
2dx/dt+dy/dt-y=exp(-t)
dx/dt+x+y=0 其中初始條件:x(0)=1.5,y(0)=0
首先求解微分方程的通解:
s=dsolve("2*Dx+Dy-y=exp(-t)","Dx+x+y=0");%求解微分方程組的通解
>> s.x
%微分方程組變數x的通解
ans =
-C1*exp((1+2^(1/2))*t)-C2*exp(-(2^(1/2)-1)*t)+1/2*C1*exp((1+2^(1/2))*t)*2^(1/2)-1/2*C2*exp(-(2^(1/2)-1)*t)*2^(1/2)-1/2*exp(-t)
>> s.y
%微分方程組變數y的通解
ans =
C1*exp((1+2^(1/2))*t)+C2*exp(-(2^(1/2)-1)*t)
然後根據初始條件,求解微分方程組的特解:
>> s=dsolve("2*Dx+Dy-y=exp(-t)","dx+x+y=0","x(0)=1.5","y(0)=0");%微分方程組在給定初始條件下的特解