pianfeng1986

鐵蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 60.9
帖子: 25
在線: 10.5小時(shí)
蟲號: 1136766
注冊: 2010-11-01
性別: GG
專業(yè): 聲學(xué)

[求助] 求教matlab線性優(yōu)化求最小值的問題

這兩天遇到這樣一個(gè)問題，目前遇到一個(gè)函數(shù) F(a,b,c,x)，如果定義a,b,c，可以畫出F在x一個(gè)區(qū)間上的曲線，但無法求得F的具體解析式。  我的目的是把F曲線與手頭上的實(shí)驗(yàn)曲線達(dá)到最好吻合，同時(shí)求出相應(yīng)的參數(shù)a,b,c。

   目前我的思路是這樣的：試驗(yàn)曲線有10000組數(shù)據(jù) （x,y)，我用最小二乘的思想，求出每一個(gè)x上對應(yīng)的F與y之差平方和，即sum((y-F)^2),  然后給出a, b , c的范圍，用fmincon來求解，但是寫了程序，發(fā)現(xiàn)問題，程序無法求出相應(yīng)解。錯(cuò)誤為  Warning: Large-scale (trust region) method does not currently solve this type of problem,
switching to medium-scale (line search).
> In fmincon at 274
  In main at 19
Optimization terminated: magnitude of directional derivative in search
direction less than 2*options.TolFun and maximum constraint violation
  is less than options.TolCon.
No active inequalities.

請問大神我該怎么辦，是函數(shù)用錯(cuò)了嗎？還是用別的方法？

主程序如下：
clear all;clc;
global nr nt apt pidt vb freq nfrq

ui=sqrt(-1);
nr=50;
nt=51;
apt=100;
pidt=10d0;
vb=4212;

%% A，b定義約束    45 < x1<55;  0.22 < x2<0.28;  0.0005 < x3<0.001;
A=[1,0,0; 0,1,0;0,0,1;-1,0,0;0,-1,0;0,0,-1];
b=[55,0.28,0.001,-45,-0.22,-0.0005]';
x0= [49.13; 0.2592; 0.0007]; % Starting guess at the solution
[x,fval] = fmincon(@myfun,x0,A,b)
%%_______________________________________
function f = myfun(x);
global nr nt apt pidt vb freq nfrq
adf=load('y_para.dat');

fs=adf(:,1);
con=adf(:,2);
freq=fs;
nfrq=length(freq);

epson= 0.3302  ;
ita=0.6052  ;
c=1.0888  ;
ym=admres(x(1),x(2),epson,ita,c,x(3));
f=sum((real(ym)-con).^2);

回復(fù)此樓

» 猜你喜歡

2026年循環(huán)經(jīng)濟(jì)功能材料國際會議（ICFMCE 2026）已經(jīng)有0人回復(fù)
2026年第五屆電氣、電子與信息工程國際會議（ISEEIE 2026）已經(jīng)有0人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費(fèi)? 已經(jīng)有120人回復(fù)
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑已經(jīng)有5人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有5人回復(fù)
求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
請問還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有25人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運(yùn)，祝面上順利！已經(jīng)有29人回復(fù)

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

matlab非線性擬合求助已經(jīng)有7人回復(fù)
用Matlab做多變量非線性擬合求參數(shù)值，求各位大俠鼎力相助！已經(jīng)有3人回復(fù)
MATLAB求解非線性方程組已經(jīng)有5人回復(fù)
求Matlab解非線性方程代碼已經(jīng)有12人回復(fù)
matlab無約束非線性方程求最小值的問題已經(jīng)有7人回復(fù)
關(guān)于matlab微分、及非線性擬合問題已經(jīng)有5人回復(fù)
matlab解非線性方程組解不開已經(jīng)有3人回復(fù)
求Matlab高手解決線性方程組的迭代求解問題已經(jīng)有17人回復(fù)
matlab求解非線性方程組，錯(cuò)誤提示怎么解決已經(jīng)有5人回復(fù)
求助用matlab擬合三元線性回歸方程及分析已經(jīng)有9人回復(fù)
求助matlab---fsolve解非線性方程組已經(jīng)有6人回復(fù)
關(guān)于matlab線性規(guī)劃背包問題的解法求救呀~~方程知道但不知道如何解答···· 已經(jīng)有4人回復(fù)
1stopt或matlab如何求解以下的非線性方程并擬合出相應(yīng)曲線？已經(jīng)有13人回復(fù)
MATLAB非線性優(yōu)化擬合怎么改才正確已經(jīng)有3人回復(fù)
求助！matlab用fsolve函數(shù)求解非線性方程組的問題！已經(jīng)有19人回復(fù)
matlab線性矩陣不等式求解做優(yōu)化問題已經(jīng)有10人回復(fù)
matlab求解非線性方程組已經(jīng)有16人回復(fù)
求回歸方程參數(shù)估計(jì)的最優(yōu)化算法matlab代碼已經(jīng)有9人回復(fù)
求助有關(guān)Matlab有約束非線性最優(yōu)化問題已經(jīng)有9人回復(fù)
求高人指點(diǎn)用matlab求解非線性方程組，解決了追加100金幣；已經(jīng)有11人回復(fù)
matlab非線性參數(shù)擬合問題已經(jīng)有7人回復(fù)
matlab的fsove 命令求解非線性方程組已經(jīng)有6人回復(fù)
【求助】求教matlab解非線性方程組已經(jīng)有9人回復(fù)

