[討論] 請教一下會使用fmincon的高手一點簡單ꨠ…已回收

看板MATLAB作者bawpee (皮)時間15年前 (2011/01/25 11:49)推噓3(3推 0噓 11→)

留言14則, 3人參與討論串1/1

噁死胎偷，新手使用這個指令..不對的話請鞭囧> 我的目標函數： --------------------------------------------------- function f = myfun3_1(x); % 目標函數 a1 = 1; a2 = 1; p = 0.7; people = 40; money = 0.0789; delta_P = p*people*money; w = 0.0002; f = (a1*w*(x(2)^3)*x(3)) + (a2*x(3)*delta_P); ---------------------------------------------------- 然後是我的限制式： function [c,ceq] = mycon3_1(x); % 限制式 V = 10*5*3; c_initial = 1000*10^-6; c_out = 500*10^-6; P = 40; AP= 0.023/3600; G = AP*P; c = []; ceq = x(1)-((c_out+(G/x(2)))*(1-exp((-x(2)/V)*x(3)))-c_initial*exp((-x(2)/V)*x(3))); ----------------------------------------------------- 執行： clear all clc x0 = [1000 0 1]; lb = [600*10^-6 0.3 1]; ub = [800*10^-6 0.5 99999]; [x,fval] = fmincon('myfun3_1',x0,[],[],[],[],lb,ub,'mycon3_1') ----------------------------------------------------- 然後出來的結果為： Warning: Large-scale (trust region) method does not currently solve this type of problem, switching to medium-scale (line search). > In fmincon at 260 In test_CO2 at 6 Optimization terminated: first-order optimality measure less than options.TolFun and maximum constraint violation is less than options.TolCon. Active inequalities (to within options.TolCon = 1e-006): lower upper ineqlin ineqnonlin 1 2 x = 0.0006 0.5000 476.2737 fval = 1.0522e+003 ------------------------------------------------------------------------------ 問題來了..新手小弟我不曉得這樣跑出來的結果到底是不是最佳解..囧因為出現了Warning以及Optimization terminated的字眼，我想請教高手是否知道我的問題在哪邊.. 囧"> 感謝~ -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 120.126.8.117

推

Gwaewluin

01/25 13:05, , 1^F

01/25 13:05, 1^F

嗯，那不好意思再請問一下，改參數是指說我設定的所有參數都有問題而解出來的值怪怪的，還是說只有部分的參數需要更改而已? 另外想請教 x0 的設置是不是也會讓答案有所不一樣? ※ 編輯: bawpee 來自: 120.126.8.117 (01/25 13:41)

推

Gwaewluin

01/25 14:55, , 2^F

01/25 14:55, 2^F

→

Gwaewluin

01/25 14:56, , 3^F