 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 907 | 回復(fù): 3 | ||
清輝12365新蟲 (初入文壇)
|
[求助]
動力學(xué)參數(shù)擬合
|
|
原料是：甲醇、CO、氧氣。主反應(yīng)：甲醇、CO、氧氣合成DMC，副反應(yīng)：CO和氧氣合成CO2. 要干的活：擬合出生成DMC和CO2的動力學(xué)參數(shù)（指前因子、活化能、級數(shù)） 動力學(xué)方程為dy1/dw=f(y1、y2)和dy2/dw=f(y1、y2),具體形式見附圖。 整理后的實驗數(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)鍵難度在：主要是兩個常微分方程、數(shù)據(jù)多不知如何進(jìn)行運算。下面是自己編的程序，出現(xiàn)問題并且解決不上來，希望大神指教 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 ];%(指前因子、活化能、級數(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ù)估計 [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)的問題如圖： 出現(xiàn)的問題.jpg 方程形式.jpg 方程具體形式2.jpg |
回復(fù)此樓» 猜你喜歡
288求調(diào)劑，一志愿華南理工大學(xué)071005
已經(jīng)有5人回復(fù)
-
一志愿中國海洋大學(xué)，生物學(xué)，301分，求調(diào)劑
已經(jīng)有5人回復(fù)
-
294求調(diào)劑材料與化工專碩
已經(jīng)有12人回復(fù)
-
0856調(diào)劑，是學(xué)校就去
已經(jīng)有5人回復(fù)
-
0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章
已經(jīng)有13人回復(fù)
-
梁成偉老師課題組歡迎你的加入
已經(jīng)有11人回復(fù)
-
復(fù)試調(diào)劑
已經(jīng)有6人回復(fù)
-
0703化學(xué)調(diào)劑 ，六級已過，有科研經(jīng)歷
已經(jīng)有11人回復(fù)
-
321求調(diào)劑
已經(jīng)有10人回復(fù)
-
東華理工大學(xué)化材專業(yè)26屆碩士博士申請
已經(jīng)有8人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 這樣的動力學(xué)方程能夠有1stopt擬合出參數(shù)來么
已經(jīng)有4人回復(fù)
- 使用MATLAB擬合動力學(xué)參數(shù)問題
已經(jīng)有18人回復(fù)
- 動力學(xué)方程組參數(shù)求解問題
已經(jīng)有9人回復(fù)
- 使用matlab最優(yōu)化方法擬合獲得多個動力學(xué)參數(shù)中的問題
已經(jīng)有4人回復(fù)
- 求大神幫忙matlab擬合函數(shù)求參數(shù)
已經(jīng)有9人回復(fù)
- 求助1stopt擬合動力學(xué)參數(shù)
已經(jīng)有4人回復(fù)
- 拜求！matlab擬合動力學(xué)參數(shù)
已經(jīng)有14人回復(fù)
- 求助動力學(xué)方程？？？
已經(jīng)有3人回復(fù)
- 求助用matlab擬合動力學(xué)方程
已經(jīng)有13人回復(fù)
- 動力學(xué)擬合求助
已經(jīng)有4人回復(fù)
- 反應(yīng)動力學(xué)的擬合有什么意義？
已經(jīng)有12人回復(fù)
- 通過實驗數(shù)據(jù)擬合，求解公式中的參數(shù)
已經(jīng)有14人回復(fù)
- 二級動力學(xué) 擬合
已經(jīng)有9人回復(fù)
- 怎樣用一組參數(shù)同時擬合兩個曲線--matlab
已經(jīng)有5人回復(fù)
- 動力學(xué)方程參數(shù)估計方法
已經(jīng)有14人回復(fù)
- 動力學(xué)參數(shù)擬合
已經(jīng)有26人回復(fù)
- matlab 擬合反應(yīng)動力學(xué)參數(shù)結(jié)果很差。大家?guī)兔匆幌?/a>
已經(jīng)有14人回復(fù)
- 為啥我的一級動力學(xué)方程擬合出來的數(shù)據(jù)R是負(fù)數(shù)呢？
已經(jīng)有16人回復(fù)
- matlab擬合擬均相動力學(xué)參數(shù)
已經(jīng)有6人回復(fù)
- 【求助】用matlab最優(yōu)化方法進(jìn)行參數(shù)擬合
已經(jīng)有17人回復(fù)
- 【求助】擬合動力學(xué)方程求助
已經(jīng)有13人回復(fù)
信彼南山
木蟲 (著名寫手)
- 應(yīng)助: 33 (小學(xué)生)
- 金幣: 4142.9
- 散金: 1221
- 紅花: 16
- 帖子: 1178
- 在線: 233.5小時
- 蟲號: 1133529
- 注冊: 2010-10-27
- 專業(yè): 導(dǎo)航、制導(dǎo)與傳感技術(shù)
2樓2015-03-05 01:28:52
清輝12365
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 23
- 帖子: 21
- 在線: 53.2小時
- 蟲號: 2804052
- 注冊: 2013-11-15
- 專業(yè): 化學(xué)反應(yīng)工程
3樓2015-03-05 08:35:46
清輝12365
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 23
- 帖子: 21
- 在線: 53.2小時
- 蟲號: 2804052
- 注冊: 2013-11-15
- 專業(yè): 化學(xué)反應(yīng)工程
|
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 ];%(指前因子、活化能、級數(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ù)估計 [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)的錯誤結(jié)果 圖1.jpg |
4樓2015-03-05 09:36:48
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 清輝12365 的主題更新
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿中國海洋大學(xué)，生物學(xué)，301分，求調(diào)劑 +5 | 1孫悟空 2026-03-17 | 5/250 |
|
|---|---|---|---|---|
|
|
[考研] 材料專碩274一志愿陜西師范大學(xué)求調(diào)劑 +8 | 薛云鵬 2026-03-13 | 8/400 |
|
|
|
[考研] 一志愿武漢理工材料工程專碩調(diào)劑 +4 | Doleres 2026-03-19 | 4/200 |
|
|
|
[考研] 一志愿 西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。 +3 | 晨昏線與星海 2026-03-19 | 3/150 |
|
|
[考研] 267一志愿南京工業(yè)大學(xué)0817化工求調(diào)劑 +10 | SUICHILD 2026-03-12 | 10/500 |
|
|
|
[考研] 085700資源與環(huán)境308求調(diào)劑 +3 | 墨墨漠 2026-03-18 | 3/150 |
|
|
|
[考研] 085410人工智能專碩317求調(diào)劑（0854都可以） +3 | xbxudjdn 2026-03-18 | 3/150 |
|
|
|
[考研]
 0854，計算機(jī)類招收調(diào)劑
+3
|
胡辣湯放糖 2026-03-15 | 6/300 |
|
|
|
[考研] 材料，紡織，生物（0856、0710），化學(xué)招生啦 +3 | Eember. 2026-03-17 | 9/450 |
|
|
|
[考研] 308求調(diào)劑 +4 | 是Lupa啊 2026-03-16 | 4/200 |
|
|
|
[考研] 材料與化工專碩調(diào)劑 +5 | heming3743 2026-03-16 | 5/250 |
|
|
|
[考研] 211本，11408一志愿中科院277分，曾在中科院自動化所實習(xí) +6 | Losir 2026-03-12 | 7/350 |
|
|
|
[考研] 283求調(diào)劑 +3 | 聽風(fēng)就是雨； 2026-03-16 | 3/150 |
|
|
|
[考研] 070300化學(xué)學(xué)碩求調(diào)劑 +6 | 太想進(jìn)步了0608 2026-03-16 | 6/300 |
|
|
|
[考研] 294求調(diào)劑 +3 | Zys010410@ 2026-03-13 | 4/200 |
|
|
|
[考研] 297一志愿上交085600求調(diào)劑 +5 | 指尖八千里 2026-03-14 | 5/250 |
|
|
|
[考研] 復(fù)試調(diào)劑 +3 | 呼呼？~+123456 2026-03-14 | 3/150 |
|
|
|
[考研] 297求調(diào)劑 +4 | 學(xué)海漂泊 2026-03-13 | 4/200 |
|
|
|
[考研] 321求調(diào)劑（食品/專碩） +3 | mxcz321 2026-03-12 | 6/300 |
|
|
|
[考研] 333求調(diào)劑 +3 | 152697 2026-03-12 | 4/200 |
|