谷根代数值分析--上机实习报告Word格式.doc
《谷根代数值分析--上机实习报告Word格式.doc》由会员分享，可在线阅读，更多相关《谷根代数值分析--上机实习报告Word格式.doc（28页珍藏版）》请在冰豆网上搜索。
的根，观察不同初试值条件下的收敛性，并给出你的结论
【解】：
2.1算法分析
（1）二分法
二分法伪代码如下所示;
输入：
，定义函数，给定精度
输出：
待解方程的根。
while
if
由于本题中恒大于0，不满足二分法的条件，故本例无法通过二分法求出方程的根。
（2）牛顿迭代法
牛顿迭代法迭代法迭代公式为：
a）.对单根条件下，牛顿迭代法平方收敛，牛顿迭代法程序框图如下所示：
b）.对重根条件下，此时迭代公式修改为
此时，牛顿迭代法至少平方收敛。
（3）割线法
割线法的迭代公式为：
]
割线法是超线性收敛，其程序流程图为
（4）斯蒂芬森迭代法
斯蒂芬森迭代法的计算格式为
斯蒂芬森迭代法程序流程图与牛顿迭代法类似。
2.2迭代结果
分别取初试值为，运用牛顿迭代法、割线法（）、，斯蒂芬森法求得根如下（误差精度，数据结够为（formatlong））
牛顿迭代法
割线法
斯蒂芬森法
未考虑重根
考虑重根
1
0.7524647887323
1.0147887323943
0.6999763241023
0.9895610856203
2
0.7954236245124
1.0000363575693
0.6999526372068
0.9999848725311
3
0.8306913539722
1.0000000002203
0.7524357432815
0.9999999999682
4
0.8597145503620
1.0000000000000
0.7784770276292
5
0.8836487289229
0.8074195637804
6
0.9034211164477
0.8305334214151
7
0.9197795580534
0.8514440083581
……
64
0.9999958165495
0.9998927397768
85
0.9999924537033
表1时，不同迭代方法迭代结果
1.8529411764705
1.1176470588235
1.9677760435023
0.9186980626830
1.7258825923114
1.0022487868386
1.9406643838244
1.6163939446397
1.0000008425224
1.5223105411680
1.0000000000001
1.4416990812267
1.3728310917166
1.5849504641737
1.3141624872713
65
1.0000085397619
1.0003773276313
68
1.0000049419954
95
1.0000075015608
表2时，不同迭代方法迭代结果
2.3分析
通过表1、表2不难分析出如下结论：
（1）牛顿法、割线法、斯蒂芬森法本题下均收敛，可以求出误差范围内的根。
（2）对于重根情况下，牛顿法的收敛速度会变慢变为线性收敛。
为解决这一问题，可以对牛顿迭代公式进行修正，此时至少能够保证迭代是平方收敛，使得求解速度大大加快。
（3）割线法的收敛速度小于牛顿法。
（4）运用斯蒂芬森法可以对牛顿迭代法进行加速。
三、本次实验的难点分析
1、本次试验是数值分析课程的第一次试验，重温了matlab的基本语法与操作。
由于本科时使用matlab主要是电力系统仿真即主要运用simlinuk模块，对于matlab中自定义函数的建立与运用还有欠缺。
这是本次试验中耗时较多的部分。
2、本次试验，原设计在牛顿迭代法程序中加入判断导数值是否是个小值（导数值为零）的过程，以防止分母出现较小的值。
但是程序运行过程中导数值始终会接近于零而跳出，得不出解。
为此不得不放弃这一思路。
四、参考文献
[1]谷根代，杨晓忠等．数值分析[M]．北京：
科学出版社，2011
[2]马东升，雷永军．数值计算方法[M]．北京：
机械工业出版社，2008
[3]周路，陈渝等译．数值方法（Matlab版）[M]．北京：
电子工业出版社，2010
实验题目插值与拟合
1掌握插值、拟合的基本概念、方法；
2理解牛顿插值法、三次样条插值原理、方法、步骤，理解最小二乘法的原理方法，会用计算机进行编程求解，懂得Runge现象及预防措施；
1区间作等距划分：
以为节点对函数进行插值逼近。
（1）分别取，用牛顿插值对进行逼近，并在同一坐标系下作出函数的图形，进行比较，写出插值函数对的逼近程度与节点个数之间的关系，并分析原因；
（2）使用三次样条插值对进行逼近，在同一坐标系下画出图形，观察样条插值函数对的逼近程度与节点个数的关系；
（3）整体插值有和局限性？
如何避免？
2已知一组数据如下，求其拟合曲线。
8
9
10
11
14
16
18
19
106.42
108.2
109.5
110
109.93
110.49
110.59
110.6
110.76
111
111.2
（1）求以上数据形如的拟合曲线及其平方误差；
（2）求以上数据形如的拟合曲线及其平方误差；
（3）通过观察
（1）、
（2）的结果，写出你对数据拟合的认识。
2.1.1插值
（1）牛顿插值法
Lagrange插值公式结构紧凑，便于理论分析。
利用插值基函数也容易到插值多项式的，Lagrange插值公式的缺点是，当插值节点增加，或其位置变化时，全部插值基函数均要随之变化，从而整个插值公式的结构也发生变化，这在实际计算中是非常不利的。
下面引入的Newton插值公式可以克服这个缺点。
Newton插值多项式可以灵活地增加插值节点进行递推计算。
该公式形式对称，结构紧凑，因而容易编写计算程序。
Newton插值多项式为：
其中：
称之为f（x）在上的k阶均差（或差商）。
实际计算中，可以由差商表计算。
，如下所示：
表1差商表
牛顿迭代法程序流程图如下所示：
（2）三次样条插值
分段低次插值光滑性较差，对于像高速飞机的机翼形线，船体放样等型值线往往要求有二阶光滑度，即有二阶连续导数。
三次样条函数
记为它是定义在区间[a,b]上的函数，满足以下两个条件：
1）.S（x）在每一个小区间上是一个三次多项式函数；
2）.在整个区间[a,b]上，其二阶导数存在且连续。
即在每个节点处的二阶导数连续。
给定函数在n+1个节点处的函数值为。
要求：
求一个三次样条函数，使其满足：
如何确定三次样条函数在每一个小区间上的三次多项式函数的系数呢？
这里只简介确定系数的思想。
对于三次样条插值呢，每一个区间上的三次函数的四个参数，在该区间上由两个端点的函数值只能够产生两个方程，仅此不足以唯一确定四个参数。
注意到三次样条函数对整体光滑性要求，其二阶导数存在且连续，从全局的角度上考虑参数个数与方程个数的关系如下：
参数：
每个小段上4个，n个小段共计4n个。
方程：
1）每个小段上由给定函数值得到2个，n个小段共计2n个；
2）光滑性要求每一个内部节点的一阶二阶导数连续，得出其左右导数相等，因此，每个节点产生2个方程，共计2（n-1）个。
现在得到了4n-2个方程，还差两个。
为此，常用的方法是对边界节点除函数值外附加要求，这就是所谓的边界条件。
需要两个，正好左右两个端点各一个。
常用如下三类边界条件：
第一边界条件
给定两个边界节点的一阶导数值：
m0,mn,即：
（x0）=m0,（xn）=mn。
第二边界条件
给定两个边界节点的二阶导数值：
M0,Mn,即：
（x0）=M0,
（xn）=Mn。
特别地，当M0和Mn都为零时，称为自然边界条件。
周期性边界条件
在两个边界的函数值，一阶导数值以及二阶导数值均相等：
即
·
（x0）=（xn）；
（x0）=（xn）。
2.1.2拟合
在科学实验或者统计研究中常常需要从一组测得的数据中去求得自变量和因变量之间的一个近似解析表达式使得误差最小,这就是给定个点进行曲线拟合的问题。
曲线拟合常用次多项式来进行拟合，即
求解下面方程就可进行曲线拟合：
对于指数形式的拟合曲线，可以令，，则有，用上式依然可以解除
2.2数据处理与问题分析
2.2.1插值
（1）运用所编程序，分别求出时，得出插值函数，并在同一坐标系下画出原函数与插值函数图像，分别如下所示：
图1n=1时原函数与插值函数图象
图2n=5时原函数与插值函数图象
图3n=10时原函数与插值函数图象
图4n=20时原函数与插值函数图象
图5n=25时原函数与插值函数图象
（2）运用所编程序，求出下三次样条插值函数，并在同一坐标系下作出原函数图象与插值函数图象。
如下图所示：
图6n=1时原函数与插值函数图象
图7n=5时原函数与插值函数图象
图8n=10时原函数与插值函数图象
2.2.2拟合
（1）形如的拟合，根据所编matlab程序，得到拟合方程为
（2）形如的拟合，根据所编matlab程序，得到拟合方程为
（3）将
（1）、
（2）与原数据散点图图象做在同一坐标系中，如下如所示：
通过此图，可以发现，指数形式的拟合函数更加和原数据走势拟合，然而不管哪一种拟合形式，都不能保证所有的数据点都在拟合曲线上。
所以，拟合可以表现出数据的真题趋势，而且，通过拟合可以合理的排除一些误差点。
2.3问题分析
2.3.1插值
（1）由2.2.1中时插值函数与原函数的图象分析可知，插值点越多，插值
升级会员 