版塊導航: 正在加載中...

登錄注冊

應《網(wǎng)絡安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網(wǎng)跟帖服務。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

zhaoshazhu

新蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 125.7
散金: 10
帖子: 51
在線: 38小時
蟲號: 1812564
注冊: 2012-05-12
專業(yè): 環(huán)境化工

[求助] 求Matlab高手指導已有1人參與

Matlab運算結(jié)果的置信區(qū)間很大，和什么有關(guān)系呢？
function KineticsEst6
clear all
clc
tspan = [0 662.25];
k0 = [0.4587 0.4971  10  12 6 8 205 0.5653];
lb = [0  0  0  0  0 0 0];
ub = [50  50  100  100 50 50 1000 50];

P0 =[0.04020 4.02010 0.18794 0 0;
   0.02228 4.45682 0.10418 0 0;
   0.01541 4.62428 0.07206 0 0;
   0.01336 4.67446 0.06244 0 0;
   0.04020 4.02010 0.18794 0 0;
   0.04020 4.02010 0.18794 0 0;
   0.04020 4.02010 0.18794 0 0
];  % 初始分壓,MPa

Pi=[0.01303 4.00374 0.96404 0.007095 0.012097;
0.00393 4.44636 0.53626 0.074475 0.005998;
0.00127 4.61707 0.37150 0.070931 0.003074;
0.00090 4.66795 0.32206 0.066794 0.002414;
0.00243 4.00519 0.96961 0.158768 0.006883;
0.00681 4.00390 0.96766 0.117792 0.009855;
0.01303 4.00374 0.96404 0.070945 0.012097
];
% 經(jīng)過Wc/F0后，各物質(zhì)分壓,MPa

