 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2510 | 回復(fù): 15 | ||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | ||
[求助]
matlab解微分方程組
|
||
|
哪位高手指點一下matlab解微分方程組！ 圖片是微分方程組。 下面是我根據(jù)網(wǎng)上查到的解法寫的程序，能運行，但是計算出來的圖像完全不對。 主函數(shù)： tspan=[0:0.01:2];%求解時間為2秒也可以 tspan=[0 2]表示自動時間步長 x0=[0;0;2.4*10^(-4);20;1.15*10^(-7);0];%初值 [t,x]=ode45(@ff,tspan,x0);%用龍格庫塔法求解 figure plot(t,x(:,1)) figure plot(t,x(:,2)) figure plot(t,x(:,3)) figure plot(t,x(:,4)) figure plot(t,x(:,5)) figure plot(t,x(:,6)) M文件 %微分方程 function dx=ff(t,x) dx=zeros(6,1); %x(1)=X1 %x(2)=X2 %x(3)=X3 %x(4)=X4 %x(5)=X5 %x(6)=X6 k1u=1.5*10^(-7) k1b=1.0*10^(-4) k2=1.3*10^(7) k3=7*10^(-1) k4=1.0*10^(10) k5=2.6*10^(9) k6=3*10^(3) Xc=0.5 MMT=0 dx(1)=2*k1u*x(3)+k1b*x(3)^2-k2*x(1)*x(5)+k3*x(2)*x(4)-2*k4*x(1)^2-k5*x(1)*x(2); dx(2)=k1b*x(3)^2+k2*x(1)*x(5)-k3*x(2)*x(4)-k5*x(1)*x(2)-2*k6*x(2); dx(3)=-k1u*x(3)-2*k1b*x(3)^2+k3*x(2)*x(4); dx(4)=-k1u*x(3)-k1b*x(3)^2-k3*x(2)*x(4); dx(5)=-k2*x(1)*x(5)+k6*x(2)^2; dx(6)=(k1u*x(3)+k1b*x(3)^2)*(1-Xc-MMT)*1.1111; end |
回復(fù)此樓» 猜你喜歡
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復(fù)
-
中德博士后2026項目5月30日截止- 國內(nèi)博士后入站、德國兩年全職
已經(jīng)有0人回復(fù)
-
物理學(xué)I論文潤色/翻譯怎么收費?
已經(jīng)有207人回復(fù)
-
面向全球申請人的中德博士后2026項目5月30日截止- 國內(nèi)博后入站德國兩年全職
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有2人回復(fù)
-
重慶交大材料學(xué)院26年碩士招生復(fù)試實施細則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。
已經(jīng)有0人回復(fù)
-
北工大宋曉艷教授團隊招收2026年9月入學(xué)博士【學(xué)博和專博名額均有】！
已經(jīng)有1人回復(fù)
-
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收10人
已經(jīng)有0人回復(fù)
-
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收調(diào)劑，10人
已經(jīng)有0人回復(fù)
-
計算與實驗科學(xué)工程國際會議(8月7-11日,香港) 二維材料分會S5與聲學(xué)換能器分會S8
已經(jīng)有0人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 用matlab求解一偏微分方程，求助！
已經(jīng)有1人回復(fù)
- 流體力學(xué)一維定常不可壓縮流體ns方程，已知壓差變化，怎么求解速度u關(guān)于x的常微分方程
已經(jīng)有0人回復(fù)
- matlab書籍
已經(jīng)有103人回復(fù)
- 求大神幫忙解題�。。�！
已經(jīng)有4人回復(fù)
- Matlab求解微分方程里面有參數(shù)帶入怎么求解，求救�。�！
已經(jīng)有0人回復(fù)
- 求助CFD軟件可以自定義模擬算法
已經(jīng)有0人回復(fù)
- 【緊急求助】 帶有微分平方項的ODE求解
已經(jīng)有11人回復(fù)
- 在控制器設(shè)計中關(guān)于大家如何給定？如果給定是方波怎么給？
已經(jīng)有2人回復(fù)
- 請教高手：PDE工具箱能解兩個偏微分方程相互迭代求解嗎
已經(jīng)有2人回復(fù)
- 【求助】求解形如y''+a*y'^2+by'+c=0 二階微分方程
已經(jīng)有14人回復(fù)
- 【資料】一些數(shù)學(xué)實驗,有程序源代碼
已經(jīng)有7人回復(fù)
3樓2011-12-02 09:47:19
buxiudehun
木蟲 (小有名氣)
- 應(yīng)助: 3 (幼兒園)
- 金幣: 1710.4
- 紅花: 1
- 帖子: 60
- 在線: 151.9小時
- 蟲號: 634805
- 注冊: 2008-10-23
- 專業(yè): 計算數(shù)學(xué)與科學(xué)工程計算
2樓2011-12-02 09:22:51
zt_chem
木蟲 (正式寫手)
- 應(yīng)助: 12 (小學(xué)生)
- 金幣: 2147.8
- 紅花: 6
- 帖子: 872
- 在線: 981.6小時
- 蟲號: 286920
- 注冊: 2006-10-15
- 性別: GG
- 專業(yè): 催化化學(xué)
4樓2011-12-03 10:26:51
5樓2011-12-03 11:26:16
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 279求調(diào)劑 +4 | 蝶舞輕繞 2026-03-29 | 4/200 |
|
|---|---|---|---|---|
|
|
[考研] 329求調(diào)劑 +7 | 星野? 2026-03-26 | 7/350 |
|
|
|
[考研] 283求調(diào)劑 +3 | A child 2026-03-28 | 3/150 |
|
|
[考研] 085600，材料與化工321分求調(diào)劑 +9 | 大饞小子 2026-03-28 | 9/450 |
|
|
|
[考研] 347求調(diào)劑 +3 | 山頂見α 2026-03-25 | 3/150 |
|
|
|
[考研] 299求調(diào)劑 +7 | 嗯嗯嗯嗯2 2026-03-27 | 7/350 |
|
|
|
[考研] 【求調(diào)劑】085601材料工程專碩 | 總分272 | +6 | 腳滑的守法公民 2026-03-27 | 6/300 |
|
|
|
[考研] 292求調(diào)劑 +14 | 鵝鵝鵝額額額額?/a> 2026-03-25 | 15/750 |
|
|
|
[考研] 315分求調(diào)劑 +7 | 26考研上岸版26 2026-03-26 | 7/350 |
|
|
|
[考研] 316求調(diào)劑 +5 | Pigcasso 2026-03-24 | 5/250 |
|
|
|
[考研] 298調(diào)劑 +3 | jiyingjie123 2026-03-27 | 3/150 |
|
|
|
[論文投稿] Journal of Mechanical Science and Technology +3 | Russ_ss 2026-03-25 | 5/250 |
|
|
|
[考研] 0703化學(xué)一志愿南京師范大學(xué)303求調(diào)劑 +3 | zzffylgg 2026-03-24 | 3/150 |
|
|
|
[考研] 求調(diào)劑323材料與化工 +7 | 1124361 2026-03-24 | 7/350 |
|
|
|
[考研] 一志愿陜師大生物學(xué)071000，298分，求調(diào)劑 +5 | SYA！ 2026-03-23 | 5/250 |
|
|
|
[考研] 調(diào)劑求收留 +7 | 果然有我 2026-03-26 | 7/350 |
|
|
|
[考研] 332求調(diào)劑 +6 | 032500 2026-03-25 | 6/300 |
|
|
|
[考研] 材料與化工328分調(diào)劑 +6 | 。，。，。，。i 2026-03-23 | 6/300 |
|
|
|
[考研] 求b區(qū)院校調(diào)劑 +4 | 周56 2026-03-24 | 5/250 |
|
|
|
[考研] 上海電力大學(xué)材料防護與新材料重點實驗室招收調(diào)劑研究生（材料、化學(xué)、電化學(xué)，環(huán)境） +4 | 我愛學(xué)電池 2026-03-23 | 4/200 |
|