版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

lxyy

銀蟲 (小有名氣)

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

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

A+B=C  A+W=R  R+B=C+W
先考慮的三個反應(yīng)都是可逆反應(yīng)，有以下幾個式子：編寫程序如下
%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 對應(yīng)時間數(shù)值

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

% 使用函數(shù)lsqnonlin()進行參數(shù)估計
[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) % 目標函數(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é)果，程序要怎么改一下初值對結(jié)果影響不大（這個程序有問題，）

回復(fù)此樓

» 猜你喜歡

2026年循環(huán)經(jīng)濟功能材料國際會議（ICFMCE 2026）已經(jīng)有0人回復(fù)
2026年第五屆電氣、電子與信息工程國際會議（ISEEIE 2026）已經(jīng)有0人回復(fù)
物理學I論文潤色/翻譯怎么收費? 已經(jīng)有251人回復(fù)
哈爾濱理工大學物理系招收物理學考研調(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ù)
請問還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復(fù)
光學工程學碩調(diào)劑信息已經(jīng)有26人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運，祝面上順利！已經(jīng)有36人回復(fù)

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

在matlab中擬合曲面的表達式已經(jīng)有29人回復(fù)
matlab 求指點動力學方程擬合過程中導(dǎo)數(shù)的獲取已經(jīng)有10人回復(fù)
matlab擬合擬均相動力學參數(shù) 已經(jīng)有6人回復(fù)
matlab-常微分方程參數(shù)估計已經(jīng)有12人回復(fù)
求教origin擬合簡單的滲流方程（又稱逾滲方程），做復(fù)合材料的應(yīng)該會遇到吧已經(jīng)有20人回復(fù)
非線性擬合估計方程參數(shù)（matlab，spss，origin都可以）已經(jīng)有4人回復(fù)
matlab非線性參數(shù)擬合問題已經(jīng)有7人回復(fù)
MATLAB求解一個比較復(fù)雜的方程已經(jīng)有14人回復(fù)
【求助】用matlab最優(yōu)化方法進行參數(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)動力學方程問題已經(jīng)有7人回復(fù)

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

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

dingd

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

計算強帖: 4
應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時
蟲號: 291104
注冊: 2006-10-28

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

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

微分方程dcBdt, dcWdt的表達式呢？

贊一下

