一、符号方程求解
在MATLAB中,求解用符号表达式表示的代数方程可由函数solve实现,其调用格式为:
solve(s)(求解符号表达式s的代数方程,求解变量为默认变量,当方程右端为0时,方程可以不标出等号和0,仅标出方程的左端)
solve(s,v)(求解符号表达式s的代数方程,求解变量为v)
solve(s1,s2,…,sn,v1,v2,…,vn)(求解符号表达式s1,s2,…,sn组成的代数方程组,求解变量分别为v1,v2,…,vn)
让我举个栗子:
求解如下方程:
代码:
x=solve('2*sin(3*x-pi/4)=1')
代码:
x= solve('x+x*exp(x)-10', 'x') %仅标出方程的左端
二、求方程f ( x ) = 0数值解的基本方法
并非所有的方程 f ( x ) = 0 都能求出精确解或解析解,存在这种解的方程就需要用数值解法求出近似解,有几种常见的数值解法:
1.二分法
2.迭代法:切线法、割线法(弦截法)
注:二分法简单方便,但收敛速度慢;
迭代法虽然收敛速度稍微快点,但需要判断能否收敛;
只要初值选取得当,切线法具有恒收敛且收敛速度快的优点,但需要求出函数的导数;
弦截法不需要求导数,特别是前面介绍的快速弦截法,收敛速度很快,但是需要知道两个近似的初始根值才能作出弦,要求的初始条件较多。
三、方程f(x) = 0数值解的MATLAB实现
MATLAB中求方程数值解的办法很多,有的是专用指令,有的是根据方程性质而借用其他专用指令求得的
这里我们就先介绍求函数零点指令fzero吧
求解方程f ( x ) = 0的实数根也就是求函数f ( x)的零点。MATLAB中设有求函数f (x)零点的指令fzero,可用它来求方程的实数根。该指令的使用格式为:fzero (fun, x0, options)
①输入参数fun为函数f (x)的字符表达式、内联函数名或M函数文件名。
②输入参数x0为函数某个零点的大概位置(不要取零)或存在的区间[xi,xj],要求函数f (x)在x0点左右变号,即f (xi)f (xj) < 0。
③输入参数options可有多种选择,若用optimset (‘disp’, ‘iter’)代替options 时,将输出寻找零点的中间数据。
④该指令无论对多项式函数还是超越函数都可以使用,但是每次只能求出函数的一个零点,因此在使用前需摸清函数零点数目和存在的大体范围。为此,一般先用绘图指令plot, fplot或ezplot画出函数f (x)的曲线,从图上估计出函数零点的位置。
让我再举个栗子:
如:求方程 x^2 + 4sin(x) = 25 的实数根(-2π<x < 2π)
解:
(1)若fun为函数f (x)的字符表达式
①首先要确定方程实数根存在的大致范围。为此,先将方程变成标准形式f(x) =x2 + 4sin(x) - 25 = 0,作f(x)的曲线图:
x=-2*pi:0.1:2*pi;
f=x.^2+4*sin(x)-25;
plot(x,f);grid on;
从曲线上可以看出,函数的零点大约在x1= - 4和x2=5附近
②直接使用指令fzero求出方程在x1≈ - 4时的根
x1= fzero ('x^2+4*sin(x)-25',-4)
若键入:
fzero ('x^2+4*sin(x)-25',-4, optimset('disp', 'iter'))
将会显示迭代过程
中间数据表明,求根过程中不断缩小探测范围,最后得出- 4附近满足精度的近似根
③求x2≈ 5的根:
x2= fzero ('x^2+4*sin(x)-25',5)
(2) 若fun为函数f (x)的M函数文件名
将方程x2 + 4sin(x) = 25编成M函数文件(实用中在函数较为复杂、而又多次重复调用时,才这样做),用fzero求解。
①在M文件编辑调试窗中键入:
function yy=fan (x)
yy= x^2+4*sin(x)-25;
以fan为文件名存盘,退出编辑调试窗,回到指令窗。
②确定根的大体位置;
③在指令窗中键入下述指令可求出 - 4 附近的根:
x1= fzero ('fan',-4)
键入下述指令可求出5附近的根:
x2= fzero ('fan',5)
我还要再举个栗子:
如:求f(x)=x-10^x+2=0在x0=0.5附近的根
解:
x=-2.5:0.01:0.5;
fx=x-10.^x+2;
plot(x,fx)
从f(x)的曲线看出曲线的零点有两个,一个在x=-2附近,另一个在x=0.5附近
①建立函数文件funx.m。(function [输出变量列表]=函数名(输入变量列表))
functionfx=funx(x)
fx=x-10.^x+2;
②调用fzero函数求根。
z=fzero('funx',0.5)