数值计算方法教学大纲
《数值计算方法》教学大纲
课程编号：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.重点、难点
重点：迭代法，牛顿法。
难点：牛顿法。
3.说明：该内容在求解方程解的算法中具有普遍的应用意义，学生必须掌握。
(六) 线性方程组解法（6学时）
具体内容：高斯消去及高斯主元消去法，高斯消去法的变形，迭代法，龙贝格求积算法简介，数值微分，理查逊外推加速方法。
1.基本要求
（1）熟练掌握高斯主元消去法和了解高斯消去法的变形。
（2）掌握几种常用的简单迭代法。
2.重点、难点
重点：高斯主元消去法及其变形。
难点：高斯消去法
3.说明：这部分内容讲述几种基本的线性方程的求解方法，为复杂度较大的方程组系统的求解打下基础。
四、教学安排及方式
五、考核方式
笔试(半开卷)。
各教学环节占总分的比例:平时测验及作业：40%，期末考试：60%
六、推荐教材与参考资料
王能超，数值分析简明教程(修订版)，华中科技大学出版社,2002
（执笔人：张义门审核人：张鹤鸣）
2005年8月20日
数值分析总复习提纲教材
数值分析总复习提纲 数值分析课程学习的内容看上去比较庞杂，不同的教程也给出了不同的概括，但总的来说无非是误差分析与算法分析、基本计算与基本算法、数值计算与数值分析三个基本内容。在实际的分析计算中，所采用的方法也无非是递推与迭代、泰勒展开、待定系数法、基函数法等几个基本方法。 一、误差分析与算法分析 误差分析与算法设计包括这样几个方面： （一）误差计算 1、截断误差的计算 截断误差根据泰勒余项进行计算。 基本的问题是 (1)1 ()(01)(1)! n n f x x n θεθ++<<<+，已知ε求n 。 例1．1：计算e 的近似值，使其误差不超过10－6。 解：令f(x)=e x ,而f （k)(x)=e x ,f （k)(0)=e 0=1。由麦克劳林公式，可知 211(01)2!!(1)! n x x n x x e e x x n n θθ+=+++++<<+ 当x=1时，1 111(01)2! !(1)! e e n n θθ=+++ ++ <<+ 故3 (1)(1)!(1)! n e R n n θ=<++。 当n ＝9时，R n (1)<10－6，符合要求。此时， e≈2.718 285。 2、绝对误差、相对误差及误差限计算 绝对误差、相对误差和误差限的计算直接利用公式即可。 基本的计算公式是： ①e(x)＝x *－x ＝△x ＝dx ② *()()()ln r e x e x dx e x d x x x x ==== ③(())()()()e f x f x dx f x e x ''== ④(())(ln ())r e f x d f x = ⑤121212121122121122((,))(,)(,)(,)()(,)()x x x x e f x x f x x dx f x x dx f x x e x f x x e x ''''=+=+ ⑥121212((,)) ((,))(,) f x x f x x f x x εδ=
北师大网络教育 数值分析 期末试卷含答案
注：1、教师命题时题目之间不留空白； 2、考生不得在试题纸上答题，教师只批阅答题册正面部分，若考北师大网络教育——数值分析——期末考试卷与答案 一．填空题（本大题共4小题，每小题4分，共16分） 1.设有节点012,,x x x ，其对应的函数()y f x =的值分别为012,,y y y ，则二次拉格朗日插值基函数0()l x 为 。 2.设()2f x x =，则()f x 关于节点0120,1,3x x x ===的二阶向前差分为 。 3.设110111011A -????=--????-??，233x ?? ??=?? ???? ，则1A ＝ ，1x = 。 4. 1n +个节点的高斯求积公式的代数精确度为 。 二．简答题（本大题共3小题，每小题8分，共24分） 1. 哪种线性方程组可用平方根法求解？为什么说平方根法计算稳定？ 2. 什么是不动点迭代法？()x ?满足什么条件才能保证不动点存在和不动点迭代序列收敛于()x ?的不动点？ 3. 设n 阶矩阵A 具有n 个特征值且满足123n λλλλ>≥≥≥ ，请简单说明求解矩阵A 的主特征值和特征向量的算法及流程。 三．求一个次数不高于3的多项式()3P x ，满足下列插值条件： i x 1 2 3 i y 2 4 12 i y ' 3 并估计误差。（10分） 四．试用1,2,4n =的牛顿-科特斯求积公式计算定积分1 01 1I dx x =+? 。（10分） 五．用Newton 法求()cos 0f x x x =-=的近似解。（10分） 六．试用Doolittle 分解法求解方程组：
注：1、教师命题时题目之间不留空白； 2、考生不得在试题纸上答题，教师只批阅答题册正面部分，若考 12325610413191963630 x x x -?????? ??????-=?????? ??????----?????? （10分） 七．请写出雅可比迭代法求解线性方程组1231231 23202324 812231530 x x x x x x x x x ++=?? ++=??-+=? 的迭代格式，并 判断其是否收敛？（10分） 八．就初值问题0(0)y y y y λ'=??=?考察欧拉显式格式的收敛性。（10分）
数值计算实验课题目
数值实验课试题 本次数值实验课结课作业，请按题目要求内容写一篇文章。按题目要求 人数自由组合，每组所选题目不得相同(有特别注明的题目除外)。试题如下： 1）解线性方程组的Gauss 消去法和列主元Gauss 消去法(2人)/*张思珍，巩艳华*/ 用C 语言将不选主元和列主元Gauss 消去法编写成通用的子程序，然后用你编写的程序求解下列84阶的方程组 ???? ?????? ? ??=??????????? ????????????? ? ?1415151515768 168 168 168 1681684 8382321 x x x x x x 参考书目： 1.《计算机数值方法》，施吉林、刘淑珍、陈桂芝编 2.《数值线性代数》，徐树方、高立、张平文编 3.《数值分析简明教程》，王能超编 2）解线性方程组的平方根法(4人)/*朱春成、黄锐奇、张重威、章杰*/ 用C 语言将平方根法和改进的平方根法编写成通用的子程序，然后用你编写的程序求解对称正定方程组b Ax =，其中 (1)b 随机的选取，系数矩阵为100阶矩阵 ?????? ???? ? ? ?101 1101 1101 1101 1101110 ； (2)系数矩阵为40阶的Hilbert 矩阵，即系数矩阵A 的第i 行第j 列元素为 1 1-+= j i a ij ，向量b 的第i 个分量为∑=-+ = n j i j i b 1 1 1. 参考书目： 1.《计算机数值方法》，施吉林、刘淑珍、陈桂芝编 2.《数值线性代数》，徐树方、高立、张平文编
3.《数值分析简明教程》，王能超编 3）三对角线方程组的追赶法(3人)/*黄佳礼、唐伟、韦锡倍*/ 用C 语言将三对角线方程组的追赶法法编写成通用的子程序，然后用你编写的程序求解如下84阶三对角线方程组 ???? ?????? ? ??=??????????? ????????????? ? ?1415151515768 168 168 168 16816 84 8382321 x x x x x x 参考书目： 1.《计算机数值方法》，施吉林、刘淑珍、陈桂芝编 2.《数值分析简明教程》，王能超编 4）线性方程组的Jacobi 迭代法(3人)/*周桂宇、杨飞、李文军*/ 用C 语言将Jacobi 迭代法编写成独立的子程序，并用此求解下列方程组， 精确到小数点后5位 ???? ? ??=????? ??????? ? ?-149012 2111221 3 2 1 x x x 参考书目： 1.《计算机数值方法》，施吉林、刘淑珍、陈桂芝编 2.《数值线性代数》，徐树方、高立、张平文编 3.《数值分析简明教程》，王能超编 5）线性方程组的Gauss-Seidel 迭代法(3人)/*张玉超、范守平、周红春*/ 用C 语言将Gauss-Seidel 迭代法编写成独立的子程序，并用此求解下列方程组，精确到小数点后5位 ???? ? ??=????? ??????? ? ?--39721 1111112 3 2 1 x x x 参考书目： 1.《计算机数值方法》，施吉林、刘淑珍、陈桂芝编 2.《数值线性代数》，徐树方、高立、张平文编 3.《数值分析简明教程》，王能超编 6）解线性方程组的最速下降法法(2人)/*赵育辉、阿热孜古丽*/ 用C 语言将最速下降法编写成通用的子程序，然后用你编写的程序求解对称
数值计算方法课程设计
重庆邮电大学 数学与应用数学 专业 《数值计算方法》课程设计 姓名： 李金徽 王莹 刘姝楠 班级： 1131001 1131002 1131002 学号： 2010213542 2010213570 2010213571 设计时间： 2012-6-4 指导教师： 朱伟
一、课程设计目的 在科学计算与工程设计中，我们常会遇到求解线性方程组的问题，对于系数矩阵为低阶稠密矩阵的线性方程组，可以用直接法进行消元，而对于系数矩阵为大型稀疏矩阵的情况，直接法就显得比较繁琐，而迭代法比较适用。比较常用的迭代法有Jacobi 迭代与Gauss － seidel 迭代。本文基于两种方法设计算法，并比较他们的优劣。 二、课程设计内容 给出Jacobi 迭代法和Gauss-Seidel 迭代法求解线性方程组的算法思想和MATLAB 程序实现，并对比分析这两种算法的优劣。 三、问题的分析（含涉及的理论知识、算法等） Jacobi 迭代法 方程组迭代法的基本思想和求根的迭代法思想类似，即对于线性 方程组Ax = b( 其中n n n R b R R A ∈?∈,），即方程组 )1(2211222221211 1212111?? ???? ?=+?++??=+?++=+?++n n nn n n n n n n b x a x a x a b x a x a x a b x a x a x a 将系数矩阵A 写为 )2(000000 21122 12122 11U L D a a a a a a a a a A n n n n nn --≡??? ?? ? ? ??---- ??????? ??----??????? ??= 若选取D M =，则U L A M N +=-=，方程组)1(转化为等价方程组 b x U L Dx ++=)(
数值分析复习提纲
数值分析（英）复习提纲 考试以基本概念为主，书上以前布置的计算机题目都不作要求。 第一章Solving equations 1.1 THE BISECTION METHOD (a) 熟练掌握二分法; (b) 对于给定解的误差精度要求能够熟练计算所需二分法步数，参考书上28页内容。 习题5,6 1.2 FIXED-POINT ITERATION (a) 熟练掌握不动点迭代方法求方程的根；掌握不动点迭代方法的线性收敛性与收敛率; 此节书后习题不作要求。 1.4 NEWTON’S METHOD （a）熟练掌握方程求根的NEWTON’S METHOD：Example 1.11, 1.12, 1.13 （b）对于重根熟练掌握Theorem 1.12, Theorem 1.13 习题2,5,7 第二章Systems of Equations 2.2 THE LU FACTORIZATION （a）掌握矩阵LU分解方法； （b）会使用LU分解方法求线性方程组的解:Example 2.5, 2.6, 2.7 2.3 SOURCES OF ERROR 本节只要掌握矩阵范数的定义，参阅90页 2.4 THE PA = LU FACTORIZATION 熟练掌握2.4.2 Permutation matrices, 2.4.3 PA = LU factorization: Example 2.16, 2.17, 2.18 习题4 2.5 ITERATIVE METHODS 熟练掌握Jacobi Method，Gauss–Seidel Method. 习题2
第三章Interpolation 3.1 DATA AND INTERPOLATING FUNCTIONS： （a）熟练掌握Lagrange interpolation （b）熟练掌握Newton’s divided differences 习题1,2,5 3.2 INTERPOLATION ERROR 熟练掌握定理3.4, Example 3.8, 习题1,2,4 第四章Least Squares 4.1.1 Inconsistent systems of equations 熟练掌握Normal equations for least squares：Example 4.1, Example 4.2 习题1,2 第五章Numerical Differentiation and Integration 5.1 NUMERICAL DIFFERENTIATION 熟练掌握一阶导数的Two-point forward-difference formula，Three-point centered-difference formula 熟练掌握二阶导数的Three-point centered-difference formula for second derivative 习题1,2,5,8,9 5.2 NEWTON–COTES FORMULAS FOR NUMERICAL INTEGRATION 熟练掌握Composite Trapezoid Rule，Example 5.8，习题1 第六章Ordinary Differential Equations 6.1.1 Euler’s Method (a) 熟练掌握Euler方法（6.7）: Example6.2 习题5 6.2.2 The explicit Trapezoid Method 熟练掌握The explicit Trapezoid Method（6.29）：Example6.10 习题1
数值分析学期期末考试试题与答案(A)
期末考试试卷（A 卷） 2007学年第二学期 考试科目： 数值分析 考试时间：120 分钟 学号 姓名 年级专业 一、判断题（每小题2分，共10分） 1. 用计算机求 1000 1000 1 1 n n =∑时，应按照n 从小到大的顺序相加。 （ ） 2. 为了减少误差,进行计算。 （ ） 3. 用数值微分公式中求导数值时，步长越小计算就越精确。 （ ） 4. 采用龙格－库塔法求解常微分方程的初值问题时，公式阶数越高，数值解越精确。（ ） 5. 用迭代法解线性方程组时，迭代能否收敛与初始向量的选择、系数矩阵及其演变方式有 关，与常数项无关。 （ ） 二、填空题（每空2分，共36分） 1. 已知数a 的有效数为0.01，则它的绝对误差限为________，相对误差限为_________. 2. 设1010021,5,1301A x -????????=-=-????????-???? 则1A =_____，2x =______，Ax ∞ =_____. 3. 已知5 3 ()245,f x x x x =+-则[1,1,0]f -= ,[3,2,1,1,2,3]f ---= . 4. 为使求积公式 1 1231 ()()(0)33 f x dx A f A f A f -≈- ++? 的代数精度尽量高，应使1A = ，2A = ，3A = ，此时公式具有 次的代数精度。 5. n 阶方阵A 的谱半径()A ρ与它的任意一种范数A 的关系是 . 6. 用迭代法解线性方程组AX B =时，使迭代公式(1) ()(0,1,2,)k k X MX N k +=+=产 生的向量序列{ }() k X 收敛的充分必要条件是 . 7. 使用消元法解线性方程组AX B =时，系数矩阵A 可以分解为下三角矩阵L 和上三角矩
曲线拟合的数值计算方法实验
曲线拟合的数值计算方法实验 【摘要】实际工作中，变量间未必都有线性关系，如服药后血药浓度与时间的关系；疾病疗效与疗程长短的关系；毒物剂量与致死率的关系等常呈曲线关系。曲线拟合（curve fitting）是指选择适当的曲线类型来拟合观测数据，并用拟合的曲线方程分析两变量间的关系。曲线直线化是曲线拟合的重要手段之一。对于某些非线性的资料可以通过简单的变量变换使之直线化，这样就可以按最小二乘法原理求出变换后变量的直线方程，在实际工作中常利用此直线方程绘制资料的标准工作曲线，同时根据需要可将此直线方程还原为曲线方程，实现对资料的曲线拟合。常用的曲线拟合有最小二乘法拟合、幂函数拟合、对数函数拟合、线性插值、三次样条插值、端点约束。 关键词曲线拟合、最小二乘法拟合、幂函数拟合、对数函数拟合、线性插值、三次样条插值、端点约束 一、实验目的 1.掌握曲线拟合方式及其常用函数指数函数、幂函数、对数函数的拟合。 2.掌握最小二乘法、线性插值、三次样条插值、端点约束等。 3.掌握实现曲线拟合的编程技巧。 二、实验原理 1.曲线拟合 曲线拟合是平面上离散点组所表示的坐标之间的函数关系的一种数据处理方法。用解析表达式逼近离散数据的一种方法。在科学实验或社会活动中，通过 实验或观测得到量x与y的一组数据对（X i ，Y i )（i=1，2，...m），其中各X i 是彼此不同的。人们希望用一类与数据的背景材料规律相适应的解析表达式，y=f(x，c）来反映量x与y之间的依赖关系，即在一定意义下“最佳”地逼近或 拟合已知数据。f(x，c)常称作拟合模型，式中c=(c 1，c 2 ，…c n )是一些待定参 数。当c在f中线性出现时，称为线性模型，否则称为非线性模型。有许多衡量拟合优度的标准，最常用的一种做法是选择参数c使得拟合模型与实际观测值在
数值计算方法课程报告
课程报告 课程名称______《数值计算》 __ 学生学院_____机电工程学院___ 专业班级_____微电子（1）班____ 学号________ 学生姓名_______________ 指导教师_____ ________ XXXX年XX月XX日
姓 名： 线 学 号 ： 订 装专 业：学院： 广东工业大学考试试卷( A ) 课程名称: 数值计算试卷满分100 分考试时间: 2015 年 12 月 26 日(第 17 周星期六) 题号一二三四五六七八九十总分 评卷得分 评卷签名 复核得分 复核签名 “数值计算”考试要求 “数值计算”考试以开卷形式进行。在“数值计算”课程考试日（2015 年12 月 19 日，第 12 周星期五）考试时间，在考试教室领取试题，在 2015 年12 月 26 日（第 17 周星期六）进行答辩。不参加答辩者将取消考试成绩。 “数值计算”考试结果要求独立在计算机上完成，可使用Matlab或 C 程序编程实现。考试结果将以报告书形式提交，内容包括对问题描述、计算程序以及算例、计算结果、分析组成。计算程序要求具有通用性，能够处理异常情况，可以输入问题、算法参数、算例及初始值，在计算过程中显示当前计算状态、计算完成后显示计算结果。上述内容将作为试卷成绩的主要评定依据。特别提醒，不得使用教师在讲课和实验时的范例作为考试结果。报告书具体格式参考毕业设计手册。 以考生学号命名的文件夹存放程序及报告书电子版，以班级为单位刻录在一张光盘中，与打印版报告书一起由班长和学习委员一起上交任课教师。 数值计算课程总成绩将由试卷成绩（70%）、平时成绩（30%）组成。
数值计算方法复习提纲
数值计算方法复习提纲 第一章 数值计算中的误差分析 1．了解误差及其主要来源，误差估计； 2．了解误差(绝对误差、相对误差)和有效数字的概念及其关系； 3．掌握算法及其稳定性，设计算法遵循的原则。 1、 误差的来源 模型误差 观测误差 截断误差 舍入误差 2误差与有效数字 绝对误差 E （x ）=x-x * 绝对误差限ε εε+≤≤-**x x x 相对误差 ***/)(/)()(x x x x x x x E r -≈-= 有效数字 m n a a a x 10.....021*?±= 若 n m x x -?≤ -102 1 *，称*x 有n 位有效数字。 有效数字与误差关系 （1） m 一定时，有效数字n 越多，绝对误差限越小； （2） *x 有n 位有效数字，则相对误差限为)1(1 1021 )(--?≤ n r a x E 。 选择算法应遵循的原则 1、 选用数值稳定的算法，控制误差传播；
例 ?= 10 1dx e x e I x n n e I nI I n n 11101 - =-=- △!n x n =△x 0 2、 简化计算步骤，减少运算次数； 3、 避免两个相近数相减，和接近零的数作分母； 避免 第二章 线性方程组的数值解法 1．了解Gauss 消元法、主元消元法基本思想及算法； 2．掌握矩阵的三角分解，并利用三角分解求解方程组； （Doolittle 分解；Crout 分解；Cholesky 分解；追赶法） 3．掌握迭代法的基本思想，Jacobi 迭代法与Gauss-Seidel 迭代法； 4．掌握向量与矩阵的范数及其性质,迭代法的收敛性及其判定 。 本章主要解决线性方程组求解问题，假设n 行n 列线性方程组有唯一解，如何得到其解 ?? ??? ? ?=+++=+++=+++n n nn n n n n n n b x a x a x a b x a x a x a b x a x a x a (22112222212111212111) 两类方法，第一是直接解法，得到其精确解； 第二是迭代解法，得到其近似解。
数值分析期末考试复习题及其答案.doc
数值分析期末考试复习题及其答案 1. 已知325413.0,325413* 2* 1==X X 都有6位有效数字，求绝对误差限。（4分） 解： 由已知可知,n=6 5.01021 ,0,6,10325413.0016*1=?= =-=?=ε绝对误差限n k k X 2分 620* 21021,6,0,10325413.0-?=-=-=?=ε绝对误差限n k k X 2分 2. 已知?????=001A 220 - ???? ?440求21,,A A A ∞ （6分） 解： {},88,4,1max 1==A 1分 {},66,6,1max ==∞A 1分 () A A A T max 2λ= 1分 ?????=001A A T 420 ?? ?? ? -420?????001 220 - ?????440=?????001 080 ???? ?3200 2分 {}3232,8,1max )(max ==A A T λ 1分 24322==A 3. 设3 2 )()(a x x f -= （6分） ① 写出f(x)=0解的Newton 迭代格式 ② 当a 为何值时，)(1k k x x ?=+ （k=0,1……）产生的序列{}k x 收敛于2 解： ①Newton 迭代格式为： x a x x x a x a x x a x x x f x f x x k k k k k k k k k k 665)(665)(6)()(')(2 2 32 1 += +=---=-=+? 3分
②时迭代收敛即当222,112 10)2(',665)('2<<-<-=-=a a x a x ?? 3分 4. 给定线性方程组Ax=b ，其中：? ??=1 3A ??? 22，??????-=13b 用迭代公式)()()()1(k k k Ax b x x -+=+α（k=0,1……）求解Ax=b ，问取什么实数α，可使迭代收 敛 （8分） 解： 所给迭代公式的迭代矩阵为?? ? --? ??--=-=ααααα21231A I B 2分 其特征方程为 0) 21(2)31(=----= -αλα ααλλB I 2分 即，解得αλαλ41,121-=-= 2分 要使其满足题意，须使1)(数值分析实验报告1
实验一误差分析 实验1.1（病态问题） 实验目的：算法有“优”与“劣”之分，问题也有“好”与“坏”之别。对数值方法的研究而言，所谓坏问题就是问题本身对扰动敏感者，反之属于好问题。通过本实验可获得一个初步体会。 数值分析的大部分研究课题中，如线性代数方程组、矩阵特征值问题、非线性方程及方程组等都存在病态的问题。病态问题要通过研究和构造特殊的算法来解决，当然一般要付出一些代价（如耗用更多的机器时间、占用更多的存储空间等）。 问题提出：考虑一个高次的代数多项式 显然该多项式的全部根为1,2,…,20共计20个，且每个根都是单重的。现考虑该多项式的一个扰动 其中ε（1.1）和（1.221,,,a a 的输出b ”和“poly ε。 （1（2 （3）写成展 关于α solve 来提高解的精确度，这需要用到将多项式转换为符号多项式的函数poly2sym,函数的具体使用方法可参考Matlab 的帮助。 实验过程： 程序： a=poly(1:20); rr=roots(a); forn=2:21 n form=1:9 ess=10^(-6-m);
ve=zeros(1,21); ve(n)=ess; r=roots(a+ve); -6-m s=max(abs(r-rr)) end end 利用符号函数：（思考题一）a=poly(1:20); y=poly2sym(a); rr=solve(y) n
很容易的得出对一个多次的代数多项式的其中某一项进行很小的扰动，对其多项式的根会有一定的扰动的，所以对于这类病态问题可以借助于MATLAB来进行问题的分析。 学号：06450210 姓名：万轩 实验二插值法
太原理工大学数值计算方法实验报告
本科实验报告 课程名称：计算机数值方法 实验项目：方程求根、线性方程组的直接解 法、线性方程组的迭代解法、代数插值和最 小二乘拟合多项式 实验地点：行勉楼 专业班级： ******** 学号： ********* 学生姓名： ******** 指导教师：李誌，崔冬华 2016年 4 月 8 日
y = x*x*x + 4 * x*x - 10; return y; } float Calculate(float a,float b) { c = (a + b) / 2; n++; if (GetY(c) == 0 || ((b - a) / 2) < 0.000005) { cout << c <<"为方程的解"<< endl; return 0; } if (GetY(a)*GetY(c) < 0) { return Calculate(a,c); } if (GetY(c)*GetY(b)< 0) { return Calculate(c,b); } } }; int main() { cout << "方程组为：f(x)=x^3+4x^2-10=0" << endl; float a, b; Text text; text.Getab(); a = text.a; b = text.b; text.Calculate(a, b); return 0; } 2.割线法： // 方程求根（割线法）.cpp : 定义控制台应用程序的入口点。// #include "stdafx.h" #include"iostream"
心得体会 使用不同的方法，可以不同程度的求得方程的解，通过二分法计算的程序实现更加了解二分法的特点，二分法过程简单，程序容易实现，但该方法收敛比较慢一般用于求根的初始近似值，不同的方法速度不同。面对一个复杂的问题，要学会简化处理步骤，分步骤一点一点的循序处理，只有这样，才能高效的解决一个复杂问题。
《数值分析》总复习题-2013年-附部分答案
工程硕士《数值分析》总复习题（2013年用） [由教材中的习题、例题和历届考试题选编而成,供教师讲解和学生复习用] 疑课上的笔记整理，如有错漏，欢迎指出；碍于本人水平有限，部分题目未有解答。祝各位考试顺利！ 一. 解答下列问题: 1）下列所取近似值有多少位有效数字( 注意根据什么? ): a) 对 e = 2.718281828459045…,取* x = 2.71828 ( 答: 6 位 (因为它是按四舍五入来的) ) b) 数学家祖冲之取 113355 作为π的近似值. ( 答: 7 位 ( 按定义式 621113355 10 1415929.31415926.3-?≤-=-ΛΛπ 推得 ) ) c) 经过四舍五入得出的近似值12345，-0.001, 90.55000, 它们的有效 数字位数分别为 5 位, 1 位, 7 位。 2) 简述下名词: a) 截断误差 (不超过60字) (见书P.5) 答：它是指在构造数值计算方法时，用有限过程代替无限过程或用容易计算 的方法代替不容易计算的方法，其计算结果所存在的误差 b) 舍入误差 (不超过60字) (见书P.6) 答：对原始数据、中间计算结果和最后计算结果，都只能取有限位数表示， 这就要求进行“舍入”，这时所产生的误差就是舍入误差。 c) 算法数值稳定性 (不超过60字) (见书P.9) 答：是指算法在执行过程中，某阶段所产生的小误差在随后的阶段中不会被 积累或放大，从而不会严重降低全部计算的精确度。 3) 试推导( 按定义或利用近似公式 ): 计算3 x 时的相对误差约等于x 的相对 误差的3倍。 (参考书P.7例1.2.3)
数值计算方法学习心得
数值计算方法学习心得 ------一个代码的方法是很重要，一个算法的思想也很重要，但 在我看来，更重要的是解决问题的方法，就像爱因斯坦说的内容比 思维本身更重要。 我上去讲的那次其实做了挺充分的准备，程序的运行，pdf文档，算法公式的推导，程序伪代码，不过有一点缺陷的地方，很多细节 没有讲的很清楚吧，下来之后也是更清楚了这个问题。 然后一学期下来，总的来说，看其他同学的分享，我也学习到 许多东西，并非只是代码的方法，更多的是章胜同学的口才，攀忠 的排版，小冯的深入挖掘…都是对我而言比算法更加值得珍惜的东西，又骄傲地回想一下，曾同为一个项目组的我们也更加感到做项 目对自己发展的巨大帮助了。 同时从这些次的实验中我发现以前学到的很多知识都非常有用。 比如说，以前做项目的时候，项目导师一直要求对于要上传的 文件尽量用pdf格式，不管是ppt还是文档，这便算是对产权的一种 保护。 再比如代码分享，最基础的要求便是——其他人拿到你的代码 也能运行出来，其次是代码分享的规范性，像我们可以用轻量级Ubuntu Pastebin，以前做过一小段时间acm，集训队里对于代码的分享都是推荐用这个，像数值计算实验我觉得用这个也差不多了，其 次项目级代码还是推荐github（被微软收购了），它的又是可能更 多在于个人代码平台的搭建，当然像readme文档及必要的一些数据 集放在上面都更方便一些。
然后在实验中，发现debug能力的重要性，对于代码错误点的 正确分析，以及一些与他人交流的“正规”途径，讨论算法可能出 错的地方以及要注意的细节等，比如acm比赛都是以三人为一小组，讨论过后，讲了一遍会发现自己对算法理解更加深刻。 然后学习算法，做项目做算法一般的正常流程是看论文，尽量 看英文文献，一般就是第一手资料，然后根据论文对算法的描述， 就是如同课上的流程一样，对算法进一步理解，然后进行复现，最 后就是尝试自己改进。比如知网查询牛顿法相关论文，会找到大量 可以参考的文献。 最后的最后，想说一下，计算机专业的同学看这个数值分析， 不一定行云流水，但肯定不至于看不懂写不出来，所以我们还是要 提高自己的核心竞争力，就是利用我们的优势，对于这种算法方面 的编程，至少比他们用的更加熟练，至少面对一个问题，我们能思 考出对应问题的最佳算法是哪一个更合适解决问题。 附记: 对课程的一些小建议： 1. debug的能力不容忽视，比如给一个关于代码实现已知错误的代码给同学们，让同学们自己思考一下，然后分享各自的debug方法，一步一步的去修改代码，最后集全班的力量完成代码的debug，这往往更能提升同学们的代码能力。 2. 课堂上的效率其实是有点低的，可能会给学生带来一些负反馈，降低学习热情。 3. 总的来说还是从这门课程中学到许多东西。 数值分析学习心得体会
数值计算方法实验5
实验报告 学院（系）名称： 主程序部分列选主元部分
实验结果： 一．列主元消去法 输入各个数据，最终使用列选主元法，得到结果为：x1=x2=x3=1二．高斯-赛德尔迭代法 输入各个数据，输出每一步迭代数据，最终结果为：x1=0.285716，附录（源程序及运行结果） 一．列主元高斯消去法 #include MATLAB与数值分析课程总结 姓名：董建伟 学号：2015020904027 一：MATLAB部分 1.处理矩阵－容易 矩阵的创建 （1）直接创建注意 a中括号里可以用空格或者逗号将矩阵元素分开 b矩阵元素可以是任何MATLAB表达式，如实数复数等 c可以调用赋值过的任何变量，变量名不要重复，否则会被覆盖 （2）用MATLAB函数创建矩阵如：a空阵［］ b rand/randn——随机矩阵 c eye——单位矩阵 d zeros ——0矩阵 e ones——1矩阵 f magic——产生n阶幻方矩阵等 向量的生成 （1）用冒号生成向量 （2）使用linspace和logspace分别生成线性等分向量和对 数等分向量 矩阵的标识和引用 （1）向量标识 （2）“0 1”逻辑向量或矩阵标识 （3）全下标，单下标，逻辑矩阵方式引用 字符串数组 （1）字符串按行向量进行储存 （2）所有字符串用单引号括起来 （3）直接进行创建 矩阵运算 （1）注意与数组点乘，除与直接乘除的区别，数组为乘方对应元素的幂 （2）左右除时斜杠底部靠近谁谁是分母 （3）其他运算如，inv矩阵求逆，det行列式的值, eig特征值，diag 对角矩阵 2.绘图－轻松 plot－绘制二维曲线 （1）plot(x)绘制以x为纵坐标的二维曲线 plot(x,y) 绘制以x为横坐标，y为纵坐标的二维曲线 x,y为向量或矩阵 （2）plot（x1,y1,x2,y2，。。。。。。）绘制多条曲线，不同字母代替不同颜色:b蓝色，y黄色，r红色，g绿色 （3）hold on后面的pl ot图像叠加在一起 hold off解除hold on命令，plot将先冲去窗口已有图形（4）在hold后面加上figure，可以绘制多幅图形 （5）subplot在同一窗口画多个子图 三维图形的绘制 (1)plot3(x,y,z,’s’) s是指定线型，色彩，数据点形的字 符串 (2)［X,Y］＝meshgrid(x,y)生成平面网格点 (3)mesh(x,y,z，c)生成三维网格点，c为颜色矩阵 (4)三维表面处理mesh命令对网格着色，surf对网格片着色 (5)contour绘制二维等高线 (6)axis(［x1,xu,y1,yu］)定义x,y的显示范围 3．编程－简洁 （1）变量命名时可以由字母，数字，下划线，但是不得包含空格和标点 （2）最常用的数据类型只有双精度型和字符型，其他数据类型只在特殊条件下使用 （3）为得到高效代码，尽量提高代码的向量化程度，避免使用循环结构 第一章绪论（1-4） 一、误差来源及分类 二、误差的基本概念 1.绝对误差及绝对误差限 2.相对误差及相对误差限 3.有效数字 三、数值计算的误差估计 1.函数值的误差估计 2.四则运算的误差估计 四、数值计算的误差分析原则 第二章插值（1.2.4-8） 一、插值问题的提法（定义）、插值条件、插值多项式的存在唯一性 二、拉格朗日插值 1.拉格朗日插值基函数的定义、性质 2.用拉格朗日基函数求拉格朗日多项式 3.拉格朗日插值余项（误差估计） 三、牛顿插值 1.插商的定义、性质 2.插商表的计算 3.学会用插商求牛顿插值多项式 四、等距节点的牛顿插值 1.差分定义、性质及计算（向前、向后和中心） 2.学会用差分求等距节点下的牛顿插值公式 五、学会求低次的hermite插值多项式 六、分段插值 1.分段线性插值 2.分段三次hermite插值 3.样条插值 第三章函数逼近与计算（1-6） 一、函数逼近与计算的提法（定义）、常用两种度量标准（一范数、二范数\平方逼近） 二、基本概念 连续函数空间、最佳一次逼近、最佳平方逼近、内积、内积空间、偏差与最小偏差、偏差点、交错点值、平方误差 三、学会用chebyshev定理求一次最佳一致逼近多项式，并估计误差（最大偏差） 四、学会在给定子空间上通过解方程组求最佳平方逼近，并估计误差（平方误差） 五、正交多项式（两种）定义、性质，并学会用chebyshev多项式性质求特殊函数的（降阶）最佳一次逼近多项式 六、函数按正交多项式展开求最佳平方逼近多项式，并估计误差 七、一般最小二乘法（多项式拟合）求线性拟合问题 第四章数值分析（1-4） 一、数值求积的基本思想及其机械求积公式 一、单项选择题（每小题3分，共15分） 1. 3.142和3.141分别作为的近似数具有（ ）和（ ）位有效数字. A ．4和3 B ．3和2 C ．3和4 D ．4和4 2. 已知求积公式 ，则＝（ ） A ． B ． C ． D ． 3. 通过点 的拉格朗日插值基函数满足（ ） A ． ＝0， B ． ＝0， C ．＝1， D ． ＝1， 4. 设求方程 的根的牛顿法收敛，则它具有（ ）敛速。 A ．超线性 B ．平方 C ．线性 D ．三次 5. 用列主元消元法解线性方程组 作第一次消元后得到的第3个方程（ ）. A ． B ． C ． D ． π()()2 1 121 1()(2)636f x dx f Af f ≈ ++? A 1613122 3()()0011,,,x y x y ()()01,l x l x ()00l x ()110l x =() 00l x ()111 l x =() 00l x ()111 l x =() 00l x ()111 l x =()0 f x =12312312 20 223332 x x x x x x x x ++=?? ++=??--=?232 x x -+=232 1.5 3.5 x x -+=2323 x x -+= 单项选择题答案 1.A 2.D 3.D 4.C 5.B 二、填空题（每小题3分，共15分） 1. 设, 则 ， . 2. 一阶均差 3. 已知时，科茨系数 ，那么 4. 因为方程 在区间 上满 足 ，所以 在区间内有根。 5. 取步长，用欧拉法解初值问题 的计算公 式 . 填空题答案 230.5 1.5 x x -=-T X )4,3,2(-==1||||X 2||||X =()01,f x x = 3n =()()() 33301213,88C C C === () 3 3C =()420 x f x x =-+=[]1,2()0 f x =0.1h =()211y y y x y ?'=+?? ?=? 数值分析心得体会 篇一：学习数值分析的经验 数值分析实验的经验、感受、收获、建议班级：计算131 学号：XX014302 姓名：曾欢欢 数值分析实验主要就是学习MATLAB的使用以及对数值分析类容的应用，可以使学生更加理解和记忆数值分析学得类容，也巩固了MATLAB的学习，有利于以后这个软件我们的使用。在做实验中，我们需要具备较好的编程能力、明白MATLAB软件的使用以及掌握数值分析的思想，才能让我们独立自主的完成该作业，如果是上述能力有限的同学，需要借助MATLAB的书以及网络来完成实验。数值分析实验对于我来说还是有一定难度，所以我课下先复习了MATLAB的使用方法以及编写程序的基本类容，借助互联网和同学老师资源完成了数值分析得实验的内容。在实验书写中，我复习了各种知识，所以我认为这门课程是有必要且是有用处的，特别是需要处理大量实验数据的人员，很有必要深入了解学习它，这样在以后的工作学习里面就减少了很多计算问题也提高了实验结果的精确度。 学习数值分析的经验、感受、收获、建议数值分析的内容包括插值与逼近，数值微分与数值积分，非线性方程与线性方程组的数值解法，矩阵的特征值与特征向量计算，常微分方程数值解等。 首先我们必须明白数值分析的用途。通常所学的其他数学类学科都是由公式定理开始，从研究他们的定义，性质再到证明与应用。但实际上，尤其是工程，物理，化学等其它具体的学科。往往我们拿到 手的只是通过实验得到的数据。如果是验证性试验，需要代回到公式 进行分析，验证。但往往更多面对的是研究性或试探性试验，无具体 公式定理可代。那就必须通过插值，拟合等计算方法进行数据处理以得到一个相对可用的一般公式。还有许多计算公式理论上非常复杂，在工程中不实用，所以必须根据实际情况把它转化成多项式近似表 示。学习数值分析，不应盲目记公式，因为公事通常很长且很乏味。其次，应从公式所面临的问题以及用途出发。比如插值方法，就 是就是把实验所得的数据看成是公式的解，由这些解反推出一个近似公式，可以具有局部一般性。再比如说拟合，在插值的基础上考虑实 验误差，通过拟合能将误差尽可能缩小，之后目的也是得到一个具有 一定条件下的一般性的公式。。建议学习本门课程要结合知识与实际，比如在物理实验里面很多MATLAB与数值分析课程总结
数值分析知识点
数值计算方法期末考试题
数值分析心得体会