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

登錄注冊

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

24小時熱門版塊排行榜

返回列表

ygqsgc

新蟲 (初入文壇)

應助: 0 (幼兒園)
金幣: 20.5
帖子: 9
在線: 31.4小時
蟲號: 5046571
注冊: 2016-09-22
專業(yè): 能源化工

[求助] MATLAB計算動力學參數(shù)時出現(xiàn)誤差大的問題已有2人參與

昨天修改了某個地方導致今天求解參數(shù)誤差極大，但我看了半天也沒找到原因，也忘記昨天修改的什么地方了
程序如下，圖1是昨天的，圖2 是今天的，差距有些大，希望大神能夠解惑

CODE:

function Kinetics2

clear all; clc

format long

k0=[0 0 0 0 0 0 0 0];

tspan=[0 2 4 6 8 10 15 20 25 30 40 50 60 75 90 105 120];             %參數(shù)初值

lb=[0 0 0 0 0 0 0 0];  ub = [1 1 1 1 1 1 1 1]*1e1;  %lb、ub：參數(shù)下限和上限

Y0=[0 0];

Kinetics1=[0        0        0

2        0.101406296        0.025750304

4        0.148139381        0.077208075

6        0.171971874        0.097848253

8        0.200101464        0.131543106

10        0.223088343        0.149258466

15        0.268349285        0.206242074

20        0.319194899        0.249534539

25        0.352139308        0.280899049

30        0.373849138        0.309100951

40        0.4145939        0.348550324

50        0.442947836        0.381294016

60        0.460584916        0.395704095

75        0.483226647        0.421979897

90        0.499672966        0.430587898

105        0.50291736        0.452012783

120        0.518449036        0.452107899];

yexp=Kinetics1(1:17,2:3);            

%使用函數(shù)fmincon()進行參數(shù)估計

[k,fval,flag]=fmincon(@ObjFunc4Fmincon,k0,[],[],[],[],lb,ub,[],[],Y0,yexp);

fprintf('\n 使用函數(shù)fmincon()估計得到的參數(shù)值為:\n')

fprintf('\tk1=%.4f\n',k(1)), fprintf('\tk2=%.4f\n',k(2)), fprintf('\tk3=%.4f\n',k(3)),fprintf('\tk4=%.4f\n',k(4)), fprintf('\tk5=%.4f\n',k(5)),fprintf('\tk6=%.4f\n',k(6)),fprintf('\tk7=%.4f\n',k(7)),fprintf('\tk8=%.4f\n',k(8));

fprintf(' The sum of the squares is:%.1e\n\n', fval);

k_fmincon=k;

ts=0:1:max(tspan);

[ts ys] = ode45(@func,ts,Y0,[],k);

yy = [Kinetics1(:,2) Kinetics1(:,3)];

plot(ts,ys(:,1),'b',tspan,yy(:,1),'bo',ts,ys(:,2),'r',tspan,yy(:,2),'ro'),

legend('Y1的計算值','Y1的實驗值','Y2的計算值','Y2的實驗值','Location','best');

function f=ObjFunc4Fmincon(k,Y0,yexp)

tspan=[0 2 4 6 8 10 15 20 25 30 40 50 60 75 90 105 120];

[t,Y]=ode45(@func,tspan,Y0,[],k);

y(:,1)=Y(:,1);y(:,2:3)=Y(:,1:2);

f=sum((y(:,2)-yexp(:,1)).^2)+sum((y(:,3)-yexp(:,2)).^2)

function dYdt=func(t,Y,k)

k1=k(1);k2=k(2);k3=k(3);k4=k(4);k5=k(5);k6=k(6);k7=k(7);k8=k(8);

