实验指导书_非线性方程(组)数值解法.docVIP

实验指导书_非线性方程（组）的数值解法 隔根区间的求法（确定根的初始近似值）：作图法，逐步搜索法等 求根的方法：二分法，迭代法，牛顿法，割线法，抛物线法，迭代法的加速等 实验平台：MATLAB软件 说明：与方程求根有关的MATLAB函数另见MATLAB学习资料 一、搜索根的方法及其MATLAB程序 求解非线性方程根的近似值时，首先需要判断方程有没有根？如果有根，有几个根？如果有根，需要搜索根所在的区间或确定根的初始近似值（简称初始值）.搜索根的近似位置的常用方法有三种:作图法、逐步搜索法和二分法等，使用这些方法的前提是高等数学中的零点定理. 1．作图法及其MATLAB程序 （1）作函数在区间 [a,b] 的图形的MATLAB程序一 x=a:h:b; % h是步长 y=f(x); plot(x,y) grid, gtext(y=f(x)) 说明： ①此程序在MATLAB的工作区输入，运行后即可出现函数的图形.此图形与轴交点的横坐标即为所要求的根的近似值. ②区间[a,b] 的两个端点的距离 b-a 和步长h的绝对值越小，图形越精确. （2）作函数在区间 [a,b] 的图形的MATLAB程序二 将化为，其中是两个相等的简单函数 x=a:h:b; y1=h(x); y2=g(x); plot(x, y1, x, y2) grid,gtext( y1=h(x),y2=g(x)) 说明：此程序在MATLAB的工作区输入，运行后即可出现函数和的图形.两图形交点的横坐标即为所要求的根的近似值. 2．逐步搜索法及其MATLAB程序 逐步搜索法也称试算法.它是求方程根的近似值位置的一种常用的方法. 逐步搜索法依赖于寻找连续函数满足与异号的区间.一旦找到区间，无论区间多大，通过某种方法总会找到一个根. MATLAB的库函数中没有逐步搜索法的程序，现根据逐步搜索法的计算步骤和它的收敛判定准则编写其MATLAB程序如下，命名为zhubuss.m. 输入输出说明：a和b的值，tol， k=(b-a)/h+1, 方程根的近似值r. function [k,r]=zhubuss(a,b,h,tol) % 输入的量： % a和b是闭区间[a,b]的左、右端点； %h是步长； % tol是预先给定的精度. % 运行后输出的量： % k是搜索点的个数； % r是方程在[a,b]上的实根的近似值，tol； X=a:h:b; Y=funs(X); %f(x)的M文件，对应名称为funs.m的函数文件 n=fix((b-a)/h)+1; m=0; X(n+1)=X(n); Y(n+1)=Y(n); for k=2:n X(k)=a+k*h; Y(k)=funs(X(k)); %程序中调用的funs.m为函数 sk=Y(k)*Y(k-1); if sk=0, m=m+1; r(m)=X(k); end xielv=(Y(k+1)-Y(k))*(Y(k)-Y(k-1)); if (abs(Y(k))tol)( xielv=0) m=m+1;r(m)=X(k); end end 例1 用逐步搜索法的MATLAB程序分别求方程和在区间上的根的近似值，要求精度是0.0001. 解? 将逐步搜索法的MATLAB程序保存名为zhubuss.m的M文件. 建立M文件funs.m function y=funs(x) ??????? ?? ?y=2.*x.^3+2.*x.^2-3.*x-3 在MATLAB工作窗口输入如下程序 [k,r]=zhubuss(-2,2,0.001,0.0001) 运行后输出的结果 k =4001 r = -1.2240?? -1.0000?? -1.0000?? -0.9990??? 1.2250 即搜索点的个数为k =4 001，其中有5个是方程 的近似根，即r = -1.224 0，-1.000 0，-1.000 0，-0.999 0，1.225 0，其精度为0.000 1. 在程序中将y=2.*x.^3+2.*x.^2-3.*x-3用y=sin(cos(2.*x.^3)) 代替，可得到方程在区间 上的根的近似值如下 r =? -1.9190?? -1.7640?? -1.5770?? -1.3300?? -0.9220??? 0.9230??? 1.3310 1.5780??? 1.7650??? 1.9200 ? ?二分法及其MATLAB程序 ? 二分法的MATLAB程序 二分法的MATLAB主程序function [k,x,wuca,yx]=erfen(a,b,abtol) a(1)=a; b(1)=b; ya=fun(a(1)); yb=fun(b(1));