 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 1 | 1/1 | 返回列表 |
| 查看: 1315 | 回復(fù): 0 | |||
| 【懸賞金幣】回答本帖問(wèn)題，作者757272131將贈(zèng)送您 10 個(gè)金幣 | |||
757272131新蟲 (初入文壇)
|
[求助]
使用matlab最小二乘法擬合求解常微分方程組未知參數(shù)
|
||
|
function k1k2k32 format long clear all clc tspan = [0 6 24 44 68 72 74 92 104 116]';%%這是時(shí)間 yexp= [3.111,3.639,3.887,4.289,4.658,5.531,6.218,6.979,7.111,7.114]';%%%這是菌落總數(shù) x0 = [3.111]; k0 = [0.1 3 45 4 37 8.1]; lb = [0 0 0 0 0 0]; ub = [100 100 100 100 100 100]; % % % opts = statset('nlinfit'); % % opts.RobustWgtFun = 'bisquare'; % % mdll = fitnlm(tspan,yexp,ObjFunc,b0,'Options') % mdll = fitnlm(tspan,yexp,ObjFunc,b0,'Options',opts) options = optimoptions(@lsqnonlin,'Algorithm','trust-region-reflective'); [k,resnorm,residual,exitflag,output,lambda,jacobian] = ... lsqnonlin(@ObjFunc,k0,lb,ub,[],tspan,x0,yexp); ci = nlparci(k,residual,jacobian); fprintf('\n\n使用函數(shù)lsqnonlin()估計(jì)得到的參數(shù)值為:\n') fprintf('\tk1 = %.9f ± %.9f\n',k(1),ci(1,2)-k(1)) fprintf('\tk2 = %.9f ± %.9f\n',k(2),ci(2,2)-k(2)) fprintf('\tk3 = %.9f ± %.9f\n',k(3),ci(3,2)-k(3)) fprintf('\tk4 = %.9f ± %.9f\n',k(4),ci(4,2)-k(4)) fprintf('\tk5 = %.9f ± %.9f\n',k(5),ci(5,2)-k(5)) fprintf('\tk6 = %.9f ± %.9f\n',k(6),ci(6,2)-k(6)) fprintf(' The sum of the squares is: %.9e\n\n',resnorm) tsa=0:0.01:max(tspan); [tsa ysa]=ode45(@KineticsEqs,tsa,x0,[],k); figure(1), plot(tsa,ysa(:,1),'b',tspan,yexp,'or'),legend('計(jì)算值','實(shí)驗(yàn)值','Location','best'); function f = ObjFunc(k,tspan,x0,yexp) % 目標(biāo)函數(shù) [tspan y] = ode45(@KineticsEqs,tspan,x0,[],k); f = (y-yexp); function dCdt = KineticsEqs(tspan ,C,k) T=4*(tspan >=0 & tspan <=30)+10*(tspan >31 & tspan <=90)+4*(tspan >90 ); % 假設(shè)這是我寫的溫度關(guān)于時(shí)間的分段函數(shù) umax=((k(2)*(T-k(3))*(T-k(4))*(T-k(4)))/(((k(5)-k(4))*(T-k(5))-(k(5)-k(3))*(k(5)+k(4)-2*T))*(k(5)-k(4)))); dCdt=(1/(1+exp(-4*(tspan -k(1)))))*umax*(1-exp(C-k(6))); 這是根據(jù)他人的代碼更改的，得到的結(jié)果是置信區(qū)間過(guò)大？想有什么辦法解決？以及想用fitnlm函數(shù)進(jìn)行擬合，達(dá)到一些參數(shù)評(píng)價(jià)t值p值等，但用fitnlm運(yùn)行不起來(lái)。 @beefly |
回復(fù)此樓» 猜你喜歡
幫我的英語(yǔ)口語(yǔ)老師找學(xué)生
已經(jīng)有4人回復(fù)
-
計(jì)算三階力常數(shù)用thirdorder擴(kuò)胞時(shí)參數(shù)如何選擇？臨近原子數(shù)怎么取值
已經(jīng)有7人回復(fù)
-
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有203人回復(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)有1人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑
已經(jīng)有2人回復(fù)
-
求調(diào)劑
已經(jīng)有0人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑有論文
已經(jīng)有0人回復(fù)
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 757272131 的主題更新
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 301求調(diào)劑 +3 | Liyouyumairs 2026-03-11 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 接受26屆調(diào)劑生 +22 | 豬豬豬毛 2026-03-06 | 23/1150 |
|
|
|
[考研] 290求調(diào)劑 +3 | 柯淮然 2026-03-10 | 8/400 |
|
|
|
[考研] 材料371求調(diào)劑 +7 | 鱷魚? 2026-03-11 | 9/450 |
|
|
|
[考研] 0856材料與化工353分求調(diào)劑 +11 | NIFFFfff 2026-03-09 | 11/550 |
|
|
|
[考研] 308求調(diào)劑 +4 | 是Lupa啊 2026-03-08 | 7/350 |
|
|
|
[考研] 311求調(diào)劑 +3 | 牛乳糖的卡卡 2026-03-10 | 3/150 |
|
|
|
[考研] 0860求調(diào)劑（272分） +3 | lllllcsjsj 2026-03-05 | 4/200 |
|
|
|
[考研] 337一志愿華南理工材料求調(diào)劑 +5 | mysdl 2026-03-07 | 5/250 |
|
|
|
[考研] 材料與化工求調(diào)劑 +11 | 與冬清寧 2026-03-07 | 14/700 |
|
|
|
[考研] 一志愿：武漢理工，材料工程，英二數(shù)二 總分314 +3 | 2202020125 2026-03-10 | 4/200 |
|
|
|
[考研] 一志愿清華深研院材料專碩294分，專業(yè)課111分，本科中南大學(xué)材料，有六級(jí)，有工作經(jīng)驗(yàn) +4 | H14528 2026-03-04 | 4/200 |
|
|
|
[考研] 化工284求調(diào)劑 +12 | 小黑想變白 2026-03-07 | 12/600 |
|
|
|
[考研] 320求調(diào)劑 +4 | 魏zy 2026-03-08 | 4/200 |
|
|
|
[考研] 安徽農(nóng)業(yè)大學(xué)材料與化學(xué)學(xué)院0856材料招收調(diào)劑 +3 | akakk47 2026-03-05 | 9/450 |
|
|
|
[考研] 085701環(huán)境工程317分求調(diào)劑 +9 | 6汆尼9 2026-03-07 | 9/450 |
|
|
|
[考研] 一志愿鄭大071000分?jǐn)?shù)282求調(diào)劑 +3 | 研研顏 2026-03-05 | 7/350 |
|
|
|
[考研] 材料277分求調(diào)劑 +13 | 飯飯星球 2026-03-04 | 14/700 |
|
|
|
[考研] 08工科求調(diào)劑 +3 | 隆LLL 2026-03-06 | 4/200 |
|
|
|
[考研] 求材料調(diào)劑 +4 | berdmond 2026-03-05 | 4/200 |
|