f1=(k1.*k2.*k3.*k4.*(1-Y(1)).*6.25.*(1-Y(2)).*3.125-k5.*k6.*k7.*k8.*Y(1).*6.25.*Y(2).*3.125)./(((k1.*k2.*k7+k1.*k2.*k4).*(1-Y(1)).*6.25+(k5.*k3.*k4+k2.*k3.*k4).*(1-Y(2)).*3.125+(k5.*k6.*k7+k1.*k6.*k7).*Y(1).*6.25+(k2.*k6.*k8+k5.*k7.*k8).*Y(2).*3.125+(k1.*k3.*k4+k1.*k2.*k3).*(1-Y(1)).*6.25.*(1-Y(2)).*3.125+(k1.*k6.*k7+k5.*k6.*k4).*(1-Y(1)).*6.25.*Y(1).*6.25+(k2.*k4.*k8+k5.*k3.*k8).*(1-Y(2)).*3.125+(k6.*k7.*k8+k5.*k6.*k8).*Y(1).*6.25.*Y(2).*3.125).*6.125);

f2=(k1.*k2.*k3.*k4.*(1-Y(1)).*6.25.*(1-Y(2)).*3.125-k5.*k6.*k7.*k8.*Y(1).*6.25.*Y(2).*3.125)./(((k1.*k2.*k7+k1.*k2.*k4).*(1-Y(1)).*6.25+(k5.*k3.*k4+k2.*k3.*k4).*(1-Y(2)).*3.125+(k5.*k6.*k7+k1.*k6.*k7).*Y(1).*6.25+(k2.*k6.*k8+k5.*k7.*k8).*Y(2).*3.125+(k1.*k3.*k4+k1.*k2.*k3).*(1-Y(1)).*6.25.*(1-Y(2)).*3.125+(k1.*k6.*k7+k5.*k6.*k4).*(1-Y(1)).*6.25.*Y(1).*6.25+(k2.*k4.*k8+k5.*k3.*k8).*(1-Y(2)).*3.125+(k6.*k7.*k8+k5.*k6.*k8).*Y(1).*6.25.*Y(2).*3.125).*3.125);

dYdt=[f1;f2];

回復此樓

» 猜你喜歡

【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有0人回復
幫我的英語口語老師找學生已經(jīng)有4人回復
物理學I論文潤色/翻譯怎么收費? 已經(jīng)有154人回復
計算三階力常數(shù)用thirdorder擴胞時參數(shù)如何選擇？臨近原子數(shù)怎么取值已經(jīng)有7人回復
2026年循環(huán)經(jīng)濟功能材料國際會議（ICFMCE 2026）已經(jīng)有0人回復
2026年第五屆電氣、電子與信息工程國際會議（ISEEIE 2026）已經(jīng)有0人回復
哈爾濱理工大學物理系招收物理學考研調劑已經(jīng)有0人回復
0702一志愿吉大B區(qū)求調劑已經(jīng)有2人回復

1樓 2021-07-17 20:38:43

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

獨孤神宇

版主 (知名作家)

應助: 490 (碩士)
貴賓: 0.008
金幣: 31016.3
散金: 802
紅花: 122
沙發(fā): 1
帖子: 5600
在線: 856.5小時
蟲號: 3522474
注冊: 2014-11-06
性別: GG
專業(yè): 機械動力學
管轄: 計算模擬

【答案】應助回帖

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

y(:,1)=Y(:,1);y(:,2:3)=Y(:,1:2);
f=sum((y(:,2)-yexp(:,1)).^2)+sum((y(:,3)-yexp(:,2)).^2)
檢查一下這里，尤其是第一行

發(fā)自小木蟲Android客戶端

贊一下

回復此樓

數(shù)值計算

2樓2021-07-17 23:47:54

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

dingd

鐵桿木蟲 (職業(yè)作家)

計算強帖: 4
應助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時
蟲號: 291104
注冊: 2006-10-28

【答案】應助回帖

★ ★ ★ ★ ★ ★ ★
感謝參與，應助指數(shù) +1
獨孤神宇: 金幣+5, 鼓勵交流 2021-07-18 11:01:36
ygqsgc: 金幣+2, ★★★很有幫助 2021-07-23 15:17:58

用1stOpt試了下，似乎是多解吧：

1：
均方差(RMSE): 0.133593203954135
殘差平方和(SSR): 0.571108612567395
相關系數(shù)(R): 0.999611248630951
相關系數(shù)之平方(R^2): 0.999222648389528
修正R平方(Adj. R^2): 0.999028310594448
確定系數(shù)(DC): 0.0502878953255795
F統(tǒng)計(F-Statistic): 1.32838081213258

