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

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

24小時(shí)熱門版塊排行榜

返回列表

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

lxyy

銀蟲(chóng) (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 237.1
散金: 480
紅花: 9
帖子: 182
在線: 210.5小時(shí)
蟲(chóng)號(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è)式子：編寫(xiě)程序如下
%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è)程序有問(wèn)題，）

回復(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)有94人回復(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)問(wèn)還有沒(méi)有用Latex寫(xiě)文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有25人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運(yùn)，祝面上順利！已經(jīng)有28人回復(fù)

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

在matlab中擬合曲面的表達(dá)式已經(jīng)有29人回復(fù)
matlab 求指點(diǎn) 動(dòng)力學(xué)方程擬合過(guò)程中導(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ù)擬合問(wèn)題已經(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ù)
【求助】最小二乘法編寫(xiě)指數(shù)擬合方程已經(jīng)有19人回復(fù)
【求助】使用Matlab擬合反應(yīng)動(dòng)力學(xué)方程問(wèn)題已經(jīng)有7人回復(fù)

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

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

lxyy

銀蟲(chóng) (小有名氣)

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

引用回帖:

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

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

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

贊一下

回復(fù)此樓

» 本帖附件資源列表

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

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

dingd

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

計(jì)算強(qiáng)帖: 4
應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時(shí)
蟲(chóng)號(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

銀蟲(chóng) (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 237.1
散金: 480
紅花: 9
帖子: 182
在線: 210.5小時(shí)
蟲(chóng)號(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á)式呢？

我寫(xiě)在附件里面了

贊一下

回復(fù)此樓

» 本帖附件資源列表

歡迎監(jiān)督和反饋：小木蟲(chóng)僅提供交流平臺(tái)，不對(duì)該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識(shí)產(chǎn)權(quán)問(wè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的回帖

dingd

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

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

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

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

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

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è)回答

普通表情龍兔虎貓

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

[考研] 301求調(diào)劑 +5	yy要上岸呀 2026-03-17	5/250	2026-03-17 20:20 by peike
[考研] 293求調(diào)劑 +7	zjl的號(hào) 2026-03-16	12/600	2026-03-17 18:22 by 重科小霸王
[考研] 070300化學(xué)319求調(diào)劑 +4	錦鯉0909 2026-03-17	4/200	2026-03-17 18:21 by 重科小霸王
[考研] 302求調(diào)劑 +9	負(fù)心者當(dāng)誅 2026-03-11	9/450	2026-03-17 17:13 by ruiyingmiao
[考研] 268求調(diào)劑 +6	好運(yùn)連綿不絕 2026-03-12	7/350	2026-03-17 14:56 by 呦呦憂郁
[碩博家園] 深圳大學(xué)碩士招生（2026秋，傳感器方向，僅錄取第一志愿） +4	xujiaoszu 2026-03-11	9/450	2026-03-17 10:29 by xujiaoszu
[考研] 278求調(diào)劑 +3	Yy7400 2026-03-13	3/150	2026-03-17 08:24 by laoshidan
[考研] 一志愿，福州大學(xué)材料專碩339分求調(diào)劑 +3	木子momo青爭(zhēng) 2026-03-15	3/150	2026-03-17 07:52 by laoshidan
[考研] 297一志愿上交085600求調(diào)劑 +5	指尖八千里 2026-03-14	5/250	2026-03-14 17:26 by a不易
[考研] 328求調(diào)劑 +3	5201314Lsy！ 2026-03-13	6/300	2026-03-14 15:31 by hyswxzs
[考研] 330求調(diào)劑 +3	?醬給調(diào)劑跪了 2026-03-13	3/150	2026-03-14 10:13 by JourneyLucky
[考研] 279求調(diào)劑 +3	抓著星星的女孩 2026-03-10	3/150	2026-03-13 23:47 by userper
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +6	邱gl 2026-03-12	7/350	2026-03-13 23:24 by 邱gl
[考研] 一志愿中科院，化學(xué)方向，295求調(diào)劑 +4	一氧二氮 2026-03-11	4/200	2026-03-13 22:35 by JourneyLucky
[考研] 26調(diào)劑/材料/英一數(shù)二/總分289/已過(guò)A區(qū)線 +6	步川酷紫123 2026-03-13	6/300	2026-03-13 21:59 by 星空星月
[考研] 求調(diào)劑 +7	18880831720 2026-03-11	7/350	2026-03-13 16:10 by JourneyLucky
[考研] 考研調(diào)劑 +4	芬達(dá)46 2026-03-12	4/200	2026-03-13 16:04 by ruiyingmiao
[考研] 土木第一志愿276求調(diào)劑，科研和技能十分豐富，求新興方向的導(dǎo)師收留 +3	土木小天才 2026-03-12	3/150	2026-03-13 15:01 by JourneyLucky
[考研] 328化工專碩求調(diào)劑 +4	。，。，。，。i 2026-03-12	4/200	2026-03-13 14:44 by JourneyLucky
[考研] 0817化學(xué)工程與技術(shù)考研312分調(diào)劑 +3	T123 tt 2026-03-12	3/150	2026-03-13 10:49 by houyaoxu