 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 1 | 1/1 | 返回列表 |
| 查看: 511 | 回復(fù): 0 | ||
zhaoshazhu新蟲 (小有名氣)
|
[求助]
求助Matlab
|
|
請(qǐng)高手幫我查看一下代碼有什么問(wèn)題，謝謝！下面是我估算出來(lái)的初值和殘差。但是殘差比實(shí)驗(yàn)值還要大，代碼有問(wèn)題可以我Matlab不懂，求大神指導(dǎo)。 function DMM170 clear all clc k0 = [0.01 0.006 0.01 8 0.6 9 0.8 1 0.3 0.007]; lb = [0 0 0 0 0 0 0 0 0 0]; ub = [inf inf inf inf inf inf inf inf inf inf]; P0 =[0.016664 4.642001 0.340144 0 0 0; 0.019859 4.576200 0.402387 0 0 0; 0.024201 4.481438 0.492568 0 0 0; 0.017000 4.641373 0.340098 0 0 0; 0.020132 4.575695 0.402343 0 0 0; 0.024499 4.480897 0.492509 0 0 0; 0.017128 4.641134 0.340081 0 0 0; 0.020678 4.574688 0.402254 0 0 0; 0.025811 4.478524 0.492248 0 0 0; 0.046146 4.056541 0.891733 0 0 0 ]; % 初始分壓,MPa Pi=[0.0052847 4.6395855 0.3479159 0.0037084 0.0032052 0.0003002; 0.0072011 4.5720394 0.4119683 0.0044129 0.0040198 0.0003585; 0.0089072 4.4774599 0.5037451 0.0048294 0.0046554 0.0004029; 0.0049652 4.6375496 0.3490244 0.0032731 0.0047272 0.0004605; 0.0063178 4.5709579 0.4129024 0.0039176 0.0053778 0.0005265; 0.0073769 4.4786796 0.5039073 0.0043511 0.0051445 0.0005405; 0.0040522 4.6367309 0.3498930 0.0030793 0.0057272 0.0005174; 0.0051847 4.5702655 0.4137574 0.0037354 0.0063233 0.0007336; 0.0062153 4.4770637 0.5052473 0.0043140 0.0064064 0.0007533; 0.0191512 4.0567809 0.9105013 0.0078227 0.0050752 0.0006686 ]; % 經(jīng)過(guò)Wc/F0后，各物質(zhì)分壓,MPa % 使用函數(shù)lsqnonlin()進(jìn)行參數(shù)估計(jì) opt=optimset('Algorithm','levenberg-marquardt'); [k,resnorm,residual,exitflag,output,lambda,jacobian] = lsqnonlin(@ObjFunc,k0,[],[],opt,P0,Pi); ci = nlparci(k,residual,jacobian); fprintf('\n\n使用函數(shù)lsqnonlin()估計(jì)得到的參數(shù)值為:\n') fprintf('\tk1 = %.4f ± %.4f\n',k(1),ci(1,2)-k(1)) fprintf('\tk2 = %.4f ± %.4f\n',k(2),ci(2,2)-k(2)) fprintf('\tk3 = %.4f ± %.4f\n',k(3),ci(3,2)-k(3)) fprintf('\tk4 = %.4f ± %.4f\n',k(4),ci(4,2)-k(4)) fprintf('\tk5 = %.4f ± %.4f\n',k(5),ci(5,2)-k(5)) fprintf('\tk6 = %.4f ± %.4f\n',k(6),ci(6,2)-k(6)) fprintf('\tk7 = %.4f ± %.4f\n',k(7),ci(7,2)-k(7)) fprintf('\tk8 = %.4f ± %.4f\n',k(8),ci(8,2)-k(8)) fprintf('\tk9 = %.4f ± %.4f\n',k(9),ci(9,2)-k(9)) fprintf('\tk10 = %.4f ± %.4f\n',k(10),ci(10,2)-k(10)) fprintf('\t殘差平方和 = %.4f\n',resnorm) fprintf('\texitflag = %.4f\n',exitflag) fprintf('\tresidual = %.4f\n',residual) % ------------------------------------------------------------------ function f = ObjFunc(k,P0,Pi) % 目標(biāo)函數(shù) [m,n] = size(P0); Pcal = zeros(m,n); tspan = [0 0.82; 0 0.97; 0 1.19; 0 1.02; 0 1.21; 0 1.48; 0 1.37; 0 1.62; 0 1.98; 0 2.16 ]; % 即Wc/F0,g.h/mol for i = 1:m [t PP] = ode45(@Euqations,tspan(i,  ,P0(i,,[],k);Pcal(i, = PP(end,;end f= Pcal-Pi; % ------------------------------------------------------------------ function dPdt = Euqations(t, P, k) % here t = Wc / F0 denom = 1+k(4)*P(1)+k(5)*P(3)+k(6)*P(4)+k(7)*P(5)+k(8)*P(6); % k(4) = KDMM,k(5) = KME ,k(6)=KHPM,k(7)=KPDO,k(8)=KNPA,k(9)=Kp1,k(10)=Kp2 theA =k(4)*(P(1)*P(2)-P(4)*P(3)/k(9)*P(2)) / denom^2; theB =k(6)*( P(4)*P(2)-P(5)*P(3)/k(10)*P(2))/ denom^2; theC =k(7)*P(5)*P(2)/denom^2; r1 = k(1)*theA; r2 = k(2)*theB; r3 = k(3)*theC; dPDMMdt = -r1; dPHdt = -2*r1-2*r2-r3; dPMEdt = r1+r2; dPHPMdt = r1-r2; dPPDOdt = r2-r3; dPNPAdt = r3; dPdt = [dPDMMdt;dPHdt;dPMEdt;dPHPMdt;dPPDOdt;dPNPAdt]; 使用函數(shù)lsqnonlin()估計(jì)得到的參數(shù)值為: k1 = 0.0157 ± 5.6278 k2 = 0.0058 ± 5.8687 k3 = 0.0098 ± 462.8444 k4 = 8.0000 ± 2964.0969 k5 = 0.6002 ± 129.8991 k6 = 9.0000 ± 8644.0718 k7 = 0.8000 ± 33677.3330 k8 = 1.0000 ± 365138.2060 k9 = 0.2990 ± 15.8150 k10 = 0.0007 ± 0.2817 殘差平方和 = 0.0026 exitflag = 4.0000 residual = 0.0080 residual = 0.0086 residual = 0.0103 residual = 0.0080 residual = 0.0090 residual = 0.0113 residual = 0.0082 residual = 0.0096 residual = 0.0125 residual = 0.0189 residual = -0.0044 residual = -0.0041 residual = -0.0061 residual = -0.0042 residual = -0.0049 residual = -0.0094 residual = -0.0054 residual = -0.0073 residual = -0.0127 residual = -0.0164 residual = -0.0044 residual = -0.0055 residual = -0.0062 residual = -0.0049 residual = -0.0057 residual = -0.0056 residual = -0.0049 residual = -0.0056 residual = -0.0059 residual = -0.0107 residual = -0.0003 residual = -0.0003 residual = 0.0002 residual = 0.0007 residual = 0.0009 residual = 0.0014 residual = 0.0018 residual = 0.0021 residual = 0.0027 residual = 0.0002 residual = -0.0032 residual = -0.0040 residual = -0.0046 residual = -0.0047 residual = -0.0054 residual = -0.0051 residual = -0.0057 residual = -0.0063 residual = -0.0064 residual = -0.0051 residual = -0.0003 residual = -0.0004 residual = -0.0004 residual = -0.0005 residual = -0.0005 residual = -0.0005 residual = -0.0005 residual = -0.0007 residual = -0.0008 residual = -0.0007 >> |
回復(fù)此樓» 猜你喜歡
調(diào)劑
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復(fù)
-
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有237人回復(fù)
-
中德博士后2026項(xiàng)目5月30日截止- 國(guó)內(nèi)博士后入站、德國(guó)兩年全職
已經(jīng)有0人回復(fù)
-
面向全球申請(qǐng)人的中德博士后2026項(xiàng)目5月30日截止- 國(guó)內(nèi)博后入站德國(guó)兩年全職
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有2人回復(fù)
-
重慶交大材料學(xué)院26年碩士招生復(fù)試實(shí)施細(xì)則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。
已經(jīng)有0人回復(fù)
-
北工大宋曉艷教授團(tuán)隊(duì)招收2026年9月入學(xué)博士【學(xué)博和專博名額均有】！
已經(jīng)有1人回復(fù)
-
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收10人
已經(jīng)有0人回復(fù)
-
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收調(diào)劑，10人
已經(jīng)有0人回復(fù)
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 zhaoshazhu 的主題更新
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 22408 359分調(diào)劑 +4 | Qshers 2026-03-27 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 一志愿北化085600材料專碩275|有文章專利｜求調(diào)劑 +7 | Micky11223 2026-03-25 | 7/350 |
|
|
|
[考研] 320分，材料與化工專業(yè)，求調(diào)劑 +9 | 一定上岸aaa 2026-03-27 | 13/650 |
|
|
|
[考研] 0703化學(xué) +10 | 妮妮ninicgb 2026-03-27 | 10/500 |
|
|
|
[考研] 299求調(diào)劑 +7 | 嗯嗯嗯嗯2 2026-03-27 | 7/350 |
|
|
|
[考研] 311求調(diào)劑 +9 | lin0039 2026-03-26 | 9/450 |
|
|
|
[考研] 調(diào)劑 +3 | 好好讀書。 2026-03-28 | 3/150 |
|
|
|
[考研] 張芳銘-中國(guó)農(nóng)業(yè)大學(xué)-環(huán)境工程專碩-298 +4 | 手機(jī)用戶 2026-03-26 | 4/200 |
|
|
|
[考研] 材料與化工085600，總分304，本科有兩篇sci參與，求調(diào)劑 +10 | 幸運(yùn)的醬醬 2026-03-22 | 12/600 |
|
|
|
[考研] 348求調(diào)劑 +4 | 小懶蟲不懶了 2026-03-27 | 5/250 |
|
|
|
[考研] 317求調(diào)劑 +7 | 蛋黃咸肉粽 2026-03-26 | 7/350 |
|
|
|
[考研] 327求調(diào)劑 +7 | prayer13 2026-03-23 | 7/350 |
|
|
|
[考研]
材料調(diào)劑
 5+4
|
想要一壺桃花水 2026-03-25 | 10/500 |
|
|
|
[考研] 求調(diào)劑 +8 | Auroracx 2026-03-22 | 8/400 |
|
|
|
[考研] 生物學(xué) 296 求調(diào)劑 +4 | 朵朵- 2026-03-26 | 6/300 |
|
|
|
[考研]
 085602 289分求調(diào)劑
+8
|
WWW西西弗斯 2026-03-24 | 8/400 |
|
|
|
[考研] 263求調(diào)劑 +6 | yqdszhdap－ 2026-03-22 | 10/500 |
|
|
|
[考研] 上海電力大學(xué)材料防護(hù)與新材料重點(diǎn)實(shí)驗(yàn)室招收調(diào)劑研究生（材料、化學(xué)、電化學(xué)，環(huán)境） +4 | 我愛學(xué)電池 2026-03-23 | 4/200 |
|
|
|
[考研] 一志愿北化315 求調(diào)劑 +3 | akrrain 2026-03-24 | 3/150 |
|
|
|
[考研] 一志愿河北工業(yè)大學(xué)0817化工278分求調(diào)劑 +7 | jhybd 2026-03-23 | 12/600 |
|