1、实验一实验一 舍入误差与数值稳定性 实验目的与要求:实验目的与要求:通过上机编程,复习巩固以前所学程序设计语言及上机操作指令 通过上机计算,了解舍入误差所引起的数值不稳定性 实验内容:实验内容:通过正反两个实例的计算,了解利用计算机进行数值计算中舍入误差所引起的数值不稳定性,深入理解初始小的舍入误差可能造成误差积累从而对计算结果的巨大影响。通过实际编程,了解运算次序对计算结果的影响,了解实数运算符合的结合律和分配律在计算机里不一定成立。概要概要舍入误差在计算方法中是一个很重要的概念。在实际计算中,如果选用了不同的算法,由于舍入误差的影响,将会得到截然不同的结果。因此,选取稳定的算法,在实际计算
2、中是十分重要的。N2211 311-12 2 N N+1jSj设,已知其精确值为()1)编制按从大到小的顺序计算SN的程序 2)编制按从小到大的顺序计算SN的程序 3)按2种顺序分别计算S1000,S10000,S30000,并指出有效位数 实验题目:实验题目:实验二 插值法 目的与要求:目的与要求:熟悉拉格朗日插值多项式和牛顿插值多项式,注意其不同特点;实验内容:实验内容:通过拉格朗日插值和牛顿插值多项式的两个实例的计算,了解两种求解方法,分析他们的优缺点。按下列数据 作五次插值,并求x=0.46,x=0.55,x=0.60时的函数近似值 iix0.300.420.500.580.660.7
3、2y1.044031.084621.118031.156031.198171.23223 实验题目:实验题目:编写一个用牛顿前插公式计算函数值的程序,要求先输出差分表,再计算x点的函数值,并应用于下面的问题:求x=21.4时的三次插值多项式的值 iix2021222324y1.301031.322221.342421.361731.38021实验三数值积分 目的与要求:目的与要求:通过实际计算体会各种方法的精确度;会编写用龙贝格算法求定积分的程序。实验内容:实验内容:通过实际计算体会各种方法的精确度并且会编写用龙贝格算法求定积分的程序 龙贝格算法龙贝格算法 用事后估计法控制精度 ),)(102
4、/12nkknnnxfhTT;/)(nabhnnkhkx)2/1(2/1)(3/122nnnnTTTS)_(15/122nnnnSSSC)(63/122nnnnCCCR|2nnRR 实验题目:实验题目:用龙贝格方法上机计算(按学号尾数mod 3 选择):dxxxxdxxdxex)sin25.01/(sin)3sin)cos2sin()2)12/022/028.002实验四 常微分方程数值解法 目的与要求:熟悉求解常微分方程初值问题的有关方法和理论,主要是欧拉法和改进欧拉法;会编制上述方法的计算程序,针对实习题编制程序,并上机计算其所需要的结果;通过对各种求解方法的计算实习,体会各种解法的功能,
5、优缺点及适用场合,会选取适当的求解方法。实验内容:实验内容:熟悉求解常微分方程初值问题的有关方法和理论,主要是欧拉法和改进欧拉法,并且通过对各种求解方法的计算实习,体会各种解法的功能,优缺点及适用场合,会选取适当的求解方法。实验题目:实验题目:202yxyy 30 x实验五方程求根 目的与要求:通过对二分法和牛顿迭代法作编程练习和上机运算,进一步体会它们在方程求根中的不同特点;比较二者的计算速度和计算精度。实验内容:通过对二分法和牛顿迭代法作编程练习和上机运算,进一步体会它们在方程求根中的不同特点 二分法二分法 算法:算法:给定区间a,b,并设与符号相反,取为根的容许误差,为的容许误差。(1)
6、令c=(a+b)/2(2)如果(c-a)或,则输出,结束;否则执行(3)(3)如果,则令;否则则令,重复(1),(2),(3)。算法:算法:给定初值 ,为根的容许误差,为 的容许误差,N为迭代次数的容许值。(1)如果 =0或迭代次数大于N,则算法失败,结束;否则执行(2)。(2)计算 =-(3)若 或 ,则输出 ,程序结束;否则执行(4)。(4)令 =,转向(1)。牛顿迭代法牛顿迭代法 x0)(xf)(xfx1x0)()(0 xxoffxx01)(1xfx1x0 x1 实验题目:实验题目:1.用二分法求方程 在区间1,1.5上的根,要求求出具有3位有效数的近似根。104)(23xxxf2.用牛
7、顿法求方程x3-3x-1=0在x=2附近的根实验六线性方程组数值解法 目的与要求:熟悉求解线性方程组的有关理论和方法;会编制列主元消去法、LU 分解法、雅可比及高斯塞德尔迭代法德程序;通过实际计算,进一步了解各种方法的优缺点,选择合适的数值方法。实验内容:会编制列主元消去法、LU 分解法、雅可比及高斯塞德尔迭代法德程序,进一步了解各种方法的优缺点。实验题目:实验题目:1.用列主元消去法解方程组 2.用LU分解法求解方程组A Ax=b b 035.3643x.5072x.1835x.2137.2623x.43712x347x.11.1833.555x2.304x0.101x321321321A=
8、,b=3812659732745812214911427 3.用雅可比迭代法解方程组 4.用高斯-塞尔德迭代法解方程组 16321382825321321321xxxxxxxxx3612363311420238321321321xxxxxxxxx实验七矩阵的特征值与特征向量的计算 目的与要求:领会求矩阵特征值及特征向量的幂法的理论及其方法;会编制幂法的计算程序,并用来计算有关问题。实验内容:编制幂法的计算程序,并用来计算有关问题 实验题目:实验题目:求矩阵求矩阵 的最大特征值和特征向的最大特征值和特征向 量,当特征值有量,当特征值有3位小数稳定时迭代终止位小数稳定时迭代终止 210121012A