《C语言程序设计教学资料》chap6s.pptVIP

其迭代关系描述为 root x1+x2 /2 x2 root; /*为下一次迭代准备*/ x1 root; 该重复过程由误差控制。 * * 例题，求x3-5x2+16x-80 0。 #include void main float x0, x1, x2, fx0, fx1, fx2; do scanf %f,%f, x1, x2 ; fx1 pow x1, 3 -5*pow x1, 2 +16*x1-80; fx2 pow x2, 3 -5*pow x2, 2 +16*x2-80; while fx1*fx2 0 ; * * x0 x1+x2 /2.0; fx0 pow x0, 3 - 5*pow x0, 2 +16*x0-80; while fabs fx0 1e-5 x0 x1+x2 /2.0; fx0 pow x0, 3 - 5*pow x0, 2 +16*x0-80; if fx0*fx1 0 x2 x0; fx2 fx0; * * else x1 x0; fx1 fx0; printf the root %f\n, x0 ; 牛顿迭代法 * * 0 Y X x* xk Pk xk+1 xk+2 Pk+1 Pk+2 xk+3 设xk是方程f x 0的精确解x*的附近的一个猜测解，过点 Pk xk, f x 作f x 的切线。该切线方程为： y f xk +f’ xk * x-xk 与X轴的交点是方程：f xk +f’ xk * x-xk 0的解。 其解为： xk＋1＝ xk－ f xk ∕f’ xk 牛顿迭代法的迭代公式。 经过有限次迭代后，可以求得符合误差要求的近似根。 * * 例题，用牛顿迭代法计算一个正实数a的平方根。 算法设计 迭代关系：x x－f x / f’ x x－ x*x－a /2x x+a/x /2 给定初值：x a，可以任意 给定误差? ：误差是控制迭代的条件，当误差大于?时，要继续迭代。 可以是：fabs x*x-a ? * * #include void main double a, x; scanf %lf, a ; x a/3.0; while fabs x*x-a 1e-6 x 0.5* x+a/x ; printf %lf\n, x ; * * 弦截法 * * 0 Y X x x2 x1 f x1 f x2 f x x* a.先取f x 0的两个粗略解x1与x2。若f x1 与f x2 符号相反，在区间 x1,x2 中至少有一个根。 b.连接f x1 和f x2 两点，此线（即弦）交X轴于x。 x点坐标可求出： x x1*f x2 -x2*f x1 / f x2 -f x1 c.求出f x ，若f x 与f x1 同符号，则根必在区间 x,x2 内，此时将x作为新的x1。若f x 与f x2 同符号，则根必在区间 x1,x 内，此时将x作为新的x2。 d.重复，直到|f x | ?为止。 x就是所求的根。 * * 例题，求x3-5x2+16x-80 0。 给定初值：x1, x2 f x1 , f x2 迭代关系：x x1*f x2 -x2*f x1 / f x2 -f x1 x1 x 或 x2 x /*为下一次迭代准备*/ 给定精度? ：精度是控制迭代的条件，当误差大于? 时， 要继续迭代。 可以是：fabs f x ? * * void main int i, j; for i 1; i 10; i++ for j 1; j i; j++ printf %d*%d %d , i, j, i*j ; printf \n ; * * 循环的嵌套 * *** ***** ******* ********* *********** ************* *************** * * 行的控制 i:1～8 ‘*’的个数ｊ与当前行的关系 j 2*i-1 ‘*’前面的空格k与行的关系: 开始时,第一行有8个空格 每多一行,少一个空格 k 8-i while i 8 for k 1;k 8-i;k++ 输出空格; for j 1;j 2*i-1 ;j++ 输出* 循环的嵌套 void main int i, j, k; i 1; while i 8 for k 1;k 8-i;k++ printf ; for j 1;j 2*i-1 ;j++ printf * ; printf \n ; i++; * * 6.7 几种循环