 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 3221 | 回復(fù): 18 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||
小小莓郵鐵蟲 (小有名氣)
|
[求助]
matlab擬合參數(shù)顯示初始點(diǎn)是一個(gè)局部最小值 已有5人參與
|
||
|
小弟最近剛學(xué)習(xí)用matlab擬合參數(shù)，但編寫完運(yùn)行顯示Initial point is a local minimum.得到的k值為k0，求給位大神給指點(diǎn)一下怎么解決這種問題！急求，新人啥也不懂啊！先在此謝過啦！ 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ù)此樓» 猜你喜歡
調(diào)劑
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復(fù)
-
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有163人回復(fù)
-
中德博士后2026項(xiàng)目5月30日截止- 國(guó)內(nèi)博士后入站、德國(guó)兩年全職
已經(jīng)有0人回復(fù)
-
面向全球申請(qǐng)人的中德博士后2026項(xiàng)目5月30日截止- 國(guó)內(nèi)博后入站德國(guó)兩年全職
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有2人回復(fù)
-
重慶交大材料學(xué)院26年碩士招生復(fù)試實(shí)施細(xì)則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。
已經(jīng)有0人回復(fù)
-
北工大宋曉艷教授團(tuán)隊(duì)招收2026年9月入學(xué)博士【學(xué)博和專博名額均有】！
已經(jīng)有1人回復(fù)
-
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收10人
已經(jīng)有0人回復(fù)
-
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收調(diào)劑，10人
已經(jīng)有0人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 拜求！matlab擬合動(dòng)力學(xué)參數(shù)
已經(jīng)有14人回復(fù)
- 關(guān)于matlab兩個(gè)自變量、一個(gè)因變量曲線擬合問題
已經(jīng)有3人回復(fù)
- 求助一個(gè)最小二乘法matlab擬合程序。
已經(jīng)有7人回復(fù)
- 求助，matlab/ 1stOpt 自定義微分方程擬合實(shí)驗(yàn)數(shù)據(jù)來求方程中的參數(shù)
已經(jīng)有4人回復(fù)
- 用Matlab做多變量非線性擬合求參數(shù)值，求各位大俠鼎力相助！
已經(jīng)有3人回復(fù)
- 怎樣用一組參數(shù)同時(shí)擬合兩個(gè)曲線--matlab
已經(jīng)有5人回復(fù)
- MATLAB微分方程參數(shù)擬合問題，求大神
已經(jīng)有7人回復(fù)
- Matlab 含有兩個(gè)參數(shù)的公式擬合
已經(jīng)有4人回復(fù)
- 用MATLAB怎么打開點(diǎn)云并把它擬合到一個(gè)平面上？
已經(jīng)有8人回復(fù)
- 模式搜索MATLAB程序
已經(jīng)有9人回復(fù)
- matlab 擬合反應(yīng)動(dò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ù)時(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)過一個(gè)點(diǎn)
已經(jīng)有8人回復(fù)
- 【求助】使用MATLAB怎么實(shí)現(xiàn)擬合力場(chǎng)參數(shù)的程序？【已解決】
已經(jīng)有15人回復(fù)
戴鋼盔的豬頭
木蟲 (知名作家)
- 應(yīng)助: 200 (高中生)
- 金幣: 5083.9
- 散金: 1310
- 紅花: 18
- 帖子: 6001
- 在線: 319.9小時(shí)
- 蟲號(hào): 2854040
- 注冊(cè): 2013-12-07
- 性別: GG
- 專業(yè): 波譜分析與成像分析
3樓2014-02-26 22:09:59
月只藍(lán)
主管區(qū)長(zhǎng) (職業(yè)作家)
-
專家經(jīng)驗(yàn): +1059 - 計(jì)算強(qiáng)帖: 8
- 應(yīng)助: 1712 (講師)
- 貴賓: 8.888
- 金幣: 68128.7
- 散金: 1938
- 紅花: 443
- 沙發(fā): 4
- 帖子: 4373
- 在線: 3291.4小時(shí)
- 蟲號(hào): 1122189
- 注冊(cè): 2010-10-14
- 專業(yè): 宇宙學(xué)
- 管轄: 計(jì)算模擬區(qū)
2樓2014-02-26 22:02:02
dingd
鐵桿木蟲 (職業(yè)作家)
- 計(jì)算強(qiáng)帖: 4
- 應(yīng)助: 1641 (講師)
- 金幣: 15037.3
- 散金: 101
- 紅花: 234
- 帖子: 3410
- 在線: 1223.7小時(shí)
- 蟲號(hào): 291104
- 注冊(cè): 2006-10-28
【答案】應(yīng)助回帖
★
感謝參與，應(yīng)助指數(shù) +1
fegg7502: 金幣+1, 鼓勵(lì)交流 2014-02-27 08:43:32
感謝參與，應(yīng)助指數(shù) +1
fegg7502: 金幣+1, 鼓勵(lì)交流 2014-02-27 08:43:32
|
x2的起始值為0（0 8.32 0）有問題吧？會(huì)導(dǎo)致計(jì)算值全部為0； 還有公式還可以化簡(jiǎn)，否則有過擬合現(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; 可簡(jiǎn)化為： 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
pdl9527
專家顧問 (小有名氣)
-
專家經(jīng)驗(yàn): +8 - 應(yīng)助: 100 (初中生)
- 金幣: 2110.3
- 紅花: 23
- 帖子: 282
- 在線: 112小時(shí)
- 蟲號(hào): 1227333
- 注冊(cè): 2011-03-09
- 性別: GG
- 專業(yè): 分離過程
- 管轄: 計(jì)算模擬
5樓2014-02-27 00:00:32
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 266分，求材料冶金能源化工等調(diào)劑 +6 | 哇呼哼呼哼 2026-03-27 | 8/400 |
|
|---|---|---|---|---|
|
|
[考研] 求調(diào)劑 +8 | 張zz111 2026-03-27 | 9/450 |
|
|
[考研] 275求調(diào)劑 +10 | Micky11223 2026-03-25 | 13/650 |
|
|
|
[考研] 一志愿南師大0703化學(xué) 275求調(diào)劑 +4 | Ripcord上岸 2026-03-27 | 4/200 |
|
|
|
[考研] 308求調(diào)劑 +7 | 墨墨漠 2026-03-25 | 7/350 |
|
|
|
[考研] 322求調(diào)劑 +4 | 我真的很想學(xué)習(xí) 2026-03-23 | 4/200 |
|
|
|
[考研] 298調(diào)劑 +3 | jiyingjie123 2026-03-27 | 3/150 |
|
|
|
[考研] 276求調(diào)劑。有半年電池和半年高分子實(shí)習(xí)經(jīng)歷 +10 | 材料學(xué)257求調(diào)劑 2026-03-23 | 11/550 |
|
|
|
[考研]
材料學(xué)碩，求調(diào)劑
 6+5
|
糖葫蘆888ll 2026-03-22 | 10/500 |
|
|
|
[考研] 一志愿北化求調(diào)劑 +3 | Jsman 2026-03-22 | 3/150 |
|
|
|
[考研] 271求調(diào)劑 +6 | 生如夏花… 2026-03-22 | 6/300 |
|
|
|
[考研]
 085602 289分求調(diào)劑
+8
|
WWW西西弗斯 2026-03-24 | 8/400 |
|
|
|
[考研] 化學(xué)調(diào)劑一志愿上海交通大學(xué)336分-本科上海211 +4 | 小魚愛有機(jī) 2026-03-25 | 4/200 |
|
|
|
[考研] 334分 一志愿武理-080500 材料求調(diào)劑 +4 | 李李不服輸 2026-03-25 | 4/200 |
|
|
|
[考研] 求調(diào)劑 +3 | 李李不服輸 2026-03-25 | 3/150 |
|
|
|
[考研]
材料考研調(diào)劑生
+3
|
黃粱一夢(mèng)千年 2026-03-24 | 3/150 |
|
|
|
[考研] 341求調(diào)劑(一志愿湖南大學(xué)070300) +5 | 番茄頭--- 2026-03-22 | 6/300 |
|
|
|
[考研] 求調(diào)劑一志愿海大，0703化學(xué)學(xué)碩304分，有大創(chuàng)項(xiàng)目，四級(jí)已過 +6 | 幸運(yùn)哩哩 2026-03-22 | 10/500 |
|
|
|
[考研] 求調(diào)劑 +5 | Zhangbod 2026-03-21 | 7/350 |
|
|
|
[考研] 材料與化工（0856）304求 B區(qū) 調(diào)劑 +3 | 邱gl 2026-03-21 | 3/150 |
|