回復(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小時
蟲號: 1131417
注冊: 2010-10-25
專業(yè): 化學反應(yīng)工程

引用回帖:

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

我寫在附件里面了

贊一下

回復(fù)此樓

» 本帖附件資源列表

歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺，不對該內(nèi)容負責。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識產(chǎn)權(quán)問題，其責任在于用戶本人，如對版權(quán)有異議，請聯(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的回帖

dingd

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

計算強帖: 4
應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時
蟲號: 291104
注冊: 2006-10-28

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

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

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

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

贊一下

回復(fù)此樓

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

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

lxyy

銀蟲 (小有名氣)

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

引用回帖:

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

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

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

時間：t0 = [5,10, 15, 30, 45, 60,90,120,160,200,240,300,360];
不同時刻對應(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)督和反饋：小木蟲僅提供交流平臺，不對該內(nèi)容負責。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識產(chǎn)權(quán)問題，其責任在于用戶本人，如對版權(quán)有異議，請聯(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的回帖

dingd

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

計算強帖: 4
應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時
蟲號: 291104
注冊: 2006-10-28

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

數(shù)據(jù)項B、W有對應(yīng)的微分方程嗎？還是自由變量？

贊一下

回復(fù)此樓

6樓2012-03-28 14:58:43

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

lxyy

銀蟲 (小有名氣)

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

引用回帖:

6樓: Originally posted by dingd at 2012-03-28 14:58:43:
數(shù)據(jù)項B、W有對應(yīng)的微分方程嗎？還是自由變量？

有的，我寫在附件里了，真是太謝謝了

贊一下

回復(fù)此樓

» 本帖附件資源列表

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

2012-03-28 15:05:19, 32.5 K

7樓2012-03-28 15:05:21

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

dingd

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

計算強帖: 4
應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時
蟲號: 291104
注冊: 2006-10-28

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

怎么感覺和擠牙膏似的��！

公式寫的對嗎？是否好好檢查過？和1樓的代碼有不少不同，比如對r2的定義，"dCB/wdt"應(yīng)該是“dCB/dt”嗎？公式（2）和（5）都給的是“dCB/dt”，那“dCC/dt”呢？

贊一下

回復(fù)此樓

8樓2012-03-28 16:30:05

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

lxyy

銀蟲 (小有名氣)

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

引用回帖:

8樓: Originally posted by dingd at 2012-03-28 16:30:05:
怎么感覺和擠牙膏似的�。�

公式寫的對嗎？是否好好檢查過？和1樓的代碼有不少不同，比如對r2的定義，"dCB/wdt"應(yīng)該是“dCB/dt”嗎？公式（2）和（5）都給的是“dCB/dt”，那“dCC/dt”呢？

dCB/wdt"應(yīng)該是“dCB/dt，沒錯
對不起我太心急了，公式已經(jīng)檢查過，我保證這次沒錯了,真的

贊一下

回復(fù)此樓

» 本帖附件資源列表

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

2012-03-28 16:38:07, 34.5 K

9樓2012-03-28 16:38:11

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

dingd

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

計算強帖: 4
應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時
蟲號: 291104
注冊: 2006-10-28

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

供參考：

CODE:

Constant keq1 = 3.94, keq3 = 1.08;

ConstStr k1minus= k1plus/keq1,k3minus= k3plus/keq3,

         r1=k1plus*cA*Cb-k1minus*cC,

         r2=k2*cA*cW,

         r3=k3plus*cB*cR-k3minus*cC*cW;

Variable t,ca,cb,cc,cr,cw;

ODEFunction cA' = -r1-r2;

            cB' = -r1-r3;

            cC' = r1+r3;

            cR' = r2-r3;

            cW' = r3-r2;

Data;

//t,cA,cB,cC,cR,cW

0        5.939601        5.203973        0        0.0282484        0.078956

5        3.742508811        4.213566482        0.990406938        0.252498405        0.090415512

10        3.045716838        3.387013322        1.816960097        0.215546741        0.081072022

15        2.504754523        2.901991595        2.301981824        0.193427102        0.099462604

30        1.65098702        2.036686653        3.167286766        0.170236582        0.095144044

45        1.271633833        1.655460403        3.548513016        0.124730781        0.10334349

60        1.129887499        1.630815394        3.573158025        0.129798695        0.097764543

90        1.01239208        1.460748979        3.74322444        0.144649391        0.073614958

120        1.007417603        1.429180292        3.774793127        0.123477124        0.094249307

160        0.990005723        1.499216939        3.704756481        0.135524067        0.076072022

200        0.970073726        1.493966669        3.71000675        0.130791851        0.076387812

240        0.96903024        1.429705429        3.77426799        0.119510017        0.101202216

300        0.980873385        1.457988124        3.745985296        0.124189274        0.097022161

360        0.98477262        1.453402965        3.750570455        0.1200936        0.0762608

均方差(RMSE): 0.384192283597719
殘差平方和(SSE): 9.59424120044197
相關(guān)系數(shù)(R): 0.998354081222552
相關(guān)系數(shù)之平方(R^2): 0.996710871493726
決定系數(shù)(DC): 0.995560449063788

參數(shù) 最佳估算
-------------------- -------------
k1plus 0.0122448530178222
k2 0.044692817673755
k3plus -0.0135416873776772

贊一下

回復(fù)此樓

10樓2012-03-28 17:10:04

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

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 lxyy 的主題更新

返回列表

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

[考研] 301求調(diào)劑 +6	yy要上岸呀 2026-03-17	6/300	2026-03-17 23:58 by 星空星月
[考研] 311求調(diào)劑 +9	冬十三 2026-03-15	10/500	2026-03-17 22:53 by lbsjt
[考研] 303求調(diào)劑 +3	睿08 2026-03-17	5/250	2026-03-17 22:01 by 睿08
[考研] 296求調(diào)劑 +5	大口吃飯身體健 2026-03-13	5/250	2026-03-17 21:05 by 不惑可樂
[考研] 326求調(diào)劑 +5	上岸的小葡 2026-03-15	6/300	2026-03-17 17:26 by ruiyingmiao
[考研] 考研化學學碩調(diào)劑，一志愿985 +4	張vvvv 2026-03-15	6/300	2026-03-17 17:15 by ruiyingmiao
[考研] 化學工程321分求調(diào)劑 +11	大米飯！ 2026-03-15	14/700	2026-03-17 17:11 by ruiyingmiao
[考研] 289求調(diào)劑 +6	步川酷紫123 2026-03-11	6/300	2026-03-17 10:23 by Sammy2
[考研] 286求調(diào)劑 +3	lemonzzn 2026-03-16	5/250	2026-03-16 20:43 by lemonzzn
[考研] 333求調(diào)劑 +3	文思客 2026-03-16	7/350	2026-03-16 18:21 by 文思客
[考研] 326求調(diào)劑 +4	諾貝爾化學獎覬?/a> 2026-03-15	7/350	2026-03-16 17:11 by 諾貝爾化學獎覬?/a>
[考研] 0703化學調(diào)劑，求各位老師收留 +8	秋有木北 2026-03-14	8/400	2026-03-16 15:21 by 哦哦123
[考研] 085600材料與化工求調(diào)劑 +13	enenenhui 2026-03-13	14/700	2026-03-16 15:19 by 了了了了。。
[考研] 327求調(diào)劑 +6	拾光任染 2026-03-15	11/550	2026-03-15 22:47 by 拾光任染
[考研] 070305求調(diào)劑 +3	mlpqaz03 2026-03-14	4/200	2026-03-15 11:04 by peike
[考研] 297求調(diào)劑 +4	學海漂泊 2026-03-13	4/200	2026-03-14 11:51 by 熱情沙漠
[考研] 招收0805（材料）調(diào)劑 +3	18595523086 2026-03-13	3/150	2026-03-14 00:33 by 123%、
[考研] 求材料調(diào)劑 085600英一數(shù)二總分302 前三科235 精通機器學習一志愿哈工大 +4	林yaxin 2026-03-12	4/200	2026-03-13 22:04 by 星空星月
[考研] 材料專碩350 求調(diào)劑 +4	王金科 2026-03-12	4/200	2026-03-13 16:02 by ruiyingmiao
[考研] 308求調(diào)劑 +3	是Lupa啊 2026-03-12	3/150	2026-03-13 14:30 by 求調(diào)劑zz

24小時熱門版塊排行榜

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

» 猜你喜歡

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

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

» 本帖附件資源列表

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

» 本帖附件資源列表

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

» 本帖附件資源列表

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

» 本帖附件資源列表

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

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