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

返回列表

自制奶酪

銅蟲(chóng) (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 3
散金: 6
帖子: 106
在線(xiàn): 53.4小時(shí)
蟲(chóng)號(hào): 3447629
注冊(cè): 2014-09-28
性別: GG
專(zhuān)業(yè): 綠色有機(jī)化學(xué)

[求助] 用matlab求動(dòng)力學(xué)參數(shù) 已有2人參與

我不知道MATLAB怎么輸入求動(dòng)力學(xué)函數(shù)的方程，抄了木蟲(chóng)大佬的內(nèi)容，請(qǐng)大家?guī)臀野阉侥苓\(yùn)行出來(lái)，動(dòng)力學(xué)參數(shù)擬合問(wèn)題，求助大家~~~

function odes_fit
format long
clear all
clc
k0=[0 0 0 0 0 0 0 0];%參數(shù)初值
lb=[0 0 0 0 0 0 0 0];ub=[+inf +inf +inf +inf +inf +inf +inf +inf];%lb、ub:參數(shù)下限和上限
x0=[0 0 0 0 0];

data=...
[0.5       1          1.5          2          3          4          5          6
      18.7 6.5 2.7 2.1 1.7 1.5 1.2 0.2;
      48.9 31.3 15.3 6.9 2.2 1.8 1.7 0.5;
      19.3 34.7 41.2 40.9 37.2 31.7 26.1 19.4;
      9.7 19.6 28.7 34.5 38.4 38.2 36.3 33.4;
      0.2 2.6 5.7 9.3 17.1 24.8 30.5 36.7;
];
x0=data(1,2:end);
tspan = [data(:,1)'];
yexp = [data(2:end,2) data(2:end,3) data(2:end,4) data(2:end,5) data(2:end,6) data(2:end,7) data(2:end,8)];
%使用函數(shù)fmincon()進(jìn)行參數(shù)估計(jì)
[k,fval,flag]=fmincon(@ObjFunc4Fmincon,k0,[],[],[],[],lb,ub,[],[],x0,yexp);
fprintf('\n\n使用函數(shù)lsqnonlin()估計(jì)得到的參數(shù)值為:\n')
fprintf('\tk1 = %.9f \n',k(1))
fprintf('\tk2 = %.9f \n',k(2))
fprintf('\tk3 = %.9f \n',k(3))
fprintf('\tk4 = %.9f \n',k(4))
fprintf('\tk5 = %.9f \n',k(5))
fprintf('\tk6 = %.9f \n',k(6))
fprintf('The sum of the squares is:%.le\n\n',fval)
k_fmincon=k;
%使用函數(shù)lsqnonlin()進(jìn)行參數(shù)估計(jì)
[k,resnorm,residual,exitflag,output,lambda,jacobian]=...
lsqnonlin(@ObjFunc4LNL,k0,lb,ub,[],x0,yexp);
ci=nlparci(k,residual,jacobian);
fprintf('\n\n使用函數(shù)lsqnonlin()估計(jì)得到的參數(shù)值:\n'),output
%以函數(shù)fmincon()估計(jì)得到的結(jié)果為初值，使用函數(shù)lsqnonlin()進(jìn)行參數(shù)估計(jì)
k0=k_fmincon;
[k,resnorm,residual,exitflag,output,lambda,jacobian]=lsqnonlin(@ObjFunc4LNL,k0,lb,ub,[],x0,yexp);
ci=nlparci(k,residual,jacobian);
fprintf('\n\n以fmincon()的結(jié)果為初值，使用函數(shù)lsqnonlin()估計(jì)得到的參數(shù)值:\n')
output
figure(1)

ts=0

(max(tspan)-min(tspan))/100):max(tspan);
[ts, ys] = ode45(@KineticsEqs,ts,x0,[],k);
yy = [data(:,2) data(:,3) data(:,4) data(:,5) data(:,6)];
figure(1)
plot(ts,ys(:,1),'b',tspan,yy(:,1),'bo');
figure(2)
plot(ts,ys(:,2),'r',tspan,yy(:,2),'ro');
figure(3)
plot(ts,ys(:,3),'k',tspan,yy(:,3),'ko');
figure(4)
plot(ts,ys(:,4),'g',tspan,yy(:,4),'ko');
figure(5)
plot(ts,ys(:,4),'g',tspan,yy(:,4),'ko');
figure(6)
plot(ts,ys(:,6),'m',tspan,yy(:,6),'mo');
%legend('C1的計(jì)算值','C1的實(shí)驗(yàn)值','C2的計(jì)算值','C2的實(shí)驗(yàn)值','C3的計(jì)算值','C3的實(shí)驗(yàn)值','C4的計(jì)算值','C4的實(shí)驗(yàn)值','C5的計(jì)算值','C5的實(shí)驗(yàn)值','C6的計(jì)算值','Location','best');

