 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 3209 | 回復(fù): 18 | |||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | |||
小小莓郵鐵蟲 (小有名氣)
|
[求助]
matlab擬合參數(shù)顯示初始點是一個局部最小值 已有5人參與
|
||
|
小弟最近剛學(xué)習(xí)用matlab擬合參數(shù)，但編寫完運行顯示Initial point is a local minimum.得到的k值為k0，求給位大神給指點一下怎么解決這種問題！急求，新人啥也不懂啊！先在此謝過啦！ function zixie clear all; clc; data=[0 8.32 0; 1 7.37 0.211; 2 6.74 0.738; 3 5.92 1.106; 4 5.68 1.296; 5 5.01 1.583; 6 6.09 1.364; 8 5.35 1.84; 10 4.99 2]; x0=[8.32 0]; k0=[10 10 10 10 10 ]; lb=[0 0 0 0 0]; ub=[+inf +inf +inf +inf +inf]; yexp=data(:,2:3); [k,resnorm,residual,exitflag,output,lambda,jacobin]=lsqnonlin(@objFunc1,k0,lb,ub,[],x0,yexp) fprintf('\tk1=%.4\n',k(1)), fprintf('\tk2=%.4\n',k(2)), fprintf('\tk3=%.4\n',k(3)), fprintf('\tk4=%.4\n',k(4)), fprintf('\tk5=%.4\n',k(5)) function f=objFunc1(k,x0,yexp) tspan=[0 1 2 3 4 5 6 8 10]; [t x]=ode45(@funceqs,tspan,x0,[],k); y(:,1)=x(:,1); y(:,2)=x(:,2); f1=y(:,1)-yexp(:,1); f2=y(:,2)-yexp(:,2); f=[f1;f2]; function dxdt=funceqs(t,x,k) dx1dt=-k(1)*k(2)*x(1)*sqrt(k(3)*x(2))/(k(1)+k(2)*x(1))-k(4)*k(5)*x(2); dx2dt=k(1)*k(2)*x(1)*sqrt(k(3)*x(2))/(k(1)+k(2)*x(1))+k(5)*x(2); dxdt=[dx1dt;dx2dt]; Initial point is a local minimum. Optimization completed because the size of the gradient at the initial point is less than the default value of the function tolerance. <stopping criteria details> k = 10 10 10 10 10 Optimization completed: The final point is the initial point. The first-order optimality measure, 0.000000e+00, is less than options.TolFun = 1.000000e-06. Optimization Metric Options relative first-order optimality = 0.00e+00 TolFun = 1e-06 (default) >> |
回復(fù)此樓» 猜你喜歡
請問還有沒有用Latex寫文章的小伙伴們？
已經(jīng)有0人回復(fù)
-
光學(xué)工程學(xué)碩調(diào)劑信息
已經(jīng)有32人回復(fù)
-
物理學(xué)I論文潤色/翻譯怎么收費?
已經(jīng)有245人回復(fù)
-
歡迎加入課題組
已經(jīng)有0人回復(fù)
-
散金幣，求好運，祝面上順利！
已經(jīng)有78人回復(fù)
-
華東師范大學(xué)芯片設(shè)計徐瓏真實水平咋樣
已經(jīng)有1人回復(fù)
-
調(diào)劑
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復(fù)
-
中德博士后2026項目5月30日截止- 國內(nèi)博士后入站、德國兩年全職
已經(jīng)有0人回復(fù)
-
面向全球申請人的中德博士后2026項目5月30日截止- 國內(nèi)博后入站德國兩年全職
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有2人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 拜求！matlab擬合動力學(xué)參數(shù)
已經(jīng)有14人回復(fù)
- 關(guān)于matlab兩個自變量、一個因變量曲線擬合問題
已經(jīng)有3人回復(fù)
- 求助一個最小二乘法matlab擬合程序。
已經(jīng)有7人回復(fù)
- 求助，matlab/ 1stOpt 自定義微分方程擬合實驗數(shù)據(jù)來求方程中的參數(shù)
已經(jīng)有4人回復(fù)
- 用Matlab做多變量非線性擬合求參數(shù)值，求各位大俠鼎力相助！
已經(jīng)有3人回復(fù)
- 怎樣用一組參數(shù)同時擬合兩個曲線--matlab
已經(jīng)有5人回復(fù)
- MATLAB微分方程參數(shù)擬合問題，求大神
已經(jīng)有7人回復(fù)
- Matlab 含有兩個參數(shù)的公式擬合
已經(jīng)有4人回復(fù)
- 用MATLAB怎么打開點云并把它擬合到一個平面上？
已經(jīng)有8人回復(fù)
- 模式搜索MATLAB程序
已經(jīng)有9人回復(fù)
- matlab 擬合反應(yīng)動力學(xué)參數(shù)結(jié)果很差。大家?guī)兔匆幌?/a>
已經(jīng)有14人回復(fù)
- matlab如何編寫共享參數(shù)擬合程序
已經(jīng)有8人回復(fù)
- matlab參數(shù)擬合
已經(jīng)有6人回復(fù)
- matlab微分方程組參數(shù)擬合，以周為單位求解，匯總后以年為單位進行數(shù)值比較
已經(jīng)有4人回復(fù)
- matlab擬合方程參數(shù)時初值的選擇
已經(jīng)有15人回復(fù)
- matlab擬合模型參數(shù)
已經(jīng)有5人回復(fù)
- matlab非線性參數(shù)擬合問題
已經(jīng)有7人回復(fù)
- 【求助】用matlab最優(yōu)化方法進行參數(shù)擬合
已經(jīng)有17人回復(fù)
- 【求助】如何使得matlab擬合的曲線強制經(jīng)過一個點
已經(jīng)有8人回復(fù)
- 【求助】使用MATLAB怎么實現(xiàn)擬合力場參數(shù)的程序？【已解決】
已經(jīng)有15人回復(fù)
小小莓郵
鐵蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 0.2
- 帖子: 52
- 在線: 48.8小時
- 蟲號: 1245029
- 注冊: 2011-03-25
- 性別: GG
- 專業(yè): 化學(xué)反應(yīng)工程
9樓2014-02-27 10:05:29
月只藍
主管區(qū)長 (職業(yè)作家)
-
專家經(jīng)驗: +1059 - 計算強帖: 8
- 應(yīng)助: 1712 (講師)
- 貴賓: 8.888
- 金幣: 68127.7
- 散金: 1938
- 紅花: 443
- 沙發(fā): 4
- 帖子: 4373
- 在線: 3291.4小時
- 蟲號: 1122189
- 注冊: 2010-10-14
- 專業(yè): 宇宙學(xué)
- 管轄: 計算模擬區(qū)
2樓2014-02-26 22:02:02
戴鋼盔的豬頭
木蟲 (知名作家)
- 應(yīng)助: 200 (高中生)
- 金幣: 5083.9
- 散金: 1310
- 紅花: 18
- 帖子: 6001
- 在線: 319.9小時
- 蟲號: 2854040
- 注冊: 2013-12-07
- 性別: GG
- 專業(yè): 波譜分析與成像分析
3樓2014-02-26 22:09:59
dingd
鐵桿木蟲 (職業(yè)作家)
- 計算強帖: 4
- 應(yīng)助: 1641 (講師)
- 金幣: 15037.3
- 散金: 101
- 紅花: 234
- 帖子: 3410
- 在線: 1223.7小時
- 蟲號: 291104
- 注冊: 2006-10-28
【答案】應(yīng)助回帖
★
感謝參與，應(yīng)助指數(shù) +1
fegg7502: 金幣+1, 鼓勵交流 2014-02-27 08:43:32
感謝參與，應(yīng)助指數(shù) +1
fegg7502: 金幣+1, 鼓勵交流 2014-02-27 08:43:32
|
x2的起始值為0（0 8.32 0）有問題吧？會導(dǎo)致計算值全部為0； 還有公式還可以化簡，否則有過擬合現(xiàn)象。 dx1dt=-k1*k2*x1*sqrt(k3*x2)/(k1+k2*x1)-k4*k5*x2; dx2dt=k1*k2*x1*sqrt(k3*x2)/(k1+k2*x1)+k5*x2; 可簡化為： dx1dt=-x1*sqrt(x2)/(k1+k2*x1)-k4*k5*x2; dx2dt=x1*sqrt(x2)/(k1+k2*x1)+k5*x2; |
4樓2014-02-26 23:10:52
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
10