1樓 2013-05-12 11:57:04

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

csgt0

榮譽(yù)版主 (著名寫手)

彩色掛圖

專家經(jīng)驗(yàn): +2
計(jì)算強(qiáng)帖: 2
應(yīng)助: 367 (碩士)
貴賓: 0.457
金幣: 7807
散金: 8897
紅花: 20
沙發(fā): 7
帖子: 2809
在線: 515.8小時(shí)
蟲號: 449046
注冊: 2007-11-02
專業(yè): 金屬非晶態(tài)、準(zhǔn)晶和納米晶
管轄: 計(jì)算模擬

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

★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
pianfeng1986: 金幣+5 2013-05-16 15:07:20

用lsqcurvefit試試

贊一下

回復(fù)此樓

showmethemoney

2樓2013-05-14 15:54:36

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

xxz903

金蟲 (小有名氣)

應(yīng)助: 33 (小學(xué)生)
金幣: 1045.6
散金: 190
紅花: 4
帖子: 253
在線: 66.7小時(shí)
蟲號: 1360336
注冊: 2011-08-03
性別: GG
專業(yè): 化工系統(tǒng)工程

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

★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
csgt0: 金幣+1, 謝謝 2013-05-14 16:19:41
pianfeng1986: 金幣+5, ★有幫助 2013-05-16 15:07:06

線性問題直接用origin就行啊~matlab優(yōu)化最小誤差也行~~~不過優(yōu)化擬合最好用的還是1stopt~有1.5破解的絕對能解決你的問題

贊一下

回復(fù)此樓

3樓2013-05-14 16:01:07

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

pianfeng1986

鐵蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 60.9
帖子: 25
在線: 10.5小時(shí)
蟲號: 1136766
注冊: 2010-11-01
性別: GG
專業(yè): 聲學(xué)

引用回帖:

3樓: Originally posted by xxz903 at 2013-05-14 16:01:07
線性問題直接用origin就行啊~matlab優(yōu)化最小誤差也行~~~不過優(yōu)化擬合最好用的還是1stopt~有1.5破解的絕對能解決你的問題

我所求的函數(shù)木有具體的解析式，1stopt不好用啊

贊一下

回復(fù)此樓

4樓2013-05-14 21:35:05

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

pianfeng1986

鐵蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 60.9
帖子: 25
在線: 10.5小時(shí)
蟲號: 1136766
注冊: 2010-11-01
性別: GG
專業(yè): 聲學(xué)

引用回帖:

2樓: Originally posted by csgt0 at 2013-05-14 15:54:36
用lsqcurvefit試試

謝謝我研究下

回復(fù)此樓

5樓2013-05-14 21:36:51

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

xxz903

金蟲 (小有名氣)

應(yīng)助: 33 (小學(xué)生)
金幣: 1045.6
散金: 190
紅花: 4
帖子: 253
在線: 66.7小時(shí)
蟲號: 1360336
注冊: 2011-08-03
性別: GG
專業(yè): 化工系統(tǒng)工程

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

1stopt可以自動(dòng)尋找公式的~~而且擬合度還蠻高的~

贊一下

回復(fù)此樓