function f = ObjFunc(k,tspan,x0,yexp) % 目標(biāo)函數(shù)
[t, Xsim] = ode45(@KineticsEqs,tspan,x0,[],k);
Xsim1=Xsim(:,1);
Xsim2=Xsim(:,2);
Xsim3=Xsim(:,3);
Xsim4=Xsim(:,4);
Xsim5=Xsim(:,5);
Xsim6=Xsim(:,6);

ysim(:,1) = Xsim1(2:end);
ysim(:,2) = Xsim2(2:end);
ysim(:,3) = Xsim3(2:end);
ysim(:,4) = Xsim4(2:end);
ysim(:,5) = Xsim5(2:end);
ysim(:,6) = Xsim6(2:end);

f = [(ysim(:,1)-yexp(:,1)) (ysim(:,2)-yexp(:,2)) (ysim(:,3)-yexp(:,3)) (ysim(:,5)-yexp(:,5)) ...
(ysim(:,6)-yexp(:,6))];

function dCdt = KineticsEqs(t,C,k) % ODE模型方程
C1=C(1);C2=C(2);C3=C(3);C4=C(4);C5=C(5);C6=C(6);
k1=k(1);k2=k(2);k3=k(3);k4=k(4);k5=k(5);k6=k(6);k7=k(7);
dC1dt = -k1*C1-k7*C1;
dC2dt = k1*C1-k2*C2-k5*C2;
dC3dt = k2*C2-k3*C3-k6*C3;
dC4dt = k3*C3-k4*C5;
dC5dt = k4*C5;
dC6dt = k5*C2+k6*C3+k7*C1;

dCdt = [dC1dt;dC2dt;dC3dt;dC4dt;dC5dt;dC6dt]; @beefly

回復(fù)此樓

» 猜你喜歡

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

簡(jiǎn)單快樂(lè)

1樓 2021-01-05 20:42:15

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

獨(dú)孤神宇

版主 (知名作家)

應(yīng)助: 490 (碩士)
貴賓: 0.008
金幣: 31016.3
散金: 802
紅花: 122
沙發(fā): 1
帖子: 5600
在線(xiàn): 856.5小時(shí)
蟲(chóng)號(hào): 3522474
注冊(cè): 2014-11-06
性別: GG
專(zhuān)業(yè): 機(jī)械動(dòng)力學(xué)
管轄: 計(jì)算模擬

【答案】應(yīng)助回帖

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

這個(gè)代碼不完整

lsqnonlin 函數(shù) ObjFunc4LNL 缺失

贊一下

回復(fù)此樓

數(shù)值計(jì)算

2樓2021-01-05 21:03:26

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

自制奶酪