參數(shù)       最佳估算
-------------------- -------------
k1 0.127201864728704
k2 0.733927119057626
k3 0.153739275177941
k4 0.802137681283057
k5 0.132858854951337
k6 0.190315517131239
k7 0.547878187848255
k8 0.303434451713526

2：
均方差(RMSE): 0.133593202692168
殘差平方和(SSR): 0.571108601777619
相關系數(shù)(R): 0.999611260251213
相關系數(shù)之平方(R^2): 0.999222671621019
修正R平方(Adj. R^2): 0.999028339634003
確定系數(shù)(DC): 0.0502879151372654
F統(tǒng)計(F-Statistic): 1.32838018712688

參數(shù)       最佳估算
-------------------- -------------
k1 0.477143075782806
k2 0.955948247341439
k3 0.221367547288537
k4 0.15762852341868
k5 0.999996228278175
k6 0.39574211416766
k7 0.175285259028897
k8 0.365239513210575

3：
均方差(RMSE): 0.133593217021487
殘差平方和(SSR): 0.571108724292805
相關系數(shù)(R): 0.999611167855669
相關系數(shù)之平方(R^2): 0.999222486901774
修正R平方(Adj. R^2): 0.999028108733572
確定系數(shù)(DC): 0.0502875041008147
F統(tǒng)計(F-Statistic): 1.32840750182457

參數(shù)       最佳估算
-------------------- -------------
k1 0.174350587941431
k2 0.96796438540332
k3 0.997777812182903
k4 0.830312095815316
k5 0.99999988343245
k6 0.757243090774864
k7 0.197815090185031
k8 0.39528880243776

贊一下

回復此樓

3樓2021-07-18 10:50:32

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

ygqsgc

新蟲 (初入文壇)

應助: 0 (幼兒園)
金幣: 20.5
帖子: 9
在線: 31.4小時
蟲號: 5046571
注冊: 2016-09-22
專業(yè): 能源化工

引用回帖:

2樓: Originally posted by 獨孤神宇 at 2021-07-17 23:47:54
y(:,1)=Y(:,1);y(:,2:3)=Y(:,1:2);
f=sum((y(:,2)-yexp(:,1)).^2)+sum((y(:,3)-yexp(:,2)).^2)
檢查一下這里，尤其是第一行

謝謝，但是我修改這里好像沒什么變化。我今天又發(fā)現(xiàn)初值變化會很大的影響結果，請問您有什么建議嗎？

贊一下

回復此樓

4樓2021-07-18 16:37:00

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

獨孤神宇

版主 (知名作家)

應助: 490 (碩士)
貴賓: 0.008
金幣: 31016.3
散金: 802
紅花: 122
沙發(fā): 1
帖子: 5600
在線: 856.5小時
蟲號: 3522474
注冊: 2014-11-06
性別: GG
專業(yè): 機械動力學
管轄: 計算模擬

【答案】應助回帖

★ ★ ★
ygqsgc: 金幣+3, ★★★很有幫助 2021-07-23 15:18:05

引用回帖:

4樓: Originally posted by ygqsgc at 2021-07-18 16:37:00
謝謝，但是我修改這里好像沒什么變化。我今天又發(fā)現(xiàn)初值變化會很大的影響結果，請問您有什么建議嗎？...

這個需要多次嘗試了，可以將結果作為初始值代入反復計算。

可以參考三樓的結果，將其作為初值代入計算

贊一下

回復此樓

數(shù)值計算

5樓2021-07-18 17:32:29

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

ygqsgc

新蟲 (初入文壇)

應助: 0 (幼兒園)
金幣: 20.5
帖子: 9
在線: 31.4小時
蟲號: 5046571
注冊: 2016-09-22
專業(yè): 能源化工

引用回帖:

5樓: Originally posted by 獨孤神宇 at 2021-07-18 17:32:29
這個需要多次嘗試了，可以將結果作為初始值代入反復計算。

可以參考三樓的結果，將其作為初值代入計算...

謝謝！

回復此樓

6樓2021-07-18 18:50:07

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

ygqsgc

新蟲 (初入文壇)

應助: 0 (幼兒園)
金幣: 20.5
帖子: 9
在線: 31.4小時
蟲號: 5046571
注冊: 2016-09-22
專業(yè): 能源化工

