程序设计基础第3章(2010级).pptVIP

* 例3.25求一元方程f(x)=0的根 用某种数学方法导出等价的形式x=g(x)，算法步骤为： （1）选一个方程的近似根，赋给变量x0； （2）x1=x0，然后计算x0 =g(x1) ； （3）当x0与x1的差的绝对值还未小于指定的精度要求时，重复步骤（2）的计算。 若方程有根，并且用上述方法计算出来的近似根序列收敛，则按上述方法求得的x0就认为是方程的近似根。上述算法用程序的形式表示为： {x0=初始近似根； do { x1=x0； x0=g(x1)； /*按特定的方程计算新的近似根*/ } while ( fabs(x0-x1)Epsilon)； printf(“方程的近似根是%f\n”，x0)； } 第3章 程序控制结构 * （例3.15）求一元方程3x3+4x2-2x+5=0的一个近似解 设: f(x)=3x3+4x2-2x+5 用牛顿迭代方法构造x=F(x)的形式，即 xk+1=xk- F(xk)/F’(xk) k=0,1,2…. 以上所产生的数列如果收敛，可以求出一个近似解。 xk+1-xk=- F(xk)/F’(xk) 当F(xk)/F’(xk) 的绝对值小于某个指定的很小数ε时，就求出近似解， 取ε=1.0e-6 方程在[-2，-3]之间有一个解。 * # include “stdio.h” # include “math.h” # define Epsiln 1.0e-6 void main( ) { double x, d; x=-2.0; do {d=(((3.0*x+4.0)*x-2.0)*x+5)/((9.0*x+8.0)*x-2.0); x=x-d; }while(fabs(d) Epsiln); printf(“The root=%f\n”,x); } * 设 f(x)是区间[ a ,b ]上的连续函数，并且f(a)*f(b)0。求方程f(x)=0的解。 二分迭代法： 第3章 程序控制结构 X1=a, x2=b root=(x1+x2)/2 fabs (f(root))1e-6 f(x1)*f(root)0 x2=root x1=root Y N root=(x1+x2)/2 输出：root * 例3.26用二分法求sin(x)=0在[3π/4，3π/2]内的解 # include “math.h” /*二分法迭代*/ # include “stdio.h” # define PI 3.14159 void main( ) { double x1, x2, root; x1=3*PI/4; x2=3*PI/2 ; /*假设a=PI*3/2,b=PI*3/2*/ root=(x1+x2)/2; while(fabs(sin(root))1e-6) /*假设f(x)=sin(x)*/ { if(sin(x1)*sin(root)0) x2=root; else x1=root; root=(x1+x2)/2; } printf(“solutions:%10.5f \n”,root); } 第3章 程序控制结构 * 例3.27 求两个非负整数m，n的最大公约数（公因子） 问题分析： 辗转相除法：假设mn0，一定存在q和r，使得 m=qn+r ⑴ 如果r≠0，m=n，n=r，重做上面的工作，存在新的q和r，使得新得m，n具有上面表达式⑴的关系。如果新r≠0，再次重复上面的工作，直到r=0，这时m，n的最大公约数就是n。 上面的工作是一种迭代，将 （m， n）=（n，r）=（r1，r2）=（r2，r3）=· · · 通过求最大公约数，可以很容易求得最小公倍数，判断两个非负整数是否互质（最大公约数是1）。 第3章 程序控制结构 * 求最大公约数程序 #include “iostream.h” void main( ) { int m, n, r; cout“