最优化问题的MATLAB求解.ppt

主要内容
代数方程的求解
■无约束最优化间题求解
■有约束最优化问题的计算机求解
整数规划问题的计算机求解
2018/10/17
代数方程的求解
■代数方程的图解法
■多项式型方程的准解析解法
般非线性方程数值解
2018/10/17
1代数方程的图解法

【例1】用图解法求c-3sin(41+2)+4-051c0(21)=
3>> ezplot('exp(-3*t)*sin(4""t+2)+
4“exp(-)“cos(2t)-’,[05])
hold on, line(lo, 5,[0,0)
sx;t=
vpa(exp (3*t)*sin(4*t+2)
4“exp(-)Cos(2“t)-)
2018/10/17

/2
【例2】用图解法求
cos(x+y=)+y
第一个方程曲线
ezplot(x 2*exp (x*y 2/2)+exp(-x/2)*sin(x y)')
第二个方程曲线
hold on
ezplot(y 2"cos(y+x 2)+x 2 exp(x+y)')
2018/10/17

【例3】用图解方法求
-y+=0
> ezplot('x 2+y 2-1);
hold on
ezplot(-y+)
2018/10/17
ean eano ..
egn
n)
%最简调用方式
Lx,.]=solve(eqn, eqn2
%直接得出根
] =solve(eqn, eqn2
ean
%同上,并指定变量
2018/10/17
【例4】用 solve函数求{x2+y2-1=0
-y+=0
syms x y
[x,y]= solve(x^2+y2-1=0,
75x3/100-y+9/10=0)
检验
管>>[eval('x."2+-1’),
eval(75*-y+9/10)]
2018/10/17
+3y3+2
【例5】求解
1/2
x2+3y+
十
/4
>>[x,y,Z]= solve(x+3y^3+2*z2=1/2’,
x^2+3xy+z^3=2’,x^3+2*z+2xy2=2/4’)
检验:
>>err=[x+3y."3+-1/2
2+3“y+
+2z+2“-2/4];
norm(double (eval(err))
2018/10/17
另一种方法:
[x,y,2]= solve(x+32y^3+2*z2=1/2’,
x^2+3*y+z^3=2’,’x^3+2*zxy^2=2/4’);
err=[x+3y.^3+2*-1/2,+3*y+-2,
3+2z.“-2/4]
norm(double (eval(err))
2018/10/17
【例6】求解
X-+X+-+
syms x y
[x,y]= solve(x^2/2+x+3/2+2/y+5/(2*y^2)+
3/x3=0”,y/2+3/(2*x)+1/x4+
)
检验
size (x)
>err=[+x+3/2+2./y+5./()+
3,y/2+3./(2x)+1./+
norm(double (eval (err))
2018/10/17