引用回帖:

3樓: Originally posted by dingd at 2021-07-18 10:50:32
用1stOpt試了下，似乎是多解吧：

1：
均方差(RMSE): 0.133593203954135
殘差平方和(SSR): 0.571108612567395
相關系數(shù)(R): 0.999611248630951
相關系數(shù)之平方(R^2): 0.999222648389528
修正R平方(Adj. R^2) ...

謝謝您！我剛下載了這個軟件，還不會用，冒昧問一句能看一下您的代碼嗎？

贊一下

回復此樓

7樓2021-07-19 15:29:27

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

相關版塊跳轉我要訂閱樓主 ygqsgc 的主題更新

返回列表

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

[考研] 新疆大學地質與礦業(yè)工程學院招生 +24	another12 2026-03-04	32/1600	2026-03-10 01:15 by 1111--
[考研] 一志愿鄭州大學，學碩，物理化學， 333求調劑 +7	李魔女斗篷 2026-03-04	8/400	2026-03-09 19:58 by L135790
[考研] 一志愿山東大學，總分327，英語二79，有論文，有競賽，已過四六級 +3	木木目目1 2026-03-09	3/150	2026-03-09 19:52 by yuningshan
[考研] 求調劑，數(shù)一英一274分 +4	小菲會努力 2026-03-08	4/200	2026-03-09 12:40 by 一定上岸喲_
[考研] 一志愿武理314求調劑 +4	（￣～￣;） 2026-03-08	5/250	2026-03-08 23:12 by L135790
[考研] 269求調劑 +3	朔朔話 2026-03-08	4/200	2026-03-08 20:39 by 熱情沙漠
[考研] 求調劑 +3	拾柒12。 2026-03-08	3/150	2026-03-08 17:33 by bingxueer79
[考研] 求調劑，一志愿華中科大0702，數(shù)一英一，293 +4	小羅露一二 2026-03-07	4/200	2026-03-08 16:36 by 星空星月
[考研] 0703化學求調劑學校 +3	一志愿京區(qū)211 2026-03-04	3/150	2026-03-07 21:35 by yinhuanshun
[考研] 0307化學求調劑 +6	0ok0k 2026-03-07	6/300	2026-03-07 20:10 by pies112
[考研] 267化工調劑求助 +7	聰少OZ 2026-03-04	7/350	2026-03-07 17:33 by 無際的草原
[考博] 2026年博士名額撿漏 +4	科研ya 2026-03-04	7/350	2026-03-06 16:05 by 科研ya
[考研] 材料與化工304求B區(qū)調劑 +4	邱gl 2026-03-06	4/200	2026-03-06 15:51 by 聰明的大松鼠
[考研] 316求調劑 +3	林小星發(fā)大財 2026-03-05	5/250	2026-03-06 09:05 by Iveryant
[考研] 282求調劑 +7	夕～日 2026-03-05	8/400	2026-03-05 21:31 by zzpnuli111
[考研] 材料調劑 +5	ws 上岸鴨 2026-03-05	5/250	2026-03-05 20:06 by houyaoxu
[考研] 求調劑，學校研究所都可以，材料與化工267分 +6	wmx1 2026-03-05	6/300	2026-03-05 20:01 by houyaoxu
[考研] 320材料一志愿西工大專碩英二數(shù)二有兩年光伏工作經(jīng)驗 +5	錘某人 2026-03-04	8/400	2026-03-05 09:24 by oxidpl
[考研] 070300化學 280 一志愿太原理工求調劑 +8	拾玖壹 2026-03-04	8/400	2026-03-05 08:52 by anny19840123
[考研] 0703化學一志愿東北大學 279分求調劑 +4	聶樂平 2026-03-03	4/200	2026-03-04 16:31 by zhukairuo

24小時熱門版塊排行榜

[求助] MATLAB計算動力學參數(shù)時出現(xiàn)誤差大的問題 已有2人參與

» 猜你喜歡

【答案】應助回帖

【答案】應助回帖

【答案】應助回帖

[求助] MATLAB計算動力學參數(shù)時出現(xiàn)誤差大的問題已有2人參與