用matlab解含分段函數(shù)的一階微分方程
瀏覽：281715
已知一階常微分方程，
g＋0.047*du/dt＋ u/6.7＝(13-u)/18
其中，當(dāng)sin(10/pi*t)>=0時，g=1.18sin(10/pi*t)
當(dāng)sin(10/pi*t)<0時， g=0
u的初值為u(0)=0，求t>=0時的解。畫出圖像，并求u的極大值。
曾經(jīng)將g寫成，0.59 sin(10/pi*t)+abs(0.59*sin(10/pi*t)),用dsolve求解，可惜出錯了。
請大家?guī)兔Γx謝！（用其他函數(shù)解也可以。）
程序如下
fun=inline(['((13-u)/18-(sin(10*t/pi)>0)*',...
'1.18*sin(10*t/pi)-u/6.7)/0.047'],'t','u');
[t,u]=ode45(fun,[0,10],[0]);
plot(t,u)
說明g這樣表示的:
gt=(sin(10*t/pi)>0)*1.18*sin(10*t/pi);
感謝蘿卜網(wǎng)友
g＋0.047*du/dt＋ u/6.7＝(13-u)/18
其中，當(dāng)sin(10/pi*t)>=0時，g=1.18sin(10/pi*t)
當(dāng)sin(10/pi*t)<0時， g=0
u的初值為u(0)=0，求t>=0時的解。畫出圖像，并求u的極大值。
曾經(jīng)將g寫成，0.59 sin(10/pi*t)+abs(0.59*sin(10/pi*t)),用dsolve求解，可惜出錯了。
請大家?guī)兔Γx謝！（用其他函數(shù)解也可以。）
程序如下
fun=inline(['((13-u)/18-(sin(10*t/pi)>0)*',...
'1.18*sin(10*t/pi)-u/6.7)/0.047'],'t','u');
[t,u]=ode45(fun,[0,10],[0]);
plot(t,u)
說明g這樣表示的:
gt=(sin(10*t/pi)>0)*1.18*sin(10*t/pi);
感謝蘿卜網(wǎng)友
技術(shù)鄰APP
工程師必備
工程師必備
- 項目客服
- 培訓(xùn)客服
- 平臺客服
TOP
