_alarrn90

無蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: -28
紅花: 1
帖子: 291
在線: 118.5小時
蟲號: 1016060
注冊: 2010-05-11
性別: GG
專業(yè): 高分子材料的加工與成型

[求助] 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ù)此樓

» 猜你喜歡

求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
物理學I論文潤色/翻譯怎么收費? 已經(jīng)有160人回復(fù)
請問還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復(fù)
光學工程學碩調(diào)劑信息已經(jīng)有32人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運，祝面上順利！已經(jīng)有76人回復(fù)
華東師范大學芯片設(shè)計徐瓏真實水平咋樣已經(jīng)有1人回復(fù)
調(diào)劑已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有0人回復(fù)
中德博士后2026項目5月30日截止- 國內(nèi)博士后入站、德國兩年全職已經(jīng)有0人回復(fù)

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

用matlab求解一偏微分方程，求助！已經(jīng)有1人回復(fù)
流體力學一維定常不可壓縮流體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ù)學實驗,有程序源代碼已經(jīng)有7人回復(fù)

1樓 2011-12-01 18:53:32

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

buxiudehun

木蟲 (小有名氣)

應(yīng)助: 3 (幼兒園)
金幣: 1710.4
紅花: 1
帖子: 60
在線: 151.9小時
蟲號: 634805
注冊: 2008-10-23
專業(yè): 計算數(shù)學與科學工程計算

【答案】應(yīng)助回帖

★
dbb627(金幣+1): 歡迎交流 2011-12-02 09:29:51
_alarrn90(金幣+2): 2011-12-10 10:44:38

可能性很多，首先你確定沒有打錯公式，其次你可以選擇別的求解方法，不過這個問題看起來非線性性不是很強，步長小一些再試試。。。

贊一下(1人)

回復(fù)此樓

2樓2011-12-02 09:22:51

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

_alarrn90

無蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: -28
紅花: 1
帖子: 291
在線: 118.5小時
蟲號: 1016060
注冊: 2010-05-11
性別: GG
專業(yè): 高分子材料的加工與成型

多謝，我試一下。

回復(fù)此樓

3樓2011-12-02 09:47:19

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

zt_chem

木蟲 (正式寫手)

應(yīng)助: 12 (小學生)
金幣: 2147.8
紅花: 6
帖子: 872
在線: 981.6小時
蟲號: 286920
注冊: 2006-10-15
性別: GG
專業(yè): 催化化學

【答案】應(yīng)助回帖

_alarrn90(金幣+2): 2011-12-10 10:44:43

為什么在程序最后加上‘end’？加上以后程序運行不下去�。苛硗�，步長減小以后結(jié)果變化很大

贊一下

回復(fù)此樓

4樓2011-12-03 10:26:51

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

_alarrn90

無蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: -28
紅花: 1
帖子: 291
在線: 118.5小時
蟲號: 1016060
注冊: 2010-05-11
性別: GG
專業(yè): 高分子材料的加工與成型

多謝，但是我感覺加不加end沒什么影響吧。去掉了end計算出來的結(jié)果也不太對，你還看出其他什么問題沒有？

贊一下

回復(fù)此樓

5樓2011-12-03 11:26:16

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

dbb627

榮譽版主 (著名寫手)

專家經(jīng)驗: +4
計算強帖: 12
應(yīng)助: 289 (大學生)
貴賓: 0.589
金幣: 24640.4
散金: 551
紅花: 61
沙發(fā): 1
帖子: 1246
在線: 1794.8小時
蟲號: 149791
注冊: 2005-12-29
性別: GG
專業(yè): 污染控制化學
管轄: 計算模擬

【答案】應(yīng)助回帖

你的方程有些問題,氧氣的方程有關(guān)于空間的擴散項,沒包涵.第四個初值是pH吧,怎么取到20了,氧氣濃度也太低了吧,你是不是把條件寫錯了

贊一下

回復(fù)此樓

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.

6樓2011-12-09 23:39:40

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

_alarrn90

無蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: -28
紅花: 1
帖子: 291
在線: 118.5小時
蟲號: 1016060
注冊: 2010-05-11
性別: GG
專業(yè): 高分子材料的加工與成型

確實是氧氣濃度的初值有問題，我忘記乘上環(huán)境壓力，但是修改后的圖趨勢沒什么變化。擴散的影響在這里可以忽略。[PH]的初值確實是20.問一下您對照下面微分方程的話，程序本身有沒有錯誤？因為計算出來的圖的趨勢與文獻明顯不同。另外我想問一下初值以及步長的設(shè)定對圖的變化趨勢影響大不大，或者還有沒有其他影響圖形趨勢的因素？非常感謝。

贊一下

回復(fù)此樓

7樓2011-12-10 10:43:54

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

dbb627

榮譽版主 (著名寫手)

