版塊導(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
在線: 53.4小時(shí)
蟲(chóng)號(hào): 3447629
注冊(cè): 2014-09-28
性別: GG
專業(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ù)此樓

» 猜你喜歡

2026年循環(huán)經(jīng)濟(jì)功能材料國(guó)際會(huì)議（ICFMCE 2026）已經(jīng)有0人回復(fù)
2026年第五屆電氣、電子與信息工程國(guó)際會(huì)議（ISEEIE 2026）已經(jīng)有0人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有114人回復(fù)
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑已經(jīng)有5人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有5人回復(fù)
求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
請(qǐng)問(wèn)還有沒(méi)有用Latex寫(xiě)文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有26人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運(yùn)，祝面上順利！已經(jīng)有34人回復(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
在線: 856.5小時(shí)
蟲(chóng)號(hào): 3522474
注冊(cè): 2014-11-06
性別: GG
專業(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
在線: 53.4小時(shí)
蟲(chóng)號(hào): 3447629
注冊(cè): 2014-09-28
性別: GG
專業(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
在線: 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
在線: 32分鐘
蟲(chóng)號(hào): 29690859
注冊(cè): 2022-06-09
專業(yè): 動(dòng)力學(xué)與控制

您好

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

回復(fù)此樓

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

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

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

返回列表

普通表情龍兔虎貓

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

[考研] 268求調(diào)劑 +7	好運(yùn)連綿不絕 2026-03-12	8/400	2026-03-17 20:28 by xilongliang
[考研] 281求調(diào)劑（0805） +3	煙汐憶海 2026-03-16	7/350	2026-03-17 20:16 by peike
[考研] 328求調(diào)劑，英語(yǔ)六級(jí)551，有科研經(jīng)歷 +3	生物工程調(diào)劑 2026-03-16	8/400	2026-03-17 19:03 by Wangjingyue
[考研] 考研化學(xué)學(xué)碩調(diào)劑，一志愿985 +4	張vvvv 2026-03-15	6/300	2026-03-17 17:15 by ruiyingmiao
[考博] 26申博 +4	八旬速覽 2026-03-16	4/200	2026-03-17 13:00 by 輕松不少隨
[考研] 機(jī)械專碩325，尋找調(diào)劑院校 +3	y9999 2026-03-15	5/250	2026-03-16 19:58 by y9999
[基金申請(qǐng)] 今年的國(guó)基金是打分制嗎？ 50+3	zhanghaozhu 2026-03-14	3/150	2026-03-16 17:07 by 北京萊茵潤(rùn)色
[考研] 304求調(diào)劑 +3	曼殊2266 2026-03-14	3/150	2026-03-16 16:39 by houyaoxu
[考研] 085600調(diào)劑 +5	漾漾123sun 2026-03-12	6/300	2026-03-16 15:58 by 漾漾123sun
[考研] 285求調(diào)劑 +6	ytter 2026-03-12	6/300	2026-03-16 15:05 by njzyff
[考研] 0856求調(diào)劑 +3	劉夢(mèng)微 2026-03-15	3/150	2026-03-16 10:00 by houyaoxu
[考博] 東華理工大學(xué)化材專業(yè)26屆碩士博士申請(qǐng) +6	zlingli 2026-03-13	6/300	2026-03-15 20:00 by ryzcf
[考研] 080500，材料學(xué)碩302分求調(diào)劑學(xué)校 +4	初識(shí)可樂(lè) 2026-03-14	5/250	2026-03-14 21:08 by peike
[考研] 材料工程327求調(diào)劑 +3	xiaohe12w 2026-03-11	3/150	2026-03-14 20:20 by ms629
[考研] 279求調(diào)劑 +3	Dizzy123@ 2026-03-10	3/150	2026-03-13 23:02 by JourneyLucky
[考研] 332求調(diào)劑 +3	zjy101327 2026-03-11	6/300	2026-03-13 22:48 by JourneyLucky
[考研] 0703化學(xué)調(diào)劑 +4	快樂(lè)的香蕉 2026-03-11	4/200	2026-03-13 22:41 by JourneyLucky
[考研] 0856材料與化工301求調(diào)劑 +5	奕束光 2026-03-13	5/250	2026-03-13 22:00 by 星空星月
[考研] 【考研調(diào)劑求收留】 +3	Ceciilia 2026-03-11	3/150	2026-03-13 20:18 by JourneyLucky
[考研] 295求調(diào)劑 +3	小匕仔汁 2026-03-12	3/150	2026-03-13 15:17 by vgtyfty

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

24小時(shí)熱門(mén)版塊排行榜

自制奶酪

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

» 猜你喜歡

獨(dú)孤神宇

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

自制奶酪

dingd

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

流浪在地球

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产高清中文字幕,99re热久久亚洲综合精品成人,熟妇一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频