数值计算方法例题--邢斌分解.doc

数值计算方法例题 分别用单点弦割法和双点弦割法求方程=的根，要求. 解： 因为，令，则由于的判别式，故没有极值点.由于=-70，，因此在（1，2）内仅有一根. （1）用单点弦割法，迭代公式为 ， =1，2，… 取，，计算结果如下： 0 1 5 1.368 898 095 1 2 6 1.368 798 076 2 1.304 347 826 7 1.368 809 227 3 1.376 053 620 8 1.368 807 983 4 1.3680 010 58 9 1.368 808 121 此时，已满足精度要求，故取=1.368808121即可. 若采用双点弦割法，迭代公式为 ， =1，2，… 仍取，，则有： 0 1 5 1.368 807 460 1 2 6 1.368 808 108 2 1.304 347 826 3 1.357 912 305 4 1.369 013 327 取=1.368808108即可. 用牛顿法求解，要求. 解：由上题知，在（1，2）内有一个根，且0，，故取，利用牛顿迭代公式： ， =1，2，… 计算结果如下： 0 2 3 1.368 869 419 1 1.6 4 1.368 808 109 2 1.383 388 704 5 1.368 808 108 所以取=1.368808108. 已知函数的函数表如下： 0.40 0.55 0.65 0.80 0.90 1.05 0.410 75 0.578 15 0.696 75 0.88 11 1.026 52 1.253 82 求四次牛顿差值多项式，并由此求的近似值. 解 构造差商表如下： 一阶差商 二阶差商 三阶差商 四阶差商 0.40 0.410 75 0.55 0.578 15 1.116 00 0.65 0.696 75 1.186 00 0.280 00 0.80 0.888 11 1.275 73 0.358 93 0.197 33 0.90 1.026 52 1.384 10 0.433 48 0.213 00 0.031 34 故四次牛顿差值多项式为 =0.41075+1.11600（-0.4）+0.28000（- 0.4）(- 0.55)+ 0.19733(- 0.4)(- 0.55)(- 0.65)+0.03134(- 0.4)(- 0.55)× (- 0.65)(- 0.80) 所以. 给出自然对数In和它的导数的数表如下： 0.40 0.50 0.70 0.80 In -0.916 291 -0.693 147 -0.356 675 -0.223 144 2.50 2.00 1.43 1.25 用拉格朗日插值公式求In0.60. 用埃尔米特插值公式求In0.60. 解 列表计算如下： 0 -0.166 667 =-15.833 333 1 0.666 667 =1.666 667 2 0.666 667 =-1.666 667 3 -0.166 667 =15.833 333 用拉格朗日插值公式得 In0.60≈== （2）用埃尔米特插值公式得 得 In0.60≈=-0.510899 求使达到最小. 解 由于题中提分达到最小，实际上是上求，使其称为的最佳平方逼近多项式，故b,a满足正规方程组 实际计算有： ， ， ， 故有 解得，. 利用正交化方法求[0，1]上带权的前三个正交多项式，，. 解 =-，利用公式 =1 = = =，=0，1 以及内积定义，得 所以 再由 得 故得 . 从地面发射一枚火箭，在最初80s内，记录其加速度如下表，试求火箭在第80s时的速度. t(s) 0 10 20 30 40 50 60 70 80 a(m·) 30.00 31.63 33.44 35.47 37.75 40.33 43.29 46.69 50.67 解 应用复化辛普森公式计算。此题中=80，9个节点，故，.由于火箭从地面向上发射，，因此 ≈·20[30.00+4×（31.63+35.47+40.33+46.69）+2×（33.44+37.75+43.29）+50.67]=3087.0