实验一应数2班201430100209见雪婷.docVIP

实验报告1 非线性方程的迭代解法 班级专业14应数2 姓名见雪婷 学号201430100209 成绩 一、实验目的 掌握方程迭代法思想，会用二分法、牛顿法、弦截法求解方程。 二、实验内容、计算过程与结果 题目 1、（P34 第1题） 用二分法求方程内的根，精确到即误差不超过. 程序 a=3; b=4; epls=0.5e-3; %counter = 0; x=(a+b)/2; while abs(b-a)epls if f_2(a)*f_2(x)0 b=x; else a=x; end x=(a+b)/2; %counter = counter + 1; %counter, end x, function y=f_2(x) y = x^3 - 2 * x^2 - 4 * x-7; 结果 x = 3.6321 题目 2、（P34 第4题） 将一元非线性方程写成收敛的迭代公式，并求其在附近的根，精确到. 程序 epsl=0.005; x(1)=0.5; n=1; x(n+1)=phy(x(n)); while abs(x(n+1)-x(n))=epls; n=n+1; x(n+1)=phy(x(n)); end x, function y=phy(x) y=log(2*cos(x)); 结果 260 至 264 列 0.5398 0.5398 0.5398 0.5398 0.5398 题目 2、（P34 第4题） 将一元非线性方程写成收敛的迭代公式，并求其在附近的根，精确到. 程序 k=1; x(k)=0.5; epls=0.5e-2; x(k+1)=f_iterative(x(k)); while abs(x(k+1)-x(k))epls k=k+1; x(k+1)=f_iterative(x(k)); end x, function y=f_iterative(x) y=2*cos(x)-exp(x); 结果 260 至 264 列 1.1129 -2.1592 -1.2255 0.3834 0.3876 题目 3、（P35 第11题）分别用下列方法求方程在附近的根，要求有三位有效数字。 （1）牛顿法，取； （2）弦截法，取； （3）快速弦截法，取. 程序 （将程序粘贴在此） (1) clear; k=1; x(k)=pi/4; epls=0.5e-3; x(k+1)=x(k)-f_iterative(x(1))/df_iterative(x(1)); while abs(x(k+1)-x(k))epls k=k+1; x(k+1)=x(k)-f_iterative(x(k))/df_iterative(x(k)); end x, function y=f_iterative(x) y=4*cos(x)-exp(x); function y=df_iterative(x) y=-4*sin(x)-exp(x); (2) clear; k=1; x(k)=pi/4; x(k+1)=pi/2; epls=0.5e-3; x(k+2)=x(k+1)-(f_iterative(x(k+1))/ (f_iterative(x(k+1))-f_iterative(x(k))))*(x(k+1)-x(k)); while abs(x(k+1)-x(k))epls k=k+1; x(k+2)=x(k+1)-(f_iterative(x(k+1))/ (f_iterative(x(k+1))-f_iterative(x(1))))*(x(k+1)-x(1)); end x, function y=f_iterative(x) y=4*cos(x)-exp(x); (3) clear; k=1; x(k)=pi/4; x(k+1)=pi/2; epls=0.5e-3; x(k+2)=x(k+1)-(f_iterative(x(k+1))/ (f_iterative(x(k+1))-f_iterative(x(k))))*(x(k+1)-x(k)); while abs(x(k+1)-x(k))epls k=k+1; x(k+2)=x(k+1)-(f_iterative(x(k+1))/ (f_iterative(x(k