245_362677_数值分析实验报告.doc

数值分析课程设计报告书 系部名称 ： 应 用 数 理 系 学生姓名 ： 专业名称 ： 信息与计算科学 班 级 ： 信息0702 时间 ： 2009年6月8日至2009年6月19日 实验一 实验内容 1.用三次样条插值的三弯矩法,编制第一与第二种边界条件的程序.已知数据如下： 0.2 0.4 0.6 0.8 1.0 0.9798652 0.9177710 0.8080348 0.6386093 0 求的三次样条插值函数满足： (1)自然边界条件 (2)第一种边界条件 要求输出用追赶法解出的弯矩向量(,,,)及(i=0,1,2,3,4,5,6,7,8)的值.并画出的图形. 实验原理简述 我们先利用一阶导数在节点（i=1,2,…n-1）上的连 性以及边界条件，列出确定二阶导数（i=0,1,2,…n-1）的线性方程组（三弯矩方程），并由此解出然后用表达是即可 （三弯矩方程） ， ， . 其中有（）个未知数，而方程只有（n-1）个，当满足第一种边界条件时，可的另两个方程 ， 如果令，将上述方程综合后的一下矩阵形式： = 可以证明此方程组满足追赶法的条件，我们用追赶法可得M的 值，将其带入公式即得 对第二种边界条件，直接的端点方程并且令，则又得三弯矩方程同理即可求得解。 实验结果 实验分析 某些实际工程中不仅要求曲线连续，而且要求曲线的曲率也连续，通过上述结果我们可以看出三次样条完全满足此要求。且在解方程中我们采用了追赶法，其也有效地解决了高斯消去法中的舍入误差积累的问题。 实验二 实验内容 编制以离散点的正交多项式为基的最小二乘拟合程序,并用于对下列数据做三次多项式最小二乘拟合. -1.0 -0.5 0.0 0.5 1.0 1.5 2.0 -4.447 -0.452 0.551 0.048 -0.447 0.549 4.552 取权1,求出拟合曲线,输出,,及平方误差,并画出的图形. 实验原理简述 由多项式拟合公式：求的a我们即的拟合系数，带入多项式即可，在一般多项式拟合中A的第k列即为的(k-1)次方，而在正交多项式拟合中，尤其性质我们有A的非对角线元素为0，所以我们有迭代公式 即的其结果。但在此题中我采用的勒让德正交多项式做基，由于勒让德的正交性是在[-1 1],而我们的离散点并不全在其内，所以我们不能再用上述公式进行简化计算了，相反由于对于勒让德正交多项式而言的A也就不在具有向一般多项式那样的性质了，所以我们必须循环赋值这一就复杂问题了，但其计算结果还是较满意的。 实验结果 一般多项式结果 计算的平方和误差为2.176190481775883e-005。 勒让德正交多项式做基 实验分析 通过结果我们可以看出一般多项式拟合与勒让德正交多项式做基的结果误差都较小，效果较好，并且进一步增加多项式得系数得知一般多项式，的拟合结果将出现较大误差，这是因为要解方程，而随着系数的增大，希尔伯特矩阵是病态的，所以解将出现大的误差拟合自然偏差较大。而勒让德正交多项式有于其在[-1 1]上的正交性，在此段上的解并没有出现较大的误差相反结果还是很理想的，对于在其它段上的我们只要做一个变换，则结果也是较好的。 实验三 实验内容 给出积分 ①, ②, ③ (1)运用龙贝格求积公式计算上述积分I的值,要求到时 结束,输出T表及I的近似值. (2)用5点高斯求积公式及复化3点高斯求积公式计算上述积分,并输出I的近似值. 3)分析比较各种计算结果. 实验原理简述 龙贝格求积公式 高斯求积 复化求积 实验结果 结果为I=0.422724. ① 五点高斯法： 0.422720775202602 复化3点高斯：0.422735763060519。 ② 五点高斯法：1.542126458733056 复化3点高斯：0.925275412602127 龙贝格求积公式：1.542125687670212 ③ 五点高斯法：0.202732641800508 复化3点高斯：-0.025024125653057 龙贝格求积公式：0.203470023940993 实验分析 有实验结果显然可以看出龙贝格的精度最高，五点高斯法次之，最后是复化3点高斯法，由于龙贝格有误差限的控制，所以精度有r控制。对于高斯法，我们知道由于其自身的构造性保证了其具有（2n+1）的精度。而复化3点高斯法我用的是2段的三点高斯，效果自然差些。而对于复化的我们也可以给定误差限，而不停地分段来增加精度，也可已达到要求的精度，只不过迭代的次数多些罢了。 实验四 实验内容 比较求一阶导数的数值方法,给出函数.利用某距离点函数值,必要时给定端点导数值,分别用中心差分,数值积分求导和理查森外推