 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 3102 | 回復(fù): 21 | |||
ybnav銀蟲 (小有名氣)
|
[交流]
【求助】向各位大俠求助matlab求解微分方程組遇到的一個(gè)問題 已有7人參與
|
|
向各位大俠求助matlab求解微分方程組遇到的一個(gè)問題： 假設(shè)有以下的微分方程組 u'(t)=-a*f1(t)*y(t); v'(t)=-b*v(t)+a*f1(t)*y(t); w'(t)=b*v(t); x'(t)=-a*f2(t)*v(t)+a*f1(t)*y(t)+a*f3(t)*y(t)+c*y(t); y'(t)= a*f2(t)*v(t) - a*f1(t)*y(t) - a*f3(t)*y(t) - c*y(t); 其中 f1(t)=0 if u(t)=0 =1 if u(t)>0; f2(t)=0 if x(t)=0 =1 if x(t)>0; f3(t)=0 if v(t)=0 =1 if v(t)>0; a=1;b=0.1;c=0.5; 且u(0)=1000;v(0)=1;w(0)=0;x(0)=200;y(0)=0.現(xiàn)在要得到各函數(shù)隨t的變化曲線，如果微分方程組中沒有f1(t)，f2(t)，f2(t)這三個(gè)方程，那么在matlab中可以直接用ode來得到，但是現(xiàn)在有這三個(gè)方程，我就沒辦法了，試了好幾種方法都沒用,想請(qǐng)教各位高手，在matlab下怎么才能求解這個(gè)方程？ [ Last edited by ybnav on 2011-3-11 at 22:27 ] |
回復(fù)此樓» 猜你喜歡
311求調(diào)劑
已經(jīng)有9人回復(fù)
-
302求調(diào)劑
已經(jīng)有4人回復(fù)
-
被我言中：新模板不強(qiáng)調(diào)格式了，假專家開始管格式了
已經(jīng)有4人回復(fù)
-
303求調(diào)劑
已經(jīng)有5人回復(fù)
-
考研求調(diào)劑
已經(jīng)有4人回復(fù)
-
085601專碩，總分342求調(diào)劑，地區(qū)不限
已經(jīng)有4人回復(fù)
-
085601材料工程專碩求調(diào)劑
已經(jīng)有5人回復(fù)
-
一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑
已經(jīng)有6人回復(fù)
-
能源材料化學(xué)課題組招收碩士研究生8-10名
已經(jīng)有6人回復(fù)
-
296求調(diào)劑
已經(jīng)有5人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- Matlab求解偏微分方程組
已經(jīng)有22人回復(fù)
- matlab數(shù)值求解邊界條件微分方程組
已經(jīng)有7人回復(fù)
- 高金求助matlab解微分方程組
已經(jīng)有12人回復(fù)
- 求助各位大俠，關(guān)于一個(gè)四氫呋喃的問題。謝謝了
已經(jīng)有17人回復(fù)
- 求matlab求解一個(gè)方程組的問題
已經(jīng)有6人回復(fù)
- 求助各位大俠一個(gè)去羥基化的問題，謝謝！
已經(jīng)有8人回復(fù)
- 【求助】matlab 求解微分方程中的未知參數(shù)
已經(jīng)有20人回復(fù)
- 【求助】matlab求解非線性方程組，并畫圖處理。要求y,z是實(shí)數(shù)解！
已經(jīng)有18人回復(fù)
- 【求助】MATLAB 有限差分法（FDM）求解偏微分方程
已經(jīng)有22人回復(fù)
- 【求助】積分微分方程matlab求解
已經(jīng)有6人回復(fù)
- 【求助】matlab怎么求解偏微分方程組啊，先謝謝了
已經(jīng)有13人回復(fù)
★
小木蟲(金幣+0.5):給個(gè)紅包，謝謝回帖交流
小木蟲(金幣+0.5):給個(gè)紅包，謝謝回帖交流
|
把你的函數(shù)貼出來，我稍改改就可以求解了。 或者參考：http://www.gaoyang168.com/bbs/viewthread.php?tid=2764636### 把微分方程寫成分段函數(shù)的形式即可。 |
2樓2011-03-11 15:42:54
ybnav
銀蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 681.3
- 帖子: 65
- 在線: 47.6小時(shí)
- 蟲號(hào): 508614
- 注冊(cè): 2008-02-21
- 性別: GG
- 專業(yè): 信息安全
|
就是這個(gè)微分方程組 u'(t)=-*f1(t)*y(t); v'(t)=-0.1*v(t)+*f1(t)*y(t); w'(t)=0.1*v(t); x'(t)=-f2(t)*v(t)+f1(t)*y(t)+f3(t)*y(t)+0.5*y(t); y'(t)= f2(t)*v(t) -f1(t)*y(t) -f3(t)*y(t) -0.5*y(t); 其中 f1(t)=0 if u(t)=0 =1 if u(t)>0; f2(t)=0 if x(t)=0 =1 if x(t)>0; f3(t)=0 if v(t)=0 =1 if v(t)>0; 且u(0)=1000;v(0)=1;w(0)=0;x(0)=200;y(0)=0. 只要求在matlab中畫出u(t)、v(t)、w(t)、x(t)、y(t)在t的某一段區(qū)間內(nèi)，比如[0,60]的曲線圖，還是沒想出來改怎么辦，同學(xué)說可能需要迭代，我也不知道怎么辦，萬分感謝版主同志！ [ Last edited by ybnav on 2011-3-11 at 22:28 ] |
3樓2011-03-11 17:31:40
ybnav
銀蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 681.3
- 帖子: 65
- 在線: 47.6小時(shí)
- 蟲號(hào): 508614
- 注冊(cè): 2008-02-21
- 性別: GG
- 專業(yè): 信息安全
|
剛方程寫錯(cuò)了，應(yīng)該是 就是這個(gè)微分方程組 u'(t)=-f1(t)*y(t); v'(t)=-0.1*v(t)+f1(t)*y(t); w'(t)=0.1*v(t); x'(t)=-f2(t)*v(t)+f1(t)*y(t)+f3(t)*y(t)+0.5*y(t); y'(t)= f2(t)*v(t) -f1(t)*y(t) -f3(t)*y(t) -0.5*y(t); 其中 f1(t)=0 if u(t)=0 =1 if u(t)>0; f2(t)=0 if x(t)=0 =1 if x(t)>0; f3(t)=0 if v(t)=0 =1 if v(t)>0; 且u(0)=1000;v(0)=1;w(0)=0;x(0)=200;y(0)=0. 只要求在matlab中畫出u(t)、v(t)、w(t)、x(t)、y(t)在t的某一段區(qū)間內(nèi)，比如[0,60]的曲線圖，還是沒想出來改怎么辦，同學(xué)說可能需要迭代，我也不知道怎么辦，萬分感謝版主同志！ [ Last edited by ybnav on 2011-3-11 at 22:31 ] |
4樓2011-03-11 17:33:08
5樓2011-03-11 18:40:53
ybnav
銀蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 681.3
- 帖子: 65
- 在線: 47.6小時(shí)
- 蟲號(hào): 508614
- 注冊(cè): 2008-02-21
- 性別: GG
- 專業(yè): 信息安全
6樓2011-03-11 18:58:21
7樓2011-03-11 19:02:59
ybnav
銀蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 681.3
- 帖子: 65
- 在線: 47.6小時(shí)
- 蟲號(hào): 508614
- 注冊(cè): 2008-02-21
- 性別: GG
- 專業(yè): 信息安全
8樓2011-03-11 19:09:39
change0618
鐵桿木蟲 (著名寫手)
方丈大師
- 應(yīng)助: 44 (小學(xué)生)
- 金幣: 17724.5
- 紅花: 17
- 帖子: 2413
- 在線: 546.7小時(shí)
- 蟲號(hào): 496517
- 注冊(cè): 2008-01-19
- 專業(yè): 化學(xué)反應(yīng)工程
9樓2011-03-11 20:03:41
zzg_gm
銀蟲 (小有名氣)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 331.4
- 散金: 21
- 帖子: 234
- 在線: 118.1小時(shí)
- 蟲號(hào): 1124750
- 注冊(cè): 2010-10-17
- 專業(yè): 基礎(chǔ)物理學(xué)
10樓2011-03-11 20:33:49
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 ybnav 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑 +6 | yangfz 2026-03-17 | 6/300 |
|
|
|
[考研] 328求調(diào)劑，英語六級(jí)551，有科研經(jīng)歷 +3 | 生物工程調(diào)劑 2026-03-16 | 8/400 |
|
|
|
[考研] 326求調(diào)劑 +5 | 上岸的小葡 2026-03-15 | 6/300 |
|
|
|
[考研] 274求調(diào)劑0856材料化工 +13 | z2839474511 2026-03-11 | 14/700 |
|
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學(xué)123 2026-03-15 | 8/400 |
|
|
|
[考研] 289求調(diào)劑 +6 | 步川酷紫123 2026-03-11 | 6/300 |
|
|
|
[考研] 283求調(diào)劑 +3 | 聽風(fēng)就是雨； 2026-03-16 | 3/150 |
|
|
|
[考研] 297一志愿上交085600求調(diào)劑 +5 | 指尖八千里 2026-03-14 | 5/250 |
|
|
|
[考研] 285 求調(diào)劑 資源與環(huán)境 一志愿北京化工大學(xué) +3 | 未名考生 2026-03-10 | 3/150 |
|
|
|
[考研] 材料專碩288分求調(diào)劑 一志愿211 +4 | 在家想你 2026-03-11 | 4/200 |
|
|
|
[考研] 332求調(diào)劑 +3 | zjy101327 2026-03-11 | 6/300 |
|
|
|
[考研] 材料工程調(diào)劑 +9 | 咪咪空空 2026-03-12 | 9/450 |
|
|
|
[考研] 26調(diào)劑/材料/英一數(shù)二/總分289/已過A區(qū)線 +6 | 步川酷紫123 2026-03-13 | 6/300 |
|
|
|
[考研] 一志愿西南交大，材料專碩317求調(diào)劑 +5 | lx8568 2026-03-11 | 5/250 |
|
|
|
[考研] 281求調(diào)劑 +9 | Koxui 2026-03-12 | 11/550 |
|
|
|
[考研] 【考研調(diào)劑求收留】 +3 | Ceciilia 2026-03-11 | 3/150 |
|
|
|
[考研] 308求調(diào)劑 +3 | 是Lupa啊 2026-03-12 | 3/150 |
|
|
|
[考研] 0856化學(xué)工程280分求調(diào)劑 +4 | shenzxsn 2026-03-11 | 4/200 |
|
|
|
[考研] 求調(diào)劑 資源與環(huán)境 285 +3 | 未名考生 2026-03-10 | 3/150 |
|