计算方法实验报告.docVIP

实验报告 一、求方程f(x)=x^3-sinx-12x+1的全部根, ε=1e-6 用一般迭代法; 用牛顿迭代法; 并比较两种迭代的收敛速度。. 其次，分析得到其根所在的区间。 令，可得到. 用一阶导数分析得到和两个函数的增减区间；再用二阶导数分析得到两个函数的拐点以及凹凸区间. ③ 在直角坐标轴上描摹出和的图，在图上可以看到他们的交点，然后估计交点所在的区间，即是所要求的根的区间。经过估计，得到根所在的区间为，和. 一般迭代法 （1）算法步骤： 设为给定的允许精度，迭代法的计算步骤为： 选定初值.由确定函数，得等价形式. 计算.由迭代公式得. 如果，则迭代结束，取为解的近似值；否则，用代替，重复步骤②和步骤③. （2）程序代码： ① 在区间内， 代码： clc x0=-3.5; %初值 iter_max=100; %迭代的最大次数 ep=1e-6; %允许精度 k=0; while k=iter_max %k从0开始到iter_max循环 x1=(sin(x0)+12*x0-1).^(1/3); %代入，算出的值 if abs(x1-x0)ep %与允许精度作比较 break; %条件成立，跳出循环 end x0=x1; %条件不成立，用代替 k=k+1; %k加1 end x_star=x1, iter=k %为解的近似值，iter为迭代次数 运行结果：x_star = -3.4101 ；iter =14 ②在区间内， 代码： clc x0=0.5; %初值 iter_max=100; %迭代的最大次数 ep=1e-6; %允许精度 k=0; while k=iter_max %k从0开始到iter_max循环 x1=(1/12)*(x0.^3-sin(x0)+1); %代入，算出的值 if abs(x1-x0)ep %与允许精度作比较 break; %条件成立，跳出循环 end x0=x1; %条件不成立，用代替 k=k+1; %k加1 end x_star=x1, iter=k %为解的近似值，iter为迭代次数 结果：x_star = 0.07696；iter =6 ③在区间内， 代码： clc x0=3.5; %初值 iter_max=100; %迭代的最大次数 ep=1e-6; %允许精度 k=0; while k=iter_max %k从0开始到iter_max循环 x1=(sin(x0)+12*x0-1).^(1/3); %代入，算出的值 if abs(x1-x0)ep %与允许精度作比较 break; %条件成立，跳出循环 end x0=x1; %条件不成立，用代替 k=k+1; %k加1 end x_star=x1, iter=k %为解的近似值，iter为迭代次数 运行结果：x_star = 3.4101 ；iter =10 牛顿迭代法 （1）算法步骤： ① 选定初值，计算，. ② 按公式迭代，得新的近似值，并计算，. 对于给定的允许精度，如果，则终止迭代，取；否则，，在转回步骤②计算. （2）程序代码： ①在区间内， clc x1=-3.5; %初值 k=0; while k=100 %k从0开始到100循环 x0=x1; %将初值赋给 f0=x0.^3-sin(x0)-12*x0+1; %计算 f1=3*x0.^2-cos(x0)-12; %计算 x1=x0-f0/f1; %计算得到新的近似值 if abs(x1-x0) 1.0e-6 %与允许精度作比较 break; %条件成立，跳出循环 end k=k+1; %条件不成立，k加1 end x_star=x1, iter=k %为解的近似值，iter为迭代次数 运行结果：x_star = -3.4911；iter =2 ②在区间内， clc x1=0.5; %初值 k=0; while k=100 %k从0开始到100循环 x0=x1; %将初值赋给 f0=x0.^3-sin(x0)-12*x0+1; %计算 f1=3*x0.^2-cos(x0)-12; %计算 x1=x0-f0/f1; %计算得到新的近似值 if abs(x1-x0) 1.0e-6 %与允许精度作比较 break; %条件成立，跳出循环 end k=k+1; %条件不成立，k加1 end x_star=x1, iter=k %为解的近似值，iter为迭代次数 运行结果：x_star =0.0