[問題] ODE多變數求解已回收
最近用matlab來解Laser Gain的微分方程
式子如下
dN1/dt=-W(N1-N2)+N4/tau1
dN2/dt=W(N1-N2)-N2/tau2
dN3/dt=-S(N3-N4)-N3/tauG+N2/tau2
dN4/dt=S(N3-N4)+N3/tauG-N4/tau1
N1+N2+N3+N4=Ng
其中Ng W tau1 tau2 tauG是常數
N和S是變數
這樣兩個變數的情況下我ODE的function怎麼寫比較好
如果我的S寫成N(5)會出問題嗎
如果不行的話 應該要怎麼寫才正確
以下附上程式碼(常數宣告的地方全部省略)
function [dN]=laser_gain(t,N) or ( [dN S]=laser_gain(t,N,S) 有這種寫法嗎@@)
dN=zeros(4,1);
dN(1)=-W.*(N(1)-N(2))+N(4)./tau41;
dN(2)=W.*(N(1)-N(2))-N(2)./tau23;
dN(3)=-S.*(N(3)-N(4))-N(3)./tauG+N(2)./tau23;
dN(4)=S.*(N(3)-N(4))+N(3)./tauG-N(4)./tau41;
Ng=N(1)+N(2)+N(3)+N(4);
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.114.19.129
→
09/29 11:47, , 1F
09/29 11:47, 1F
→
09/29 12:29, , 2F
09/29 12:29, 2F