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

返回列表

jwj216

金蟲 (正式寫手)

應(yīng)助: 3 (幼兒園)
金幣: 348.5
散金: 1823
紅花: 6
帖子: 399
在線: 163.3小時(shí)
蟲號(hào): 3325146
注冊(cè): 2014-07-16
性別: MM
專業(yè): 數(shù)理統(tǒng)計(jì)

[求助] 常微分方程的參數(shù)擬合已有1人參與

參數(shù)擬合一直出錯(cuò)，但是找不出錯(cuò)誤。哪位能幫忙給解決一下？先謝謝了！

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);%利用估計(jì)參數(shù)重新運(yùn)行模型
figure(1)
plot(T1,X1(:,1),'k-');
hold on
plot(T1,X1(:,2),'k-',0:14,ID,'ko')

%定義最小二乘的目標(biāo)函數(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)];

回復(fù)此樓

» 猜你喜歡

求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費(fèi)? 已經(jīng)有165人回復(fù)
請(qǐng)問還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有32人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運(yùn)，祝面上順利！已經(jīng)有76人回復(fù)
華東師范大學(xué)芯片設(shè)計(jì)徐瓏真實(shí)水平咋樣已經(jīng)有1人回復(fù)
調(diào)劑已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有0人回復(fù)
中德博士后2026項(xiàng)目5月30日截止- 國內(nèi)博士后入站、德國兩年全職已經(jīng)有0人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

請(qǐng)教一下各位大神怎么求解 Runge–Kutta method的微分方程的各個(gè)系數(shù) 已經(jīng)有12人回復(fù)
1stopt軟件擬合非線性方程，計(jì)算結(jié)果全為0，這是為什么已經(jīng)有3人回復(fù)
求大神幫你擬合一組數(shù)據(jù)~萬分感激。。。已經(jīng)有6人回復(fù)
使用matlab最優(yōu)化方法擬合獲得多個(gè)動(dòng)力學(xué)參數(shù)中的問題已經(jīng)有4人回復(fù)
求大神幫忙擬合一個(gè)非線性方程，求出模型參數(shù) 已經(jīng)有15人回復(fù)
用1stopt軟件做微分方程組的參數(shù)估計(jì) 已經(jīng)有13人回復(fù)
求助用matlab擬合動(dòng)力學(xué)方程已經(jīng)有9人回復(fù)
求助，matlab自定義微分方程擬合實(shí)驗(yàn)數(shù)據(jù)來求方程中的參數(shù) 已經(jīng)有4人回復(fù)
通過實(shí)驗(yàn)數(shù)據(jù)擬合，求解公式中的參數(shù) 已經(jīng)有14人回復(fù)
緊急求助，利用Matlab對(duì)實(shí)驗(yàn)數(shù)據(jù)進(jìn)行擬合求解參數(shù)。已經(jīng)有27人回復(fù)
微分方程組求參數(shù)問題，求高人指教，感謝已經(jīng)有11人回復(fù)
微分方程組參數(shù)擬合的問題（多參數(shù)動(dòng)力系統(tǒng)參數(shù)估計(jì)）已經(jīng)有11人回復(fù)
微分方程組參數(shù)擬合的問題已經(jīng)有11人回復(fù)
關(guān)于matlab微分、及非線性擬合問題已經(jīng)有5人回復(fù)
MATLAB微分方程參數(shù)擬合問題，求大神已經(jīng)有7人回復(fù)
回歸和擬合有啥區(qū)別？已經(jīng)有13人回復(fù)
matlab 擬合反應(yīng)動(dòng)力學(xué)參數(shù)結(jié)果很差。大家?guī)兔匆幌?/a> 已經(jīng)有14人回復(fù)

matlab擬合方程參數(shù)時(shí)初值的選擇已經(jīng)有15人回復(fù)
求解四階微分方程組已經(jīng)有14人回復(fù)
matlab-常微分方程參數(shù)估計(jì) 已經(jīng)有12人回復(fù)
lsqnonlin函數(shù)擬合微分方程組參數(shù)擬合問題已經(jīng)有10人回復(fù)
【求助】多元方程組求解及擬合已經(jīng)有9人回復(fù)
重金加人民幣求擬合wilson方程和NRTL方程參數(shù)的程序已經(jīng)有17人回復(fù)
【求助】常微分方程組的解析解已經(jīng)有3人回復(fù)
【求助】最小二乘法求參數(shù) 已經(jīng)有6人回復(fù)

青春就是用來折騰的！

1樓 2014-08-16 15:18:40

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

月只藍(lán)

主管區(qū)長 (職業(yè)作家)

專家經(jīng)驗(yàn): +1059
計(jì)算強(qiáng)帖: 8
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68127.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時(shí)
蟲號(hào): 1122189
注冊(cè): 2010-10-14
專業(yè): 宇宙學(xué)
管轄: 計(jì)算模擬區(qū)

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

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

內(nèi)容已刪除

贊一下

回復(fù)此樓

MATLAB、MS小問題、普通問題請(qǐng)發(fā)帖求助！時(shí)間精力有限，恕不接受無償私信求助。

2樓2014-08-16 16:59:14

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

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