專家經(jīng)驗: +4
計算強帖: 12
應(yīng)助: 289 (大學生)
貴賓: 0.589
金幣: 24640.4
散金: 551
紅花: 61
沙發(fā): 1
帖子: 1246
在線: 1794.8小時
蟲號: 149791
注冊: 2005-12-29
性別: GG
專業(yè): 污染控制化學
管轄: 計算模擬

【答案】應(yīng)助回帖

引用回帖:

7樓: Originally posted by _alarrn90 at 2011-12-10 10:43:54:
確實是氧氣濃度的初值有問題，我忘記乘上環(huán)境壓力，但是修改后的圖趨勢沒什么變化。擴散的影響在這里可以忽略。[PH]的初值確實是20.問一下您對照下面微分方程的話，程序本身有沒有錯誤？因為計算出來的圖的趨勢與 ...

程序本身沒看出來啥問題，初值和步長對圖形趨勢很有影響的，這個步長取太大，有些極值，拐點就跳過了

贊一下

回復(fù)此樓

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.

8樓2011-12-10 11:08:07

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

_alarrn90

無蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: -28
紅花: 1
帖子: 291
在線: 118.5小時
蟲號: 1016060
注冊: 2010-05-11
性別: GG
專業(yè): 高分子材料的加工與成型

這是模擬出的圖和文獻中圖的對比

贊一下

回復(fù)此樓

9樓2011-12-10 11:09:40

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

_alarrn90

無蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: -28
紅花: 1
帖子: 291
在線: 118.5小時
蟲號: 1016060
注冊: 2010-05-11
性別: GG
專業(yè): 高分子材料的加工與成型

那您覺得怎么調(diào)整可能會好一些？

贊一下

回復(fù)此樓

10樓2011-12-10 11:14:33

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 307求調(diào)劑 +3	余意卿 2026-03-21	3/150	2026-03-23 10:32 by Iveryant
[考研] 352求調(diào)劑 +3	大米飯！ 2026-03-22	3/150	2026-03-22 23:28 by king123！
[考研] 310求調(diào)劑 +4	baibai1314 2026-03-16	4/200	2026-03-22 20:19 by edmund7
[考研] 資源與環(huán)境調(diào)劑申請(333分) +5	holy J 2026-03-21	5/250	2026-03-21 22:42 by Catalysis25
[考研] 求調(diào)劑 +4	要好好無聊 2026-03-21	4/200	2026-03-21 18:57 by 學員8dgXkO
[考研] 0805 316求調(diào)劑 +3	大雪深藏 2026-03-18	3/150	2026-03-21 18:55 by 學員8dgXkO
[考研] 297求調(diào)劑 +3	喜歡還是不甘心 2026-03-20	3/150	2026-03-21 18:33 by 學員8dgXkO
[考研] 297求調(diào)劑 +11	戲精丹丹丹 2026-03-17	12/600	2026-03-21 17:47 by ColorlessPI
[考研] 336求調(diào)劑 +5	rmc8866 2026-03-21	5/250	2026-03-21 17:24 by 學員8dgXkO
[考研] 266求調(diào)劑 +3	哇呼哼呼哼 2026-03-20	3/150	2026-03-21 16:46 by barlinike
[考研] 機械專碩299求調(diào)劑至材料 +3	kkcoco25 2026-03-16	4/200	2026-03-21 03:52 by JourneyLucky
[考研] 307求調(diào)劑 +3	wyyyqx 2026-03-17	3/150	2026-03-21 03:20 by JourneyLucky
[考研] 材料專業(yè)求調(diào)劑 +6	hanamiko 2026-03-18	6/300	2026-03-21 00:24 by JourneyLucky
[考研] 296求調(diào)劑 +6	www_q 2026-03-18	10/500	2026-03-20 23:56 by JourneyLucky
[考研] 321求調(diào)劑 +9	何潤采123 2026-03-18	11/550	2026-03-20 23:19 by JourneyLucky
[考研] 295復(fù)試調(diào)劑 +8	簡木ChuFront 2026-03-19	8/400	2026-03-20 20:44 by zhukairuo
[考研] 一志愿南京航空航天大學大學，080500材料科學與工程學碩 +5	@taotao 2026-03-20	5/250	2026-03-20 20:16 by JourneyLucky
[考研] 288求調(diào)劑，一志愿華南理工大學071005 +5	ioodiiij 2026-03-17	5/250	2026-03-19 18:22 by zcl123
[考研] 【同濟軟件】軟件（085405）考研求調(diào)劑 +3	2026eternal 2026-03-18	3/150	2026-03-18 19:09 by 搏擊518
[考研] 收復(fù)試調(diào)劑生 +4	雨后秋荷 2026-03-18	4/200	2026-03-18 14:16 by elevennnne