% 使用函數(shù)lsqnonlin()進行參數(shù)估計
options=optimset('MaxFunEvals',1000000,'MaxIter',400000)
[k,resnorm,residual,exitflag,output,lambda,jacobian] = lsqnonlin(@ObjFunc,k0,lb,ub,[],P0,Pi);
ci = nlparci(k,residual,jacobian);
fprintf('\n\n使用函數(shù)lsqnonlin()估計得到的參數(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))
% ------------------------------------------------------------------
function f = ObjFunc(k,P0,Pi)          % 目標函數(shù)
[m,n] = size(P0);
Pcal = zeros(m,n);
tspan =[0 662.25];       % 即Wc/F0,g.h/mol
for i = 1:m
[t PP] = ode45(@Euqations,tspan,P0(i,

,[],k);
Pcal(i,

= PP(end,

;
end
f= Pcal-Pi;

% ------------------------------------------------------------------
function dPdt = Euqations(t, P, k) % here t = Wc / F0
denom = 1+k(3)*P(1)+k(5)*P(4)+k(5)*P(3)+k(6)*P(5); % k(3) = KDMM, k(4) = KME ,k(5)=KHPM,k(6)=KPDO,k(7)=Kp1,k(8)=Kp2
theA =k(3)*P(1)*P(2)*(1-P(4)*P(3)/k(7)*P(1)*P(2)^2) / denom;
theB =k(5)* P(4)*P(2)*(1-P(5)*P(3)/k(8)*P(4)*P(2)^2)/ denom;
r1 = k(1)*theA;
r2 = k(2)*theB;

dPDMMdt = -r1;
dPHdt = -2*r1-2*r2;
dPMEdt = r1+r2;
dPHPMdt = r1-r2;
dPPDOdt = r2;

dPdt = [dPDMMdt;dPHdt;dPMEdt;dPHPMdt;dPPDOdt];
使用函數(shù)lsqnonlin()估計得到的參數(shù)值為:
k1 = 0.4599 ± 89509.6910
k2 = 0.4971 ± 101889.0949
k3 = 10.0001 ± 1946900.8019
k4 = 12.0000 ± 8639183.4320
k5 = 5.9997 ± 420897.9438
k6 = 7.9999 ± 2837488.5903
k7 = 205.0000 ± 1472315356.6507
k8 = 0.5653 ± 3696268.6136
結(jié)果和我設的初值一樣，是不是就沒有計算程序。

回復此樓

» 猜你喜歡

光學工程學碩調(diào)劑信息已經(jīng)有32人回復
歡迎加入課題組已經(jīng)有0人回復
物理學I論文潤色/翻譯怎么收費? 已經(jīng)有214人回復
散金幣，求好運，祝面上順利！已經(jīng)有78人回復
華東師范大學芯片設計徐瓏真實水平咋樣已經(jīng)有1人回復
調(diào)劑已經(jīng)有0人回復
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有0人回復
中德博士后2026項目5月30日截止- 國內(nèi)博士后入站、德國兩年全職已經(jīng)有0人回復
面向全球申請人的中德博士后2026項目5月30日截止- 國內(nèi)博后入站德國兩年全職已經(jīng)有0人回復
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有2人回復
重慶交大材料學院26年碩士招生復試實施細則已出。歡迎加入光子學微結(jié)構(gòu)與器件課題組。已經(jīng)有0人回復

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

求大師指點（焦急等待指點）已經(jīng)有5人回復
origin8.0 非線性自定義函數(shù)擬合，求指導已經(jīng)有5人回復
跪求高手幫忙看看這組數(shù)據(jù)如何進行擬合已經(jīng)有7人回復
ELK簡單計算失敗，求高手解答！已經(jīng)有4人回復
matlab里擬合后，知道y坐標求x坐標已經(jīng)有4人回復
求在matlab中進行分類的程序已經(jīng)有5人回復
求大神用matlab給我處理一組數(shù)據(jù)，有人沒有��？已經(jīng)有15人回復
求Origin或Matlab高手非線性擬合公式已經(jīng)有5人回復
求助高手�。。。。。。。。。�！matlab最小約束的多變量尋優(yōu)函數(shù)求解方程已經(jīng)有5人回復
尋求matlab程序高手已經(jīng)有3人回復
求高手解答Matlab求解微分方程的初始值和邊界值問題已經(jīng)有5人回復
matlab 求解非線性方程解出虛數(shù)解已經(jīng)有7人回復
學習matlab算法編程，推薦一些有用的資料或者書籍吧已經(jīng)有21人回復
求Matlab高手解決線性方程組的迭代求解問題已經(jīng)有17人回復
matlab回歸三參數(shù) 求高手解救已經(jīng)有7人回復
【請教】 2D Gaussian Fitting based on Matlab 已經(jīng)有4人回復
求SWMM軟件高手指導已經(jīng)有15人回復
matlab 氣泡邊緣提取求助高手��！已經(jīng)有5人回復
matlab中dydt的編寫要求已經(jīng)有12人回復

1樓 2014-10-24 10:01:14

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

回帖支持 ( 顯示支持度最高的前 50 名 )

月只藍

主管區(qū)長 (職業(yè)作家)

專家經(jīng)驗: +1059
計算強帖: 8
應助: 1712 (講師)
貴賓: 8.888
金幣: 68127.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時
蟲號: 1122189
注冊: 2010-10-14
專業(yè): 宇宙學
管轄: 計算模擬區(qū)

【答案】應助回帖

感謝參與，應助指數(shù) +1

這算是很普遍的問題了，初值給得不好。建議添加相關(guān)系數(shù)或者決定系數(shù)定量地考察擬合效果，只看輸出的參數(shù)數(shù)值沒什么意義。

贊一下(1人)

回復此樓

MATLAB、MS小問題、普通問題請發(fā)帖求助！時間精力有限，恕不接受無償私信求助。

2樓2014-10-24 16:25:38

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

普通回帖

月只藍

主管區(qū)長 (職業(yè)作家)

專家經(jīng)驗: +1059
計算強帖: 8
應助: 1712 (講師)
貴賓: 8.888
金幣: 68127.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時
蟲號: 1122189
注冊: 2010-10-14
專業(yè): 宇宙學
管轄: 計算模擬區(qū)

【答案】應助回帖

樓上說初值給得不好是因為 MATLAB 的lsqnonlin 對初值的依賴性很大

贊一下

回復此樓

MATLAB、MS小問題、普通問題請發(fā)帖求助！時間精力有限，恕不接受無償私信求助。

3樓2014-10-24 16:27:49

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

zhaoshazhu

新蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 125.7
散金: 10
帖子: 51
在線: 38小時
蟲號: 1812564
注冊: 2012-05-12
專業(yè): 環(huán)境化工

引用回帖:

3樓: Originally posted by 月只藍 at 2014-10-24 16:27:49
樓上說初值給得不好是因為 MATLAB 的lsqnonlin 對初值的依賴性很大

您可以幫我設計一下嗎，因為我不是搞模擬的，這只是我論文中的一部分，我是初學者有點不入門。程序里需要加一些殘差平方和，相關(guān)系數(shù)或者其他的，您可以幫我嗎？謝謝了

贊一下

回復此樓

4樓2014-10-24 21:14:08

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

zhaoshazhu

新蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 125.7
散金: 10
帖子: 51
在線: 38小時
蟲號: 1812564
注冊: 2012-05-12
專業(yè): 環(huán)境化工

引用回帖:

2樓: Originally posted by 月只藍 at 2014-10-24 16:25:38
這算是很普遍的問題了，初值給得不好。建議添加相關(guān)系數(shù)或者決定系數(shù)定量地考察擬合效果，只看輸出的參數(shù)數(shù)值沒什么意義。

我現(xiàn)在改初值就會運行的特別慢，計算不出結(jié)果，是因為我的程序?qū)懙牟粚�？您可以幫我看一下嗎，我是初學者，對Matlab不怎么入門，幫我在程序中相應的殘差平方和，殘差，相關(guān)系數(shù)吧。。。謝謝啦

贊一下

回復此樓

5樓2014-10-24 21:17:26

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

zhaoshazhu

新蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 125.7
散金: 10
帖子: 51
在線: 38小時
蟲號: 1812564
注冊: 2012-05-12
專業(yè): 環(huán)境化工

引用回帖:

3樓: Originally posted by 月只藍 at 2014-10-24 16:27:49
樓上說初值給得不好是因為 MATLAB 的lsqnonlin 對初值的依賴性很大

非常感謝您能幫我解答。。

贊一下

回復此樓

6樓2014-10-24 21:19:19

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

月只藍

主管區(qū)長 (職業(yè)作家)

專家經(jīng)驗: +1059
計算強帖: 8
應助: 1712 (講師)
貴賓: 8.888
金幣: 68127.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時
蟲號: 1122189
注冊: 2010-10-14
專業(yè): 宇宙學
管轄: 計算模擬區(qū)

【答案】應助回帖

引用回帖:

5樓: Originally posted by zhaoshazhu at 2014-10-24 21:17:26
我現(xiàn)在改初值就會運行的特別慢，計算不出結(jié)果，是因為我的程序?qū)懙牟粚幔磕梢詭臀铱匆幌聠�，我是初學者，對Matlab不怎么入門，幫我在程序中相應的殘差平方和，殘差，相關(guān)系數(shù)吧。。。謝謝啦...

參見：http://www.gaoyang168.com/bbs/viewthread.php?tid=7603645

贊一下

回復此樓

MATLAB、MS小問題、普通問題請發(fā)帖求助！時間精力有限，恕不接受無償私信求助。

7樓2014-10-25 09:01:06

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

zhaoshazhu

新蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 125.7
散金: 10
帖子: 51
在線: 38小時
蟲號: 1812564
注冊: 2012-05-12
專業(yè): 環(huán)境化工

引用回帖:

7樓: Originally posted by 月只藍 at 2014-10-25 09:01:06
參見：http://www.gaoyang168.com/bbs/viewthread.php?tid=7603645...

謝謝您的指導

回復此樓

8樓2014-10-25 09:54:31

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

zhaoshazhu

新蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 125.7
散金: 10
帖子: 51
在線: 38小時
蟲號: 1812564
注冊: 2012-05-12
專業(yè): 環(huán)境化工

引用回帖:

7樓: Originally posted by 月只藍 at 2014-10-25 09:01:06
參見：http://www.gaoyang168.com/bbs/viewthread.php?tid=7603645...

我之前也嘗試過1stop但是我需要龍格庫塔法計算所以不能使用。我現(xiàn)在也不知道初值怎么選取，是一個很頭疼的問題，有時候根本運算不出結(jié)果，是因為我的程序不對的原因嗎？我想讓您多幫我一下。

贊一下

回復此樓

9樓2014-10-25 10:00:27

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

月只藍

主管區(qū)長 (職業(yè)作家)

專家經(jīng)驗: +1059
計算強帖: 8
應助: 1712 (講師)
貴賓: 8.888
金幣: 68127.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時
蟲號: 1122189
注冊: 2010-10-14
專業(yè): 宇宙學
管轄: 計算模擬區(qū)

【答案】應助回帖

引用回帖:

9樓: Originally posted by zhaoshazhu at 2014-10-25 10:00:27
我之前也嘗試過1stop但是我需要龍格庫塔法計算所以不能使用。我現(xiàn)在也不知道初值怎么選取，是一個很頭疼的問題，有時候根本運算不出結(jié)果，是因為我的程序不對的原因嗎？我想讓您多幫我一下。...

程序能跑通就是對的。
出不來結(jié)果和初值關(guān)系很大，用MATLAB，目前我也沒有很好的辦法能解決這個問題。

我推薦你聯(lián)系一下 dingd，他有高版本的 1stopt軟件，他的軟件幫你做這個問題，就很簡單，而且結(jié)果比MATLAB的好。

贊一下

回復此樓

MATLAB、MS小問題、普通問題請發(fā)帖求助！時間精力有限，恕不接受無償私信求助。

10樓2014-10-25 10:57:21

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

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

返回列表

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

[考研] 293求調(diào)劑 +5	加一一九 2026-03-24	5/250	2026-03-24 11:54 by syl20081243
[考研] 材料調(diào)劑 +5	匹克i 2026-03-23	5/250	2026-03-24 08:50 by dick_runner
[考研] 298-一志愿中國農(nóng)業(yè)大學-求調(diào)劑 +11	手機用戶 2026-03-17	12/600	2026-03-23 23:51 by 熱情沙漠
[考研] 341求調(diào)劑(一志愿湖南大學070300) +5	番茄頭--- 2026-03-22	6/300	2026-03-23 23:45 by Txy@872106
[考研] 材料專業(yè)求調(diào)劑 +11	hanamiko 2026-03-18	11/550	2026-03-23 23:12 by peike
[考研] 265求調(diào)劑 +10	梁梁校校 2026-03-17	10/500	2026-03-23 21:17 by 一切OK
[考研] 350求調(diào)劑 +6	weudhdk 2026-03-19	6/300	2026-03-23 15:47 by tangyuan0840221
[考研] 求調(diào)劑材料學碩080500，總分289分 5+3	@taotao 2026-03-19	21/1050	2026-03-23 10:17 by 冠c哥
[考研] 306求調(diào)劑 +5	來好運來來來 2026-03-22	5/250	2026-03-22 16:17 by BruceLiu320
[考研] 289求調(diào)劑 +7	懷瑾握瑜l 2026-03-20	7/350	2026-03-22 15:57 by ColorlessPI
[考研] 085600材料與化工306 +4	z1z2z3879 2026-03-21	4/200	2026-03-21 23:44 by ms629
[考研] 286分人工智能專業(yè)請求調(diào)劑愿意跨考！ +4	lemonzzn 2026-03-17	8/400	2026-03-21 22:49 by lemonzzn
[考研] 初試 317 +7	半拉月丙 2026-03-20	7/350	2026-03-21 22:26 by peike
[考研] 材料工程專碩 348分求調(diào)劑 +3	冬辭. 2026-03-17	5/250	2026-03-21 18:47 by 學員8dgXkO
[考研] 一志愿深大，0703化學，總分302，求調(diào)劑 +4	七月-七七 2026-03-21	4/200	2026-03-21 18:20 by 學員8dgXkO
[考研] 0805材料320求調(diào)劑 +3	深海物語 2026-03-20	3/150	2026-03-21 15:46 by 無際的草原
[考研] 南昌大學材料專碩311分求調(diào)劑 +6	77chaselx 2026-03-20	6/300	2026-03-21 07:24 by JourneyLucky
[考研] 330求調(diào)劑 +4	小材化本科 2026-03-18	4/200	2026-03-20 23:13 by JourneyLucky
[考研] 求調(diào)劑一志愿南京航空航天大學289分 +3	@taotao 2026-03-19	3/150	2026-03-20 21:34 by JourneyLucky
[考研] 【同濟軟件】軟件（085405）考研求調(diào)劑 +3	2026eternal 2026-03-18	3/150	2026-03-18 19:09 by 搏擊518

24小時熱門版塊排行榜

[求助] 求Matlab高手指導 已有1人參與

» 猜你喜歡

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

【答案】應助回帖

【答案】應助回帖

【答案】應助回帖

【答案】應助回帖

[求助] 求Matlab高手指導已有1人參與

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