数值计算初步 - 计算机科学与技术.doc

数值计算初步 晋中学院 宫建平编 目 录 第1章 误差和MARLAB的计算精度 1 1.1 误差 1 1.1.1 误差的来源 1 1.1.2 有关误差的一些概念 2 1.2 MATLAB中的数值计算精度 3 1.2.1 浮点数及其运算特点 3 1.2.2 MATLAB中的数值计算精度 4 1.3 设计算法的若干原则 6 1.3.1 算法的数值稳定性 6 1.3.2 设计算法的若干原则 6 第2章 求解非线性方程f(x)=0 9 2.1 求解f(x)=0的MATLAB符号法 9 2.2 求方程f(x)=0数值解的基本方法 11 2.2.1 求实根的二分法原理 11 2.2.2 迭代法 12 2.2.3 切线法 13 2.2.4 割线法(弦截法) 15 2.3 方程f(x)=0数值解的M ATLAB实现 15 2.3.1 代数方程的求根指令roots 15 2.3.2 求函数零点指令fzero 16 2.4 求解非线性方程组数值解的迭代法 18 2.5 求方程组数值解的指令 19 第3章 求解线性代数方程组的直接法 21 3.1 线性代数方程组求解概论 21 3.1.1 线性代数方程组的矩阵表示 21 3.1.2 线性代数方程组解的性质 21 3.2 恰定线性代数方程组求解 22 3.2.1 克莱姆法则 22 3.2.2 高斯消去法 23 3.3 矩阵的三角分解 24 3.3.1 高斯消去法和三角矩阵 24 3.3.2 矩阵的三角分解 25 3.4 线性代数方程组数值解和矩阵三角分解的MATLAB实现 26 3.4.1 齐次线性代数方程组求解指令 26 3.4.2 求解非齐次线性代数方程组的MATLAB方法 27 3.4.3 矩阵分解指令 31 第4章 求解线性代数方程组和计算矩阵特征值的迭代法 35 4.1 求解线性代数方程组的迭代法 35 4.1.1 迭代法的基本原理 35 4.1.2 雅可比迭代法 35 4.1.3 赛德尔迭代法 36 4.1.4 迭代法的敛散性 37 4.2 方阵特征值和特征向量的计算 40 4.2.1 方阵特征方程的求解 40 4.2.2 计算特征值和特征向量的迭代法 40 4.3 矩阵一些特征参数的MATLAB求算 42 4.3.1 求方阵特征值的有关指令 42 4.3.2 矩阵的正交三角分解指令 45 4.3.3 计算范数和矩阵谱半径的指令 46 第5章 插值法和数据拟合 49 5.1 多项式插值 49 5.1.1 代数多项式插值的基本原理 49 5.1.2 两种常见插值法 50 5.1.3 插值多项式的误差估计 53 5.2 分段三次插值和三次样条函数 53 5.2.1 分段三次Hermte插值 54 5.2.2 三次样条插值的基本原理 54 5.2.3 三次样条插值函数的一种具体推导方法 55 5.3 插值法在MATALAB中的实现 57 5.3.1 一元函数的插值(查表)指令 57 5.3.2 三次插值和三次样条插值指令 57 5.4 数据的曲线拟合 59 5.4.1 数据拟合的最小二乘法 59 5.4.2 超定方程组的最小二乘解 60 5.5 多项式运算在MATLAB中的实现 61 5.5.1 多项式及其系数向量 61 5.5.2 多项式运算 62 5.6 曲线拟合在MATLAB中的实现 64 5.6.1 数据的多项式曲线拟合 64 5.6.2多项式数据拟合应用的扩充 65 第6章 数 值 积 分 67 6.1 计算积分的M ATLAB符号法 67 6.1.1常微分方程的MATLAB符号表示法 67 6.2 牛顿_柯特斯求积公式 70 6.2.1 牛顿_柯特斯求积公式推导 70 6.2.2 牛顿_柯特斯求积公式的误差估计 72 6.3 几个低次牛顿_柯特斯求积公式 73 6.3.1 矩形求积公式 73 6.3.2 梯形求积公式 74 6.3.3 抛物线求积公式 74 6.4 复合求积公式及其M ATLAB 实现 75 6.4.1 复合矩形求积法及其MATLAB实现 75 8.4.2 复合梯形求积法及其MATLAB实现 76 6.5 变步长复合求积及其M ATLAB 实现 78 6.5.1 复合抛物线求积公式 78 6.5.2 变步长复合抛物线求积公式 78 6.5.3 求数值积分的指令quad和quadl 79 第7章 常微分方程初值问题的数值解 83 7.1求解常微分方程的MATLAB符号法 83 7.1.1常微分方程的MATLAB符号表示法 83 7.1.2 求解常微分方程的符号法指令dsolve 84 7.2 常微分方程数值解的基本原理 86 7.2.1 求常微分方程数值解的基本原理 86 7.2.2 泰勒展开法 87 7.2.3 龙格_库塔法 88 7.3 常微分方程初值问