非线性方程求根.pptVIP

非线性方程求根第1页，共27页，星期日，2025年，2月5日

MATLAB命令例1：解方程解：在MATLAB里输入命令：X=solve(‘8*x^9+17*x^3-3*x=-1’)可解出9个解。例2：解方程解：X=solve(‘sin(x)=0’),很可惜运行后只输出一个根x=0.缺点：用solve命令不能求出周期函数所对应的全部根。第2页，共27页，星期日，2025年，2月5日

例3：解非线性方程组解：在MATLAB命令窗口输入命令：E1=sym(x^x-4=0);E2=sym(2*x*y+x=1);[x,y]=solve( E1,E2)x1=double(x),y1=double(y)出来的结果为：x=log(4)/lambertw(log(4))y=-1/2*(log(4)-lambertw(log(4)))/log(4)第3页，共27页，星期日，2025年，2月5日

x1=2y1=-0.2500注：MATLAB系统只能做数值运算，并没有符号运算功能，符号运算工具箱（symbolicmathtoolbox)则扩充了MATLAB这方面的功能，它是由Maple的核心来完成。第4页，共27页，星期日，2025年，2月5日

例4：解方程解：将方程化为，在matlab窗口输入命令fa=[8,0,0,0,0,0,17,0,-3,1];xk=roots(fa)运行后得所有根。缺点：命令roots只能求为多项式时方程的根。第5页，共27页，星期日，2025年，2月5日

一、图解法求方程的正根，设分别画出的图形，两条曲线的交点即为原方程的根，从图中观察，根大约为0.38。第6页，共27页，星期日，2025年，2月5日

二、二分法对于求解给定区间的根，二分法是一种既简单又稳健的方法，可以与图解法结合使用。第7页，共27页，星期日，2025年，2月5日

最终解经过n步迭代后，区间长度变为：即为可能的最大误差第8页，共27页，星期日，2025年，2月5日

当给定容许误差时，所需最小迭代步数为：用于二分法计算的函数bisec_n.functionbisec_n(f_name,a,c)tolerance=0.000001;it_limit=30;fprintf(It.abcf(a));fprintf(f(b)f(c)\n);it=0;Ya=feval(f_name,a);Yc=feval(f_name,c);if(Ya*Yc0)fprintf(\n\nStoppedbecausef(a)f(c)0\n);第9页，共27页，星期日，2025年，2月5日

elsewhile1it=it+1;b=(a+c)/2;Yb=feval(f_name,b);fprintf(%3.0f%10.6f,%10.6f,it,a,b);fprintf(%10.6f,%10.6f,%10.6f,%10.6f\n,c,Ya,Yb,Yc);if(abs(c-a)=tolerance)fprintf(Toleranceissatisfied.\n);breakendif(itit_limit)fprintf(Iterationlimitexceeded.\n);breakendif(Ya*Yb=0)c=b;Yc=Yb;elsea=b;Ya=Yb;endendfprintf(Finalresult:Root=%12.6f\n,b);end第10页，共27页，星期日，2025年，2月5日

例1求函数