matlab 代数方程与最优化问题的求解PPT课件.pptx

代数方程的求解 代数方程的图解法
一元方程的图解法
例：
>> ezplot('exp(-*t)…
*sin(*t+)+*exp…
(-.*t)*cos(*t)-…
.',[ ])
>> hold on,
>> line([,],[,])
% 同时绘制横轴
第页/共页
验证：
>> syms t ; t=.; vpa(exp(-*t)*sin(*t+)+…
*exp(-.*t)*cos(*t)-.)
ans =
-.e-
第页/共页
二元方程的图解法
例：
>> ezplot('x^*exp(-x*y^/)+exp(-x/)*sin(x*y)')
% 第一个方程曲线
>> hold on
% 保护当前坐标系
>> ezplot(‘x^ *…
cos(x+y^) +…
y^*exp(x+y)')
第页/共页
方程的图解法
仅适用于一元、
二元方程的求根
问题。
第页/共页
多项式型方程的准解析解法
例：
>> ezplot('x^+y^-'); hold on % 绘制第一方程的曲线
>> ezplot(‘.*x^-y+.’) % 绘制第二方程
为关于x的次多项式方程
应有对根。图解法只能
显示求解方程的实根。
第页/共页
一般多项式方程的根可为实数，也可为复数。
可用MATLAB符号工具箱中的solve( )函数。
最简调用格式：
S=solve(eqn,eqn,…,eqnn)
（返回一个结构题型变量S，。）
直接得出根： （变量返回到MATLAB工作空间）
[x,…]=solve(eqn,eqn,…,eqnn)
同上，并指定变量
[x,…]=solve(eqn,eqn,…,eqnn，’x,…’)
第页/共页
例：
>> syms x y;
>> [x,y]=solve('x^+y^-=','*x^/-y+/=')
x =
[ -.]
[ -.-.*i]
[ -.+.*i]
[ .]
[ .-.*i]
[ .+.*i]
y =
[ .]
[ .-.*i]
[ .+.*i]
[ .]
[ -.-.*i]
[ -.+.*i]
第页/共页
验证
>> [eval('x.^+y.^-') eval('*x.^/-y+/')]
ans =
[ ,