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

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

返回列表

當(dāng)前只顯示滿足指定條件的回帖，點(diǎ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ù)此樓

» 猜你喜歡

071000生物學(xué)求調(diào)劑，初試成績(jī)343 已經(jīng)有7人回復(fù)
0703化學(xué)調(diào)劑，求導(dǎo)師收已經(jīng)有8人回復(fù)
0856材料化工調(diào)劑總分330 已經(jīng)有11人回復(fù)
本科新能源科學(xué)與工程，一志愿華理能動(dòng)285求調(diào)劑已經(jīng)有9人回復(fù)
332求調(diào)劑已經(jīng)有4人回復(fù)
085600，專業(yè)課化工原理，321分求調(diào)劑已經(jīng)有3人回復(fù)
材料工程294求調(diào)劑已經(jīng)有3人回復(fù)
295材料工程專碩求調(diào)劑已經(jīng)有7人回復(fù)
305求調(diào)劑已經(jīng)有5人回復(fù)
食品工程專碩一志愿中海洋309求調(diào)劑已經(jīng)有8人回復(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
專業(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è)回答

信彼南山

木蟲(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的回帖

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

普通表情龍兔虎貓

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

[考研] 085600，專業(yè)課化工原理，321分求調(diào)劑 +3	大饞小子 2026-03-28	3/150	2026-03-28 20:10 by wxiongid
[考研] 材料與化工272求調(diào)劑 +9	阿斯蒂芬2004 2026-03-28	9/450	2026-03-28 15:21 by 1018329917
[考研] 320分，材料與化工專業(yè)，求調(diào)劑 +9	一定上岸aaa 2026-03-27	13/650	2026-03-28 15:00 by 神馬都不懂
[考研] 求調(diào)劑 +3	QiMing7 2026-03-25	4/200	2026-03-28 14:30 by QiMing7
[考研] 394求調(diào)劑 +3	好事多磨靜候佳?/a> 2026-03-26	5/250	2026-03-28 14:24 by 唐沐兒
[考研] 317求調(diào)劑 +6	十閑wx 2026-03-24	6/300	2026-03-28 13:27 by Iveryant
[考研] 0703本科鄭州大學(xué)求調(diào)劑 +3	nhj_ 2026-03-25	3/150	2026-03-28 13:24 by Iveryant
[考研] 0703化學(xué)/290求調(diào)劑/本科經(jīng)歷豐富/工科也可 +9	丹青奶蓋 2026-03-26	10/500	2026-03-28 07:45 by barnett0632
[考研] 一志愿上海理工能源動(dòng)力（085800）310分求調(diào)劑 +3	zhangmingc 2026-03-27	4/200	2026-03-27 19:01 by 給你你注意休息
[考研] 一志愿西北大學(xué) 總分282 英語(yǔ)一62 求調(diào)劑 +7	18419759900 2026-03-25	8/400	2026-03-27 16:38 by 18419759900
[考研] 材料與化工085600，總分304，本科有兩篇sci參與，求調(diào)劑 +10	幸運(yùn)的醬醬 2026-03-22	12/600	2026-03-27 16:08 by muchong357
[考研] 305求調(diào)劑 +5	哇盧卡庫(kù) 2026-03-26	5/250	2026-03-27 14:01 by laoshidan
[考研] 一志愿華東理工大學(xué)081700，初試分?jǐn)?shù)271 +6	kotoko_ik 2026-03-23	7/350	2026-03-27 12:29 by 惠州彭于晏
[考研] 315調(diào)劑 +4	0860求調(diào)劑 2026-03-26	5/250	2026-03-27 11:23 by wangjy2002
[考研] 0703化學(xué)338求調(diào)劑！ +6	Zuhui0306 2026-03-26	7/350	2026-03-27 10:35 by shangxh
[考研] 342求調(diào)劑 +3	加油a李zs 2026-03-26	3/150	2026-03-27 00:29 by wxiongid
[考研] 求b區(qū)院校調(diào)劑 +4	周56 2026-03-24	5/250	2026-03-25 17:12 by yishunmin
[考研] 277分求調(diào)劑，跨調(diào)材料 +3	考研調(diào)劑lxh 2026-03-24	3/150	2026-03-24 13:52 by JourneyLucky
[考研] 085404電子信息284分求調(diào)劑 +4	13659058978 2026-03-24	4/200	2026-03-24 12:15 by syl20081243
[考研] 344求調(diào)劑 +3	desto 2026-03-24	3/150	2026-03-24 10:09 by 搏擊518

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

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

清輝12365

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

» 猜你喜歡

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

清輝12365

信彼南山

清輝12365

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

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