當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

lxyy

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 237.1
散金: 480
紅花: 9
帖子: 182
在線: 210.5小時(shí)
蟲號(hào): 1131417
注冊(cè): 2010-10-25
專業(yè): 化學(xué)反應(yīng)工程

[求助] matlab擬合方程參數(shù)時(shí)初值的選擇

A+B=C  A+W=R  R+B=C+W
先考慮的三個(gè)反應(yīng)都是可逆反應(yīng)，有以下幾個(gè)式子：編寫程序如下
%30℃
clear all
clc
global keq1 keq2 t0 c0

keq1 = 3.94;       % 平衡常數(shù)k1
keq3 = 1.08; % 平衡常數(shù)k3

t0 = [5,10, 15, 30, 45, 60,90,120,160,200,240,300,360];
c0 = [5.939601,5.203973,0,0.022484,0.078956]       % A B C R W 組成的初始濃度
ci=[3.742508811 4.213566482 0.990406938 0.252498405 0.090415512;
3.045716838 3.387013322 1.816960097 0.215546741 0.081072022;
2.504754523 2.901991595 2.301981824 0.193427102 0.099462604;
1.65098702 2.036686653 3.167286766 0.170236582 0.095144044;
1.271633833 1.655460403 3.548513016 0.124730781 0.10334349;
1.129887499 1.630815394 3.573158025 0.129798695 0.097764543;
1.01239208 1.460748979 3.74322444 0.144649391 0.073614958;
1.007417603 1.429180292 3.774793127 0.123477124 0.094249307;
0.990005723 1.499216939 3.704756481 0.135524067 0.076072022;
0.970073726 1.493966669 3.71000675 0.130791851 0.076387812;
0.96903024 1.429705429 3.77426799 0.119510017 0.101202216;
0.980873385 1.457988124 3.745985296 0.124189274 0.097022161;
0.98477262 1.453402965 3.750570455 0.1200936 0.0762608];

% cA cC cR 對(duì)應(yīng)時(shí)間數(shù)值

k0=[1 1 1]; % k1+, k2, k3+ 初始向值
lb = [0 0 0];
ub = [1000  1000 500];  % 上下限

% 使用函數(shù)lsqnonlin()進(jìn)行參數(shù)估計(jì)
[k,resnorm,residual,exitflag] = lsqnonlin(@ObjFunc,k0,lb,ub,[],ci);

k1plus=k(1);
k1minus= k1plus/keq1;
k2=k(2);
k3plus=k(3);
k3minus= k3plus/keq3;

% ------------------------------------------------------------------
function f = ObjFunc(k,ci) % 目標(biāo)函數(shù)
global t0  c0

[t,c_cal] = ode45(@Euqations,t0,c0,[],k);
f= sum((c_cal-ci)^2);

% ------------------------------------------------------------------
function dcdt = Euqations(t,c,k)
global keq1 keq2

k1plus=k(1);
k1minus= k1plus/keq1;
k2=k(2);
k3plus=k(3);
k3minus= k3plus/keq3;

cA=c(:,1);
cB=c(:,2);
cC=c(:,3);
cR=c(:,4);
cW=c(:,5);

dcAdt = k1plus*cA*cB-k1minus*cC+k2*cA*cW;
dcCdt = -k1plus*cA*cB+k1minus*cC+k3plus*cB*cR-k3minus*cC*cW;
dcRdt = k2*cA*cW-k3plus*cB*cR+k3minus*cC*cW;
dcdt = [dcAdt; dcCdt;dcRdt];
但是k的初值一直影響結(jié)果，程序要怎么改一下初值對(duì)結(jié)果影響不大（這個(gè)程序有問題，）

回復(fù)此樓

» 猜你喜歡

2026年循環(huán)經(jīng)濟(jì)功能材料國(guó)際會(huì)議（ICFMCE 2026）已經(jīng)有0人回復(fù)
2026年第五屆電氣、電子與信息工程國(guó)際會(huì)議（ISEEIE 2026）已經(jīng)有0人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有227人回復(fù)
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑已經(jīng)有5人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有5人回復(fù)
求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
請(qǐng)問還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有24人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運(yùn)，祝面上順利！已經(jīng)有17人回復(fù)

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