銅蟲(chóng) (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 3
散金: 6
帖子: 106
在線(xiàn): 53.4小時(shí)
蟲(chóng)號(hào): 3447629
注冊(cè): 2014-09-28
性別: GG
專(zhuān)業(yè): 綠色有機(jī)化學(xué)

引用回帖:

2樓: Originally posted by 獨(dú)孤神宇 at 2021-01-05 21:03:26
這個(gè)代碼不完整

lsqnonlin 函數(shù) ObjFunc4LNL 缺失

您好，這個(gè)要怎么修改啊，求助

贊一下

回復(fù)此樓

簡(jiǎn)單快樂(lè)

3樓2021-01-05 21:07:08

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

dingd

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

計(jì)算強(qiáng)帖: 4
應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線(xiàn): 1223.7小時(shí)
蟲(chóng)號(hào): 291104
注冊(cè): 2006-10-28

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
自制奶酪: 金幣+20, ★有幫助 2021-01-05 21:39:36
獨(dú)孤神宇: 金幣+5, 鼓勵(lì)交流 2021-01-06 17:08:32

參考下1stOpt的計(jì)算結(jié)果：
Weighted Root of Mean Square Error (RMSE): 3.52215359203778
Weighted Sum of Squared Residual: 434.194807406662
Correlation Coef. (R): 0.967264353331916
R-Square: 0.935600329226609
Determination Coef. (DC): 0.698480780131114
F-Statistic: -0.866225958773538

Parameter Best Estimate
-------------------- -------------
k1 1.18206935637483
k2 1.23804230026633
k3 0.392064236557302
k4 0.237399839174299
k5 2.4090372987211E-18
k6 1.44475738382174E-24
k7 0.719127412563461
f Initial Value 19.2237430015159

贊一下

回復(fù)此樓

4樓2021-01-05 21:34:03

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

流浪在地球

新蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 80.8
散金: 8
帖子: 33
在線(xiàn): 32分鐘
蟲(chóng)號(hào): 29690859
注冊(cè): 2022-06-09
專(zhuān)業(yè): 動(dòng)力學(xué)與控制

您好

發(fā)自小木蟲(chóng)Android客戶(hù)端

回復(fù)此樓

5樓2022-11-17 13:58:31

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

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

返回列表

普通表情龍兔虎貓

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

[考研] 調(diào)劑的同學(xué)，走過(guò)路過(guò)，不要錯(cuò)過(guò) +9	likeihood 2026-03-06	14/700	2026-03-08 14:37 by Wh2wyx1314.
[考研] 337一志愿華南理工材料求調(diào)劑 +4	mysdl 2026-03-07	4/200	2026-03-08 11:00 by 無(wú)際的草原
[考研] 085701環(huán)境工程317分求調(diào)劑 +9	6汆尼9 2026-03-07	9/450	2026-03-08 06:41 by 劉兵
[考研] 334求調(diào)劑 +7	Trying] 2026-03-06	7/350	2026-03-08 03:55 by lfhuang
[考研] 求調(diào)劑，一志愿華中科大0702，數(shù)一英一，293 +3	小羅露一二 2026-03-07	3/150	2026-03-07 23:58 by 1710470214
[考研] 一志愿鄭大071000分?jǐn)?shù)282求調(diào)劑 +3	研研顏 2026-03-05	7/350	2026-03-07 22:33 by 帆船哥
[考研] 08工科 +4	li李樂(lè)成 2026-03-06	4/200	2026-03-07 21:46 by 韓應(yīng)許
[考研] 一志愿211 化學(xué)305分求調(diào)劑 +3	0703楊?lèi)?05分 2026-03-05	3/150	2026-03-07 21:15 by yinhuanshun
[考研] 材料專(zhuān)碩調(diào)劑 +5	慕辰123 2026-03-05	8/400	2026-03-07 15:54 by xuesongli73
[考研] 289求調(diào)劑 +3	馳光絕景 2026-03-04	4/200	2026-03-06 12:34 by 馳光絕景
[考研] 320材料一志愿西工大專(zhuān)碩英二數(shù)二有兩年光伏工作經(jīng)驗(yàn) +5	錘某人 2026-03-04	8/400	2026-03-05 09:24 by oxidpl
[考研] 調(diào)劑材料學(xué)碩 +6	詞凝Y 2026-03-02	6/300	2026-03-05 07:10 by 昆工985
[考研] 0703化學(xué)求調(diào)劑，總分320分，一志愿中南民族大學(xué) +3	19286712440 2026-03-04	3/150	2026-03-04 18:30 by houyaoxu
[考研] 本科太原理工采礦工程，求調(diào)劑 +3	onlx 2026-03-01	3/150	2026-03-04 15:57 by Stephen_ym
[考研] 085600 英一數(shù)二272求調(diào)劑 5+6	vida_a 2026-03-01	47/2350	2026-03-04 13:35 by vida_a
[考研] 331求調(diào)劑 +3	zzZ&zZ 2026-03-03	3/150	2026-03-04 10:37 by aaa顯卡批發(fā)王總
[考研] 化工270求調(diào)劑 +10	什么名字qwq 2026-03-02	10/500	2026-03-03 13:54 by NUAAZXWS
[考研] 278求調(diào)劑 +3	滿(mǎn)天星11_22 2026-03-02	3/150	2026-03-03 13:51 by Iveryant
[考研] 26考研報(bào)考西工大材料308分求調(diào)劑 +4	weizhong123 2026-03-01	5/250	2026-03-03 12:22 by weizhong123
[考研] 【2026 碩士調(diào)劑】課題組招收調(diào)劑生 +3	考研版棒棒 2026-03-02	5/250	2026-03-03 01:45 by kkky.