数值分析课程设计研究报告(分).docVIP

数值分析课程设计报告 设计题1、2、3、5 学院、系： 专 业： 姓 名： 学 号： 任课教师： 提交日期： 电子邮箱： [设计题一] 3 1.1问题分析与设计思路 3 1.2程序清单 4 1.4 结果分析 7 1.5设计总结 7 [设计题二] 8 2.1问题分析与设计思路 8 2.2程序清单 8 2.3 运行结果 10 2.4结果分析与设计总结 10 [设计题三] 11 3.1问题分析与设计思路 11 3.2程序清单 11 3.3 运行结果 13 3.4结果分析与设计总结 13 [设计题五] 14 4.1问题分析与设计思路 14 4.2程序清单 15 4.3 运行结果 20 4.4结果分析 21 【数值分析课程设计总结】 22 [设计题一] 设计实验验证Hilbert矩阵的病态性。 1.1问题分析与设计思路 在求解任何反问题的过程中通常会遇到病态矩阵问题，而且病态矩阵问题还未有很好的解决方法，尤其是长方形、大型矩阵。目前主要有Tikhonov、奇异值截断、奇异值修正等方法。　 求解方程组时对数据的小扰动很敏感的矩阵就是病态矩阵。解线性方程组Ax=b时，若对于系数矩阵A及右端项b的小扰动δA、δb,方程组(A+δA)χ=b+δb的解χ与原方程组Ax=b的解差别很大，则称矩阵A为病态矩阵。方程组的近似解χ一般都不可能恰好使剩余r=b-Aχ为零，这时χ亦可看作小扰动问题Aχ=b-r(即δA=0，δb=-r)的解,所以当A为病态时,即使剩余很小,仍可能得到一个与真解相差很大的近似解。 因此，设计思路如下： 令x0=（1,1…..1），计算出b=Hx0，求出b，然后再用高斯消去法球解Hx=b，得到近似解x，然后利用标准差： 比较x与x0之间的误差。截图是取了几个n去计算，看一下随着n的增大误差的变化情况。 1.2程序清单 qm1.m gauss_liezhu1.m 1.4 结果分析 按照N的递增顺序取了9个误差数据，制成散点折线图如上所示。 由此可以看出，此矩阵求解方程组时对数据的小扰动很敏感实验验证Hilbert矩阵的病态性令x0=（1,1…..1），计算出b=Hx0，求出b，然后再用高斯消去法球解Hx=b，得到近似解x，然后利用标准差 比较x与x0之间的误差。 [设计题二] 1225年，达芬奇研究了方程并得到它的一个近似根。没有人知道他用什么方法得到它。设计两种方法去计算，并比较这两种方法。 2.1问题分析与设计思路 。第二种方法是用牛顿迭代法。 二分法 优点：条件和方法简单(只要求f(x)连续即可)，方法收敛；缺点：收敛速度慢，不易求偶数重根（如图）. Newton 迭代法迭代公式 2.2程序清单 2.3 运行结果 二分法： Newton迭代法 2.4结果分析与设计总结 通过二分法与Newton迭代法得出的答案相同。 确定求方程近似根的三种方法 翻书了解编程步骤 总结本章重点知识： 1.熟悉区间二分法； 2.熟悉迭代法的建立，会使用收敛定理； 3.熟悉Newton迭代法及其几何意义; 4.熟悉收敛阶的定义； 5.熟悉Newton迭代法的收敛阶; [设计题三] 某飞机头部的光滑外形曲线的型值点坐标由下表给出: 0 1 2 3 4 5 6 7 8 9 10 0 70 130 210 337 578 776 1012 1142 1462 1841 0 57 78 103 135 182 214 244 256 272 275 试建立其合适的模拟曲线（未必是用拟合方法），并求在点x＝100，250，400，500，800处的函数值y及一阶、二阶导数值y’，y”。绘出模拟曲线的图形。 3.1问题分析与设计思路 曲线拟合的一般步骤 绘制散点图，选择合适的曲线类型　　一般根据资料性质结合专业知识便可确定资料的曲线类型，不能确定时，可在方格坐标纸上绘制散点图，根据散点的分布，选择接近的、合适的曲线类型。　　 进行变量变换使变换后的两个变量呈直线关系。　　 按最小二乘法原理求线性方程和方差分析　 （四）将直线化方程转换为关于原变量X、Y的函数表达式　三次定义:函数S(x)∈C2[a,b] ，且在每个小区间[ xj,xj+1 ]上是三次多项式，其中　a =x0 x1..