版塊導航: 正在加載中...

登錄注冊

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

24小時熱門版塊排行榜

返回列表

change0618

鐵桿木蟲 (著名寫手)

方丈大師

計算強帖: 1
應助: 44 (小學生)
金幣: 17724.5
紅花: 17
帖子: 2413
在線: 546.7小時
蟲號: 496517
注冊: 2008-01-19
專業(yè): 化學反應工程

cB=????; % cB物料衡算表達
cW=????; % cW物料衡算表達

原來我給你寫過，這地方用物料衡算的計算式，就是將B，W的濃度用其它組份的濃度表示出來。你不理解啥意思嗎？

贊一下

回復此樓

11樓2012-03-28 18:20:48

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

lxyy

銀蟲 (小有名氣)

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

引用回帖:

11樓: Originally posted by change0618 at 2012-03-28 18:20:48:
cB=????; % cB物料衡算表達
cW=????; % cW物料衡算表達

原來我給你寫過，這地方用物料衡算的計算式，就是將B，W的濃度用其它組份的濃度表示出來。你不理解啥意思嗎？

這個我有試過，但是程序還是有問題，運行不了

贊一下

回復此樓

12樓2012-03-28 20:06:23

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

lxyy

銀蟲 (小有名氣)

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

引用回帖:

10樓: Originally posted by dingd at 2012-03-28 17:10:04:
供參考：

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 ...

您好，我想再請教下，用1stopt模擬的話，如果不知道k1+、k1-、k3+、k3-、之間的關系，這樣要求的參數(shù)里面就多了k1-、k3-、兩個參數(shù)，這樣擬合的結果會好嗎

贊一下

回復此樓

13樓2012-03-28 20:26:03

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

zl030091126

金蟲 (正式寫手)

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

【答案】應助回帖

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

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

贊一下

回復此樓

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

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

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

lxyy

銀蟲 (小有名氣)

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

引用回帖:

14樓: Originally posted by zl030091126 at 2012-03-28 22:16:11:
要說從程序上調整來使得對初值不依賴，那么就必須換掉調用的優(yōu)化函數(shù)了，你用的是這個函數(shù)lsqnonlin，可以試試別的函數(shù)，或者自己編個優(yōu)化函數(shù)。

這個我不是很懂，謝謝了

贊一下

回復此樓

15樓2012-03-28 22:22:22

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

fielding

木蟲 (正式寫手)

應助: 5 (幼兒園)
金幣: 3544.5
散金: 135
紅花: 4
帖子: 672
在線: 219.3小時
蟲號: 822444
注冊: 2009-08-05
性別: GG
專業(yè): 能源化工

【答案】應助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應助指數(shù) +1
lxyy: 金幣+10, ★★★很有幫助 2012-03-30 12:30:33

function KineticsEst
clear all
clc

x0 = [4.639975 5.289842 0 0.039974 0.063806];       % A B C R W 組成的初始濃
k0=[1 1 1]; % k1, k2, k3 初始向值
lb = [0 0 0];
ub = [1  +inf +inf];  % 上下限
KineticsData3;
yexp = Kinetics(:,2:6);
% 使用函數(shù)lsqnonlin()進行參數(shù)估計
[k,resnorm,residual,exitflag,output,lambda,jacobian]=...
lsqnonlin(@ObjFunc,k0,lb,ub,[],x0,yexp);

fprintf('\n\n使用函數(shù)lsqnonlin()估計得到的參數(shù)值為:\n')
fprintf('\tk1 = %.4f ± %.4f\n',k(1))
fprintf('\tk2 = %.4f ± %.4f\n',k(2))
fprintf('\tk3 = %.4f ± %.4f\n',k(3))
% ------------------------------------------------------------------
function f = ObjFunc(k,x0,yexp) % 目標函數(shù)
tspan = [0;5;10;15;30;46;60;90;120;160;200;240;300;360];
[t x]  = ode45(@Euqations,tspan,x0,[],k);
y(:,1) = x(:,1);
y(:,2:5) = x(:,2:5);
f1 = y(:,1) - yexp(:,1);
f2 = y(:,2) - yexp(:,2);
f3 = y(:,3) - yexp(:,3);
f4 = y(:,4) - yexp(:,4);
f5 = y(:,5) - yexp(:,5);
f = [f1; f2; f3; f4;f5];
% ------------------------------------------------------------------
function dcdt = Euqations(t,x,k)

r1 = k(1)*x(1)*x(2)-k(1)*x(3)/3.94;
r2 = k(2)*x(1)*x(5);
r3 = k(3)*x(2)*x(4)-k(3)*x(3)*x(5)/1.08;
dcAdt = -r1-r2;
dcBdt = -r1-r3;
dcCdt = r1+r3;
dcRdt = r2-r3;
dcWdt = -r2+r3;
dcdt = [dcAdt;dcBdt; dcCdt;dcRdt;dcWdt;];

贊一下

回復此樓

16樓2012-03-29 18:30:31

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

相關版塊跳轉我要訂閱樓主 lxyy 的主題更新

返回列表

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

[考研] 材料工程專碩調劑 +5	204818@lcx 2026-03-17	5/250	2026-03-17 17:27 by Little-xue
[考研] 材料專碩326求調劑 +6	墨煜姒莘 2026-03-15	7/350	2026-03-17 17:10 by ruiyingmiao
[考研] 312求調劑 +4	陌宸希 2026-03-16	5/250	2026-03-17 17:09 by ruiyingmiao
[考研] 274求調劑0856材料化工 +13	z2839474511 2026-03-11	14/700	2026-03-17 16:51 by share_joy
[考研] 26考研求調劑 +6	丶宏Sir 2026-03-13	6/300	2026-03-17 16:13 by 醉在風里
[考研] 材料與化工304求B區(qū)調劑 +7	邱gl 2026-03-11	8/400	2026-03-17 09:36 by 努力學習賺彩禮
[論文投稿] 有沒有大佬發(fā)小論文能帶我個二作 +3	增銳漏人 2026-03-17	4/200	2026-03-17 09:26 by xs74101122
[考研] 278求調劑 +3	Yy7400 2026-03-13	3/150	2026-03-17 08:24 by laoshidan
[考研] 304求調劑 +3	曼殊2266 2026-03-14	3/150	2026-03-16 16:39 by houyaoxu
[考研] 0703一志愿211 285分求調劑 +5	ly3471z 2026-03-13	5/250	2026-03-16 16:16 by 哦哦123
[考研] 285求調劑 +6	ytter 2026-03-12	6/300	2026-03-16 15:05 by njzyff
[考研] 26考研一志愿中國石油大學(華東)305分求調劑 +3	嘉年新程 2026-03-15	3/150	2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 085601材料工程315分求調劑 +3	yang_0104 2026-03-15	3/150	2026-03-15 10:58 by peike
[考研] 本科南京大學一志愿川大藥學327 +3	麥田耕者 2026-03-14	3/150	2026-03-14 20:04 by 外星文明
[考研] 材料371求調劑 +9	鱷魚? 2026-03-11	11/550	2026-03-13 22:53 by JourneyLucky
[考研] 材料專碩288分求調劑一志愿211 +4	在家想你 2026-03-11	4/200	2026-03-13 22:49 by JourneyLucky
[考研] 337一志愿華南理工0805材料求調劑 +7	mysdl 2026-03-11	9/450	2026-03-13 22:43 by JourneyLucky
[考研] 311求調劑 +3	冬十三 2026-03-13	3/150	2026-03-13 20:41 by JourneyLucky
[考研] 295求調劑 +3	小匕仔汁 2026-03-12	3/150	2026-03-13 15:17 by vgtyfty
[考研] 0817化學工程與技術考研312分調劑 +3	T123 tt 2026-03-12	3/150	2026-03-13 10:49 by houyaoxu

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

24小時熱門版塊排行榜

change0618

lxyy

lxyy

zl030091126

【答案】應助回帖

lxyy

fielding

【答案】應助回帖

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