matlab求解二元函数极值
依然是机房中的R2010a版本
命令:
1、x=fminsearch(fun,x0)或x=fminunc(fun,x0)求极小值点x,初值选为x0
2、[x,fmin]=fminsearch(fun,x0)或[x,fmin]=fminunc(fun,x0)
3、fminsearch采用单纯形法,fminunc采用牛顿法
除了fminsearch和fminunc这两种命令外,建立函数还可以用不同的方法:
建立函数的方法
以p191task2_2为例子,采用字符串建立函数
% p191task2_2
%求min(f(x))=(4*x(1)^2+2*x(2)^2+4*x(1)*x(2)+2*x(2)+1)*exp(x(1))
clc,clear
f='(4*x(1)^2+2*x(2)^2+4*x(1)*x(2)+2*x(2)+1)*exp(x(1))';
digits 15
x0=[0,0];
[x,fmin]=fminsearch(f,x0);
vpa(x,15)
vpa(fmin,15)
然而此题的问题在于,无法找到合适的初值,初值每变动一次,fmin也随之变动。
为此,画出函数图像:
[x,y]=meshgrid(-30:0.3:-10,-30:0.3:-10);
f=(4*x.^2+ 2*y.^2+4.*x.*y+2.*y+1).*exp(x);
figure(1)
surf(x,y,f)
figure(2)
contour(x,y,f,100)
补充:contour为画等高线命令
contour(x,y,z,n) n为等高线条数;
由图可以大致得到,函数在所选区域内全部为正,最小值接近0,在另外选择的若干区域中,结果类似,看出最小值大概是0。
以p191task2_3为例子,用.m文件建立函数
先创建.m函数
%p192task2_3 fun
%fun2_3.m
function f=fun2_3(x)
f=4*x(1)^2+5*x(1)*x(2)+2*x(2)^2;
再在同一个目录下写
%p191task2_3.m
clc,clear
[x,y]=meshgrid(-10:0.3:10,-10:0.3:10);
z=4*x.^2+5*x.*y+2*y.^2;
figure(1)
surf(x,y,z)
figure(2)
contour(x,y,z)
x0=[-3,-3];
% [x,fmin]=fminsearch(@fun2_3,x0)
[x,fmin]=fminunc(@fun2_3,x0)
结果: