MATLAB 非线性方程求解与函数极值计算

MATLAB 非线性方程求解与函数极值计算

单变量非线性方程求解:

函数调用 x=fzero(filename,x0)

eg:求f(x)=x-1/x+5在x0=-5和x0=1作为迭代初值时的根。

f=@(x)x-1./x+5;
x1=fzero(f,-5)
x2=fzero(f,1)
x3=fzero(f,0.1)
fplot(f)

eg:求x^2-1=0的根。

f=@(x)x.^2-1;
x=[];
x0=-0.25:0.001:0.25;
for x00=x0
    x=[x,fzero(f,x00)];
end
plot(x0,x,'-o')
xlabel('初值');
ylabel('方程的根');
axis([-0.25,0.25,-1,1])

非线性方程组的求解:

函数调用格式： x=fsolve(filename,x0,option)

eg:

f=@(x)x-1./x+5;
x1=fsolve(f,-5,optimset('Display','off'))
x2=fsolve(f,1,optimset('Display','off'))
x3=fsolve(f,0.1,optimset('Display','off'))

eg:

f=@(x)[sin(x(1))+x(2)+x(3)^2*exp(x(1)),x(1)+x(2)+x(3),x(1)*x(2)*x(3)];
f([1,1,1])
x=fsolve(f,[1,1,1],optimset('Display','off'))
f(x)

函数极值的计算:

无约束最优化问题:

eg:求函数f(x)=x-1/x+5,在区间（-10，-1）和（1，10）上的最小值点。

f=@(x)x-1./x+5;
[xmin,fmin]=fminbnd(f,-10,-1)
[xmin,fmin]=fminbnd(f,1,10)

有约束最优化问题:

约束条件可细化为：

1.线性不等式约束
2.线性等式约束
3.非线性不等式约束
4.非线性等式约束
5.X的上界和下界

eg:

f=@(x)0.4*x(2)+x(1)^2+x(2)^2-x(1)*x(2)+1/30*x(1)^3;
x0=[0.5;0.5];
A=[-1,-0.5;-0.5,-1];
b=[-0.4;-0.5];
c=[0;0];
option=optimset('Display','off');
[xmin,fmin]=fmincon(f,x0,A,b,[],[],c,[],[],option)

eg:

a=[10,30,16.667,0.555,22.2221];
b=[10,50,29,29.888,49.988];
c=[10,18,20,14,25];
f=@(x)sum(c.*sqrt((x(1)-a).^2+(x(2)-b).^2));
[xmin,fmin]=fminsearch(f,[15,30])