版塊導(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í)熱門(mén)版塊排行榜

返回列表

清輝12365

新蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 23
帖子: 21
在線: 53.2小時(shí)
蟲(chóng)號(hào): 2804052
注冊(cè): 2013-11-15
專業(yè): 化學(xué)反應(yīng)工程

[求助] 動(dòng)力學(xué)參數(shù)擬合

原料是：甲醇、CO、氧氣。主反應(yīng)：甲醇、CO、氧氣合成DMC，副反應(yīng)：CO和氧氣合成CO2.
要干的活：擬合出生成DMC和CO2的動(dòng)力學(xué)參數(shù)（指前因子、活化能、級(jí)數(shù)）
動(dòng)力學(xué)方程為dy1/dw=f(y1、y2)和dy2/dw=f(y1、y2),具體形式見(jiàn)附圖。
整理后的實(shí)驗(yàn)數(shù)據(jù)：N0(原料初始的摩爾數(shù))、P（壓力）、t(溫度)、x1(甲醇初始的摩爾分?jǐn)?shù))、x2(CO初始的摩爾分?jǐn)?shù))、x3(O2初始的摩爾分?jǐn)?shù))、y1(DMC產(chǎn)物的摩爾分?jǐn)?shù))、y2(O2產(chǎn)物的摩爾分?jǐn)?shù))、W=0.375(催化劑的質(zhì)量)。
關(guān)鍵難度在：主要是兩個(gè)常微分方程、數(shù)據(jù)多不知如何進(jìn)行運(yùn)算。下面是自己編的程序，出現(xiàn)問(wèn)題并且解決不上來(lái)，希望大神指教
function YYHJIA
clear all;clc
global t T P N0 X1 X2 X3 ;
t=[140 140 140.2 143.4 150 150.1 150.2 151.2 153.7 159.2 160 161.3 140.2 141.1 149.7 150.3 150.8 159.7 161.2 138 150.3 160 160.2 162.3];
%溫度
T=t+273.15;
P=[0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.3 0.3 0.3 0.3 0.3];
%壓力
N0=[0.3091 0.3260 0.2796 0.2856 0.2801 0.3260 0.3091 0.2853 0.3190 0.2853 0.3962 0.3190 0.3085 0.2793 0.2804 0.2856 0.3083 0.3085 0.2788 0.2368 0.26 0.2597 0.3207 0.2344];
%原料初始摩爾
x1=[0.624 0.5 0.5310 0.5710 0.53 0.5 0.624 0.572 0.511 0.572 0.524 0.511 0.625 0.531 0.529 0.571 0.625 0.625 0.532 0.501 0.57 0.571 0.555 0.506];
%原料甲醇占的摩爾分?jǐn)?shù)
x2=[0.25 0.334 0.349 0.287 0.348 0.334 0.25 0.286 0.345 0.286 0.352 0.345 0.25 0.349 0.348 0.286 0.25 0.25 0.35 0.334 0.287 0.288 0.334 0.337];
%原料CO摩爾分?jǐn)?shù)
x3=[0.127 0.166 0.121 0.142 0.122 0.166 0.127 0.142 0.144 0.142 0.124 0.144 0.125 0.12 0.123 0.143 0.124 0.125 0.118 0.165 0.142 0.141 0.111 0.157];
%原料O2的摩爾分?jǐn)?shù)
Y1=[0.015 0.016 0.015 0.017 0.02 0.021 0.021 0.021 0.021 0.025 0.023 0.025 0.018 0.018 0.025 0.026 0.026 0.03 0.03 0.023 0.028 0.036 0.034 0.038];
%出口DMC的摩爾數(shù)
Y2=[0.018 0.02 0.018 0.033 0.04 0.041 0.042 0.042 0.045 0.054 0.053 0.054 0.031 0.032 0.05 0.051 0.05 0.059 0.058 0.043 0.051 0.063 0.061 0.066];
%出口CO的摩爾數(shù)
k0=[8e4 2e4 -0.5 0.5 1 4e5 3e4 0.5 0.5 ];%(指前因子、活化能、級(jí)數(shù))初始量
lb=[0 0 -3 -3 -3 0 0 -3 -3 ];%下限
ub=[+inf +inf +inf +inf +inf +inf +inf +inf +inf];%上限
%使用lsqnonlin函數(shù)進(jìn)行參數(shù)估計(jì)
[k,resnorm,residual,exitflag,output,lambda,jacobian]=lsqnonlin(@ObjFunc,k0,lb,ub,[],T,P,N0,X1,X2,X3,Y1,Y2);
ci=nlparci(k,residual,jacobian);
function f=ObjFunc(k,T,P,N0,X1,X2,X3,Y1,Y2)
wspan=[0 0.375];
Y0=[Y1' Y2']
[m,n]=size(Y0);
Ycal = zeros(m,n);
for i=1:m
[w yy]=ode45(@Equations,wspan,Y0(i,

,[],k);
Ycal(i,

=yy(end,

;
end
f=Ycal-Y0;
function dydw=Equations(w,yy,k)
%如何表達(dá)呢？這么多數(shù)字
r1=k(1)*exp*(-k(2)./T).*P.^(k(3)+k(4)+k(5))*...
(yy(1).*(1.5*x1-2)+0.5*yy(2).*x1).^(k(3))*...
(yy(1).*(1.5*x2-1)+yy(2).*(0.5*x2-1)+x2).^(k(4))*...
(yy(1).*(1.5*x3-0.5)+yy(2).*(0.3*x3-0.5)+x3).^k(5);
r2=k(6)*exp*(-k(7)./T).*P.^(k(8)+k(9))*...
(yy(1).*(1.5*x2-1)+yy(2).*(0.5*x2-1)+x2).^k(8)*...
(yy(1).*(1.5*x3-0.5)+yy(2).*(0.3*x3-0.5)+x3).^k(9);
dy1dw=(1+1.5*yy(1)+0.5*yy(2)).*(r1.*(1.5*yy(1)+1)+0.5*yy(1).*r2)./N0;
dy2dw=(1+1.5*yy(1)+0.5*yy(2)).*(r2.*(1+0.5*yy(2))+1.5*r1.*yy(2))./N0;
dydw=[dy1dw;dy2dw];
出現(xiàn)的問(wèn)題如圖：

出現(xiàn)的問(wèn)題.jpg

方程形式.jpg

方程具體形式2.jpg

回復(fù)此樓

» 猜你喜歡

307求調(diào)劑已經(jīng)有9人回復(fù)
288求調(diào)劑已經(jīng)有15人回復(fù)
085700資源與環(huán)境308求調(diào)劑已經(jīng)有5人回復(fù)
294求調(diào)劑材料與化工專碩已經(jīng)有14人回復(fù)
0817 化學(xué)工程 299分求調(diào)劑有科研經(jīng)歷有二區(qū)文章已經(jīng)有18人回復(fù)
265求調(diào)劑已經(jīng)有3人回復(fù)
296求調(diào)劑已經(jīng)有6人回復(fù)
0703化學(xué)調(diào)劑已經(jīng)有13人回復(fù)
085600材料與化工求調(diào)劑已經(jīng)有15人回復(fù)
生物學(xué)調(diào)劑招人�。�！已經(jīng)有4人回復(fù)

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

這樣的動(dòng)力學(xué)方程能夠有1stopt擬合出參數(shù)來(lái)么已經(jīng)有4人回復(fù)
使用MATLAB擬合動(dòng)力學(xué)參數(shù)問(wèn)題已經(jīng)有18人回復(fù)
動(dòng)力學(xué)方程組參數(shù)求解問(wèn)題已經(jīng)有9人回復(fù)
使用matlab最優(yōu)化方法擬合獲得多個(gè)動(dòng)力學(xué)參數(shù)中的問(wèn)題已經(jīng)有4人回復(fù)
求大神幫忙matlab擬合函數(shù)求參數(shù) 已經(jīng)有9人回復(fù)
求助1stopt擬合動(dòng)力學(xué)參數(shù) 已經(jīng)有4人回復(fù)
拜求！matlab擬合動(dòng)力學(xué)參數(shù) 已經(jīng)有14人回復(fù)
求助動(dòng)力學(xué)方程？？？已經(jīng)有3人回復(fù)
求助用matlab擬合動(dòng)力學(xué)方程已經(jīng)有13人回復(fù)
動(dòng)力學(xué)擬合求助已經(jīng)有4人回復(fù)
反應(yīng)動(dòng)力學(xué)的擬合有什么意義？已經(jīng)有12人回復(fù)
通過(guò)實(shí)驗(yàn)數(shù)據(jù)擬合，求解公式中的參數(shù) 已經(jīng)有14人回復(fù)
二級(jí)動(dòng)力學(xué) 擬合已經(jīng)有9人回復(fù)
怎樣用一組參數(shù)同時(shí)擬合兩個(gè)曲線--matlab 已經(jīng)有5人回復(fù)
動(dòng)力學(xué)方程參數(shù)估計(jì)方法已經(jīng)有14人回復(fù)
動(dòng)力學(xué)參數(shù)擬合已經(jīng)有26人回復(fù)
matlab 擬合反應(yīng)動(dòng)力學(xué)參數(shù)結(jié)果很差。大家?guī)兔匆幌?/a> 已經(jīng)有14人回復(fù)

為啥我的一級(jí)動(dòng)力學(xué)方程擬合出來(lái)的數(shù)據(jù)R是負(fù)數(shù)呢？已經(jīng)有16人回復(fù)
matlab擬合擬均相動(dòng)力學(xué)參數(shù) 已經(jīng)有6人回復(fù)
【求助】用matlab最優(yōu)化方法進(jìn)行參數(shù)擬合已經(jīng)有17人回復(fù)
【求助】擬合動(dòng)力學(xué)方程求助已經(jīng)有13人回復(fù)

1樓 2015-03-04 16:31:29

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

信彼南山

木蟲(chóng) (著名寫(xiě)手)

應(yīng)助: 33 (小學(xué)生)
金幣: 4142.9
散金: 1221
紅花: 16
帖子: 1178
在線: 233.5小時(shí)
蟲(chóng)號(hào): 1133529
注冊(cè): 2010-10-27
專業(yè): 導(dǎo)航、制導(dǎo)與傳感技術(shù)

exp后面是不要*的
其余的太長(zhǎng)，沒(méi)看完

贊一下(1人)

回復(fù)此樓

2樓2015-03-05 01:28:52

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

清輝12365

新蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 23
帖子: 21
在線: 53.2小時(shí)
蟲(chóng)號(hào): 2804052
注冊(cè): 2013-11-15
專業(yè): 化學(xué)反應(yīng)工程

引用回帖:

2樓: Originally posted by 信彼南山 at 2015-03-05 01:28:52
exp后面是不要*的
其余的太長(zhǎng)，沒(méi)看完

我也看到了，請(qǐng)問(wèn)一下，還有什么錯(cuò)誤嗎？？麻煩了~~

贊一下

回復(fù)此樓

3樓2015-03-05 08:35:46

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

清輝12365

新蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 23
帖子: 21
在線: 53.2小時(shí)
蟲(chóng)號(hào): 2804052
注冊(cè): 2013-11-15
專業(yè): 化學(xué)反應(yīng)工程

引用回帖:

2樓: Originally posted by 信彼南山 at 2015-03-05 01:28:52
exp后面是不要*的
其余的太長(zhǎng)，沒(méi)看完

function YYHJIA
clear all;clc
global T P N0 x1 x2 x3;
t=[140 140 140.2 143.4 150 150.1 150.2 151.2 153.7 159.2 160 161.3 140.2 141.1 149.7 150.3 150.8 159.7 161.2 138 150.3 160 160.2 162.3];
%溫度
T=t+273.15;
P=[0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.3 0.3 0.3 0.3 0.3];
%壓力
N0=[0.3091 0.3260 0.2796 0.2856 0.2801 0.3260 0.3091 0.2853 0.3190 0.2853 0.3962 0.3190 0.3085 0.2793 0.2804 0.2856 0.3083 0.3085 0.2788 0.2368 0.26 0.2597 0.3207 0.2344];
%原料初始摩爾
x1=[0.624 0.5 0.5310 0.5710 0.53 0.5 0.624 0.572 0.511 0.572 0.524 0.511 0.625 0.531 0.529 0.571 0.625 0.625 0.532 0.501 0.57 0.571 0.555 0.506];
%原料甲醇占的摩爾分?jǐn)?shù)
x2=[0.25 0.334 0.349 0.287 0.348 0.334 0.25 0.286 0.345 0.286 0.352 0.345 0.25 0.349 0.348 0.286 0.25 0.25 0.35 0.334 0.287 0.288 0.334 0.337];
%原料CO摩爾分?jǐn)?shù)
x3=[0.127 0.166 0.121 0.142 0.122 0.166 0.127 0.142 0.144 0.142 0.124 0.144 0.125 0.12 0.123 0.143 0.124 0.125 0.118 0.165 0.142 0.141 0.111 0.157];
%原料O2的摩爾分?jǐn)?shù)
Y1=[0.015 0.016 0.015 0.017 0.02 0.021 0.021 0.021 0.021 0.025 0.023 0.025 0.018 0.018 0.025 0.026 0.026 0.03 0.03 0.023 0.028 0.036 0.034 0.038];
%出口DMC的摩爾數(shù)
Y2=[0.018 0.02 0.018 0.033 0.04 0.041 0.042 0.042 0.045 0.054 0.053 0.054 0.031 0.032 0.05 0.051 0.05 0.059 0.058 0.043 0.051 0.063 0.061 0.066];
%出口CO的摩爾數(shù)
k0=[8e4 2e4 -0.5 0.5 1 4e5 3e4 0.5 0.5 ];%(指前因子、活化能、級(jí)數(shù))初始量
lb=[0 0 -3 -3 -3 0 0 -3 -3 ];%下限
ub=[+inf +inf +inf +inf +inf +inf +inf +inf +inf];%上限
%使用lsqnonlin函數(shù)進(jìn)行參數(shù)估計(jì)
[k,resnorm,residual,exitflag,output,lambda,jacobian]=lsqnonlin(@ObjFunc,k0,lb,ub,[],T,P,N0,x1,x2,x3,Y1,Y2);
ci=nlparci(k,residual,jacobian);
function f=ObjFunc(k,T,P,N0,x1,x2,x3,Y1,Y2)
wspan=[0 0.375];
Y0=[Y1' Y2'];
X=[N0' P' T' x1' x2' x3' Y1' Y2']
[m,n]=size(Y0);
Ycal = zeros(m,n);
for i=1:m
[w yy]=ode45(@Equations,wspan,Y0(i,

,[],k);
Ycal(i,

=yy(end,

;
end
f=Ycal-Y0;
function dydw=Equations(w,yy,k,N0,P,T,x1,x2,x3)
global T P N0 x1 x2 x3;
%如何表達(dá)呢？這么多數(shù)字
r1=k(1)*exp(-k(2)./T).*P.^(k(3)+k(4)+k(5))*...
(yy(1).*(1.5*x1-2)+0.5*yy(2).*x1).^(k(3))*...
(yy(1).*(1.5*x2-1)+yy(2).*(0.5*x2-1)+x2).^(k(4))*...
(yy(1).*(1.5*x3-0.5)+yy(2).*(0.3*x3-0.5)+x3).^k(5);
r2=k(6)*exp(-k(7)./T).*P.^(k(8)+k(9))*...
(yy(1).*(1.5*x2-1)+yy(2).*(0.5*x2-1)+x2).^k(8)*...
(yy(1).*(1.5*x3-0.5)+yy(2).*(0.3*x3-0.5)+x3).^k(9);
dy1dw=(1+1.5*yy(1)+0.5*yy(2)).*(r1.*(1.5*yy(1)+1)+0.5*yy(1).*r2)./N0;
dy2dw=(1+1.5*yy(1)+0.5*yy(2)).*(r2.*(1+0.5*yy(2))+1.5*r1.*yy(2))./N0;
dydw=[dy1dw;dy2dw];
圖像為出現(xiàn)的錯(cuò)誤結(jié)果

圖1.jpg

贊一下

回復(fù)此樓

4樓2015-03-05 09:36:48

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

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

返回列表

普通表情龍兔虎貓

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

[考研] 085600材料與化工求調(diào)劑 +14	enenenhui 2026-03-13	15/750	2026-03-19 21:53 by allen-yin
[考研] 一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑 +3	葵梓衛(wèi)隊(duì) 2026-03-18	5/250	2026-03-19 19:35 by 給你你注意休息
[考博] 東華理工大學(xué)化材專業(yè)26屆碩士博士申請(qǐng) +8	zlingli 2026-03-13	8/400	2026-03-19 16:32 by 輕松不少隨
[考研] 0703化學(xué)調(diào)劑 +4	18889395102 2026-03-18	4/200	2026-03-19 16:13 by 30660438
[考研] 271材料工程求調(diào)劑 +6	.6lL 2026-03-18	6/300	2026-03-19 15:41 by haoshis
[考研] 085600材料與化工求調(diào)劑 +6	緒幸與子 2026-03-17	6/300	2026-03-19 13:27 by houyaoxu
[考研] 材料考研調(diào)劑 +3	xwt。 2026-03-19	3/150	2026-03-19 11:22 by w沐陽(yáng)w
[考研] 328求調(diào)劑，英語(yǔ)六級(jí)551，有科研經(jīng)歷 +4	生物工程調(diào)劑 2026-03-16	12/600	2026-03-19 11:10 by 生物工程調(diào)劑
[考研] 材料專碩英一數(shù)二306 +5	z1z2z3879 2026-03-18	5/250	2026-03-19 07:43 by BruceLiu320
[考研] 08工科 320總分求調(diào)劑 +5	梨花珞晚風(fēng) 2026-03-17	5/250	2026-03-18 14:49 by haxia
[考研] 297求調(diào)劑 +8	戲精丹丹丹 2026-03-17	8/400	2026-03-18 14:30 by laoshidan
[考研] 工科材料085601 279求調(diào)劑 +6	困于星晨 2026-03-17	6/300	2026-03-18 10:21 by kkcoco25
[考博] 26申博 +4	八6八68 2026-03-16	4/200	2026-03-17 13:00 by 輕松不少隨
[考研] 機(jī)械專碩325，尋找調(diào)劑院校 +3	y9999 2026-03-15	5/250	2026-03-16 19:58 by y9999
[基金申請(qǐng)] 今年的國(guó)基金是打分制嗎？ 50+3	zhanghaozhu 2026-03-14	3/150	2026-03-16 17:07 by 北京萊茵潤(rùn)色
[考研] 304求調(diào)劑 +3	曼殊2266 2026-03-14	3/150	2026-03-16 16:39 by houyaoxu
[考研] 0856專碩279求調(diào)劑 +5	加油加油！? 2026-03-15	5/250	2026-03-15 11:58 by 2020015
[考研] 22408總分284求調(diào)劑 +3	InAspic 2026-03-13	3/150	2026-03-15 11:10 by zhq0425
[考研] 本科南京大學(xué)一志愿川大藥學(xué)327 +3	麥田耕者 2026-03-14	3/150	2026-03-14 20:04 by 外星文明
[基金申請(qǐng)] 現(xiàn)在如何回避去年的某一個(gè)專家，不知道名字 +3	zk200107 2026-03-12	6/300	2026-03-14 17:13 by zk200107

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

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

清輝12365

[求助] 動(dòng)力學(xué)參數(shù)擬合

» 猜你喜歡

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

信彼南山

清輝12365

清輝12365

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