| 8 | 1/1 | 返回列表 |
| 查看: 2179 | 回復(fù): 7 | |||
pianfeng1986鐵蟲 (初入文壇)
|
[求助]
求教matlab線性優(yōu)化求最小值的問題
|
|
這兩天遇到這樣一個(gè)問題, 目前遇到一個(gè)函數(shù) F(a,b,c,x), 如果定義a,b,c,可以畫出F在x一個(gè)區(qū)間上的曲線,但無(wú)法求得F的具體解析式。 我的目的是把F曲線與手頭上的實(shí)驗(yàn)曲線達(dá)到最好吻合,同時(shí)求出相應(yīng)的參數(shù)a,b,c。 目前我的思路是這樣的:試驗(yàn)曲線有10000組數(shù)據(jù) (x,y),我用最小二乘的思想,求出每一個(gè)x上對(duì)應(yīng)的F與y之差平方和,即sum((y-F)^2), 然后給出a, b , c的范圍,用fmincon來(lái)求解,但是寫了程序,發(fā)現(xiàn)問題,程序無(wú)法求出相應(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. 請(qǐng)問大神我該怎么辦,是函數(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); |
榮譽(yù)版主 (著名寫手)
彩色掛圖
![]() |
專家經(jīng)驗(yàn): +2 |

金蟲 (小有名氣)
鐵蟲 (初入文壇)
鐵蟲 (初入文壇)
金蟲 (小有名氣)
鐵蟲 (初入文壇)
金蟲 (小有名氣)
| 8 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料專碩306英一數(shù)二 +5 | z1z2z3879 2026-03-16 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6 | 薛云鵬 2026-03-15 | 6/300 |
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學(xué)123 2026-03-15 | 8/400 |
|
|
[考研] 268求調(diào)劑 +5 | 好運(yùn)連綿不絕 2026-03-12 | 5/250 |
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +7 | 邱gl 2026-03-11 | 8/400 |
|
|
[考研] 274求調(diào)劑 +5 | 時(shí)間點(diǎn) 2026-03-13 | 5/250 |
|
|
[考研] 材料與化工一志愿南昌大學(xué)327求調(diào)劑推薦 +7 | Ncdx123456 2026-03-13 | 8/400 |
|
|
[教師之家] 焦慮 +7 | 水冰月月野兔 2026-03-13 | 9/450 |
|
|
[考研] 材料工程327求調(diào)劑 +3 | xiaohe12w 2026-03-11 | 3/150 |
|
|
[基金申請(qǐng)] 現(xiàn)在如何回避去年的某一個(gè)專家,不知道名字 +3 | zk200107 2026-03-12 | 6/300 |
|
|
[考研] 復(fù)試調(diào)劑 +4 | z1z2z3879 2026-03-14 | 5/250 |
|
|
[考研] 一志愿中科院,化學(xué)方向,295求調(diào)劑 +4 | 一氧二氮 2026-03-11 | 4/200 |
|
|
[考研] 0703化學(xué)一志愿211 總分320求調(diào)劑 +5 | 瑪卡巴卡啊哈 2026-03-11 | 5/250 |
|
|
[考研] 311求調(diào)劑 +3 | 冬十三 2026-03-13 | 3/150 |
|
|
[考研] 求調(diào)劑 +5 | 一定有學(xué)上- 2026-03-12 | 5/250 |
|
|
[考研] 材料301分求調(diào)劑 +5 | Liyouyumairs 2026-03-12 | 5/250 |
|
|
[考研] 0856化學(xué)工程280分求調(diào)劑 +4 | shenzxsn 2026-03-11 | 4/200 |
|
|
[考研] 08食品或輕工求調(diào)劑,本科發(fā)表3篇sci一區(qū)top論文,一志愿南師大食品科學(xué)與工程 +3 | 我是一個(gè)兵, 2026-03-10 | 3/150 |
|
|
[考博] 26讀博 +4 | Rui135246 2026-03-12 | 10/500 |
|
|
[考研] 333求調(diào)劑 +3 | 152697 2026-03-12 | 4/200 |
|