| 2 | 1/1 | 返回列表 |
| 查看: 1494 | 回復: 1 | |||
[求助]
常微分方程的參數(shù)擬合 已有1人參與
|
|
參數(shù)擬合一直出錯,但是找不出錯誤。哪位能幫忙給解決一下?先謝謝了! function Par_SI clear clc ID=[1 3 25 72 222 282 256 233 189 123 70 25 11 4]; X0=[762,1];%初始條件 lb=[0,0];%約束條件,參數(shù)值的下界 ub=[0.1,1];%約束條件,參數(shù)值的下界 par1guess=[0.01 0.1];%參數(shù)值的初始 options=optimset('Display','final','MaxIter',2000,'MaxFunEvals',2000);%求最優(yōu)解的優(yōu)化條件 [par1,fval]=fmincon(LSmin,par1guess,[],[],[],[],lb,ub,[],options,ID,X0); %利用優(yōu)化算法得到的參數(shù)擬合數(shù)據(jù) [T1,X1]=ode45(SImodel,[0 14],X0,[],par1);%利用估計參數(shù)重新運行模型 figure(1) plot(T1,X1(:,1),'k-'); hold on plot(T1,X1(:,2),'k-',0:14,ID,'ko') %定義最小二乘的目標函數(shù) function d=LSmin(par1,ID,X0) [T,x]=ode45(SImodel,[0:1:14],X0,[],par1) d=norm(x(:,2)'-ID)^2; %定義微分方程 function dy=SImodel(t,y,par1) dy=zeros(2,1); beta=par1(1); gamma=par1(2); dy=[-beta*y(1)*y(2);beta*y(1)*y(2)-gamma*y(2)]; |

主管區(qū)長 (職業(yè)作家)
![]() |
專家經(jīng)驗: +1059 |

| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 263求調劑 +6 | yqdszhdap- 2026-03-22 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 307求調劑 +3 | 余意卿 2026-03-21 | 3/150 |
|
|
[考研] 276求調劑 +3 | YNRYG 2026-03-21 | 4/200 |
|
|
[考研] 一志愿中南大學化學學碩0703總分337求調劑 +4 | niko- 2026-03-22 | 4/200 |
|
|
[考研] 一志愿 西北大學 ,070300化學學碩,總分287,雙非一本,求調劑。 +3 | 晨昏線與星海 2026-03-20 | 3/150 |
|
|
[考研] 269專碩求調劑 +6 | 金恩貝 2026-03-21 | 6/300 |
|
|
[考研] 384求調劑 +3 | 子系博 2026-03-22 | 4/200 |
|
|
[考研] 廣西大學材料導師推薦 +3 | 夏夏夏小正 2026-03-17 | 5/250 |
|
|
[考研] 0805 316求調劑 +3 | 大雪深藏 2026-03-18 | 3/150 |
|
|
[考研] 求助 +5 | 夢里的無言 2026-03-21 | 6/300 |
|
|
[考研] 材料 271求調劑 +5 | 展信悅_ 2026-03-21 | 5/250 |
|
|
[考研] 22 350 本科985求調劑,求老登收留 +3 | 李軼男003 2026-03-20 | 3/150 |
|
|
[考研] 材料學學碩080502 337求調劑-一志愿華中科技大學 +4 | 順順順mr 2026-03-18 | 5/250 |
|
|
[考研] 265求調劑 +9 | 梁梁校校 2026-03-17 | 9/450 |
|
|
[考研] 295材料求調劑,一志愿武漢理工085601專碩 +5 | Charlieyq 2026-03-19 | 5/250 |
|
|
[考研] 一志愿南理工085701環(huán)境302求調劑院校 +3 | 葵梓衛(wèi)隊 2026-03-20 | 3/150 |
|
|
[考研] 一志愿中國海洋大學,生物學,301分,求調劑 +5 | 1孫悟空 2026-03-17 | 6/300 |
|
|
[考研] 288求調劑,一志愿華南理工大學071005 +5 | ioodiiij 2026-03-17 | 5/250 |
|
|
[考研] 收復試調劑生 +4 | 雨后秋荷 2026-03-18 | 4/200 |
|
|
[考研] 0703化學調劑 +3 | 妮妮ninicgb 2026-03-17 | 3/150 |
|