6樓2013-05-14 21:41:15

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

pianfeng1986

鐵蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 60.9
帖子: 25
在線: 10.5小時(shí)
蟲號: 1136766
注冊: 2010-11-01
性別: GG
專業(yè): 聲學(xué)

引用回帖:

6樓: Originally posted by xxz903 at 2013-05-14 21:41:15
1stopt可以自動(dòng)尋找公式的~~而且擬合度還蠻高的~

我這個(gè)其實(shí)非擬合問題了，優(yōu)化求最小值的問題~~ 不知道您想法是如何的？

贊一下

回復(fù)此樓

7樓2013-05-16 15:06:56

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

xxz903

金蟲 (小有名氣)

應(yīng)助: 33 (小學(xué)生)
金幣: 1045.6
散金: 190
紅花: 4
帖子: 253
在線: 66.7小時(shí)
蟲號: 1360336
注冊: 2011-08-03
性別: GG
專業(yè): 化工系統(tǒng)工程

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

★
csgt0: 金幣+1, 謝謝 2013-05-16 17:14:33

最小值問題可以使用matlab的優(yōu)化工具箱，也可以使用1stopt的優(yōu)化功能（min 或者max）都可以的問題不大

贊一下

回復(fù)此樓

8樓2013-05-16 16:10:41

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

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 材料，紡織，生物（0856、0710），化學(xué)招生啦 +3	Eember. 2026-03-17	7/350	2026-03-17 20:20 by 花125533
[考研] 材料與化工求調(diào)劑 +6	為學(xué)666 2026-03-16	6/300	2026-03-17 20:15 by peike
[考研] 344求調(diào)劑 +4	knight344 2026-03-16	4/200	2026-03-17 17:27 by ruiyingmiao
[考研] 326求調(diào)劑 +5	上岸的小葡 2026-03-15	6/300	2026-03-17 17:26 by ruiyingmiao
[考研] 085601求調(diào)劑 +4	Du.11 2026-03-16	4/200	2026-03-17 17:08 by ruiyingmiao
[考研] 286求調(diào)劑 +3	lemonzzn 2026-03-16	5/250	2026-03-16 20:43 by lemonzzn
[文學(xué)芳草園] 伙伴們，祝我生日快樂吧 +17	myrtle 2026-03-10	26/1300	2026-03-16 18:32 by 青橙Ln
[考研] 0854控制工程 359求調(diào)劑可跨專業(yè) +3	626776879 2026-03-14	9/450	2026-03-16 17:42 by 626776879
[考研] 環(huán)境工程調(diào)劑 +6	大可digkids 2026-03-16	6/300	2026-03-16 17:16 by barlinike
[考研] 0703化學(xué)調(diào)劑 290分有科研經(jīng)歷，論文在投 +7	膩膩gk 2026-03-14	7/350	2026-03-16 10:12 by houyaoxu
[考研] 本科南京大學(xué)一志愿川大藥學(xué)327 +3	麥田耕者 2026-03-14	3/150	2026-03-14 20:04 by 外星文明
[考研] 材料080500調(diào)劑求收留 +3	一顆meteor 2026-03-13	3/150	2026-03-14 10:54 by peike
[考研] 四川大學(xué)085601材料工程專碩初試294求調(diào)劑 +4	祝我們好在冬天 2026-03-11	4/200	2026-03-13 21:39 by peike
[考研] 315求調(diào)劑 +9	小羊小羊_ 2026-03-11	10/500	2026-03-13 21:13 by SXNU李老師
[考研] 329求調(diào)劑 +3	miaodesi 2026-03-12	4/200	2026-03-13 20:53 by 18595523086
[考研] 281求調(diào)劑 +9	Koxui 2026-03-12	11/550	2026-03-13 20:50 by Koxui
[考研] 求調(diào)劑 +5	一定有學(xué)上- 2026-03-12	5/250	2026-03-13 18:31 by ms629
[考研] 0817化學(xué)工程與技術(shù)考研312分調(diào)劑 +3	T123 tt 2026-03-12	3/150	2026-03-13 10:49 by houyaoxu
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3	崔wj 2026-03-12	4/200	2026-03-12 19:33 by 求調(diào)劑zz
[考博] 2026年博士申請 +3	QwQwQW10 2026-03-11	3/150	2026-03-12 17:58 by gxch43