版塊導(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)版塊排行榜

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

返回列表

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

清輝12365

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

應(yīng)助: 0 (幼兒園)
金幣: 23
帖子: 21
在線: 53.2小時(shí)
蟲(chóng)號(hào): 2804052
注冊(cè): 2013-11-15
專(zhuān)業(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ù)此樓

» 猜你喜歡

本科新能源科學(xué)與工程，一志愿華理能動(dòng)285求調(diào)劑已經(jīng)有5人回復(fù)
一志愿華北電力大學(xué)能動(dòng)專(zhuān)碩，293，求調(diào)劑已經(jīng)有5人回復(fù)
調(diào)劑310 已經(jīng)有7人回復(fù)
食品工程專(zhuān)碩一志愿中海洋309求調(diào)劑已經(jīng)有6人回復(fù)
275求調(diào)劑已經(jīng)有14人回復(fù)
283求調(diào)劑已經(jīng)有3人回復(fù)
085602 化工專(zhuān)碩 338分求調(diào)劑已經(jīng)有12人回復(fù)
312，生物學(xué)求調(diào)劑已經(jīng)有3人回復(fù)
化學(xué)0703 調(diào)劑 306分一志愿211 已經(jīng)有4人回復(fù)
311（085601）求調(diào)劑已經(jīng)有3人回復(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的回帖

清輝12365

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

應(yīng)助: 0 (幼兒園)
金幣: 23
帖子: 21
在線: 53.2小時(shí)
蟲(chóng)號(hào): 2804052
注冊(cè): 2013-11-15
專(zhuān)業(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的回帖

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

信彼南山

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

應(yīng)助: 33 (小學(xué)生)
金幣: 4142.9
散金: 1221
紅花: 16
帖子: 1178
在線: 233.5小時(shí)
蟲(chóng)號(hào): 1133529
注冊(cè): 2010-10-27
專(zhuān)業(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
專(zhuān)業(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的回帖

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

普通表情龍兔虎貓

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

[考研] 調(diào)劑310 +6	溫柔的晚安 2026-03-25	7/350	2026-03-28 16:04 by qq2000
[考研] 復(fù)試調(diào)劑 +3	raojunqi0129 2026-03-28	3/150	2026-03-28 15:27 by 落�？伤�
[考研] 一志愿廈門(mén)大學(xué)化學(xué)學(xué)碩307求調(diào)劑 +10	y7czhao 2026-03-26	10/500	2026-03-28 14:23 by 唐沐兒
[考研] 本科新能源科學(xué)與工程，一志愿華理能動(dòng)285求調(diào)劑 +3	AZMK 2026-03-28	5/250	2026-03-28 13:55 by AZMK
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +8	邱gl 2026-03-27	8/400	2026-03-28 12:42 by 唐沐兒
[考研] 求調(diào)劑推薦材料 304 +15	荷包蛋hyj 2026-03-26	15/750	2026-03-28 04:13 by fmesaito
[考研] 291求調(diào)劑 +7	孅華 2026-03-22	7/350	2026-03-28 04:02 by fmesaito
[考研] 285求調(diào)劑 +4	AZMK 2026-03-27	7/350	2026-03-27 20:59 by AZMK
[考研] 一志愿西北大學(xué) 總分282 英語(yǔ)一62 求調(diào)劑 +7	18419759900 2026-03-25	8/400	2026-03-27 16:38 by 18419759900
[考研] 07化學(xué)280分求調(diào)劑 +10	722865 2026-03-23	10/500	2026-03-27 15:51 by Plutoqq
[考研] 安徽大學(xué)專(zhuān)碩生物與醫(yī)藥專(zhuān)業(yè)(086000)324分，英語(yǔ)已過(guò)四六級(jí)，六級(jí)521，求調(diào)劑 +4	美味可樂(lè)雞翅 2026-03-26	4/200	2026-03-27 15:27 by 星空星月
[考研] 279 分求調(diào)劑 +4	睡個(gè)好覺(jué)_16 2026-03-24	4/200	2026-03-27 15:05 by 醉在風(fēng)里
[考研] 考研調(diào)劑 +10	呼呼？~+123456 2026-03-24	10/500	2026-03-27 11:46 by wangjy2002
[考研] 材料學(xué)碩，求調(diào)劑 6+5	糖葫蘆888ll 2026-03-22	10/500	2026-03-27 08:18 by hypershenger
[考研] 284求調(diào)劑 +11	junqihahaha 2026-03-26	12/600	2026-03-27 04:37 by wxiongid
[考研] 349求調(diào)劑 +4	李木子啊哈哈 2026-03-25	4/200	2026-03-26 22:49 by fmesaito
[考研] 340求調(diào)劑 +5	話(huà)梅糖111 2026-03-24	5/250	2026-03-25 06:53 by ilovexiaobin
[考研] 一志愿吉大化學(xué)322求調(diào)劑 +4	17501029541 2026-03-23	6/300	2026-03-24 10:21 by 戴圍脖的小蚊子
[考研] 336化工調(diào)劑 +4	王大坦1 2026-03-23	5/250	2026-03-23 18:32 by allen-yin
[考研] 275求調(diào)劑 +6	shansx 2026-03-22	8/400	2026-03-22 15:27 by barlinike