在matlab中擬合曲面的表達(dá)式已經(jīng)有29人回復(fù)
matlab 求指點(diǎn) 動(dòng)力學(xué)方程擬合過程中導(dǎo)數(shù)的獲取已經(jīng)有10人回復(fù)
matlab擬合擬均相動(dòng)力學(xué)參數(shù) 已經(jīng)有6人回復(fù)
matlab-常微分方程參數(shù)估計(jì) 已經(jīng)有12人回復(fù)
求教origin擬合簡(jiǎn)單的滲流方程（又稱逾滲方程），做復(fù)合材料的應(yīng)該會(huì)遇到吧已經(jīng)有20人回復(fù)
非線性擬合估計(jì)方程參數(shù)（matlab，spss，origin都可以）已經(jīng)有4人回復(fù)
matlab非線性參數(shù)擬合問題已經(jīng)有7人回復(fù)
MATLAB求解一個(gè)比較復(fù)雜的方程已經(jīng)有14人回復(fù)
【求助】用matlab最優(yōu)化方法進(jìn)行參數(shù)擬合已經(jīng)有17人回復(fù)
【求助】使用MATLAB擬合出按照公式的系數(shù) 已經(jīng)有19人回復(fù)
如何用MATLAB擬合指數(shù)函數(shù) 已經(jīng)有3人回復(fù)
【求助】最小二乘法編寫指數(shù)擬合方程已經(jīng)有19人回復(fù)
【求助】使用Matlab擬合反應(yīng)動(dòng)力學(xué)方程問題已經(jīng)有7人回復(fù)

1樓 2012-03-28 10:24:19

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

dingd

鐵桿木蟲 (職業(yè)作家)

計(jì)算強(qiáng)帖: 4
應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時(shí)
蟲號(hào): 291104
注冊(cè): 2006-10-28

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

★
zhangguangping: 金幣+1, 謝謝指導(dǎo)！ 2012-03-28 14:16:22

最好把你的專業(yè)問題用通俗的數(shù)學(xué)語言來描述，非化工專業(yè)的很難理解的；
你的問題應(yīng)該就是一微分方程擬合問題，僅需描述清楚：

1：微分方程組表達(dá)式
2：對(duì)應(yīng)的數(shù)據(jù)
3：要求的參數(shù)

贊一下

回復(fù)此樓

4樓2012-03-28 11:07:18

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

查看全部 16 個(gè)回答

dingd

鐵桿木蟲 (職業(yè)作家)

計(jì)算強(qiáng)帖: 4
應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時(shí)
蟲號(hào): 291104
注冊(cè): 2006-10-28

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

感謝參與，應(yīng)助指數(shù) +1

微分方程dcBdt, dcWdt的表達(dá)式呢？

贊一下

回復(fù)此樓

2樓2012-03-28 10:37:30

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

lxyy

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 237.1
散金: 480
紅花: 9
帖子: 182
在線: 210.5小時(shí)
蟲號(hào): 1131417
注冊(cè): 2010-10-25
專業(yè): 化學(xué)反應(yīng)工程

引用回帖:

2樓: Originally posted by dingd at 2012-03-28 10:37:30:
微分方程dcBdt, dcWdt的表達(dá)式呢？

我寫在附件里面了

贊一下

回復(fù)此樓

» 本帖附件資源列表

歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺(tái)，不對(duì)該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識(shí)產(chǎn)權(quán)問題，其責(zé)任在于用戶本人，如對(duì)版權(quán)有異議，請(qǐng)聯(lián)系郵箱：xiaomuchong@tal.com
附件 1 : 太殘忍doc.doc

2012-03-28 10:52:17, 30.5 K

3樓2012-03-28 10:52:39

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

lxyy

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 237.1
散金: 480
紅花: 9
帖子: 182
在線: 210.5小時(shí)
蟲號(hào): 1131417
注冊(cè): 2010-10-25
專業(yè): 化學(xué)反應(yīng)工程

引用回帖:

4樓: Originally posted by dingd at 2012-03-28 11:07:18:
最好把你的專業(yè)問題用通俗的數(shù)學(xué)語言來描述，非化工專業(yè)的很難理解的；
你的問題應(yīng)該就是一微分方程擬合問題，僅需描述清楚：

1：微分方程組表達(dá)式
2：對(duì)應(yīng)的數(shù)據(jù)
3：要求的參數(shù)

不好意思，我的微分方程組就是附件里面標(biāo)（1）（2）（3）的三個(gè)式子，對(duì)應(yīng)的數(shù)據(jù)為

