数值计算方法上机作业1.pdfVIP

地球科学与技术学院 地球探测与信息技术 研11—6班 孙殿强 _______________________________________________________________________________ 2.1 2.1 实验22..11 迭代法、初始值与收敛性 实验目的： 初步认识非线性问题的迭代法与线性问题的迭代法的差别，探讨迭代法及初 始值与迭代收敛性的关系。 问题提出： 迭代法是求解非线性方程（组）的基本方法，与线性方程的情况一样，其构 造方法多种多样，但关键是怎样才能使迭代收敛且有较快的收敛速度。 实验内容： 考虑一个简单的代数方程 2 x − x−1= 0 2 1 针对上述方程，可以构造多种迭代法，如x = x −1,x =1+ ,x = x +1等。 n+1 n n+1 n+1 n x n 在实轴上取初始值x ，分别用以上迭代法做实验，记录各算法的迭代过程。 0 实验要求： （1）取定某个初始值，按如上迭代格式进行计算，它们的收敛性如何？重 复选取不同的初始值，反复实验。请读者自行设计一种比较形象的记录方式（如 利用Matlab 的图形功能），分析3种迭代法的收敛性与初始值选取的关系。 （2）对3种迭代法中的某一个，取不同的初始值进行迭代，结果如何？试 分析不同初始值的迭代结果是否有差异。 迭代法程序如下（matlab 编程处理）： function [x]=diedai(g,x0,tol,nmax) %g表示给定的迭代函数 %x0表示给定的初始值 %tol表示允许的误差 %nmax表示最大允许迭代次数 %X是输出的迭代值 X(1)=x0; for i=2:nmax X(i)=feval(g,X(i-1)); I; if(abs(X(i)-X(i-1))tol) break; end 地球科学与技术学院 地球探测与信息技术 研11—6班 孙殿强 _______________________________________________________________________________ if(i==nmax) disp(The method failed after Nmax iterations); end end X 分别编辑如下格式迭代函数： （ ） 1 function y=g(x) y=x^2-1; （2）function y=g(x) y=1+1./x; （3）function y=g(x) y=sqrt(x+1); 在matlab中运行如下代码： fplot([x.^2-x-1],[-2 3]);grid 可得如下图形： 从图中可以得知方程的两个解分别位于【-1，-0.5】及【1.5，2】这两个区间内。 实验时分别选取初始值为1、1.5、5代入 （ ）迭代格式 ， 1 1 A.选取处置为1. diedai(g,1,10^(-5),20) The method failed after Nmax iterations X = Columns 1 through 16 地球科学与技术学院 地球探测与信息技术 研11—6班 孙殿强 _______________________________________________________________________________