MATLAB二次规划和整数规划
二次規(guī)劃
若某非線性規(guī)劃的目標(biāo)函數(shù)為自變量 x 的二次函數(shù),約束條件又全是線性的,就稱這種規(guī)劃為二次規(guī)劃。
min12xTHx+fTxs.t.Ax≤bmin \frac{1}{2}x^THx+f^Tx\\ s.t.\quad Ax \leq bmin21?xTHx+fTxs.t.Ax≤b
quadprog函數(shù)
[x,fval]=quadprog(H,f,A,b,Aeq,beq,LB,UB,X0)[x,fval]=quadprog(H,f,A,b,Aeq,beq,LB,UB,X0)[x,fval]=quadprog(H,f,A,b,Aeq,beq,LB,UB,X0)
- H、f 即目標(biāo)函數(shù)中的 H、f
- 其他參數(shù)與線性規(guī)劃相同,參見(jiàn)線性規(guī)劃
- 二次規(guī)劃所求的結(jié)果為精確最優(yōu)解
實(shí)例
minf(x)=2x12?4x1x2+4x22?6x1?3x2s.t.{x1+x2≤34x1+x2≤9x1,x2≥0minf(x)=2x_1^2-4x_1x_2+4x_2^2-6x_1-3x_2\\ s.t.\begin{cases} x_1+x_2 \leq3\\ 4x_1+x_2\leq 9\\ x_1,x_2\geq 0 \end{cases}minf(x)=2x12??4x1?x2?+4x22??6x1??3x2?s.t.??????x1?+x2?≤34x1?+x2?≤9x1?,x2?≥0?
h = [4,-4;-4,8]; f = [-6;-3]; a = [1,1;4,1]; b = [3;9]; [x,fval] = quadprog(h,f,a,b,[],[],zeros(2,1)); disp(x) disp(fval)整數(shù)規(guī)劃
規(guī)劃中的變量限制為整數(shù)時(shí),稱為整數(shù)規(guī)劃。
intlinprog函數(shù)
[x,fval]=intlinpog(c,intcon,A,b,Aeq,beq,LB,UB)[x,fval]=intlinpog(c,intcon,A,b,Aeq,beq,LB,UB)[x,fval]=intlinpog(c,intcon,A,b,Aeq,beq,LB,UB)
- intcon 是只能取整數(shù)的自變量序號(hào)
- 沒(méi)有被 intcon 限定的自變量可以不必是整數(shù);被 intcon 限定的自變量也可以有上下界
- 其他參數(shù)與線性規(guī)劃相同,參見(jiàn)線性規(guī)劃
實(shí)例
minz=?3x1?2x2?x3s.t.{x3intx1,x2≥0x1+x2+x3≤74x1+2x2+x3=12minz=-3x_1-2x_2-x_3\\ s.t.\begin{cases} x_3 \quad int\\ x_1,x_2 \geq 0\\ x_1+x_2+x_3 \leq 7\\ 4x_1+2x_2+x_3 = 12 \end{cases}minz=?3x1??2x2??x3?s.t.??????????x3?intx1?,x2?≥0x1?+x2?+x3?≤74x1?+2x2?+x3?=12?
clear c = [-3;-2;-1]; intcon = 3; A = [1 1 1]; b = 7; Aeq = [4 2 1]; beq = 12; [x,fval] = intlinprog(c,intcon,A,b,Aeq,beq,[0;0;0]); disp(x) disp(fval);總結(jié)
以上是生活随笔為你收集整理的MATLAB二次规划和整数规划的全部?jī)?nèi)容,希望文章能夠幫你解決所遇到的問(wèn)題。
- 上一篇: C语言嵌入式系统编程修炼之软件架构篇
- 下一篇: This relative module