時(shí)間：t0 = [5,10, 15, 30, 45, 60,90,120,160,200,240,300,360];
不同時(shí)刻對(duì)應(yīng)的數(shù)據(jù)
c0 = [5.939601,5.203973,0,0.022484,0.078956]       % A B C R W 組成的初始濃度
ci=[3.742508811       4.213566482       0.990406938       0.252498405       0.090415512;
3.045716838       3.387013322       1.816960097       0.215546741       0.081072022;
2.504754523       2.901991595       2.301981824       0.193427102       0.099462604;
1.65098702       2.036686653       3.167286766       0.170236582       0.095144044;
1.271633833       1.655460403       3.548513016       0.124730781       0.10334349;
1.129887499       1.630815394       3.573158025       0.129798695       0.097764543;
1.01239208       1.460748979       3.74322444       0.144649391       0.073614958;
1.007417603       1.429180292       3.774793127       0.123477124       0.094249307;
0.990005723       1.499216939       3.704756481       0.135524067       0.076072022;
0.970073726       1.493966669       3.71000675       0.130791851       0.076387812;
0.96903024       1.429705429       3.77426799       0.119510017       0.101202216;
0.980873385       1.457988124       3.745985296       0.124189274       0.097022161;
0.98477262       1.453402965       3.750570455       0.1200936       0.0762608];
要求的是ki+,k2,k3+,
另外方程中的ki-=ki+/3.14,k3-=k3+/1.08

贊一下

回復(fù)此樓

» 本帖附件資源列表

歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺(tái)，不對(duì)該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識(shí)產(chǎn)權(quán)問題，其責(zé)任在于用戶本人，如對(duì)版權(quán)有異議，請(qǐng)聯(lián)系郵箱：xiaomuchong@tal.com
附件 1 : 太殘忍doc (1)1.doc

2012-03-28 12:18:21, 26 K

5樓2012-03-28 12:18:31

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

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

[考研] 一志愿蘇州大學(xué)材料工程（085601）專碩有科研經(jīng)歷三項(xiàng)國(guó)獎(jiǎng)兩個(gè)實(shí)用型專利一項(xiàng)省級(jí)立項(xiàng) +6	大火山小火山 2026-03-16	8/400	2026-03-17 15:05 by 無懈可擊111
[考研] 08工科 320總分求調(diào)劑 +4	梨花珞晚風(fēng) 2026-03-17	4/200	2026-03-17 13:38 by houyaoxu
[考研] 一志愿南京大學(xué)，080500材料科學(xué)與工程，調(diào)劑 +4	Jy? 2026-03-16	4/200	2026-03-17 11:02 by gaoqiong
[考研] 考研調(diào)劑 +3	淇ya_~ 2026-03-17	5/250	2026-03-17 09:25 by Winj1e
[考研] 一志愿，福州大學(xué)材料專碩339分求調(diào)劑 +3	木子momo青爭(zhēng) 2026-03-15	3/150	2026-03-17 07:52 by laoshidan
[考研] 321求調(diào)劑 +5	大米飯！ 2026-03-15	5/250	2026-03-16 16:33 by houyaoxu
[考研] 283求調(diào)劑 +10	小樓。 2026-03-12	14/700	2026-03-16 16:08 by 13811244083
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +5	d如愿上岸 2026-03-12	8/400	2026-03-16 15:19 by peike
[考研] 0703 物理化學(xué)調(diào)劑 +3	我可以上岸的對(duì)?/a> 2026-03-13	5/250	2026-03-16 10:50 by 我可以上岸的對(duì)?/a>
[教師之家] 焦慮 +7	水冰月月野兔 2026-03-13	9/450	2026-03-16 10:00 by Quakerbird
[考研] 327求調(diào)劑 +6	拾光任染 2026-03-15	11/550	2026-03-15 22:47 by 拾光任染
[基金申請(qǐng)] 現(xiàn)在如何回避去年的某一個(gè)專家，不知道名字 +3	zk200107 2026-03-12	6/300	2026-03-14 17:13 by zk200107
[考研] 復(fù)試調(diào)劑 +4	z1z2z3879 2026-03-14	5/250	2026-03-14 16:30 by JourneyLucky
[考研] 341求調(diào)劑 +3	番茄頭--- 2026-03-10	3/150	2026-03-13 23:07 by JourneyLucky
[考研] 材料工程調(diào)劑 +9	咪咪空空 2026-03-12	9/450	2026-03-13 22:05 by 星空星月
[考研] 工科，求調(diào)劑 +3	我887 2026-03-11	3/150	2026-03-13 21:39 by JourneyLucky
[考研] 材料工程調(diào)劑 +4	咪咪空空 2026-03-11	4/200	2026-03-13 19:57 by JourneyLucky
[考研] 求b區(qū)學(xué)校調(diào)劑 +3	周56 2026-03-11	3/150	2026-03-13 16:20 by JourneyLucky
[考研] 工科278分求調(diào)劑 +5	周慢熱啊 2026-03-12	7/350	2026-03-13 15:49 by JourneyLucky
[考研] 333求調(diào)劑 +3	152697 2026-03-12	4/200	2026-03-13 07:08 by Iveryant