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

登錄注冊

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

24小時熱門版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖

lxyy

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 237.1
散金: 480
紅花: 9
帖子: 182
在線: 210.5小時
蟲號: 1131417
注冊: 2010-10-25
專業(yè): 化學(xué)反應(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ù)此樓

» 猜你喜歡

【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有2人回復(fù)
重慶交大材料學(xué)院26年碩士招生復(fù)試實施細則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。已經(jīng)有0人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費? 已經(jīng)有267人回復(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ù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有4人回復(fù)
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來能源學(xué)院已經(jīng)有10人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑本科期間發(fā)表一篇Sci 已經(jīng)有4人回復(fù)

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

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

zl030091126

金蟲 (正式寫手)

應(yīng)助: 15 (小學(xué)生)
金幣: 303.8
散金: 605
紅花: 2
帖子: 524
在線: 141.5小時
蟲號: 1444718
注冊: 2011-10-16
性別: GG
專業(yè): 化學(xué)反應(yīng)工程

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

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

要說從程序上調(diào)整來使得對初值不依賴，那么就必須換掉調(diào)用的優(yōu)化函數(shù)了，你用的是這個函數(shù)lsqnonlin，可以試試別的函數(shù)，或者自己編個優(yōu)化函數(shù)。

贊一下

回復(fù)此樓

A--B,虛線變成實線則是你的軌跡

14樓2012-03-28 22:16:11

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

查看全部 16 個回答

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è): 化學(xué)反應(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ù)學(xué)語言來描述，非化工專業(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的回帖

查看全部 16 個回答

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

[有機交流] 甲基亞磺磺酸鈉和甲基磺酸酯反應(yīng)機理 10+3	kaobao456 2026-03-29	4/200	2026-03-30 23:16 by nBu鋰
[考研] 085701環(huán)境工程求調(diào)劑 +11	多久上課 2026-03-27	12/600	2026-03-30 21:21 by 研究僧導(dǎo)導(dǎo)
[考研] 一志愿華東師范大學(xué)有機化學(xué)專業(yè)，初試351分，復(fù)試被刷求調(diào)劑! +7	真名有冰 2026-03-29	8/400	2026-03-30 21:08 by maaj02
[考研] 277跪求調(diào)劑 +8	1915668 2026-03-27	12/600	2026-03-30 21:01 by dophin1985
[考研] 一志愿北京化工大學(xué)材料與化工（085600）296求調(diào)劑 +25	稻妻小編 2026-03-26	25/1250	2026-03-30 20:11 by 滴滴上岸呀
[考研] 材料科學(xué)與工程 317求調(diào)劑 +7	JKSOIID 2026-03-26	7/350	2026-03-30 20:10 by 滴滴上岸呀
[考研] 367求調(diào)劑 +5	芋泥啵�！� 2026-03-28	5/250	2026-03-30 19:56 by 無際的草原
[考研] 08工科，295，接受跨專業(yè)調(diào)劑 +3	lmnlzy 2026-03-30	3/150	2026-03-30 17:49 by wangjy2002
[考研] 材料專碩調(diào)劑 +11	椰椰。 2026-03-29	11/550	2026-03-30 16:21 by wangjy2002
[考研] 070300化學(xué)354求調(diào)劑 +7	101次希望 2026-03-28	7/350	2026-03-30 13:22 by Iveryant
[考研] 311求調(diào)劑 +10	lin0039 2026-03-26	10/500	2026-03-30 10:26 by herarysara
[考研] 343求調(diào)劑085601 +3	要努力學(xué)習(xí)x 2026-03-29	3/150	2026-03-29 18:35 by wxiongid
[考研] 復(fù)試調(diào)劑 +3	raojunqi0129 2026-03-28	3/150	2026-03-28 15:27 by 落睿可思
[考研] 347求調(diào)劑 +3	山頂見α 2026-03-25	3/150	2026-03-28 14:13 by 唐沐兒
[考研] 286求調(diào)劑 +4	丟掉懶惰 2026-03-27	7/350	2026-03-28 08:07 by baoball
[考研] 0703化學(xué)求調(diào)劑，各位老師看看我�。。� +5	祁祺祺 2026-03-25	5/250	2026-03-27 21:44 by 東方豬豬
[考研] 一志愿西北大學(xué) 總分282 英語一62 求調(diào)劑 +7	18419759900 2026-03-25	8/400	2026-03-27 16:38 by 18419759900
[考研] 274求調(diào)劑 +17	顧九笙要謙虛 2026-03-24	23/1150	2026-03-27 15:16 by caszguilin
[考研] 一志愿吉大071010，316分求調(diào)劑 +3	xgbiknn 2026-03-27	3/150	2026-03-27 10:36 by guoweigw
[考研] 286求調(diào)劑 +4	lim0922 2026-03-26	4/200	2026-03-27 10:28 by guoweigw

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

24小時熱門版塊排行榜

lxyy

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

» 猜你喜歡

zl030091126

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

dingd

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

lxyy

» 本帖附件資源列表

dingd

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

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产高清中文字幕,99re热久久亚洲综合精品成人,熟妇一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频