1
习题 三
1. 用辛普森求积公式计算积分dx e x ∫?1
0，并估计误差。 2. 给定数据表
x 1.8 2.0 2.2 2.4 2.6
)(x f 3.12014 4.042569 6.04241 8.03014 10.46675
分别用复化梯形求积公式和复化辛普森求积公式计算积分
.)(6
.28.1dx x f ∫
3. 分别用变步长求积方法和龙贝格求积方法计算下列积分，并估计误差：
（1）;sin 4
0dx x
x ∫π
（2）;1)1ln(102dx x x ∫++ （3）;)1ln(110dx x x ∫+ （4）.110∫+x
dx 4. 确定下列求积公式的待定参数，使其代数精度尽量高，并指出其代数精度的准确次数：
（1）);()0()()(210h f A f A h f A dx x f h
h ++?≈∫? （2）);()0()()(21022h f A f A h f A dx x f h
h ++?≈∫?
（3）).0()1()0()(2101
0f A f A f A dx x f ′++≈∫ 5. 证明求积公式 )]()([12
)]()([2)(012101
0x f x f h x f x f h dx x f x x ′?′++≈∫ 具有3次代数精度，其中01x x h ?=.
6. 利用高斯-勒让德公式计算积分
)(141
02π=+∫dx x 的近似值. 7. 利用高斯-切比雪夫求积公式计算积分
dx x x ∫???112211 的近似值。
8. 已知函数)(x f y =的如下数据：
x 0.8 0.9 1 1.1 1.2 1.4
)(x f 0.8365 0.9095 1 1.1105 1.2446 1.6017
利用李查逊外推法计算)1(y ′.
数值分析试题 一、 填空题（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 )
硕士研究生课程报告 题目顺层高边坡稳定性影响因素 及工程灾害防治 姓名曾义 专业班级岩土13级 任课教师阳军生张学民 中南大学土木工程学院
引言 近年来，随着铁路公路建设步伐加快，铁路公路等级不断提高，边坡防护建设工程中所遇到的岩土边坡安全稳定性问题也相应增多，并成为岩土工程中比较常见的技术难题。由于工程建设的需要，往往在一定程度上破坏或扰动原来较为稳定的岩土体而形成新的人工边坡，因而普遍存在着边坡稳定的问题需要解决。国家实施西部大开发战略以来，西部山区高等级公路得到迅速发展。在山区修建高等级公路不可避免会遇到大量的深挖高填路基，就目前建设的高速公路情况看:一般情况下，100km长的山区高等级公路，挖填方路基段落长度占路线总长度的60%以上。已建高速公路最高的填方已达到50多米,最高的挖方边坡高度已超过100m。尽管山区高等级公路的建设越来越倡导环境保护，尽量避免深挖高填，但路基作为公路的主要结构，其边坡稳定问题不可避免。在山区复杂多变的地质条件下建设高等级公路，其边坡稳定性问题必将受到人们的普遍关注，高边坡岩土安全状况直接关系到公路交通运输安全。 虽然计算理论方法、地质探测技术、现代监测技术、边坡加固技术及施工技术不断的在进步，但顺层边坡稳定性问题和高边坡稳定性问题，时至今日依然是国内外学者研究的热点问题，并逐步涌现出许多的新的研究方向。 1、顺倾高边坡稳定性研究现状 随着人类工程活动的发展，对边坡问题的研究也在不断深入，归纳前人对边坡问题的研究大致可分为以下几个阶段： 人们对边坡稳定性的关注和研究最早是从滑坡现象开始的（张倬元等，2001）。19世纪末和20世纪初期，伴随着欧美资本主义国家的工业化而兴起的大规模土木工程建设（如修筑铁路、公路，露天采矿，天然建材开采等），出现了较多的人工边坡，诱发了大量滑坡和崩塌，造成了很大的损失。这时，人们才开始重视边坡失稳给人类造成的危害，并开始借用一般材料分析中的工程力学理论对滑坡进行半经验、半理论的研究。 20世纪50年代，我国学者引进苏联工程地质的体系，继承和发展了“地质历史分析”法，并将其应用于滑坡的分析和研究中，对边坡稳定性研究起到了推动作用（张倬元等，1994）。该阶段学者们着重边坡地质条件的描述和边坡类型的划分，采用工程地质类比法评价边坡稳定性。 20世纪60年代，世界上几起灾难性的边坡失稳事件的发生（如意大利的瓦依昂滑坡造成近3000人死亡和巨大的经济损失）（张倬元等，1994）,使人们逐渐认识到了结构面对边坡稳定性的控制作用以及边坡失稳的时效特征，初步形
二 1 求A的LU分解，并利用分解结果求 解由紧凑格式 故 从而 故 2求证：非奇异矩阵不一定有LU分解 证明设非奇异，要说明A不一定能做LU分解，只需举出一个反例即可。现考虑矩阵，显然A为非奇异矩阵。若A有LU分解，则 故，而，显然不能同时成立。这矛盾说明A不能做LU分解，故只假定A非奇异并不能保证A能做LU分解，只有在A的前阶顺序主子式时才能保证A一定有LU分解。 3用追赶法求解如下的三对角方程组 解设有分解 由公式 其中分别是系数矩阵的主对角线元素及其下边和上边的次对角线元素，故有 从而有 故，，， 故，，， 4设A是任一阶对称正定矩阵，证明是一种向量范数 证明（1）因A正定对称，故当时，，而当时， （2）对任何实数，有 （3）因A正定，故有分解，则 故对任意向量和，总有 综上可知，是一种向量范数。 5 设，，已知方程组的精确解为 （1）计算条件数； （2）若近似解，计算剩余； （3）利用事后误差估计式计算不等式右端，并与不等式左边比较，此结果说明了什么？解（1） （2）
（3）由事后误差估计式，右端为 而左端 这表明当A为病态矩阵时，尽管剩余很小，误差估计仍然较大。因此，当A病态时，用大小作为检验解的准确度是不可靠的。 6矩阵第一行乘以一数成为，证明当时，有最小值 证明设，则 又 故 从而当时，即时，有最小值，且 7 讨论用雅可比法和高斯-赛德尔法解方程组时的收敛性。如果收敛，比较哪一种方法收敛较快，其中 解对雅可比方法，迭代矩阵 ， 故雅可比法收敛。 对高斯-赛德尔法，迭代矩阵 ，故高斯-赛德尔法收敛。 因=故高斯-赛德尔法较雅可比法收敛快。 8设，求解方程组，求雅可比迭代法与高斯-赛德尔迭代法收敛的充要条件。 解雅可比法的迭代矩阵 ， 故雅可比法收敛的充要条件是。 高斯-赛德尔法的迭代矩阵 ， 故高斯-赛德尔法收敛的充要条件是。 9 设求解方程组的雅可比迭代格式为，其中，求证：若，则相应的高斯-赛德尔法收敛。证明由于是雅可比法的迭代矩阵，故 又，故， 即，故故系数矩阵A按行严格对角占优，从而高斯-赛德尔法收敛。 10设A为对称正定矩阵，考虑迭代格式 求证：（1）对任意初始向量，收敛； （2）收敛到的解。 证明（1）所给格式可化为 这里存在是因为，由A对称正定，，故也对称正定。 设迭代矩阵的特征值为，为相应的特征向量，则与做内积，有 因正定，故，从而，格式收敛。
1、（本题5分）试确定7 22 作为π的近似值具有几位有效数字，并确定其相对误差限。 解 因为 7 22 =3.142857…=1103142857 .0-? π=3.141592… 所以 312102 11021005.0001264.0722--?=?=<=- π （2分） 这里，3,21,0=-=+-=n n m m 由有效数字的定义可知7 22 作为π的近似值具有3位有效数字。 （1分） 而相对误差限 3102 1 0005.00004138.0001264.07 22-?= <≈= -= π π πε r （2分） 2、（本题6分）用改进平方根法解方程组：??? ?? ??=????? ??????? ??--654131*********x x x ； 解 设???? ? ??????? ? ?????? ??===????? ??--11111 1 131321112323121 32 132 31 21 l l l d d d l l l LDL A T 由矩阵乘法得： 5 7,21,215 27 ,25,2323121321- ==-== -==l l l d d d （3分） 由y D x L b Ly T 1 ,-==解得 T T x y )9 23 ,97,910(,)563, 7,4(== （3分） 3、（本题6分）给定线性方程组???????=++-=+-+=-+-=-+17 7222382311387 510432143213 21431x x x x x x x x x x x x x x 1）写出Jacoib 迭代格式和Gauss-Seidel 迭代格式； 2）考查Jacoib 迭代格式和Gauss-Seidel 迭代格式的敛散性； 解 1）Jacoib 迭代格式为
自动化工程训练 —基于MATLAB的电力电子系统仿真 学院：信息科学与工程学院 仿真内容：三相桥式整流电路 班级姓名：自动化0801 肖娉 学号：0909080320 指导老师：桂武鸣老师 日期：2011.08.29--2011.09.09
电力电子技术综合了电子电路、电机拖动、计算机控制等多学科知识，是一门实践性和应用性很强的课程。由于电力电子器件自身的开关非线性，给电力电子电路的分析带来了一定的复杂性和困难，一般常用波形分析的方法来研究。仿真技术为电力电子电路的分析提供了崭新的方法。 本次工程训练的目的是初步掌握在MA TLAB/Simulink环境下电力电子系统的仿真。通过为期两周的学习，掌握一些MA TLAB的基础、Simulink环境和模型库、电力电子器件模型、变压器和电动机模型等。 MATLAB是一种科学计算软件，它是一种以矩阵为基础的交互式程序计算语言。SIMULINK是基于框图的仿真平台，它挂接在MATLAB环境上，以MATLAB的强大计算功能为基础，以直观的模块框图进行仿真和计算。 本文主要以MATLAB/SIMULINK仿真软件为基础，完成了对三相桥式整流电路带电阻、阻感、反电动势、直流电机负载的建模与仿真，并且给出了仿真结果波形，同时根据仿真结果进行了分析。证实了该方法的简便直观、高效快捷和真实准确性。
前言 第一章MATLAB/Simulink仿真的目的与意义 (1) 第二章MATLAB/Simulink的基础知识 (2) 2.1 MATLAB基础 (2) 2.1.1 MATLAB语言的功能 (2) 2.2.2 MATLAB集成环境 (3) 2.2 Simulink仿真基础 (5) 2.2.1 Simulink的模块库介绍 (6) 2.2.2 SimPowerSystems的介绍 (6) 2.2.3 Simulink部分模型介绍 (7) 2.2.4 Simulink仿真运行 (8) 第三章三相桥式可控整流电路的仿真 (10) 3.1 三相桥式整流电路 (10) 3.1 电阻、阻感和反电动势负载 (11) 3.2 直流电机负载 (16) 3.2.1 整流状态 (16) 3.2.2 有源逆变状态 (18) 第四章心得体会 (21) 参考文献 (23)
二 1求A的LU分解，并利用分解结果求 解由紧凑格式 故 从而 故 2求证：非奇异矩阵不一定有LU分解 证明设非奇异，要说明A不一定能做LU分解，只需举出一个反例即可。现考虑矩阵，显然A为非奇异矩阵。若A有LU分解，则 故，而，显然不能同时成立。这矛盾说明A不能做LU分解，故只假定A非奇异并不能保证A能做LU分解，只有在A的前阶顺序主子式 时才能保证A一定有LU分解。
3用追赶法求解如下的三对角方程组 解设有分解 由公式 其中分别是系数矩阵的主对角线元素及其下边和上边的次对角线元素，故有 从而有 故，，， 故，，，
4设A是任一阶对称正定矩阵，证明是一种向量范数 证明（1）因A正定对称，故当时，，而当时， （2）对任何实数，有 （3）因A正定，故有分解，则 故对任意向量和，总有 综上可知，是一种向量范数。 5 设，，已知方程组的精确解为 （1）计算条件数； （2）若近似解，计算剩余； （3）利用事后误差估计式计算不等式右端，并与不等式左边比较，此结果说明了什么？解（1） （2） （3）由事后误差估计式，右端为 而左端
这表明当A为病态矩阵时，尽管剩余很小，误差估计仍然较大。因此，当A病态时，用大小作为检验解的准确度是不可靠的。 6矩阵第一行乘以一数成为，证明当时，有最小值 证明设，则 又 故 从而当时，即时，有最小值，且 7讨论用雅可比法和高斯-赛德尔法解方程组时的收敛性。如果收敛，比较哪一种方 法收敛较快，其中 解对雅可比方法，迭代矩阵 ， 故雅可比法收敛。 对高斯-赛德尔法，迭代矩阵
，故高斯-赛德尔法收敛。 因=故高斯-赛德尔法较雅可比法收敛快。 8设，求解方程组，求雅可比迭代法与高斯-赛德尔迭代法收敛的充要条件。 解雅可比法的迭代矩阵 ， 故雅可比法收敛的充要条件是。 高斯-赛德尔法的迭代矩阵 ，
数值分析整理版试题及答案
例1、 已知函数表 x -1 1 2 ()f x -3 0 4 求()f x 的Lagrange 二次插值多项式和Newton 二次插值多项式。 解： （1）k x -1 1 2 k y -3 0 4 插值基函数分别为 ()()()()()()()()()() 1200102121()1211126 x x x x x x l x x x x x x x ----= ==-------- ()()()()()()()() ()()021******* ()1211122x x x x x x l x x x x x x x --+-= ==-+---+- ()()()()()()()()()()0122021111 ()1121213 x x x x x x l x x x x x x x --+-= ==-+--+- 故所求二次拉格朗日插值多项式为 () ()()()()()()()()()()2 20 2()11131201241162314 121123537623k k k L x y l x x x x x x x x x x x x x ==?? =-? --+?-+-+?+-????=---++-=+-∑ （2）一阶均差、二阶均差分别为
[]()()[]()()[][][]010********* 011201202303 ,11204 ,412 3 4,,5 2,,126 f x f x f x x x x f x f x f x x x x f x x f x x f x x x x x ---===-----= = =----=== --- k x ()k f x 一阶 二阶 -1 -3 1 0 3/ 2 2 4 4 5/6 故所求Newton 二次插值多项式为 ()()[]()[]()() ()()()20010012012,,,35 311126537623P x f x f x x x x f x x x x x x x x x x x x =+-+--=-+ +++-=+- 例2、 设2 ()32f x x x =++，[0,1]x ∈，试求()f x 在[0, 1]上关于()1x ρ=，{} span 1,x Φ=的最佳平方逼近多项式。 解： 若{}span 1,x Φ=，则0()1x ?=，1()x x ?=，且()1x ρ=，这样，有
硕士研究生培养方案（科学学位） 一、学科概况 中南大学机械工程学科创建于1955年，1960年招收研究生，1982年获得硕士学位授予权，1986年获博士学位授予权，1998年设立“机械工程”博士后科研流动站，2000年获得一级学科博士授予权，覆盖了机械制造及自动化、机械设计及理论、机械电子工程和车辆工程等4个二级学科和数字装备与计算制造、信息器件制造技术与装备等2个自主设置的二级学科，其中“机械设计及理论”与“机械制造及其自动化”学科为国家重点学科，“机械制造及其自动化”与“机械电子工程”学科为湖南省重点学科，机械工程一级学科于2007年被批准为一级学科国家重点学科。设有“高性能复杂制造”国家重点实验室，“现代复杂装备设计与极端制造”教育部重点实验室，“铝合金强流变技术与装备”教育部工程研究中心，湖南省“岩土设备设计与控制”工程研究中心，以及“金属塑性加工摩擦与润滑”、“设备测试与故障诊断中心”等1个国家重点实验室和5个省部级重点实验室、工程中心，以及国家高技术研究发展计划成果产业化基地、与国外ASM公司共建的“微电子封装技术实验室”等。 本学科致力于机械基础理论与技术集成、先进制造理论与技术等的研究，并围绕国民经济中起支柱作用以及国防和空天运载等关键技术与装备进行研究和
设计开发，在高性能材料制备与装备、信息器件制造、齿轮数字化制造、深海资源开发、车辆与工程装备、特种机器人等研究方向具有特色和优势。 二、培养目标 学位获得者应拥护中国共产党的领导，拥护社会主义制度，热爱祖国，掌握辩证唯物主义和历史唯物主义的基本原理；具有良好的科研作风、科学道德和合作精神，品行优秀，身心健康；掌握机械工程学科坚实的基础理论、系统的专门知识，掌握一定的生产实践及试验方面的知识和技能，熟练掌握一门外语，了解本学科前沿发展动态和方向，有严谨求实的工作作风和独力工作能力。成为既能从事机械工程领域的科学研究与设计工作，又可承担相关领域的教学和管理工作的高层次、高素质的科技人才。 三、学科专业主要研究方向
1. 已知325413.0,325413*2*1==X X 都有6位有效数字，求绝对误差限。（4分） 解： 由已知可知6 5.0102 1 ,0,6,10325413.0016*1=?= =-=?=ε绝对误差限n k k X 2分 620*2102 1 ,6,0,10325413.0-?= -=-=?=ε绝对误差限n k k X 2分 2. 已知?? ???=0 01 A 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分 ?????=0 1 A A T 4 2 ???? ? -420?????0 01 2 20 - ???? ?440= ?????0 01 80 ???? ?3200 2分 {}32 32,8,1max )(max ==A A T λ
1分 24322==A 3. 设32)()(a x x f -= （6分） ① 写出f(x)=0解的迭代格式 ② 当a 为何值时，)(1k k x x ?=+ （0,1……）产生的序列{}k x 收敛于 2 解： ①迭代格式为： 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. 给定线性方程组，其中：?? ?=13A ?? ?2 2，?? ? ???-=13b 用迭代公式 )()()()1(k k k Ax b x x -+=+α（0,1……）求解，问取什么实数α ，可使 迭代收敛 （8分） 解： 所给迭代公式的迭代矩阵为?? ? --???--=-=ααααα21231A I B 2分
070102 计算数学 计算数学也叫做数值计算方法或数值分析。主要内容包括代数方程、线性代数方程组、微分方程的数值数值逼近问题，矩阵特征值的求法，最优化计算问题，概率统计计算问题等等，还包括解的存在性、唯一性差分析等理论问题。我们知道五次及五次以上的代数方程不存在求根公式，因此，要求出五次以上的高次代一般只能求它的近似解，求近似解的方法就是数值分析的方法。对于一般的超越方程，如对数方程、三角方采用数值分析的办法。怎样找出比较简洁、误差比较小、花费时间比较少的计算方法是数值分析的主要课题的办法中，常用的办法之一是迭代法，也叫做逐次逼近法。迭代法的计算是比较简单的，是比较容易进行的以用来求解线性方程组的解。求方程组的近似解也要选择适当的迭代公式，使得收敛速度快，近似误差小。 在线性代数方程组的解法中，常用的有塞德尔迭代法、共轭斜量法、超松弛迭代法等等。此外，一些比消去法，如高斯法、追赶法等等，在利用计算机的条件下也可以得到广泛的应用。在计算方法中，数值逼近本方法。数值逼近也叫近似代替，就是用简单的函数去代替比较复杂的函数，或者代替不能用解析表达式表值逼近的基本方法是插值法。 初等数学里的三角函数表，对数表中的修正值，就是根据插值法制成的。在遇到求微分和积分的时候，的函数去近似代替所给的函数，以便容易求到和求积分，也是计算方法的一个主要内容。微分方程的数值解法。常微分方程的数值解法由欧拉法、预测校正法等。偏微分方程的初值问题或边值问题，目前常用的是有限元素法等。有限差分法的基本思想是用离散的、只含有限个未知数的差分方程去代替连续变量的微分方程求出差分方程的解法作为求偏微分方程的近似解。有限元素法是近代才发展起来的，它是以变分原理和剖分的方法。在解决椭圆形方程边值问题上得到了广泛的应用。目前，有许多人正在研究用有限元素法来解双曲方程。计算数学的内容十分丰富，它在科学技术中正发挥着越来越大的作用。 排名学校名称等级 1 北京大学A+ 2 浙江大学 A+ 3 吉林大学A+ 4 大连理工大学A+ 5 西安交通大学A 北京大学：http:https://www.doczj.com/doc/f54963347.html,/NewsSpecialDetailsInfo.aspx?SID=4 浙江大学：http:https://www.doczj.com/doc/f54963347.html,/NewsSpecialDetailsInfo.aspx?SID=21847 吉林大学：http:https://www.doczj.com/doc/f54963347.html,/NewsSpecialDetailsInfo.aspx?SID=5506 大连理工大学：http:https://www.doczj.com/doc/f54963347.html,/NewsSpecialDetailsInfo.aspx?SID=4388 西安交通大学：http:https://www.doczj.com/doc/f54963347.html,/NewsSpecialDetailsInfo.aspx?SID=18285
1、（本题5分）试确定7 22作为π的近似值具有几位有效数字，并确定其相对误差限。 解 因为 7 22=3.142857…=1103142857.0-? π=3.141592… 所以 3 12 10 2 110 21005.0001264.07 22--?= ?= <=- π （2分） 这里，3,21,0=-=+-=n n m m 由有效数字的定义可知7 22作为π的近似值具有3位有效数字。 （1分） 而相对误差限 3 10 2 10005.00004138.0001264.07 22-?= <≈= -= π π πε r （2分） 2、（本题6分）用改进平方根法解方程组：???? ? ??=????? ??????? ??--654131321 112321x x x ； 解 设???? ? ? ?????? ? ?????? ??===????? ? ?--11 1 11113 1321 11232312132 1 32 31 21 l l l d d d l l l LDL A T 由矩阵乘法得： 5 7,21,21527,25,2323121321- == - == -==l l l d d d （3分） 由y D x L b Ly T 1 ,-==解得 T T x y )9 23,97,910( ,)5 63, 7,4(== （3分） 3、（本题6分）给定线性方程组??? ? ? ??=++-=+-+=-+-=-+17722238231138751043214321 321431x x x x x x x x x x x x x x 1）写出Jacoib 迭代格式和Gauss-Seidel 迭代格式； 2）考查Jacoib 迭代格式和Gauss-Seidel 迭代格式的敛散性； 解 1）Jacoib 迭代格式为
例1、 已知函数表 求()f x 的Lagrange 二次插值多项式和Newton 二次插值多项式。 解： （1） 故所求二次拉格朗日插值多项式为 （2）一阶均差、二阶均差分别为 例2、 设2 ()32f x x x =++，[0,1]x ∈，试求()f x 在[0,1]上关于()1x ρ=，{}span 1,x Φ=的最佳平 方逼近多项式。 解： 若{}span 1,x Φ=，则0()1x ?=，1()x x ?=，且()1x ρ=，这样，有 所以，法方程为
011231261192 34a a ??????????=?????????? ?????????? ，经过消元得012311 62110123a a ??? ???????=???????????????????? 再回代解该方程，得到14a =，011 6 a = 故，所求最佳平方逼近多项式为* 111()46S x x =+ 例3、 设()x f x e =，[0,1]x ∈，试求()f x 在[0,1]上关于()1x ρ=，{}span 1,x Φ=的最佳平方逼近 多项式。 解： 若{}span 1,x Φ=，则0()1x ?=，1()x x ?=，这样，有 所以，法方程为 解法方程，得到00.8732a =，1 1.6902a =， 故，所求最佳平方逼近多项式为 例4、 用4n = 的复合梯形和复合辛普森公式计算积分1 ? 。 解： （1）用4n =的复合梯形公式 由于2h =，( )f x =()121,2,3k x k k =+=，所以，有 （2）用4n =的复合辛普森公式 由于2h =，( )f x =()121,2,3k x k k =+=，()12 220,1,2,3k x k k + =+=，所以，有 例5、 用列主元消去法求解下列线性方程组的解。 解：先消元 再回代，得到33x =，22x =，11x = 所以，线性方程组的解为11x =，22x =，33x = 例6、 用直接三角分解法求下列线性方程组的解。 解： 设 则由A LU =的对应元素相等，有 1114u = ，1215u =，1316u =， 2111211433l u l =?=，3111311 22 l u l =?=， 2112222211460l u u u +=?=-，2113232311 545l u u u +=?=-，
一、单项选择题（每小题3分，共15分） 1. 3.142和3.141分别作为π的近似数具有（ ）和（ ）位有效数字. A ．4和3 B ．3和2 C ．3和4 D ．4和4 2. 已知求积公式 ()()2 1 121 1()(2)636f x dx f Af f ≈ ++? ，则A ＝（ ） A ． 16 B ．13 C ．12 D ．2 3 3. 通过点 ()()0011,,,x y x y 的拉格朗日插值基函数()()01,l x l x 满足（ ） A ． ()00l x ＝0， ()110l x = B ． ()00l x ＝0， ()111l x = C ．() 00l x ＝1，()111 l x = D ． () 00l x ＝1，()111 l x = 4. 设求方程 ()0 f x =的根的牛顿法收敛，则它具有（ ）敛速。 A ．超线性 B ．平方 C ．线性 D ．三次 5. 用列主元消元法解线性方程组 1231231 220223332 x x x x x x x x ++=?? ++=??--=? 作第一次消元后得到的第3个方程（ ）. A ． 232 x x -+= B ．232 1.5 3.5 x x -+= C ． 2323 x x -+= D ． 230.5 1.5 x x -=- 单项选择题答案 1.A 2.D 3.D 4.C 5.B 得 分 评卷人 二、填空题（每小题3分，共15分）
1. 设T X )4,3,2(-=, 则=1||||X ，2||||X = . 2. 一阶均差 ()01,f x x = 3. 已知3n =时，科茨系数()()() 33301213,88C C C ===，那么 () 33C = 4. 因为方程()420 x f x x =-+=在区间 []1,2上满足 ，所以()0f x =在区间 内有根。 5. 取步长0.1h =，用欧拉法解初值问题 ()211y y y x y ?'=+?? ?=? 的计算公式 . 填空题答案 1. 9和29 2. ()() 0101 f x f x x x -- 3. 1 8 4. ()()120 f f < 5. ()12 00.1 1.1,0,1,210.11k k y y k k y +???? ?=+? ?=+???? =??L 得 分 评卷人 三、计算题（每题15分，共60分） 1. 已知函数 21 1y x = +的一组数据： 求分 段线性插值函数，并计算 () 1.5f 的近似值. 计算题1.答案 1. 解 []0,1x ∈， ()1010.510.50110x x L x x --=?+?=---% []1,2x ∈，()210.50.20.30.81221x x L x x --=?+?=-+--%
第一章 绪论 一、填空题 1、 已知 71828.2e =，求x 的近似值a 的有效数位和相对误差： 题号 精确数x x 的近似数a a 的有效数位 a 的相对误差 ⑴ e 2.7 ⑵ e 2.718 ⑶ e/100 0.027 ⑷ e/100 0.02718 2、 设原始数据x 1，x 2，x 3和x 4的近似值（每位均为有效数字）如下： a 1=1.1021，a 2=0.031，a 3=385.6，a 4=56.430 则 ⑴ a 1+a 2+a 4= ，相对误差界为 ； ⑵ a 1a 2a 3= ，相对误差界为 ； ⑶ a 2/a 4= ，相对误差界为 。 二、为使20的近似值的相对误差小于0.01%，问应取多少位有效数字？ 三、当x 接近于0时，怎样计算 x x sin cos 1-以及当x 充分大时，怎样计算 x x -+1，才会使其结果的有效数字不会严重损失。 四、在数值计算中，为了减小误差，应该尽量避免的问题有哪些？并举出相 应的实例. 五、对于序列 ,1,0,9991 =+=? n dx x x I n n ，试构造两种递推算法计算 10I ，在你构造的算法中，那一种是稳定的，说明你的理由；
第二章 插值法 １、在互异的n+1个点处满足插值条件P(x i )=y i ,(i=0,1,…n)的次数不高于n 的 多项式是( )的 (Ａ)存在且唯一 (Ｂ)存在 (Ｃ)不存在 (Ｄ)不唯一 2、当f(x)是次数不超过n 的多项式时，f(x)的插值多项式是 ( ) (Ａ)不确定 (Ｂ)次数为n (Ｃ)f(x)自身 （D ）次数超过n 3、 插值基函数的和 ∑=n j j x l )(= ( ) (Ａ)０ (Ｂ)１ (Ｃ)２ (Ｄ)不确定 4、 设f(x)=x 3-x+5，则f[20,21,22,23]= ( )； f[20,21,22,23,24]= ( ) (Ａ)０ (Ｂ)１ (Ｃ)２ (Ｄ)不确定 5、( )插值方法具有公式整齐、程序容易实现的优点，而( )插值方法 计算灵活，如果节点个数变化时，不需要重新构造多项式，它们都是( )的方法 (Ａ)构造性 (Ｂ)解方程组 (Ｃ)拉格朗日 (Ｄ)牛顿 ６、一般地，内插公式比外推公式( ),高次插值比低次插值( )，但 当插值多项式的次数高于七、八次时，最好利用( )插值公式 (Ａ)粗糙 (Ｂ)精确 (Ｃ)分段低次 (Ｄ)高次 ７、整体光滑度高，收敛性良好，且在外型设计、数值计算中应用广泛的分 段插值方法为（ ）. （Ａ）分段线性插值 （Ｂ）分段抛物插值 （Ｃ）分段三次埃尔米特插值 （Ｄ）三次样条插值。 ８、差商与差分的关系式为 f[x 0,x 1,…,x k ]=（ ）,f[x n ,x n-1,…,x n-k ]=（ ）。 （Ａ）k n k h k f !? （Ｂ）k k h k f !0? （Ｃ）k n k h k f !? （Ｄ）k k h k f !0 ?
《数值计算方法》复习试题 一、填空题： 1、 ?? ??? ?????----=410141014A ，则A 的LU 分解为 A ? ???????? ???=????????? ?? ?。 答案： ?? ????????--??????????--=1556141501 4115401411A 2、已知3.1)3(,2.1)2(, 0.1)1(===f f f ，则用辛普生（辛卜生）公式计算求 得?≈3 1 _________ )(dx x f ,用三点式求得≈')1(f 。 答案：2.367，0.25 3、1)3(,2)2(, 1)1(==-=f f f ，则过这三点的二次插值多项式中2x 的系数 为 ，拉格朗日插值多项式为 。 答案：-1， )2)(1(21 )3)(1(2)3)(2(21)(2--------= x x x x x x x L 4、近似值*0.231x =关于真值229.0=x 有( 2 )位有效数字； 5、设)(x f 可微,求方程)(x f x =的牛顿迭代格式是( )； 答案 )(1)(1n n n n n x f x f x x x '--- =+ 6、对 1)(3++=x x x f ,差商=]3,2,1,0[f ( 1 ),=]4,3,2,1,0[f ( 0 )； 7、计算方法主要研究( 截断 )误差和( 舍入 )误差； 8、用二分法求非线性方程f (x )=0在区间(a ,b )内的根时，二分n 次后的误差限为( 1 2+-n a b )； 9、求解一阶常微分方程初值问题y '= f (x ,y )，y (x 0)=y 0的改进的欧拉公
数值分析试题及答案 一、单项选择题（每小题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． 单项选择题答案 1.A 2.D 3.D 4.C 5.B 得分评卷 人 二、填空题（每小题3分，共15分） 1. 设, 则， . 2. 一阶均差 3. 已知时，科茨系数，那么 4. 因为方程在区间上满足，所以在区间内有根。 5. 取步长，用欧拉法解初值问题的计算公式.填空题答案
1. 9和 2. 3. 4. 5. 得分评卷 人 三、计算题（每题15分，共60分） 1. 已知函数的一组数据：求分段线性插值函数，并计算的近似值. 计算题1.答案 1. 解， ， 所以分段线性插值函数为 2. 已知线性方程组 （1）写出雅可比迭代公式、高斯－塞德尔迭代公式； （2）对于初始值，应用雅可比迭代公式、高斯－塞德尔迭代公式分别计算（保留小数点后五位数字）. 计算题2.答案 1.解原方程组同解变形为 雅可比迭代公式为 高斯－塞德尔迭代法公式 用雅可比迭代公式得 用高斯－塞德尔迭代公式得 3. 用牛顿法求方程在之间的近似根 （1）请指出为什么初值应取2？ （2）请用牛顿法求出近似根，精确到0.0001. 计算题3.答案
第一章 绪论 姓名 学号 班级 习题主要考察点：有效数字的计算、计算方法的比较选择、误差和误差限的计算。 1 若误差限为5105.0-?,那么近似数有几位有效数字（有效数字的计算） 解：2*103400.0-?=x ，325*102 1102 1---?=?≤-x x 故具有3位有效数字。 2 14159.3=π具有4位有效数字的近似值是多少（有效数字的计算） 解：10314159.0?= π，欲使其近似值*π具有4位有效数字，必需 41*102 1 -?≤-ππ，3*3102 1102 1--?+≤≤?-πππ，即14209.314109.3*≤≤π 3 已知2031.1=a ，978.0=b 是经过四舍五入后得到的近似值，问b a +， b a ?有几位有效数字（有效数字的计算） 解：3*1021 -?≤-a a ，2*102 1-?≤-b b ，而1811.2=+b a ，1766.1=?b a 2123****102 1 10211021)()(---?≤?+?≤ -+-≤+-+b b a a b a b a 故b a +至少具有2位有效数字。 2123*****102 1 0065.01022031.1102978.0)()(---?≤=?+?≤ -+-≤-b b a a a b b a ab
故b a ?至少具有2位有效数字。 4 设0>x ，x 的相对误差为δ，求x ln 的误差和相对误差（误差的计算） 解：已知δ=-* *x x x ，则误差为 δ=-= -* **ln ln x x x x x 则相对误差为 * * ** * * ln ln 1ln ln ln x x x x x x x x δ = -= - 5测得某圆柱体高度h 的值为cm h 20*=，底面半径r 的值为cm r 5*=， 已知cm h h 2.0||*≤-，cm r r 1.0||*≤-，求圆柱体体积h r v 2π=的绝对误差 限与相对误差限。（误差限的计算） 解：*2******2),(),(h h r r r h r r h v r h v -+-≤-ππ 绝对误差限为 πππ252.051.02052)5,20(),(2=??+????≤-v r h v 相对误差限为 %420 1 20525) 5,20() 5,20(),(2 ==??≤ -ππv v r h v 6 设x 的相对误差为%a ,求n x y =的相对误差。（函数误差的计算） 解：%* *a x x x =-， )%(* **** *na x x x n x x x y y y n n n =-≤-= - 7计算球的体积，为了使体积的相对误差限为%1，问度量半径r 时允许的相对误差限为多大（函数误差的计算）
例1、已知函数表 求() f x的Lagrange二次插值多项式和Newton二次插值多项式。 解： （1）由题可知 插值基函数分别为 故所求二次拉格朗日插值多项式为 （2）一阶均差、二阶均差分别为 均差表为
故所求Newton 二次插值多项式为 例2、 设2()32f x x x =++，[0,1]x ∈，试求()f x 在[0, 1]上关于()1x ρ=，{}span 1,x Φ=的 最佳平方逼近多项式。 解： 若{}span 1,x Φ=，则0()1x ?=，1()x x ?=，且()1x ρ=，这样，有 所以，法方程为 011231261192 34a a ??????????=?????????? ?????????? ，经过消元得012311 62110123a a ??? ???????=???????????????????? 再回代解该方程，得到14a =，011 6 a = 故，所求最佳平方逼近多项式为* 111 ()46 S x x = +
例3、 设()x f x e =，[0,1]x ∈，试求()f x 在[0, 1]上关于()1x ρ=，{}span 1,x Φ=的最佳平 方逼近多项式。 解： 若{}span 1,x Φ=，则0()1x ?=，1()x x ?=，这样，有 所以，法方程为 解法方程，得到00.8732a =，1 1.6902a =， 故，所求最佳平方逼近多项式为 例4、 用4n =的复合梯形和复合辛普森公式计算积分1?。 解： （1）用4n =的复合梯形公式 由于 2h =，()f x =，()121,2,3k x k k =+=，所以，有 （2）用4n =的复合辛普森公式 由于2h =，()f x =，()121,2,3k x k k =+=，()1 2 220,1,2,3k x k k +=+=，所以，有 例5、 用列主元消去法求解下列线性方程组的解。 解：先消元 再回代，得到33x =，22x =，11x =
“数值分析”课程 第一次小论文 郑维珍2015210459 制研15班（精密仪器系）内容：数值分析在你所在研究领域的应用。 要求：1）字数2500以上；2）要有摘要和参考文献；3）截至10.17，网络学堂提交，过期不能提交！ 数值分析在微流控芯片研究领域的应用 摘要： 作者在硕士期间即将参与的课题是微流控芯片的研制。当前，微流控芯片发展十分迅猛，而其中涉及到诸多材料学、电子学、光学、流体力学等领域的问题，加上微纳尺度上的尺寸效应，理论研究和数值计算都显得困难重重。发展该领域的数值计算，成为重中之重。本文从微流体力学、微传热学、微电磁学、微结构力学等分支入手，简要分析一下数值分析方法在该领域的应用。 微流控芯片(Microfluidic Chip)通常又称芯片实验室(Lab-On-a-Chip )，它是20世纪90年代初由瑞士的Manz和Widmer提出的[1-2]，它通过微细加工技术，将微管道、微泵、微阀、微电极、微检测元件等功能元件集成在芯片材料(基片)上，完成整个生化实验室的分析功能，具有减少样品的消耗量、节省反应和分析的时间、高通量和便携性等优点。 通常一个微流控芯片系统都会执行一个到多个微流体功能，如泵、混合、热循环、扩散和分离等，精确地操纵这些流体过程是微流控芯片的关键。因此它的研究不仅需要生命科学、MEMS、材料学、电子学、光学、流体力学等多学科领域的基础理论的支持，还需要很多数学计算。
1）微流体力学计算[3]： 对微管里的流体动力的研究主要包含了以下几个方面：(1)微管内流体的粘滞力的研究；(2)微管内气流液流的传热活动；(3)在绝热或传热的微管内两相流的流动和能量转换。这三方面的研究涵盖了在绝热、传热和多相转换条件下，可压缩和不可压缩流体在规则或不规则的微管内的流动特性研究。 由此，再结合不同的初值条件和边界条件，我们可以得到各种常微分方程或偏微分方程，而求解这些方程，就是需要很多数值分析的知识。例如，文献[4]里就针对特定的初值和边界条件，由软件求解了Navier-Stodes方程： 文献[4]专门有一章节讨论了该方程的离散化和数值求解。 微流体力学主要向两个方面发展：一方面是研究流动非定常稳定特性、分叉解及微尺寸效应下的湍流流动的机理，更为复杂的非定常、多尺度的流动特征，高精度、高分辨率的计算方法和并行算法；另一方面是将宏观流体力学的基本模型，结合微纳效应，直接用于模拟各种实际流动，解决微纳芯片生产制造中提出来的各种问题。 2）微传热方程计算： 常微分、偏微分方程的数值求解应用较为广泛的另一问题就是微流体传热问题。由传热学的相关知识，我们可以达到如下的传热学基本方程： 该方程在二维情况下经过简化和离散，可以得到如教材第三章所讲的“五点差分格式”的方程组，从而采取数值方法求解[5]。 除此之外，微结构芯片在加工和制造过程中也会有很多热学方面的问题，例如文献[6]所反映的注塑成型工艺中，就有大量的类似问题的解决。 3）微电磁学计算： 由于外加电场的作用，电渗流道中会产生焦耳热效应。许多研究者对电渗流道中的焦耳热效应进行了数值模拟研究。新加坡南洋理工大学的G. Y. Tang等在电渗流模型的基础上，考虑了与温度有关的物理系数，在固一液祸合区域内利用