《数值分析》教学大纲（新）
（Numerical Analysis）
适用专业：全校工科专业课程性质：学位课
学时数：48学时学分数：3学分
课程号：开课学期：第1学期
大纲执笔人：杨帆大纲审核人：欧志英
一、课程的地位和教学目标
数值分析是全校研究生的一门重要学位课，工科研究生应该掌握数值分析的基本知识和方法，
主要介绍用计算机解决数学问题的数值计算方法及其理论。内容新颖，起点较高，并加强了数值试验和程序设计环节。通过本课程的学习，使学生熟练掌握各种常用的数值算法的构造原理和过程分析，提高算法设计和理论分析能力，并且能够根据数学模型，提出相应的数值计算方法编制程序在计算机上算出结果。力求使学生掌握应用数值计算方法解决实际问题的常用技巧。
二、课程教学内容和基本要求
（一）数值分析引论（Introduction）（3学时）
教学重点、难点：
重点：理解有效数字和误差的概念，理解离散化方法与递推化方法。
难点：有效数字和误差的关系，递推化方法的编程实现。
教学内容和基本要求
1、数值分析课程的特点
了解数值分析的四个特点。即：面向计算机，可靠的数学理论分析，良好的计算复杂性，可进行数值实验。
2、有效数字；
理解有效数字概念，掌握有效数字确定的方法；
3、误差的概念，和误差的求解。
理解误差的基本概念，了解误差的各种来源，会误差估计；基本要求：
（二）非线性方程求根(Solving Nonlinear Equations)（5学时）
教学重点、难点：
重点：二分法，迭代法，牛顿法和割线法求解非线性方程。
难点：迭代法的收敛性和误差分析以及迭代加速。
教学内容和基本要求
1、二分法
掌握二分法及其收敛性。
2、迭代法的基本理论
迭代法的基本思想和收敛性判别定理以及误差分析。
3、迭代加速技术。
掌握迭代加速技术的基本思想和Aitken加速公式。
4、牛顿法及其割线法
掌握牛顿法和割线法及其收敛性。
（三）线性方程组的直接解法(Direction Solving Linear Systems) (5学时）
教学重点、难点：
重点：线性方程组的高斯消元法以及LU分解法。
难点：利用向量和矩阵的范数进行误差分析。
教学内容和基本要求
1、高斯消元法
掌握顺序高斯消元法和选列主元高斯消元法
2、LU分解法
掌握LU分解法。
3、误差分析
掌握向量和矩阵的范数，利用这两类范数进行误差分析。
（四）线性方程组的迭代解法（Iterative Method for Solving Linear Equations）(5学时）
教学重点、难点：
重点：三种经典迭代法的构造。
难点：迭代法的收敛性和误差分析
教学内容和基本要求
1、迭代法的一般理论
掌握迭代公式的构造和迭代法的收敛性和误差分析。
2、三种经典迭代法
掌握雅可比和高斯-塞德尔以及逐次超松子这三类迭代法求解线性方程组的迭代
格式。
（五）插值法与最小二乘拟合(Interpolation, Curve Fitting and Polynomial Approximation)（10学时）
教学重点、难点：
重点: 拉格朗日和牛顿插值法以及最小二乘法。
难点：三次样条插值
教学内容和基本要求
1、插值法的基本理论
了解插值多项式的概念，插值基函数以及插值误差。
2、拉格朗日插值法
掌握利用拉格朗日基函数构造拉格朗日插值多项式以及误差分析。
3、牛顿插值法
掌握差商的概念，利用差商构造牛顿插值公式。
4、Hermit插值及分段插值
掌握构造Hermit插值的方法。
5、三次样条插值
掌握三次样条插值函数的构造方法。
6、曲线拟合的最小二乘法
掌握最小二乘方法和最佳平方逼近方法。
（六）数值积分和数值微分(Numerical Differentiation and Integration)（10学时）
教学重点、难点：
重点：复化的求积公式和变步长求积公式。
难点：龙贝格求积公式
教学内容和基本要求
1、几个常用的求积公式
掌握插值型求积公式以及代数精度。
2、复化求积公式
掌握复化梯形公式和复化的新普生公式。
3、高斯求积公式
掌握高斯求积公式。
4、数值微分
掌握两点和三点公式。
（七）常微分方程数值解法(Numerical Solutions of Ordinary Differential Equations)（10学时）
教学重点、难点：
重点：欧拉差分格式和龙格-库塔方法。
难点：差分格式的收敛性与稳定性分析。
1、欧拉方法及其改进
掌握欧拉公式和改进的欧拉公式。
2、龙格-库塔公式
掌握龙格-库塔方法的基本思想以及龙格库塔公式。
3、收敛性与稳定性
掌握差分格式的收敛性与稳定性。
四、本课程教学建议
（一）学生应先修高等数学，线性代数，Matlab语言及程序设计。
（二）为使学生能做到对所学知识融会贯通，一方面应加强练习尤其是matlab编程方面的练习，另一方面应在课程的后半段，适时引进一些实际问题让学生建立数学模型，然后用数值方法求解。
五、本课程评价方式
闭卷考试。
六、建议教材和教学参考书
[1] 蔡大用，白峰杉《高等数值分析》，清华大学出版社，1996.
[2] 马昌风，《现代数值分析》. 国防工业出版社, 2015.
[3] Jaan Kiusalaas. Numerical Methods in Engineering with Matlab. Cambridge: Cambridge University Press, 2005.
[4] Curtis F Gerald, Patrick O Wheatley. Applied Numerical Analysis. Sixth Edition. New York: Addison-Wesley, 1999.
[5] Huang Fang Lun. Numerical Methods Using Matlab. Fourth Edition. Bei Jing: Publishing House of Electronics Industry, 2009.
[6] 朱晓临. Numerical Analysis. 合肥: 合肥工业大学出版社, 2010.
《数值分析课程设计》教学大纲 课程编号：1512110303 课程名称： 数值分析课程设计 周数/学分：3/3 先修课程：《数值分析》 适用专业： 信息与计算科学 开课教研室：应用数学教研室 一、目的与要求： 《数值分析课程设计》是实践性教学内容之一，是《数值分析》课程的辅助教学过程，是信息与计算科学专业的必修课。通过设计，使学生深化对所学理论知识的理解，掌握数值计算方法的程序设计能力，初步具备解决实际数值计算问题的能力。 二、课程设计内容： 1．掌握数值分析的基本内容。误差的基本概念，插值与拟合，数值积分，线性代数方程组的解法，非线性方程求根，常微分方程初值问题的数值解法。 2．对每部分内容设计一定难度的问题，要求学生对问题进行分析，确定解决方案。 3．进行模拟与仿真，进行结果分析，编写课程设计报告 三、课程设计步骤与方法 1．教师向学生讲解课程设计目的和要求，补充相关基本知识，布置课程设计任务。 2．学生查找资料，编程、调试程序。本步骤是课程设计的核心内容之一，要求学生分析算法，写出相应程序，并对结果进行解释 3．撰写课程设计报告。 四、课程设计的基本要求 1．算法说明正确无误，图表符合技术规范要求。 2．毎生一台计算机，要求学生使用Matlab软件或Mathematica软件编写相关程序。 3．按要求完成一篇的课程设计报告。 4．课程设计的方式：以集中学习为主；独立完成课程设计阶段规定的全部工作任务。 五、课程设计进度表 序号 内 容 所用时间 1 教师讲解，布置任务 1天 2 学生编写程序并撰写设计报告 11天
3 教师反馈意见，学生修改设计报告 3天 合计 15天 六、课程设计考核方式 平时设计环节中的表现占总成绩30％，课程设计报告和软件运行情况占总成绩70％。 执笔：赵国喜 审定：朱耀生 梁桂珍
数值分析试题 一、 填空题（2 0×2′） 1. ?? ????-=? ?????-=32,1223X A 设x =是精确值x *=的近似值，则x 有 2 位 有效数字。 2. 若f (x )=x 7－x 3＋1，则f [20,21,22,23,24,25,26,27]= 1 ， f [20,21,22,23,24,25,26,27,28]= 0 。 3. 设，‖A ‖∞＝___5 ____，‖X ‖∞＝__ 3_____， ‖AX ‖∞≤_15_ __。 4. 非线性方程f (x )=0的迭代函数x =?(x )在有解区间满足 |?’(x )| <1 ，则使用该迭代 函数的迭代解法一定是局部收敛的。 5. 区间[a ,b ]上的三次样条插值函数S (x )在[a ,b ]上具有直到 2 阶的连续导数。 6. 当插值节点为等距分布时，若所求节点靠近首节点，应该选用等距节点下牛顿差商 公式的 前插公式 ，若所求节点靠近尾节点，应该选用等距节点下牛顿差商公式的 后插公式 ；如果要估计结果的舍入误差，应该选用插值公式中的 拉格朗日插值公式 。 7. 拉格朗日插值公式中f (x i )的系数a i (x )的特点是：=∑=n i i x a 0)( 1 ；所以当 系数a i (x )满足 a i (x )>1 ，计算时不会放大f (x i )的误差。 8. 要使 20的近似值的相对误差小于%，至少要取 4 位有效数字。 9. 对任意初始向量X (0)及任意向量g ，线性方程组的迭代公式x (k +1)=Bx (k )+g (k =0,1,…)收 敛于方程组的精确解x *的充分必要条件是 ?(B)<1 。 10. 由下列数据所确定的插值多项式的次数最高是 5 。 11. 牛顿下山法的下山条件为 |f(xn+1)|<|f(xn)| 。 12. 线性方程组的松弛迭代法是通过逐渐减少残差r i (i =0,1,…,n )来实现的，其中的残差 r i ＝ (b i -a i1x 1-a i2x 2-…-a in x n )/a ii ，(i =0,1,…,n )。 13. 在非线性方程f (x )=0使用各种切线法迭代求解时，若在迭代区间存在唯一解，且f (x )
上海交通大学致远学院 《常微分方程和偏微分方程的数值解法》教学大纲 一、课程基本信息 课程名称（中文）：常微分方程和偏微分方程的数值解法 课程名称（英文）：Numerical Methods for Ordinary and Partial Differential Equations 课程代码：MA300 学分 / 学时：4学分 / 68学时 适用专业：致远学院与数学系相关专业 先修课程：偏微分方程，数值分析 后续课程：相关课程 开课单位：理学院数学系计算与运筹教研室 Office hours: 每周二19：00—21：00，地点：数学楼1204 二、课程性质和任务 本课程是致远学院和数学系应用数学和计算数学方向的一门重要专业基础课程，其主要任务是通过数学建模、算法设计、理论分析和上机实算“四位一体”的教学方法，使学生掌握常微分方程与偏微分方程数值解的基本方法、基本原理和基本理论，进一步提升同学们利用计算机解决实际问题的能力。在常微分方程部分，将着重介绍常微分方程初值问题的单步法，含各类Euler方法和Runge-Kutta方法，以及线性多步法。将简介常微分方程组和高阶常微分方程的数值方法。在偏微分方程部分，将系统介绍求解椭圆、双曲、抛物型方程的差分方法的构造方法和理论分析技巧，对于椭圆型方程的边值问题将介绍相应变分原理与有限元方法。将在课堂上实时演示讲授的核心算法的计算效果，以强调其直观效果与应用性。本课程重视实践环节建设，学生要做一定数量的大作业。 三、教学内容和基本要求 第一部分：常微分方程数值解法 1 引论 1.1回顾：一阶常微分方程初值问题及解的存在唯一性定理
《数值计算方法》教学大纲 课程编号：MI3321048 课程名称：数值计算方法英文名称：Numerical and Computational Methods 学时： 30 学分：2 课程类型：任选课程性质：任选课 适用专业：微电子学先修课程：高等数学,线性代数 集成电路设计与集成系统 开课学期：Y3开课院系：微电子学院 一、课程的教学目标与任务 目标：学习数值计算的基本理论和方法，掌握求解工程或物理中数学问题的数值计算基本方法。 任务：掌握数值计算的基本概念和基本原理，基本算法，培养数值计算能力。 二、本课程与其它课程的联系和分工 本课程以高等数学，线性代数，高级语言编程作为先修课程，为求解复杂数学方程的数值解打下良好基础。 三、课程内容及基本要求 (一) 引论（2学时） 具体内容：数值计算方法的内容和意义，误差产生的原因和误差的传播，误差的基本概念，算法的稳定性与收敛性。 1.基本要求 （1）了解算法基本概念。 （2）了解误差基本概念，了解误差分析基本意义。 2.重点、难点 重点：误差产生的原因和误差的传播。 难点：算法的稳定性与收敛性。 3.说明：使学生建立工程中和计算中的数值误差概念。 (二) 函数插值与最小二乘拟合（8学时） 具体内容：插值概念，拉格朗日插值，牛顿插值，分段插值，曲线拟合的最小二乘法。 1.基本要求 （1）了解插值概念。 （2）熟练掌握拉格朗日插值公式，会用余项估计误差。 （3）掌握牛顿插值公式。 （4）掌握分段低次插值的意义及方法。
（5）掌握曲线拟合的最小二乘法。 2.重点、难点 重点：拉格朗日插值, 余项，最小二乘法。 难点：拉格朗日插值, 余项。 3.说明：插值与拟合是数值计算中的常用方法，也是后续学习内容的基础。 (三) 第三章数值积分与微分（5学时） 具体内容：数值求积的基本思想，代数精度的概念，划分节点求积公式（梯形辛普生及其复化求积公式），高斯求积公式，数值微分。 1.基本要求 （1）了解数值求积的基本思想，代数精度的概念。 （2）熟练掌握梯形，辛普生及其复化求积公式。 （3）掌握高斯求积公式的用法。 （4）掌握几个数值微分计算公式。 2.重点、难点 重点：数值求积基本思想，等距节点求积公式，梯形法，辛普生法，数值微分。 难点：数值求积和数值微分。 3.说明：积分和微分的数值计算，是进一步的各种数值计算的基础。 (四) 常微分方程数值解法（5学时） 具体内容：尤拉法与改进尤拉法，梯形方法，龙格—库塔法，收敛性与稳定性。 1.基本要求 （1）掌握数值求解一阶方程的尤拉法，改进尤拉法，梯形法及龙格—库塔法。 （2）了解局部截断误差，方法阶等基本概念。 （3）了解收敛性与稳定性问题及其影响因素。 2.重点、难点 重点：尤拉法,龙格－库塔法，收敛性与稳定性。 难点：收敛性与稳定性问题。 3.说明：该内容是常用的几种常微分方程数值计算方法，是工程计算的重要基础。 (五) 方程求根的迭代法（4学时） 具体内容：二分法，解一元方程的迭代法，牛顿法，弦截法。 1.基本要求 （1）了解方程求根的对分法和迭代法的求解过程。 （2）熟练掌握牛顿法。 （3）掌握弦截法。 2.重点、难点 重点：迭代法，牛顿法。
数学与计算科学学院实验报告 实验项目名称常微分方程数值解 所属课程名称数值方法B 实验类型验证 实验日期 2013.11.11 班级 学号 姓名 成绩
【实验过程】（实验步骤、记录、数据、分析） 注：以下图形是通过Excel 表格处理数据得出，并未通过MATLAB 编程序所得！ 1、1(0)1dy y x dx y ?=-++???=? 由题可知精确解为：x y x e -=+，当x=0时，y(x)=0。 h=0.1 表1 h=0.1时三个方法与精确值的真值表 图1 h=0.1时三个方法走势图 步长 Euler 法 预估校正法 经典四阶库 精确值 0.1 1.010000 1.005000 1.004838 1.249080 0.2 1.029000 1.019025 1.018731 1.055455 0.3 1.056100 1.041218 1.040818 1.091217 0.4 1.090490 1.070802 1.070320 1.131803 0.5 1.131441 1.107076 1.106531 1.176851 0.6 1.178297 1.149404 1.148812 1.226025 0.7 1.230467 1.197211 1.196586 1.279016 0.8 1.287421 1.249975 1.249329 1.335536 0.9 1.348678 1.307228 1.306570 1.395322 1.0 1.413811 1.368541 1.367880 1.458127
h=0.05（此时将源程序中i的围进行扩大，即for(i=0;i<20;i++)） 表2 h=0.05时三个方法与精确值的真值表步长Euler法预估校正法经典四阶库精确值 0.05 1.002500 1.001250 1.001229 1.011721 0.10 1.007375 1.004877 1.004837 1.024908 0.15 1.014506 1.010764 1.010708 1.039504 0.20 1.023781 1.018802 1.018731 1.055455 0.25 1.035092 1.028885 1.028801 1.072710 0.30 1.048337 1.040915 1.040818 1.091217 0.35 1.063421 1.054795 1.054688 1.110931 0.40 1.080250 1.070436 1.070320 1.131801 0.45 1.098737 1.087752 1.087628 1.153791 0.50 1.118800 1.106662 1.106531 1.176851 0.55 1.140360 1.127087 1.126950 1.200942 0.60 1.163342 1.148954 1.148812 1.226025 0.65 1.187675 1.172193 1.172046 1.252062 0.70 1.213291 1.196736 1.196585 1.279016 0.75 1.240127 1.222520 1.222367 1.306852 0.80 1.268121 1.249485 1.249329 1.335536 0.85 1.297215 1.277572 1.277415 1.365037 0.90 1.327354 1.306728 1.306570 1.395322 0.95 1.358486 1.336900 1.336741 1.426362 1.00 1.390562 1.368039 1.367880 1.458127 图2 h=0.05时三个方法走势图
华南理工大学研究生课程考试 《数值分析》试卷C （2015年1月9日） 1. 考前请将密封线内各项信息填写清楚； 所有答案请按要求填写在本试卷上； 课程代码：S0003004； 4. 考试形式：闭卷； 5. 考生类别：硕士研究生； 本试卷共八大题，满分100分，考试时间为150分钟。 一、（12分）解答下列问题： 1）设近似值0x >，x 的相对误差为δ，试证明ln x 的绝对误差近似为δ。 2）利用秦九韶算法求多项式 542()681p x x x x x =-+-+ 在3x =时的值（须写出计算形式），并统计乘法次数。 （12分）解答下列问题： 1）设()235f x x =+，求[]0,1,2f 和[]0,1,2,3f 。 2）利用插值方法推导出恒等式： 33220,0[]j j i i x j i x i j =≠=-=-∑∏ 。
（1）设{}∞ =0)(k k x q 是区间[]1,0上带权1=ρ而最高次项系数为1的正交多项式族，其中1)(0=x q ，求1()q x 和2()q x 。 （2）求形如2y a bx =+的经验公式，使它与下列数据拟合： 四、（14分）对积分()10I f x dx = ?，试 （1）构造一个以012113,,424 x x x ===为节点的插值型求积公式； （2）指出所构造公式的代数精度； （3）用所得数值求积公式计算积分1 203x dx ?的精确值； （4）指出所得公式与一般的Newton-Cotes 型公式在形式上的重要区别。
（1）设?? ????=4321A ,计算1A 、()Cond A ∞和()A ρ。 （2）用列主元Gauss 消去法解方程组： 12312315410030.112x x x ????????????=????????????-?????? 六、（13分）对2阶线性方程组 11112212112222 a x a x b a x a x b +=??+=? （11220a a ≠ ） （1）证明求解此方程组的Jacobi 迭代与Gauss-Seidel 迭代同时收敛或同时发散； （2）当同时收敛时，试比较它们的收敛速度。
《安全工程数值分析》课程教学大纲 课程编号: 适用专业: 建筑安全工程专业 计划学时: 40学时计划学分: 2.0学分 一．本课程的性质和任务 安全工程数值分析是高等工科院校安全工程专业的一门重要专业选修课，并在许多领域中有着广泛的应用。本课程的任务是使学生熟悉用于数值分析的数学和力学基础知识，初步掌握利用计算机技术分析和解决工程问题的基本数值原理和方法，为学习以后专业课程创造条件。 二、课程内容及基本要求 第一章绪论 了解数课程的任务及学习方法 第二章计算机数学语言概述——MatLab 2.1 数学问题计算机求解概述 2.1.1 学习计算技术学语言的目的 2.1.2 数学问题的解析解与数值解 2.1.3 软件包的作用 2.1.4 MatLab语言的优势 2.2 MatLab语言程序设计基础 2.2.1 MatLab语言程序设计基础 2.2.2 基本数学运算 2.2.3 MatLab语言流程控制 2.2.4 MatLab函数的编写 2.2.5 二维图形绘制 2.2.6 三维图形绘制 第三章数值分析引论 3.1 数值算法的研究对象 3.1.1 了解计算方法基本理念 3.1.2 了解数值算法的特点
3.1.3 了解三类计算机算法的定义 3.2 误差分析的概念 3.2.1 了解误差和有效数字的关系 3.2.2 了解截断误差与收敛性的关系 3.2.3 了解舍入误差与数值稳定性的关系 3.2.4 了解数据误差与病态问题的关系 3.3 数值算法设计的要点 了解数值算法设计的要点 第四章数值代数 4.1 Gauss消去法 4.2 直接三角分解法 4.3 范数和误差分析 第五章插值法 5.1 Lagrange插值法 5.1.1 基本理论 5.1.2 Lagrange插值法在结构力学中的应用 5.2 Hermite插值法 5.2.1 基本理论 5.2.2 Hermite插值法在结构力学中的应用 第六章拟合 6.1 基本概念 6.2 最佳平方逼近 6.3 最小二乘法 第七章位移法 7.1 基本理论 7.2 实例分析 第八章有限单元法基本知识 8.1 变分原理 8.2 虚位移原理 8.3 势能原理 8.4 弹性力学基本方程 第九章结构有限单元法 9.1 平面拉压杆单元的有限单元分析 9.2 平面梁单元的有限单元分析 9.3 常应变三角形单元 9.4 矩形双线性单元 9.5 有限元分析应注意的问题和结果整理 三、使用大纲说明
学生实验报告实验课程名称 开课实验室 学院年级专业班 学生姓名学号 开课时间至学年学期
if(A(m,k)~=0) if(m~=k) A([k m],:)=A([m k],:); %换行 end A(k+1:n, k:c)=A(k+1:n, k:c)-(A(k+1:n,k)/ A(k,k))*A(k, k:c); %消去end end x=zeros(length(b),1); %回代求解 x(n)=A(n,c)/A(n,n); for k=n-1:-1:1 x(k)=(A(k,c)-A(k,k+1:n)*x(k+1:n))/A(k,k); end y=x; format short;%设置为默认格式显示,显示5位 （2）建立MATLAB界面 利用MA TLAB的GUI建立如下界面求解线性方程组： 详见程序。 五、计算实例、数据、结果、分析 下面我们对以上的结果进行测试，求解：
? ? ? ? ? ? ? ? ? ? ? ? - = ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? - - - - - - 7 2 5 10 13 9 14 4 4 3 2 1 13 12 4 3 3 10 2 4 3 2 1 x x x x 输入数据后点击和，得到如下结果： 更改以上数据进行测试，求解如下方程组： 1 2 3 4 43211 34321 23431 12341 x x x x ?? ???? ?? ???? ?? ???? = ?? ???? - ?? ???? - ???? ?? 得到如下结果：
研究生《数值分析》教学大纲 课程名称：数值分析 课程编号：S061005 课程学时：64 学时 课程学分： 4 适用专业：工科硕士生 课程性质：学位课 先修课程：高等数学，线性代数，计算方法，Matlab语言及程序设计 一、课程目的与要求 “数值分析”课是理工科各专业硕士研究生的学位课程。主要介绍用计算机解决数学问题的数值计算方法及其理论。内容新颖，起点较高，并加强了数值试验和程序设计环节。通过本课程的学习，使学生熟练掌握各种常用的数值算法的构造原理和过程分析，提高算法设计和理论分析能力，并且能够根据数学模型，提出相应的数值计算方法编制程序在计算机上算出结果。力求使学生掌握应用数值计算方法解决实际问题的常用技巧。 二、教学内容、重点和难点及学时安排： 第一章? 数值计算与误差分析( 4学时) 介绍数值分析的研究对象与特点，算法分析与误差分析的主要内容。 第一节数值问题与数值方法 第二节数值计算的误差分析 第三节数学软件工具----MATLAB 语言简介 重点：误差分析 第二章? 矩阵分析基础( 10学时) 建立线性空间、赋范线性空间、内积空间的概念，为学习以后各章打好基础。矩阵分解是解决数值代数问题的常用方法，掌握矩阵的三角分解、正交分解、奇异值分解，并能够编写算法程序。 第一节? 矩阵代数基础
第二节? 线性空间 第三节? 赋范线性空间 第四节? 内积空间和内积空间中的正交系 第五节矩阵的三角分解 第六节矩阵的正交分解 第七节矩阵的奇异值分解 难点：内积空间中的正交系。矩阵的正交分解。 重点：范数，施密特(Schmidt) 正交化过程，正交多项式，矩阵的三角分解, 矩阵的正交分解。 第三章? 线性代数方程组的数值方法( 12学时) 了解研究求解线性代数方程组的数值方法分类及直接法的应用范围。高斯消元法是解线性代数方程组的最常用的直接法，也是其它类型直接法的基础。在此方法基础上加以改进，可得选主元的高斯消元法、按比例增减的高斯消元法，其数值稳定性更高。掌握用列主元高斯消元法解线性方程组及计算矩阵的行列式及逆，并且能编写算法程序。掌握矩阵的直接三角分解法：列主元LU 分解，Cholesky分解。了解三对角方程组的追赶法的分解形式及数值稳定性的充分条件。掌握矩阵条件数的定义，并能利用条件数判别方程组是否病态以及对方程组的直接方法的误差进行估计。 迭代解法是求解大型稀疏方程组的常用解法。熟练掌握雅可比迭代法、高斯- 塞德尔迭代法及SOR 方法的计算分量形式、矩阵形式，并能在计算机上编出三种方法的程序用于解决实际问题。了解极小化方法：最速下降法、共轭斜量法。迭代法的收敛性分析是研究解线性代数方程组的迭代法时必须考虑的问题。对于上述常用的迭代法，须掌握其收敛的条件。而对一般的迭代法，掌握其收敛性分析的基本方法和主要结果有助于进一步探究新的迭代法。 第一节求解线性代数方程组的基本定理 第二节高斯消元法及其计算机实现 第三节矩阵分解法求解线性代数方程组 第三节? 误差分析和解的精度改进 第四节? 大型稀疏方程组的迭代法 第五节? 极小化方法 难点：列主元高斯消元法，直接矩阵三角分解。迭代法的收敛性，雅可比迭代法，高斯-塞德尔迭代法，SOR 迭代法。
机电工程学院 机械工程 陈星星 6720150109 《数值分析》课程设计实验报告 实验一 函数插值方法 一、问题提出 对于给定的一元函数)(x f y =的n+1个节点值(),0,1,,j j y f x j n ==。试用Lagrange 公式求其插值多项式或分段二次Lagrange 插值多项式。 数据如下： （1 求五次Lagrange 多项式5L ()x ，计算(0.596)f ,(0.99)f 的值。（提示：结果为(0.596)0.625732f ≈, (0.99) 1.05423f ≈） 实验步骤： 第一步：先在matlab 中定义lagran 的M 文件为拉格朗日函数 代码为： function[c,l]=lagran(x,y) w=length(x); n=w-1; l=zeros(w,w); for k=1:n+1 v=1; for j=1:n+1 if(k~=j) v=conv(v,poly(x(j)))/(x(k)-x(j)); end end l(k,:)=v; end c=y*l; end
第二步：然后在matlab命令窗口输入： >>>> x=[0.4 0.55 0.65 0.80,0.95 1.05];y=[0.41075 0.57815 0.69675 0.90 1.00 1.25382]; >>p = lagran(x,y) 回车得到： P = 121.6264 -422.7503 572.5667 -377.2549 121.9718 -15.0845 由此得出所求拉格朗日多项式为 p（x）=121.6264x5-422.7503x4+572.5667x3-377.2549x2+121.9718x-15.0845 第三步：在编辑窗口输入如下命令： >> x=[0.4 0.55 0.65 0.80,0.95 1.05]; >> y=121.6264*x.^5-422.7503*x.^4+572.5667*x.^3-377.2549*x.^2+121.9718 *x-15.0845; >> plot(x,y) 命令执行后得到如下图所示图形，然后 >> x=0.596; >> y=121.6264*x.^5-422.7503*x.^4+572.5667*x.^3-377.2549*x.^2+121.9718 *x-15.084 y =0.6257 得到f（0.596）=0.6257 同理得到f（0.99）=1.0542
数值分析习题集 （适合课程《数值方法A 》和《数值方法B 》） 长沙理工大学 第一章 绪 论 1. 设x >0,x 的相对误差为δ,求ln x 的误差. 2. 设x 的相对误差为2％,求n x 的相对误差. 3. 下列各数都是经过四舍五入得到的近似数,即误差限不超过最后一位的半个单位,试指 出它们是几位有效数字: *****123451.1021,0.031,385.6,56.430,7 1.0.x x x x x =====? 4. 利用公式(3.3)求下列各近似值的误差限: ********12412324(),(),()/,i x x x ii x x x iii x x ++其中**** 1234 ,,,x x x x 均为第3题所给的数. 5. 计算球体积要使相对误差限为1％,问度量半径R 时允许的相对误差限是多少? 6. 设028,Y =按递推公式 1n n Y Y -=…) 计算到100Y .27.982(五位有效数字),试问计算100Y 将有多大误差? 7. 求方程2 5610x x -+=的两个根,使它至少具有四位有效数字27.982). 8. 当N 充分大时,怎样求2 1 1N dx x +∞+?? 9. 正方形的边长大约为100㎝,应怎样测量才能使其面积误差不超过1㎝2 ? 10. 设 212S gt = 假定g 是准确的,而对t 的测量有±0.1秒的误差,证明当t 增加时S 的绝对 误差增加,而相对误差却减小. 11. 序列 {}n y 满足递推关系1101n n y y -=-(n=1,2,…),若0 1.41y =≈(三位有效数字), 计算到 10y 时误差有多大?这个计算过程稳定吗? 12. 计算6 1)f =, 1.4≈,利用下列等式计算,哪一个得到的结果最好? 3 -- 13. ()ln(f x x =,求f (30)的值.若开平方用六位函数表,问求对数时误差有多大?
偏微分方程数值解 一．教学目的 大量科学技术问题的数值计算都归结为偏微分方程的数值解法，应用数学专业计算方向的学生应该掌握偏微分方程数值解的基本知识和方法，重点介绍当今流行的偏微分方程数值解的两类主要方法，即有限差分法和有限元法。二．教学内容及学时分配 总学时为48学时 1、抛物型方程的有限差分法(9学时) 差分逼近的基本概念，抛物型方程的几种古典差分格式，差分格式的收敛性和稳定性概念， Lax等价性定理，研究稳定性的直接法和分离变量法，变系数方程与非线性方程的差分方法，多维问题交替方向法及分裂格式。 2、双曲型方程的差分方法(9学时) 一阶线性双曲型方程(组)的差分格式及稳定性分析，二阶线性双曲型方程的差分方法，拟线性双曲型方程(组)特征差分格式，守恒型方程的差分方法。 3、椭圆型方程差分方法(6学时) 二维poisson方程差分方程的建立，极坐标系下的差分格式，边界条件的处理，极值原理及先验估计，差分格式的收敛性。 4、变分原理与广义解(7学时) 引言，泛函的变分与泛函的极值，两点边值问题的变分原理，二阶椭圆边值问题的变分原理，Sobo1ev空间简介与微分方程广义解，古典Ritz—Galerkin 方法。 5、有限元离散方法(7学时) 两点边值问题的有限元法，二维边值问题的有限元法，有限元法解题的一般步骤。 6、形状函数与有限元空间(6学时) 一维高次元，二维矩形剖分的形状函数，三角形单元的形状函数，等参数单元，三维情形。 7、有限元解的收敛性与误差估计(4学时) Sobolev空间中的插值理论，有限元方法的收敛性与误差估计。 三.教学对象及先修课程
本课程为计算数学方向本科生 先修课程：数学分析，高等代数，数理方程，数值分析，泛函分析四．教材及主要参考书 偏微分方程数值解，陆金甫，关浩，清华大学出版社，1987 微分方程数值方法，胡建伟，胡建伟，科学出版社，1999
习题一 1.1 求下列各数的具有四位有效数字的近似值, 并指出其绝对误差限和相对误差限 )1.0ln(,121,101 1,1014321== = = x x x x 1.2 下列各数都是对准确值进行四舍五入得到的近似值, 指出它们的绝对误差限、相对误差限和有效数字的位 数。 3 * 5* 4* 3* 2* 1100.5,5000,50.31,3015.0,0315.0?=====x x x x x 1.3 为了使 3 1的近似值的相对误差不超过0.1%, 问应取几位有效数字? 1.4 怎样计算下列各题才能使得结果比较精确? (1) x x sin )sin(-+ε,其中ε充分小 (2) ? ++1 2 1N N x dx ,其中N 是充分大的正数 (3) x x sin cos 1-,其中x 充分小 (4) o 1cos 1- (5) 1001.0-e (6) )11010ln(84-- 1.5 求方程01562=+-x x 的两个根, 使至少具有四位有效数字。 习题二 2.1 证明方程043 =-+x x 在区间[1,2]内有且仅有一个根。如果用二分法求它具有五位有效数字的根，试问需对 分多少次？(不必求根) 2.2 用二分法求方程0134 =+-x x 在[0.3, 0.4]内的一个根, 精度要求2 10 2 1-?= ε。 2.3 找出下列方程的有根区间，选择适当的初始点用二分法求方程的根，精度要求2 10 -=ε。 (1) 02 =--x x ； (2) 06cos 2 =-++-x e x x ； (3) 01tan =--x x ； (4) 0sin 2=--x e x 。 2.4 考虑方程032 =-x e x ，将其改写为3 x e x ± =，取00=x ，用两种迭代公式迭代，分别收敛到1.0和-0.5附 近的两个根(取精度要求3 10-=ε)。
《计算方法》课程教学大纲 课程编号： 学时：54 学分：3 适用对象：教育技术学专业 先修课程：高等数学、线性代数 考核方式：本课程考试以笔试为主70%，兼顾学生的平时成绩30%。 使用教材及主要参考书： 使用教材： 李庆扬.《数值分析(第四版)》, 清华大学出版,2014年。 主要参考书： 1.朱建新,李有法.《高等学校教材:数值计算方法(第3版)》,高等教育出版社,2012。 2.徐萃薇,孙绳武.《计算方法引论(第4版)》,高等教育出版社,2015。 一课程的性质和任务 计算方法是教育技术学专业学生的一门专业选修课。作为计算数学的一个重要分支，它是数学科学与计算机技术结合的一门应用性很强的学科，本课程重点介绍计算机上常用的基本计算方法的原理和使用；同时对计算方法作适当的分析。 教学任务：通过本课程的学习，要使学生具有现代数学的观点和方法，并初步掌握处理计算机常用数值分析的构造思想和计算方法。同时，也要培养学生抽象思维和慎密概括的能力，使学生具有良好的开拓专业理论的素质和使用所学知识分析和解决实际问题的能力。 二教学目的与要求 教学目的：通过学习使学生了解数值计算方法的基本原理。了解计算机与数学结合的作用及课程的应用性。为今后使用计算机解决实际问题中的数值计算问题打下基础。 通过理论教学达到如下基本要求。 1．了解误差的概念 2．掌握常用的解非线性方程根的方法 3．熟练掌握线性代数方法组的解法 4．熟练掌握插值与拟合的常用方法 5．掌握数值积分方法 6．了解常微分方程初值问题的数值方法 三学时分配
四教学中应注意的问题 本课程是一门理论性较强、内容较抽象的综合课程，因此面授辅导或自学，将是不可缺少的辅助教学手段，教师在教学的过程中一定要注意理论结合实际，课堂教学并辅助上机实验，必须通过做练习题和上机实践来加深对概念的理解和掌握，熟悉公式的运用，从而达到消化、掌握所学知识的目的。同时应注重面授辅导或答疑，及时解答学生的疑难问题。 五教学内容 第一章绪论（误差） 基本内容: 第一节数值分析研究的对象和特点 第二节数值计算的误差 1.误差的来源与分类 2.误差与有效数字 3.数值运算的误差估计 第三节误差的定性分析与避免误差的危害 1.病态问题与条件数 2.算法的数值稳定性 3.避免误差危害的若干原则 教学重点难点: 重点：数值运算的误差估计。 难点：误差的定性分析与避免误差的危害。
《数值分析》教学大纲 课程名称：数值分析 课程编号：0811010001 课程学时：54学时 课程学分：3 适用专业：控制理论与控制工程、计算机软件与理论 课程性质：专业基础课 先修课程：《高等数学》、《线性代数》 大纲执笔人： 编写时间：2009年8月 一、课程性质、地位和作用 《数值分析》是计算机软件与理论、控制理论与控制工程专业的一门重要专业方向课，属必修课。其任务在于研究用计算机求解各种数学问题的数值方法及其理论，是程序设计和对数值结果进行分析的依据。本课程理论严谨，实用性强。为学生毕业后从事科学计算等相关行业的工作提供一定的基础。 二、课程教学对象、目的和要求 本课程适用于计算机软件与理论、控制理论与控制工程等相关硕士研究生专业。 课程教学目的和要求： 1、从内容上，以现代化的计算机和数学软件为工具，以数学模型为基础进行模拟研究。要求学生牢固掌握数学分析、高等代数等基础数学中常用的、行之有效的数值计算方法。 2、从能力方面，要求学生掌握从实际问题出发，建立数学模型，将数学模型问题转换成数值问题，进而研究求解数值问题的数值方法，并设计出相应的数值算法。 3、从教学方法上，注重理论联系实际，做到重概念，重方法，重应用，重能力的培养。 三、课程内容及学时分配 总学时：54学时
（一）误差：3学时 1.1 误差的来源与分类 1.2 误差与有效数字 1.3 函数的误差估计 1.4 近似数的四则运算及数值计算中需注意的几个问题 要求学生了解数值计算方法的对象和特点。理解绝对误差、相对误差和有效数字的概念及其对数值计算的影响。掌握绝对误差、相对误差和有效数字的计算方法。 （二）非线性方程求根（3学时） 2.1 二分法（分半法） 2.2 迭代法 2.3 牛顿法 2.4 牛顿法的改进 2.5 迭代法的收敛阶 2.6 劈因子法 要求学生了解二分法的基本思想，会用二分法求根，并估计误差。理解迭代法的基本思想、能熟练地建立迭代公式，并判断其收敛性。熟练掌握Newton 迭代法的原理及计算。 （三）线性代数方程组的直接法（6学时） 3.1 高斯消元法 3.2 三角分解法 要求学生了解Gauss消去法原理。掌握道立特（Doolittle）分解法和科路特（Cholesky）分解法求解方程组。 （四）解线性方程组的迭代法（6学时） 4.1 向量和矩阵的范数 4.2 线性方程组的误差分析 4.3 雅可比（Jacobi）方法和高斯赛德尔（Gauss-Seidel）方法
数值分析试卷1 一、填空题（每空2分，共30分） 1. 近似数231.0=*x 关于真值229.0=x 有____________位有效数字； 2. 设)(x f 可微，求方程)(x f x =根的牛顿迭代格式是_______________________________________________； 3. 对1)(3++=x x x f ，差商=]3,2,1,0[f _________________； =]4,3,2,1,0[f ________； 4. 已知??? ? ??-='-=1223,)3,2(A x ，则=∞||||Ax ________________，=)(1A Cond ______________________ ； 5. 求解线性方程组?????=+=+045 11532121x x x x 的高斯—赛德尔迭代格式为_______________________________________；该迭代格式迭代矩阵的谱半径=)(G ρ_______________； 二、（12分）(1)设LU A =，其中L 为下三角阵，U 为单位上三角阵。已知 ?????? ? ??------=2100121001210012A ，求L ，U 。 (2)设A 为66?矩阵，将A 进行三角分解：LU A =，L 为单位下三角阵，U 为上三角阵，试写出L 中的元素65l 和U 中的元素56u 的计算公式。 三、给定数据表如下 x 0.20.40.60.81 1.2f(x)212523202124 （1） 用三次插值多项式计算f ( 0.7 ) 的近似值； （2） 用二次插值多项式计算f ( 0.95 ) 的近似值： （3） 用分段二次插值计算 f ( x ) )2.12.0(≤≤x 的近似值能保证有几位有
数值计算方法教学大纲(本) 本着“崇术重用、服务地方”的办学理念和我校“高素质应用型人才”的培养目标，特制定了适合我校工科专业本科生的新教学大纲。 一、课程计划 课程名称：数值计算方法Numerical Calculation Method 课程定位：数学基础课 开课单位：理学院 课程类型：专业选修课 开设学期：第七学期 讲授学时：共15周，每周4学时，共60学时 学时安排：课堂教学40学时+实验教学20学时 适用专业：计算机、电科、机械等工科专业本科生 教学方式：讲授（多媒体为主）+上机 考核方式：考试60%+上机实验30%+平时成绩10% 学分：3学分 与其它课程的联系 预修课程：线性代数、微积分、常微分方程、计算机高级语言等。 后继课程：偏微分方程数值解及其它专业课程。 二、课程介绍 数值计算方法也称为数值分析，是研究用计算机求解各种数学问题的数值方法及其理论的一门学科。随着计算科学与技术的进步和发展，科学计算已经与理论研究、科学实验并列成为进行科学活动的三大基本手段，作为一门综合性的新科学，科学计算已经成为了人们进行科学活动必不可少的科学方法和工具。 数值计算方法是科学计算的核心内容，它既有纯数学高度抽象性与严密科学性的特点，又有应用的广泛性与实际实验的高度技术性的特点，是一门与计算机使用密切结合的实用性很强的数学课程.主要介绍插值法、函数逼近与曲线拟合、线性方程组迭代解法、数值积分与数值微分、非线性方程组解法、常微分方程数值解以及矩阵特征值与特征向量数值计算，并特别加强实验环节的训练以提高学生动手能力。通过本课程的学习，不仅能使学生初步掌握数值计算方法的基本理论知识，了解算法设计及数学建模思想，而且能使学生具备一定的科学计算能力和分析与解决问题的能力，不仅为学习后继课程打下良好的理论基础，也为将来从事科学计算、计算机应用和科学研究等工作奠定必要的数学基础。 科学计算是21世纪高层次人才知识结构中不可缺少的一部分，它潜移默化地影响着人们的思维方式和思想方法，并提升一个人的综合素质。
1. 正方形的边长大约为100cm ，应怎样测量才能使面积误差不超过1cm 2? 2. 已测得某场地长l 的值为110=*l m ，宽d 的值为80=*d m ，已知 2.0≤-*l l m, 1.0≤-*d d m, 试求面积ld s =的绝对误差限与相对误差限.
3.为使π的相对误差小于0.001%,至少应取几位有效数字？ 4.设x的相对误差界为δ，求n x的相对误差界. 5.设有3个近似数a=2.31，b=1.93，c=2.24，它们都有3位有效数字，试计算 p=a+bc的误差界和相对误差界，并问p的计算结果能有几位有效数字？
6. 已知33348 7.034.0sin ,314567.032.0sin ==，请用线性插值计算3367.0sin 的值，并估计截断误差. 7. 已知sin0.32=0.314567, sin0.34=0.333487, sin0.36= 0.352274,用抛物插值计算sin0.3367的值, 并估计误差. 8. 已知 1 6243sin ,sin π ππ== =请用抛物插值求sin50的值,并估计误差
9. . .6,8,7,4,1)(,5,4,3,2,1求四次牛顿插值多项式时设当==i i x f x 10. 已知4)2(,3)1(,0)1(=-=-=f f f ， 求函数)(x f 过这3点的2次牛顿插 值多项式 . 11. 设x x f =)(，并已知483240.1)2.2(,449138.1)1.2(,414214.1)0.2(===f f f ，
试用二次牛顿插值多项式计算(2.15)f 的近似值，并讨论其误差 12. 设],[)(b a x f 在上有四阶连续导数，试求满足条件)2,1,0()()(==i x f x P i i 及 )()(11x f x P '='的插值多项式及其余项表达式. 13. 给定3201219(),,1,,44f x x x x x ====试求()f x 在1944?? ???? ，上的三次埃尔米特
实验报告 实验项目名称数值积分与数值微分实验室数学实验室 所属课程名称数值逼近 实验类型算法设计 实验日期 班级 学号 姓名 成绩
实验概述： 【实验目的及要求】 本次实验的目的是熟练《数值分析》第四章“数值积分与数值微分”的相关内容，掌握复合梯形求积公式、复合辛普森求积公式、龙贝格求积公式以及高斯-勒让德公式。 本次试验要求编写复合梯形求积公式、复合辛普森求积公式、龙贝格求积公式以及高斯-勒让德公式的程序编码，并在MATLAB软件中去实现。 【实验原理】 《数值分析》第四章“数值积分与数值微分”的相关内容，包括：复合梯形求积公式、复合辛普森求积公式、龙贝格求积公式以及高斯-勒让德公式的相应算法和相关性质。 【实验环境】（使用的软硬件） 软件： MATLAB 2012a 硬件： 电脑型号：联想 Lenovo 昭阳E46A笔记本电脑 操作系统：Windows 8 专业版 处理器：Intel（R）Core（TM）i3 CPU M 350 @2.27GHz 2.27GHz 实验内容： 【实验方案设计】 第一步，将书上关于复合梯形求积公式、复合辛普森求积公式、龙贝格求积公式以及高斯-勒让德公式的内容转化成程序语言，用MATLAB实现；第二步，分别用以上求积公式的程序编码求解不同的问题。 【实验过程】（实验步骤、记录、数据、分析） 实验的主要步骤是：首先分析问题，根据分析设计MATLAB程序，利用程序算出问题答案，分析所得答案结果，再得出最后结论。 实验：用不同数值方法计算积分 (1) 取不同的步长h.分别用复合梯形及复合辛普森求积计算积分，给出误差中关于h的函数，并与积分精确值比较两个公式的精度，是否存在一个最小的h，使得精度不能再被改善？ (2) 用龙贝格求积计算完成问题（1）。 (3)用勒让德多项式确定零点，再代入计算高斯公式，使其精度达到10-4 （1）在MATLAB的Editor中建立一个M-文件，输入程序代码，实现复合梯形求积公式的程序代码如下：