matlab数学建模-非线性规划(无约束规划、有约束规划)

目录

二次规划,沃尔夫法

无约束规划

有约束规划


非线性规划的标准形式:

minf(x)

g_i{(x)}\leq 0,i=1,2,...,m

h_j{(x)}= 0,j=1,2,...,r

gi(x), hj(x)是约束条件,gi(x), hj(x)和f(x)中至少有一个是非线性函数。

非线性模型按照约束条件分:

1)无约束非线性规划模型

2)等式约束非线性规划模型

3)不等式约束非线性规划模型


二次规划,沃尔夫法

min\frac{1}{2}x^{T}Hx+f^{T}x

案例: matlab数学建模-非线性规划(无约束规划、有约束规划)_第1张图片

H和A是矩阵,f和b是列向量

代码如下

clear all
clc
H=[1 -1;-1 2];
f=[-2;-6];
A=[1 1;-1 2;2 1];
b=[2 2 3];
lb=zeros(2,1);
[x,fval,exitflag]=quadprog(H,f,A,b,[],[],lb)

matlab数学建模-非线性规划(无约束规划、有约束规划)_第2张图片

无约束规划

三个功能函数  fminbnd、fminsearch、fminunc 

fminbnd:

求固定区间,单变量函数最小值(一元函数最小值问题)

clear all
clc
[x,y_min]=fminbnd('sin(2*x)',0,2*pi)


matlab数学建模-非线性规划(无约束规划、有约束规划)_第3张图片

 fminsearch:

clear all
clc
f='3*x(1)^3+3*x(1)*x(2)^3-7*x(1)*x(2)+2*x(2)^2';
x0=[0,0];
[x,f_min]=fminsearch(f,x0)

matlab数学建模-非线性规划(无约束规划、有约束规划)_第4张图片

fminunc :

函数matlab数学建模-非线性规划(无约束规划、有约束规划)_第5张图片

 脚本调用

matlab数学建模-非线性规划(无约束规划、有约束规划)_第6张图片

 那段英文我搜了一下,Optimization completed because the size of the gradient is less than
the value of the optimality tolerance.(优化完成,因为梯度的大小小于最优容差的值。)
matlab数学建模-非线性规划(无约束规划、有约束规划)_第7张图片

 x,fval值工作区直接看,也可以命令行输入x,fval查看


有约束规划

fmincon函数,求多变量,有约束,非线性函数的最小值

案例1

matlab数学建模-非线性规划(无约束规划、有约束规划)_第8张图片

 

代码:

目标函数

function f=ex10_6a(x)
f=x(1)*x(1)*(x(2)+2)*x(3);

条件函数

function [c,ceq]=ex10_6b(x)
c(1)=350-163*x(1)^(-2.86)*x(3)^(0.86);
c(2)=10-4*10^(-3)*x(1)^(-4)*x(2)*x(3)^3;
c(3)=x(1)*(x(2)+1.5)+0.0044*x(1)^(-4)*x(2)*x(3)^3-3.7*x(3);
c(4)=375-3.56*10^5*x(1)*x(2)^(-1)*x(3)^(-2);
c(5)=4-x(3)/x(1);
ceq=0;

函数求解程序

clear all
clc
x0=[2 25 20]';
lb=[1 4.5 10]';
ub=[4 50 30]';
[x,fval,exitflag]=fmincon(@ex10_6a,x0,[],[],[],[],lb,ub,@ex10_6b)

结果

matlab数学建模-非线性规划(无约束规划、有约束规划)_第9张图片

 案例2

matlab数学建模-非线性规划(无约束规划、有约束规划)_第10张图片

 

代码

clear all
clc
%初值x0
x0=[10;10;10];
A=[-1 -2 -2;1 2 2];
b=[0;72];
[x,fval]=fmincon('-x(1)*x(2)*x(3)',x0,A,b);

matlab数学建模-非线性规划(无约束规划、有约束规划)_第11张图片

 

你可能感兴趣的:(matlab数学建模,matlab,算法,开发语言)