第二章非线性方程组()的数值解法.docVIP

2.1 方程（组）的根及其MATLAB命令 2.1.2 求解方程（组）的solve命令 求方程f(x)=q(x)的根可以用MATLAB命令: x=solve(方程f(x)=q(x),待求符号变量x) 求方程组fi(x1,…,xn)=qi(x1,…,xn) (i=1,2,…,n)的根可以用MATLAB命令: E1=sym(方程f1(x1,…,xn)=q1(x1,…,xn)); ……………………………………………………. En=sym(方程fn(x1,…,xn)=qn(x1,…,xn)); [x1,x2,…,xn]=solve(E1,E2,…,En, x1,…,xn) 2.1.4 求解方程（组）的fsolve命令 fsolve的调用格式： X=fsolve(F,X0) 2.2 搜索根的方法及其MATLAB程序 2.2.1 作图法及其MATLAB程序 作函数在区间 [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的绝对值越小，图形越精确. 作函数在区间 [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.2.2 逐步搜索法及其MATLAB程序 逐步搜索法的MATLAB主程序 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);n=(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 例2.2.4 用逐步搜索法的MATLAB程序分别求方程和在区间上的根的近似值，要求精度是0.000 1. 解 将逐步搜索法的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 2.3 二分法及其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)); %程序中调用的fun.m 为函数 if ya* yb0, disp(注意：ya*yb0,请重新调整区间端点a和b.), return end max1=-1+ceil((log(b-a)- log(abto