数值分析实验报告.doc

实验报告:函数逼近插值多项式补充 问题/：对于给函数f(x)=——~，取点母=cos + I兀、e取0, / 。取力 l+25x2 2/?+ 2 或％。试画出拟合曲线并打印出方程，与第二章计算实习题2的结果进行比较。 问题￡：对于给函数/? = 一 在区间，处取絵=十02(口％￡?…,力)，试求M l+25x 次曲线拟合，试画出拟合曲线并打印出方程，与第二章计算实习题Z的结果进行比较。 实验目的：通过编程实现牛顿插值方法和函数逼近，加深对多项式插值的理解。应用所 编程序解决实际算例。 实验要求： /. 认真分析问题，深刻理解相关理论知识并能熟练应用； 2 编写相关程序并进行实验： M 调试程序，得到最终结果： 分析解释实验结果： 按照要求完成实验报告。 实验原理： 详见《数值分析第S版》第二章、第三章相关内容。 实验内容： (/)问题 这里我们可以沿用实验报告一的代码，对其进行少量修改即可。 当—力时，代码为： cte Eg 心皿？的〃⑵g， ^-LKl+25.\t^2)； /二“ “力 炉 /二“ -1*j 仪)二的).化?〃〃厶〃?“a尹〃力 炉/二2歼 夕0二夕a沪倂““ ■〃” 4艸4 P ： Q10 二岭a （卯 zl（Q）. 5）： ■100014： 力（P%? iO）: 员二皿（U（2冶小2）.魚xO）； §0、%O. $2） 血加（￡） 如“刃 由 此 我 们 可 以 得 到 P°Z=Y6? 63^kM0+3.0962c -冷◎乡十 30 7ftAf-7.2714c -127 ? 133.各?十7.1777c - l^S+61.密 Z -1?5宓-件小3 -12.47TZ - /. 62他-16 Z0 并可以得到牛顿插值多项式在厂人 处的图形，并和原函数进行对比，得京八 // // j / / / \ \\ -?- . / ...-?- *?? 、 ■021 I I I I I I I I I ?1 418 05 -0.4 42 0 Q2 0.4 Q6 M 知牛顿插值多项式（5）函数和原函数图形 当心）、将上述代码中的 gg “改为“QG%「即可。 由 此 我 们 可 以 得 到 P^（k）=6466.夕於％十￡0%%-E*?V乡■羽％￡（T小?M503%?忆47七777S4?（T46? 彌00. M4+2.725 3c- ?存診十7^236. （T^12-39333. CM0P2636.4,6566c- 1E7?2537.3^^6+306.6^4-21.76彳小2 P.O 并可以得到牛顿插值多项式在广人少上的图形，并和原函数进行对比，得夕乡冼 今牛顿插值多项式(-20)函数和原函数图形 回顾一下实验一的结果(见殉.E),我们不难发现，仅仅是改变了乡的取值，结果发生 了很大的变化。实验一中，插值多项式与原函数产生了很大的偏差，并且随着分的段数的增 加，英误差不断变大，但是在本次实验中，我们不难发现，虽然多项式依旧存在震荡现象, 但是误差小了很多，而且随着分的段数的增加，插值多项式曲线与原函数曲线已经十分接近 实验一结果 这个例子说明：采用切比雪夫节点替代等距节点可以消除龙格现象。 (￡)问题￡： 分析问题，发现在这个问题中，我们已经知道了原函数，同时它也告诉我们所需取的〃 个点的值，所以这里可以用两种方法进行函数逼近得到拟合曲线。 首先釆用最小二乘法来考虑这个问题，编写代码如下(这里没有直接调用0乡俗函数)： 如 ait de “二/,/〃+羽垃吻 a 妙 s g ? mN 炉d卄* )二皿恥(—)?制)： 4二尹?? 久ann cm( sphere ⑷=100014 *. 9)： 员二皿(U(1+2^込2).魚3 fit(9?乡 0. vO.汹 如打,刃 我们可以得到一个三次多项式：S戶」6皈■卅小3?0.5751^^2 - 9.4553c-ir% + 0.4^12.同 时我们也可以得到它与原函数的图形，如图7初盗 f Y 1 一 —? — 、 / ? .一 、、 、— / ?1 a 8 刀启 -0.4 a 2 a 2 0.4 0 6 x 彩?“最小二乘法的结果 我们发现得到的结果和原函数产生了巨大的误差。首先观察得到的多项式.我们发现它 的￡次项系数非常小，原因是原函数是一个偶函数，这将导致奇次项系数基本为久 这里我们调整“，对结果进行观察，取日、6、f, 10, 20。我们可以得到殉?，忿?乡。 S尸4畑 ZU3703e -1^^3-2.060^^2-/. 176% -妙2 十 0 65522 S尸?4,63才小6+4.07浙r?10小5七9 476^25.冰-仲小3YS96牛小2七1.3229c?￥上十0 7^461 S尸%初於并M浙7■図小7YM60T46+6?免％，7?十殳Z SZ7%Y?