计算方法实验报告.docxVIP

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