| 5 | 1/1 | 返回列表 |
| 查看: 2195 | 回復(fù): 7 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
pianfeng1986鐵蟲 (初入文壇)
|
[求助]
求教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上對(duì)應(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. 請(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 |

金蟲 (小有名氣)
鐵蟲 (初入文壇)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0703 化學(xué) 求調(diào)劑,一志愿山東大學(xué) 342 分 +7 | Shern—- 2026-03-28 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 求調(diào)劑323材料與化工 +10 | 1124361 2026-03-24 | 10/500 |
|
|
[考研] 307求調(diào)劑 +15 | 超級(jí)伊昂大王 2026-03-24 | 16/800 |
|
|
[考研] 303求調(diào)劑 +7 | DLkz1314. 2026-03-30 | 7/350 |
|
|
[考研] 284求調(diào)劑 +14 | junqihahaha 2026-03-26 | 15/750 |
|
|
[考研] 070300一志愿211,312分求調(diào)劑院校 +4 | 小黃鴨寶 2026-03-30 | 4/200 |
|
|
[考研] 一志愿211,335分,0856,求調(diào)劑院校和導(dǎo)師 +7 | 傾____蕭 2026-03-27 | 8/400 |
|
|
[考研] 求調(diào)劑 +10 | 張zz111 2026-03-27 | 11/550 |
|
|
[考研] 327求調(diào)劑 +4 | 小卡不卡. 2026-03-29 | 4/200 |
|
|
[考研] 各位老師好,我的一志愿為北京科技大學(xué)085601材料專碩 +9 | Koxui 2026-03-28 | 9/450 |
|
|
[考研] 343求調(diào)劑085601 +3 | 要努力學(xué)習(xí)x 2026-03-29 | 3/150 |
|
|
[考研] 299求調(diào)劑 +10 | 15188958825 2026-03-25 | 10/500 |
|
|
[考研] 291求調(diào)劑 +5 | Y-cap 2026-03-29 | 6/300 |
|
|
[考研] 316求調(diào)劑 +7 | 江辭666 2026-03-26 | 7/350 |
|
|
[考研] 壓國(guó)家一區(qū)線,求導(dǎo)師收留,有恩必謝! +7 | 迷人的哈哈 2026-03-28 | 7/350 |
|
|
[考研] 299求調(diào)劑 +7 | 嗯嗯嗯嗯2 2026-03-27 | 7/350 |
|
|
[考研] 272求調(diào)劑 +7 | 腳滑的守法公民 2026-03-27 | 7/350 |
|
|
[考研] 求調(diào)劑 一志愿 本科 北科大 化學(xué) 343 +6 | 13831862839 2026-03-24 | 7/350 |
|
|
[考研] 289求調(diào)劑 +17 | 碩星赴 2026-03-23 | 17/850 |
|
|
[考研] 求調(diào)劑 +3 | 李李不服輸 2026-03-25 | 3/150 |
|