| 5 | 1/1 | 返回列表 |
| 查看: 2194 | 回復: 7 | |||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | |||
pianfeng1986鐵蟲 (初入文壇)
|
[求助]
求教matlab線性優(yōu)化求最小值的問題
|
||
|
這兩天遇到這樣一個問題, 目前遇到一個函數(shù) F(a,b,c,x), 如果定義a,b,c,可以畫出F在x一個區(qū)間上的曲線,但無法求得F的具體解析式。 我的目的是把F曲線與手頭上的實驗曲線達到最好吻合,同時求出相應(yīng)的參數(shù)a,b,c。 目前我的思路是這樣的:試驗曲線有10000組數(shù)據(jù) (x,y),我用最小二乘的思想,求出每一個x上對應(yīng)的F與y之差平方和,即sum((y-F)^2), 然后給出a, b , c的范圍,用fmincon來求解,但是寫了程序,發(fā)現(xiàn)問題,程序無法求出相應(yīng)解。錯誤為 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ù)用錯了嗎? 還是用別的方法? 主程序如下: 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); |
鐵蟲 (初入文壇)
榮譽版主 (著名寫手)
彩色掛圖
![]() |
專家經(jīng)驗: +2 |

金蟲 (小有名氣)
鐵蟲 (初入文壇)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 081200-11408-276學碩求調(diào)劑 +4 | 崔wj 2026-03-26 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 求化學調(diào)劑 +4 | wulanna 2026-03-28 | 4/200 |
|
|
[考研] 0856求調(diào)劑 +11 | zhn03 2026-03-25 | 12/600 |
|
|
[考研] 一志愿華理,數(shù)一英一285求A區(qū)調(diào)劑 +8 | AZMK 2026-03-25 | 11/550 |
|
|
[考研] 求調(diào)劑 +6 | 蘆lty 2026-03-25 | 7/350 |
|
|
[考研] 299求調(diào)劑 +7 | 嗯嗯嗯嗯2 2026-03-27 | 7/350 |
|
|
[考研] 289求調(diào)劑 +10 | 新時代材料 2026-03-27 | 10/500 |
|
|
[考研]
|
孅華 2026-03-22 | 17/850 |
|
|
[材料工程] 一志愿C9材料與化工專業(yè)總分300求調(diào)劑 +8 | 曼111 2026-03-24 | 9/450 |
|
|
[考研] 311求調(diào)劑 +3 | 希望上岸阿小楊 2026-03-23 | 3/150 |
|
|
[考研] 315分求調(diào)劑 +7 | 26考研上岸版26 2026-03-26 | 7/350 |
|
|
[考研] 085600,材料與化工321分,求調(diào)劑 +9 | 大饞小子 2026-03-27 | 9/450 |
|
|
[考研] 中國科學院深圳先進技術(shù)研究院-光纖傳感課題組招生-中國科學院大學、深圳理工大學聯(lián)培 +5 | YangTyu1 2026-03-26 | 5/250 |
|
|
[考研] 環(huán)境專碩324分求調(diào)劑推薦 +5 | 軒小寧—— 2026-03-26 | 5/250 |
|
|
[考研] 【2026考研調(diào)劑】制藥工程 284分 求相關(guān)專業(yè)調(diào)劑名額 +4 | 袁奐奐 2026-03-25 | 8/400 |
|
|
[考研] 292求調(diào)劑 +4 | 鵝鵝鵝額額額額?/a> 2026-03-24 | 4/200 |
|
|
[考研] 305分求調(diào)劑(食品工程) +5 | Sxy112 2026-03-21 | 7/350 |
|
|
[考研]
|
13659058978 2026-03-24 | 4/200 |
|
|
[考研] 環(huán)境學碩288求調(diào)劑 +8 | 皮皮皮123456 2026-03-22 | 8/400 |
|
|
[考研] 315分,誠求調(diào)劑,材料與化工085600 +3 | 13756423260 2026-03-22 | 3/150 |
|