 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 3229 | 回復(fù): 18 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | |||
小小莓郵鐵蟲 (小有名氣)
|
[求助]
matlab擬合參數(shù)顯示初始點是一個局部最小值 已有5人參與
|
||
|
小弟最近剛學(xué)習(xí)用matlab擬合參數(shù)，但編寫完運行顯示Initial point is a local minimum.得到的k值為k0，求給位大神給指點一下怎么解決這種問題！急求，新人啥也不懂�。∠仍诖酥x過啦！ 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ù)此樓» 猜你喜歡
重慶交大材料學(xué)院26年碩士招生復(fù)試實施細(xì)則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。
已經(jīng)有0人回復(fù)
-
北工大宋曉艷教授團(tuán)隊招收2026年9月入學(xué)博士【學(xué)博和專博名額均有】！
已經(jīng)有1人回復(fù)
-
物理學(xué)I論文潤色/翻譯怎么收費?
已經(jīng)有125人回復(fù)
-
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收10人
已經(jīng)有0人回復(fù)
-
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收調(diào)劑，10人
已經(jīng)有0人回復(fù)
-
計算與實驗科學(xué)工程國際會議(8月7-11日,香港) 二維材料分會S5與聲學(xué)換能器分會S8
已經(jīng)有0人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑有論文
已經(jīng)有4人回復(fù)
-
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來能源學(xué)院
已經(jīng)有28人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑 本科期間發(fā)表一篇Sci
已經(jīng)有4人回復(fù)
-
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來能源學(xué)院
已經(jīng)有21人回復(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ù)擬合，以周為單位求解，匯總后以年為單位進(jìn)行數(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)化方法進(jìn)行參數(shù)擬合
已經(jīng)有17人回復(fù)
- 【求助】如何使得matlab擬合的曲線強(qiáng)制經(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)工程
7樓2014-02-27 09:49:39
月只藍(lán)
主管區(qū)長 (職業(yè)作家)
-
專家經(jīng)驗: +1059 - 計算強(qiáng)帖: 8
- 應(yīng)助: 1712 (講師)
- 貴賓: 8.888
- 金幣: 68128.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è)作家)
- 計算強(qiáng)帖: 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ā)表 | |
|---|---|---|---|---|
|
[考研] 324分 085600材料與化工 +6 | 呆鵝oor 2026-03-27 | 6/300 |
|
|---|---|---|---|---|
|
|
[考研] 254材料與化工求調(diào)劑 +3 | 翰冬林楠 2026-03-30 | 4/200 |
|
|
|
[考研] 318求調(diào)劑 +10 | 陳晨79 2026-03-30 | 10/500 |
|
|
|
[考研] 生物學(xué)296求調(diào)劑 +8 | 湯圓包 2026-03-29 | 12/600 |
|
|
[考研] 機(jī)械學(xué)碩總分317求調(diào)劑�。。�！ +6 | Acaciad 2026-03-25 | 6/300 |
|
|
|
[考研] 322求調(diào)劑：一志愿湖南大學(xué) 材料與化工（085600），已過六級。 +10 | XX小鄧 2026-03-29 | 10/500 |
|
|
|
[考研] 一志愿中海洋材料357 +4 | 麥恩莉. 2026-03-30 | 4/200 |
|
|
|
[考研] 總分322求生物學(xué)/生化與分子/生物信息學(xué)相關(guān)調(diào)劑 +6 | 星沉uu 2026-03-26 | 7/350 |
|
|
|
[考研] 274求調(diào)劑 +6 | xiao愛同學(xué) 2026-03-30 | 6/300 |
|
|
|
[考研] 0817化工學(xué)碩調(diào)劑 +7 | 努力上岸中！ 2026-03-31 | 7/350 |
|
|
|
[考研]
 一志愿大連理工大學(xué)材料求調(diào)劑
+6
|
Gymno 2026-03-30 | 6/300 |
|
|
|
[考研] 11408總分309，一志愿東南大學(xué)求調(diào)劑，不挑專業(yè) +5 | 天賦帶到THU 2026-03-29 | 6/300 |
|
|
|
[考研] 085601材料工程找調(diào)劑 +17 | oatmealR 2026-03-29 | 18/900 |
|
|
|
[考研] 290求調(diào)劑 +3 | dfffsar 2026-03-29 | 3/150 |
|
|
|
[考研] 285求調(diào)劑 +4 | AZMK 2026-03-27 | 7/350 |
|
|
|
[考研] 315調(diào)劑 +4 | 0860求調(diào)劑 2026-03-26 | 5/250 |
|
|
|
[考研] 材料求調(diào)劑 +5 | .m.. 2026-03-25 | 5/250 |
|
|
|
[考研] 調(diào)劑求收留 +7 | 果然有我 2026-03-26 | 7/350 |
|
|
|
[考研] 中國科學(xué)院深圳先進(jìn)技術(shù)研究院-光纖傳感課題組招生-中國科學(xué)院大學(xué)、深圳理工大學(xué)聯(lián)培 +5 | YangTyu1 2026-03-26 | 5/250 |
|
|
|
[考研] 281求調(diào)劑 +6 | Koxui 2026-03-24 | 7/350 |
|