西工大计算方法作业答案
西北工业大学数值分析(附答案)
西北工业大学数值分析习题集第一章 绪 论1. 设x >0,x 的相对误差为δ,求ln x 的误差.2. 设x 的相对误差为2％,求nx 的相对误差.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 -=( n=1,2,…)计算到100Y .27.982(五位有效数字),试问计算100Y 将有多大误差?7. 求方程25610x x -+=的两个根,使它至少具有四位有效数字27.982).8. 当N 充分大时,怎样求211N 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. 计算61)f =, 1.4≈,利用下列等式计算,哪一个得到的结果最好?3--13.()ln(f x x =,求f (30)的值.若开平方用六位函数表,问求对数时误差有多大?若改用另一等价公式ln(ln(x x =-计算,求对数时误差有多大?14. 试用消元法解方程组{101012121010;2.x x x x +=+=假定只用三位数计算,问结果是否可靠?15. 已知三角形面积1sin ,2s ab c =其中c 为弧度,02c π<<,且测量a ,b ,c 的误差分别为,,.a b c ∆∆∆证明面积的误差s ∆满足.s a b cs a b c ∆∆∆∆≤++第二章 插值法1. 根据(2.2)定义的范德蒙行列式,令2000011211121()(,,,,)11n n n n n n n n n x x x V x V x x x x x x x xx x ----==证明()n V x 是n 次多项式,它的根是01,,n x x -,且 101101()(,,,)()()n n n n V x V x x x x x x x ---=--.2. 当x = 1 , -1 , 2 时, f (x)= 0 , -3 , 4 ,求f (x )的二次插值多项式.3.4. 给出cos x ,0°≤x ≤90°的函数表,步长h =1′=(1/60)°,若函数表具有5位有效数字,研究用线性插值求cos x 近似值时的总误差界.5. 设0k x x kh =+,k =0,1,2,3,求032max ()x x x l x ≤≤.6. 设jx 为互异节点(j =0,1,…,n ),求证:i) 0()(0,1,,);nk k j j j x l x x k n =≡=∑ii) 0()()1,2,,).n k jj j xx l x k n =-≡0(=∑7. 设[]2(),f x C a b ∈且()()0f a f b ==,求证21()()().8max max a x ba xb f x b a f x ≤≤≤≤≤-"8. 在44x -≤≤上给出()xf x e =的等距节点函数表,若用二次插值求xe 的近似值,要使截断误差不超过610-,问使用函数表的步长h 应取多少?9. 若2n n y =,求4n y ∆及4n y δ.10. 如果()f x 是m 次多项式,记()()()f x f x h f x ∆=+-,证明()f x 的k 阶差分()(0)k f x k m ∆≤≤是m k -次多项式,并且()0(m l f x l +∆=为正整数).11. 证明1()k k k k k k f g f g g f +∆=∆+∆.12. 证明110010.n n kkn n k k k k f gf g f g g f --+==∆=--∆∑∑13. 证明1200.n jn j yy y -=∆=∆-∆∑14. 若1011()n n n n f x a a x a x a x --=++++有n 个不同实根12,,,n x x x ,证明{10,02;, 1.1()n k njk n a k n j jx f x -≤≤-=-=='∑15. 证明n 阶均差有下列性质: i)若()()F x cf x =,则[][]0101,,,,,,n n F x x x cf x x x =;ii) 若()()()F x f x g x =+,则[][][]010101,,,,,,,,,n n n F x x x f x x x g x x x =+.16. 74()31f x x x x =+++,求0172,2,,2f ⎡⎤⎣⎦及0182,2,,2f ⎡⎤⎣⎦.17. 证明两点三次埃尔米特插值余项是(4)22311()()()()/4!,(,)k k k k R x f x x x x x x ++=ξ--ξ∈并由此求出分段三次埃尔米特插值的误差限.18. 求一个次数不高于4次的多项式()P x ,使它满足(0)(1)P P k =-+并由此求出分段三次埃尔米特插值的误差限. 19. 试求出一个最高次数不高于4次的函数多项式()P x ,以便使它能够满足以下边界条件(0)(0)0P P ='=,(1)(1)1P P ='=,(2)1P =.20. 设[](),f x C a b ∈,把[],a b 分为n 等分,试构造一个台阶形的零次分段插值函数()n x ϕ并证明当n →∞时,()n x ϕ在[],a b 上一致收敛到()f x .21. 设2()1/(1)f x x =+,在55x -≤≤上取10n =,按等距节点求分段线性插值函数()h I x ,计算各节点间中点处的()h I x 与()f x 的值,并估计误差. 22. 求2()f x x =在[],a b 上的分段线性插值函数()h I x ,并估计误差. 23. 求4()f x x =在[],a b 上的分段埃尔米特插值,并估计误差.试求三次样条插值并满足条件 i) (0.25) 1.0000,(0.53)0.6868;S S '='= ii)(0.25)(0.53)0.S S "="=25. 若[]2(),f x C a b ∈,()S x 是三次样条函数,证明 i)[][][][]222()()()()2()()()bbbbaaaaf x dx S x dx f x S x dx S x f x S x dx"-"="-"+""-"⎰⎰⎰⎰;若()()(0,1,,)i i f x S x i n ==,式中i x 为插值节点,且01n a x x x b =<<<=,则[][][]()()()()()()()()()baS x f x S x dx S b f b S b S a f a S a ""-"="'-'-"'-'⎰.26. 编出计算三次样条函数()S x 系数及其在插值节点中点的值的程序框图(()S x 可用(8.7)式的表达式).第三章 函数逼近与计算1. (a)利用区间变换推出区间为[],a b 的伯恩斯坦多项式.(b)对()sin f x x =在[]0,/2π上求1次和三次伯恩斯坦多项式并画出图形,并与相应的马克劳林级数部分和误差做比较. 2. 求证:(a)当()m f x M ≤≤时,(,)n m B f x M ≤≤. (b)当()f x x =时,(,)n B f x x =.3. 在次数不超过6的多项式中,求()sin 4f x x =在[]0,2π的最佳一致逼近多项式.4. 假设()f x 在[],a b 上连续,求()f x 的零次最佳一致逼近多项式.5. 选取常数a ,使301max x x ax≤≤-达到极小,又问这个解是否唯一?6. 求()sin f x x =在[]0,/2π上的最佳一次逼近多项式,并估计误差.7. 求()xf x e =在[]0,1上的最佳一次逼近多项式.8. 如何选取r ,使2()p x x r =+在[]1,1-上与零偏差最小?r 是否唯一?9. 设43()31f x x x =+-,在[]0,1上求三次最佳逼近多项式.10. 令[]()(21),0,1n n T x T x x =-∈,求***0123(),(),(),()T x T x T x T x .11. 试证{}*()nT x 是在[]0,1上带权ρ=的正交多项式.12. 在[]1,1-上利用插值极小化求11()f x tg x -=的三次近似最佳逼近多项式. 13. 设()xf x e =在[]1,1-上的插值极小化近似最佳逼近多项式为()n L x ,若n f L ∞-有界,证明对任何1n ≥,存在常数n α、n β,使11()()()()(11).n n n n n T x f x L x T x x ++α≤-≤β-≤≤14. 设在[]1,1-上234511315165()128243843840x x x x x x ϕ=-----,试将()x ϕ降低到3次多项式并估计误差. 15. 在[]1,1-上利用幂级数项数求()sin f x x =的3次逼近多项式,使误差不超过0.005.16. ()f x 是[],a a -上的连续奇(偶)函数,证明不管n 是奇数或偶数,()f x 的最佳逼近多项式*()n n F x H ∈也是奇(偶)函数.17. 求a 、b 使[]220sin ax b x dx π+-⎰为最小.并与1题及6题的一次逼近多项式误差作比较.1g x C a b∈(),f x、[],定义18.()()(,)()();()(,)()()()();bbaaa f g f x g x dxb f g f x g x dx f a g a =''=''+⎰⎰问它们是否构成内积?19. 用许瓦兹不等式(4.5)估计6101x dx x +⎰的上界,并用积分中值定理估计同一积分的上下界,并比较其结果.20. 选择a ,使下列积分取得最小值:1122211(),x ax dx x ax dx----⎰⎰.21. 设空间{}{}10010121,,,span x span x x 1ϕ=ϕ=,分别在1ϕ、2ϕ上求出一个元素,使得其为[]20,1x C ∈的最佳平方逼近,并比较其结果.22. ()f x x =在[]1,1-上,求在{}2411,,span x x ϕ=上的最佳平方逼近.23.sin (1)arccos ()n n x u x +=是第二类切比雪夫多项式,证明它有递推关系()()()112n n n u x xu x u x +-=-.24. 将1()sin 2f x x=在[]1,1-上按勒让德多项式及切比雪夫多项式展开,求三次最佳平方逼近多项式并画出误差图形,再计算均方误差.25. 把()arccos f x x =在[]1,1-上展成切比雪夫级数.26.2y a bx =+.27.用最小二乘拟合求.29. 编出用正交多项式做最小二乘拟合的程序框图. 30. 编出改进FFT 算法的程序框图. 31. 现给出一张记录{}{}4,3,2,1,0,1,2,3k x =,试用改进FFT 算法求出序列{}k x 的离散频谱{}k C (0,1,,7).k =第四章 数值积分与数值微分1. 确定下列求积公式中的待定参数,使其代数精度尽量高,并指明所构造出的求积公式所具有的代数精度:(1)101()()(0)()hh f x dx A f h A f A f h --≈-++⎰; (2)21012()()(0)()hh f x dx A f h A f A f h --≈-++⎰;(3)[]1121()(1)2()3()/3f x dx f f x f x -≈-++⎰;(4)[][]20()(0)()/1(0)()hf x dx h f f h ah f f h ≈++'-'⎰.2. 分别用梯形公式和辛普森公式计算下列积分:(1)120,84xdx n x =+⎰; (2)1210(1),10x e dx n x --=⎰;(3)1,4n =⎰;(4),6n =.3. 直接验证柯特斯公式(2.4)具有5次代数精度.4. 用辛普森公式求积分10xedx-⎰并计算误差.5. 推导下列三种矩形求积公式:(1)2()()()()()2ba f f x dxb a f a b a 'η=-+-⎰; (2)2()()()()()2ba f f x dxb a f b b a 'η=---⎰;(3)3()()()()()224baa b f f x dx b a f b a +"η=-+-⎰. 6. 证明梯形公式(2.9)和辛普森公式(2.11)当n →∞时收敛到积分()baf x dx⎰.7. 用复化梯形公式求积分()baf x dx⎰,问要将积分区间[],a b 分成多少等分,才能保证误差不超过ε(设不计舍入误差)?8.1x e dx-,要求误差不超过510-.9. 卫星轨道是一个椭圆,椭圆周长的计算公式是S a =θ,这里a 是椭圆的半长轴,c 是地球中心与轨道中心(椭圆中心)的距离,记h 为近地点距离,H 为远地点距离,6371R =公里为地球半径,则(2)/2,()/2a R H h c H h =++=-.我国第一颗人造卫星近地点距离439h =公里,远地点距离2384H =公里,试求卫星轨道的周长. 10. 证明等式3524sin3!5!n nn n ππππ=-+-试依据sin(/)(3,6,12)n n n π=的值,用外推算法求π的近似值.11. 用下列方法计算积分31dyy ⎰并比较结果.(1) 龙贝格方法;(2) 三点及五点高斯公式;(3) 将积分区间分为四等分,用复化两点高斯公式.12. 用三点公式和五点公式分别求21()(1)f x x =+在x =1.0,1.1和1.2处的导数值,并估计误()f x第五章 常微分方程数值解法1. 就初值问题0)0(,=+='y b ax y 分别导出尤拉方法和改进的尤拉方法的近似解的表达式，并与准确解bx ax y +=221相比较。
西北工业大学计算方法试题
x ( k +1)
=
x(k)
−
ω
A(
x
(
k
+1
)
+ 2
x(k)
)
−
b
ω >0 , k = 0,1,2,⋯
对任意初始向量 x (0) ， x (k+1) 是否收敛到方程组 Ax = b 的解？为什么？
西北工业大学考试试题（卷）-计算方法二
1 填空 1）. 近似数 x* = 0.0142 关于真值 x = 0.0139 有__为有效数字。
0
试求满足插值条件的四次多项式 p(x).
6 设有如下的常微分方程初值问题
dy dx
=
x ,1 < y
x ≤ 1.4
y(1) = 1
1）写出每步用欧拉法预估，用梯形法进行一次校正的计算格式。 2）取步长 0.2 用上述格式求解。
∫ 7 设有积分 I = 0.6 e x2 dx 0
1）取 7 个等距节点（包括端点），列出被积函数在这些点出的值（保留到小数 点后 4 位） 2）用复化 simpson 公式求该积分的近似值。
（4） 取 3 ≈ 1.732 ，迭代过程 yn+1 = yn + 0.1 3 是否稳定？______（是或否）；
∫ （5） 求积公式 3 f ( x)dx ≈ 2 f (2) 有______次代数精度。 1
2．取初值 x0 = 1.6 ，用牛顿迭代法求 3.1 的近似值 xn+1 ，要求先论证收敛性，当
xn+1 − xn ≤ 10−5 时停止迭代。
3．用最小二乘法确定 y = a 1 + bx 2 中的常数 a 和 b ，使该函数曲线拟合 x
西工大2020年4月《计算机操作系统》作业机考参考答案
C.由于程序性错误引起的循环等待
D.由于进程占有资源后不被剥夺
正确答案:
12.（）是一旦有处理请求和要求处理的数据时，CPU就应该立即处理该数据并将结果及时送回。
A.实时操作系统
B.多道批处理系统
C.分时操作系统
D.微机操作系统
正确答案:
13.当用户程序执行访管指令时，中断装置将使中央处理器（）工作。
C.内存空间利用率得以改善
D.内存空间利用率不影响
正确答案:
8.进程和程序是两个既有联系又有区别的概念，下面描述中，（）是错误的。
A.进程是动态的，程序是静态的
B.一个程序可对应多个进程
C.进程有生命周期，而程序相对是永久的
D.程序是可以并发执行的
正确答案:
9.外存（如磁盘）上存放的程序和数据（）。
A.灵活性和可适应性
B.系统效率和吞吐量
C.交互性和响应时间
D.实时性和可靠性
正确答案:
41.进程从运行状态变为阻塞状态的原因是（）。
A.输入或输出事件发生
B.时间片到
C.输入或输出事件完成
D.某个进程被唤醒
正确答案:
42.主要由于（）原因，使UNIX易于移植。
A.UNIX是由机器指令编写的
B.UNIX大部分用汇编少部分用C语言编写
正确答案:C
5.操作系统中采用缓冲技术的目的是为了增强系统（）的能力。
A.串行操作
B.控制操作
C.重执操作
D.并行操作
正确答案:
6.进程间的基本关系为（）。
A.相互独立与相互制约
B.同步与互斥
C.并行执行与资源共享
D.信息传递与信息缓冲
计算方法 课后习题答案
得到方程组
3。举例说明一个非奇异矩阵不一定存在LU分解。
例如：设
与题设相矛盾，所以一个非奇异矩阵不一定存在LU分解。
4。下列矩阵能否分解为LU（其中L为单位下三角矩阵，U为上三角矩阵）？若能分解，那么分解是否唯一？
解：
设 B可以进行LU分解，则B=
计算得
其中。 。
解：(1)由题意，可设 ，由Lagrange插值余项公式得
(2)由(1)式可知，
15．给定数据表：
1
0
2
3
构造出函数 的差商表，并写出它的三次 插值多项式.
解：利用Newton插值公式：
先作出差商表
一阶差商
二阶差商
三阶差商
0
1
3
1
3/2
13/4
1/2
2
0
3
1/6
1/3
3
2
5/3
-2/3
-5/3
证明：据题4可知，
令 ，则有 。注意到
（证明见王能超数值简明教程145页题6）
令 即有 。
9.已知 ，求差商 和 。
解：根据差商与微商的关系，有
10.已知 互异，求 。其中 。（此题有误。）（见王能超《教程》P149－题2）
解：因为 ，则
由差商性质 可知，
11.设首项系数为1的n次式 有n个互异的零点 ，证明
解：1）用梯形公式有：
事实上，
2）Simpson公式
事实上，
3）由Cotes公式有：
事实上，
2．证明Simpson公式 具有三次代数精度。
证明：
而当 时
左侧：
右侧：
计算方法的课后答案
《计算方法》习题答案 第一章 数值计算中的误差 1．什么是计算方法？（狭义解释）答：计算方法就是将所求的的数学问题简化为一系列的算术运算和逻辑运算，以便在计算机上编程上机，求出问题的数值解，并对算法的收敛性、稳定性和误差进行分析、计算。
2．一个实际问题利用计算机解决所采取的五个步骤是什么？答：一个实际问题当利用计算机来解决时，应采取以下五个步骤：实际问题→建立数学模型→构造数值算法→编程上机→获得近似结果4．利用秦九韶算法计算多项式4)(53-+-=x x x x P 在3-=x 处的值，并编程获得解。
解：400)(2345-+⋅+-⋅+=x x x x x x P ，从而 1 0 -1 0 1 -4 -3 -3 9 -24 72 -219 1 -3 8 -24 73 -223所以，多项式4)(53-+-=x x x x P 在3-=x 处的值223)3(-=-P 。
5．叙述误差的种类及来源。
答：误差的种类及来源有如下四个方面： （1）模型误差：数学模型是对实际问题进行抽象，忽略一些次要因素简化得到的，它是原始问题的近似，即使数学模型能求出准确解，也与实际问题的真解不同，我们把数学模型与实际问题之间存在的误差称为模型误差。
（2）观测误差：在建模和具体运算过程中所用的一些原始数据往往都是通过观测、实验得来的，由于仪器的精密性，实验手段的局限性，周围环境的变化以及人们的工作态度和能力等因素，而使数据必然带有误差，这种误差称为观测误差。
（3）截断误差：理论上的精确值往往要求用无限次的运算才能得到，而实际运算时只能用有限次运算的结果来近似，这样引起的误差称为截断误差（或方法误差）。
（4）舍入误差：在数值计算过程中还会用到一些无穷小数，而计算机受机器字长的限制，它所能表示的数据只能是一定的有限数位，需要把数据按四舍五入成一定位数的近似的有理数来代替。
这样引起的误差称为舍入误差。
6．掌握绝对误差（限）和相对误差（限）的定义公式。
西工大计算方法试题参考(完整版)
2002-2003第一学期一．计算及推导（5*8）1．已知* 3.141,x x π==，试确定*x 近似x 的有效数字位数。
2．有效数***1233.105,0.001,0.100x x x =-==，试确定***123x x x ++的相对误差限。
3．已知3()0.50.12f x x x =++，试计算差商[]0,1,2,3f 4．给出拟合三点(0,1),(1,0)A B ==和(1,1)C =的直线方程。
5．推导中矩形求积公式''31()()()()()224b aa b f x dx b a f f b a η+=-+-⎰ 6．试证明插值型求积公式()()nbi i ai f x dx A f x =≈∑⎰的代数精确度至少是n 次。
7．已知非线性方程()x f x =在区间[],a b内有一实根，试写出该实根的牛顿迭代公式。
8．用三角分解法求解线性方程组123121022331302x x x ⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥--⎣⎦⎣⎦⎣⎦要用二次插值多项式计算(0.63891)f 的近似值，试选择合适的插值节点进行计算，并说明所选用节点依据。
（保留5位有效数字）（12分） 三． 已知方程ln 0x x +=在(0,1)内有一实根α（1）给出求该实根的一个迭代公式，试之对任意的初始近似0(0,1)x ∈迭代法都收敛，并证明其收敛性。
（2）00.5x =试用构造的迭代公式计算α的近似值n x ，要求3110n n x x ---≤。
四． 设有方程组112233131232a x b a x b a x b ⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦⎣⎦当参数a 满足什么条件时，雅可比方法对任意的初始向量都收敛。
写出与雅可比方法对应的高斯赛德尔迭代公式。
（12分） 五．用欧拉预估校正法求解初值问题 '2 (00.2)(0)1x y y x y y ⎧=-≤≤⎪⎨⎪=⎩ 取h=0.1，小数点后保留5位。
西北工业大学数值分析(附答案)
西北工业大学数值分析习题集第一章 绪 论1. 设x >0,x 的相对误差为δ,求ln x 的误差.2. 设x 的相对误差为2％,求nx 的相对误差.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 -= ( n=1,2,…)计算到100Y .27.982(五位有效数字),试问计算100Y 将有多大误差?7. 求方程25610x x -+=的两个根,使它至少具有四位有效数字≈27.982).8. 当N 充分大时,怎样求211Ndx 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. 计算61)f =, 1.4≈,利用下列等式计算,哪一个得到的结果最好?3--13.()ln(f x x =,求f (30)的值.若开平方用六位函数表,问求对数时误差有多大?若改用另一等价公式ln(ln(x x =-计算,求对数时误差有多大?14. 试用消元法解方程组{101012121010;2.x x x x +=+=假定只用三位数计算,问结果是否可靠?15. 已知三角形面积1sin ,2s ab c =其中c 为弧度,02c π<<,且测量a ,b ,c 的误差分别为,,.a b c ∆∆∆证明面积的误差s ∆满足.s a b c s a b c ∆∆∆∆≤++第二章 插值法1. 根据(2.2)定义的范德蒙行列式,令2000011211121()(,,,,)11n n n n n n n n n x x x V x V x x x x x x x xxx ----==证明()n V x 是n 次多项式,它的根是01,,n x x -,且 101101()(,,,)()()n n n n V x V x x x x x x x ---=--.2. 当x = 1 , -1 , 2 时, f (x)= 0 , -3 , 4 ,求f (x )的二次插值多项式.3.4. 给出cos x ,0°≤x ≤90°的函数表,步长h =1′=(1/60)°,若函数表具有5位有效数字,研究用线性插值求cos x 近似值时的总误差界.5. 设0k x x kh =+,k =0,1,2,3,求032max ()x x x l x ≤≤.6. 设jx 为互异节点(j =0,1,…,n ),求证:i) 0()(0,1,,);nk kj jj x l x x k n =≡=∑ii) 0()()1,2,,).nk jj j xx l x k n =-≡0(=∑1. 设[]2(),f x C a b ∈且()()0f a f b ==,求证21()()().8max max a x ba xb f x b a f x ≤≤≤≤≤-"2. 在44x -≤≤上给出()xf x e =的等距节点函数表,若用二次插值求xe 的近似值,要使截断误差不超过610-,问使用函数表的步长h 应取多少? 3. 若2nn y =,求4n y ∆及4n y δ. 4. 如果()f x 是m 次多项式,记()()()f x f x h f x ∆=+-,证明()f x 的k 阶差分()(0)k f x k m ∆≤≤是m k -次多项式,并且()0(m l f x l +∆=为正整数).5. 证明1()k k k k k k f g f g g f +∆=∆+∆.6. 证明110010.n n kkn n k k k k f gf g f g g f --+==∆=--∆∑∑7. 证明1200.n j n j y y y -=∆=∆-∆∑8. 若1011()n n n n f x a a x a x a x --=++++有n 个不同实根12,,,n x x x ,证明{10,02;, 1.1()n k njk n a k n j jx f x -≤≤-=-=='∑9. 证明n 阶均差有下列性质: i)若()()F x cf x =,则[][]0101,,,,,,n n F x x x cf x x x =;ii) 若()()()F x f x g x =+,则[][][]010101,,,,,,,,,n n n F x x x f x x x g x x x =+.10. 74()31f x x x x =+++,求0172,2,,2f ⎡⎤⎣⎦及0182,2,,2f ⎡⎤⎣⎦.11. 证明两点三次埃尔米特插值余项是(4)22311()()()()/4!,(,)k k k k R x f x x x x x x ++=ξ--ξ∈并由此求出分段三次埃尔米特插值的误差限.12. 求一个次数不高于4次的多项式()P x ,使它满足(0)(1)P P k =-+并由此求出分段三次埃尔米特插值的误差限. 13. 试求出一个最高次数不高于4次的函数多项式()P x ,以便使它能够满足以下边界条件(0)(0)0P P ='=,(1)(1)1P P ='=,(2)1P =.14. 设[](),f x C a b ∈,把[],a b 分为n 等分,试构造一个台阶形的零次分段插值函数()n x ϕ并证明当n →∞时,()n x ϕ在[],a b 上一致收敛到()f x .15. 设2()1/(1)f x x =+,在55x -≤≤上取10n =,按等距节点求分段线性插值函数()h I x ,计算各节点间中点处的()h I x 与()f x 的值,并估计误差. 16. 求2()f x x =在[],a b 上的分段线性插值函数()h I x ,并估计误差. 17. 求4()f x x =在[],a b 上的分段埃尔米特插值,并估计误差.试求三次样条插值并满足条件i) (0.25) 1.0000,(0.53)0.6868;S S '='= ii)(0.25)(0.53)0.S S "="=19. 若[]2(),f x C a b ∈,()S x 是三次样条函数,证明 i)[][][][]222()()()()2()()()bbbbaaaaf x dx S x dx f x S x dx S x f x S x dx"-"="-"+""-"⎰⎰⎰⎰;ii) 若()()(0,1,,)i i f x S x i n ==,式中i x 为插值节点,且01n a x x x b =<<<=,则[][][]()()()()()()()()()baS x f x S x dx S b f b S b S a f a S a ""-"="'-'-"'-'⎰.20. 编出计算三次样条函数()S x 系数及其在插值节点中点的值的程序框图(()S x 可用(8.7)式的表达式).第三章 函数逼近与计算1. (a)利用区间变换推出区间为[],a b 的伯恩斯坦多项式.(b)对()sin f x x =在[]0,/2π上求1次和三次伯恩斯坦多项式并画出图形,并与相应的马克劳林级数部分和误差做比较. 2. 求证:(a)当()m f x M ≤≤时,(,)n m B f x M ≤≤. (b)当()f x x =时,(,)n B f x x =.3. 在次数不超过6的多项式中,求()sin 4f x x =在[]0,2π的最佳一致逼近多项式.4. 假设()f x 在[],a b 上连续,求()f x 的零次最佳一致逼近多项式.5. 选取常数a ,使301max x x ax≤≤-达到极小,又问这个解是否唯一?6. 求()sin f x x =在[]0,/2π上的最佳一次逼近多项式,并估计误差.7. 求()x f x e =在[]0,1上的最佳一次逼近多项式.8. 如何选取r ,使2()p x x r =+在[]1,1-上与零偏差最小?r 是否唯一? 9. 设43()31f x x x =+-,在[]0,1上求三次最佳逼近多项式.10. 令[]()(21),0,1n n T x T x x =-∈,求***0123(),(),(),()T x T x T x T x .11. 试证{}*()nTx 是在[]0,1上带权ρ=的正交多项式.12. 在[]1,1-上利用插值极小化求11()f x tg x -=的三次近似最佳逼近多项式. 13. 设()xf x e =在[]1,1-上的插值极小化近似最佳逼近多项式为()n L x ,若n f L ∞-有界,证明对任何1n ≥,存在常数n α、n β,使11()()()()(11).n n n n n T x f x L x T x x ++α≤-≤β-≤≤14. 设在[]1,1-上234511315165()128243843840x x x x x x ϕ=-----,试将()x ϕ降低到3次多项式并估计误差.15. 在[]1,1-上利用幂级数项数求()sin f x x =的3次逼近多项式,使误差不超过0.005.16. ()f x 是[],a a -上的连续奇(偶)函数,证明不管n 是奇数或偶数,()f x 的最佳逼近多项式*()n n F x H ∈也是奇(偶)函数.17. 求a 、b 使[]220sin ax b x dx π+-⎰为最小.并与1题及6题的一次逼近多项式误差作比较.18. ()f x 、[]1(),g x C a b ∈,定义()(,)()();()(,)()()()();b baaa f g f x g x dxb f g f x g x dx f a g a =''=''+⎰⎰问它们是否构成内积?19. 用许瓦兹不等式(4.5)估计6101x dx x +⎰的上界,并用积分中值定理估计同一积分的上下界,并比较其结果.20. 选择a ,使下列积分取得最小值:1122211(),x ax dx x ax dx----⎰⎰.21. 设空间{}{}10010121,,,span x span x x 1ϕ=ϕ=,分别在1ϕ、2ϕ上求出一个元素,使得其为[]20,1x C ∈的最佳平方逼近,并比较其结果.22.()f x x =在[]1,1-上,求在{}2411,,span x x ϕ=上的最佳平方逼近.23.sin (1)arccos ()n n x u x +=是第二类切比雪夫多项式,证明它有递推关系()()()112n n n u x xu x u x +-=-.24. 将1()sin2f x x=在[]1,1-上按勒让德多项式及切比雪夫多项式展开,求三次最佳平方逼近多项式并画出误差图形,再计算均方误差.25. 把()arccos f x x =在[]1,1-上展成切比雪夫级数.26.2y a bx =+.27.用最小二乘拟合求.29. 编出用正交多项式做最小二乘拟合的程序框图. 30. 编出改进FFT 算法的程序框图. 31. 现给出一张记录{}{}4,3,2,1,0,1,2,3k x =,试用改进FFT 算法求出序列{}k x 的离散频谱{}k C (0,1,,7).k =第四章 数值积分与数值微分1. 确定下列求积公式中的待定参数,使其代数精度尽量高,并指明所构造出的求积公式所具有的代数精度:(1)101()()(0)()hh f x dx A f h A f A f h --≈-++⎰; (2)21012()()(0)()hh f x dx A f h A f A f h --≈-++⎰;(3)[]1121()(1)2()3()/3f x dx f f x f x -≈-++⎰;(4)[][]20()(0)()/1(0)()hf x dx h f f h ah f f h ≈++'-'⎰.2. 分别用梯形公式和辛普森公式计算下列积分:(1)120,84xdx n x =+⎰; (2)1210(1),10x e dx n x --=⎰;(3)1,4n =⎰;(4),6n =.3. 直接验证柯特斯公式(2.4)具有5次代数精度.4. 用辛普森公式求积分10xedx-⎰并计算误差.5. 推导下列三种矩形求积公式:(1)2()()()()()2ba f f x dxb a f a b a 'η=-+-⎰; (2)2()()()()()2ba f f x dxb a f b b a 'η=---⎰;(3)3()()()()()224baa b f f x dx b a f b a +"η=-+-⎰. 6. 证明梯形公式(2.9)和辛普森公式(2.11)当n →∞时收敛到积分()baf x dx⎰.7. 用复化梯形公式求积分()baf x dx⎰,问要将积分区间[],a b 分成多少等分,才能保证误差不超过ε(设不计舍入误差)?8.1x e dx-,要求误差不超过510-.9. 卫星轨道是一个椭圆,椭圆周长的计算公式是S a =θ,这里a 是椭圆的半长轴,c 是地球中心与轨道中心(椭圆中心)的距离,记h 为近地点距离,H 为远地点距离,6371R =公里为地球半径,则(2)/2,()/2a R H h c H h =++=-.我国第一颗人造卫星近地点距离439h =公里,远地点距离2384H =公里,试求卫星轨道的周长. 10. 证明等式3524sin3!5!n nn n ππππ=-+-试依据sin(/)(3,6,12)n n n π=的值,用外推算法求π的近似值.11. 用下列方法计算积分31dyy ⎰并比较结果.(1) 龙贝格方法;(2) 三点及五点高斯公式;(3) 将积分区间分为四等分,用复化两点高斯公式.1. 用三点公式和五点公式分别求21()(1)f x x =+在x =1.0,1.1和1.2处的导数值,并估计误()f x第五章 常微分方程数值解法1. 就初值问题0)0(,=+='y b ax y 分别导出尤拉方法和改进的尤拉方法的近似解的表达式，并与准确解bx ax y +=221相比较。
西工大20年10月机考《计算方法》作业参考答案非免费
西工大20年10月机考计算方法作业试卷总分:100 得分:96要答an：网叫福到（这四个字的拼音）一、单选题 (共 30 道试题,共 60 分)1.舍入误差是( )产生的误差。
A.只取有限位数B.模型准确值与用数值方法求得的准确值C.观察与测量D.数学模型准确值与实际值正确答案:2. {A.2B.3C.4D.5正确答案:3.用 1+x近似表示ex所产生的误差是( )误差。
A.模型B.观测C.截断D.舍入正确答案:4.解线性方程组的主元素消去法中选择主元的目的是( )。
A.控制舍入误差B.减小方法误差C.防止计算时溢出D.简化计算正确答案:5.舍入误差是(?? ?)产生的误差。
A.只取有限位数B.模型准确值与用数值方法求得的准确值C.观察与测量D.数学模型准确值与实际值正确答案:6. {A.{<img ">B.{<img g">C.0D.1正确答案:7.( )是解方程组Ax=b的迭代格式x(k+1)=Mx(k)+f收敛的一个充分条件；A.{<img ">B.{<img ">C.{<img ">D.{<img >正确答案:8.-324.7500是舍入得到的近似值，它有( )位有效数字。
A.5B.6C.7D.8正确答案:9. {A.舍入B.观测C.模型D.截断正确答案:10. {A.-1B.1C.{<img ">D.0正确答案:11. {A.{<img ">B.{<img >C.{<img >D.0正确答案:12. {A.1B.2C.4D.3正确答案:13. {A.A的各阶顺序主子式不为零B.{<img ">C.{<img ">D.{<img pg">正确答案:14. {A.0B.1C.2D.{<img ">正确答案:15. {A.0B.{<img ">C.2D.1正确答案:16. {A.0B.1C.{<img s>D.{<img s>正确答案:17. 三点的高斯型求积公式的代数精度为()。
(完整word版)西工大计算方法试题参考(完整版)
2002-2003第一学期一．计算及推导（5*8）1．已知* 3.141,x x π==，试确定*x 近似x 的有效数字位数。
2．有效数***1233.105,0.001,0.100x x x =-==，试确定***123x x x ++的相对误差限。
3．已知3()0.50.12f x x x =++，试计算差商[]0,1,2,3f 4．给出拟合三点(0,1),(1,0)A B ==和(1,1)C =的直线方程。
5．推导中矩形求积公式''31()()()()()224b aa b f x dx b a f f b a η+=-+-⎰ 6．试证明插值型求积公式()()nbi i ai f x dx A f x =≈∑⎰的代数精确度至少是n 次。
7．已知非线性方程()x f x =在区间[],a b内有一实根，试写出该实根的牛顿迭代公式。
8．用三角分解法求解线性方程组123121022331302x x x ⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥--⎣⎦⎣⎦⎣⎦要用二次插值多项式计算(0.63891)f 的近似值，试选择合适的插值节点进行计算，并说明所选用节点依据。
（保留5位有效数字）（12分） 三． 已知方程ln 0x x +=在(0,1)内有一实根α（1）给出求该实根的一个迭代公式，试之对任意的初始近似0(0,1)x ∈迭代法都收敛，并证明其收敛性。
（2）00.5x =试用构造的迭代公式计算α的近似值n x ，要求3110n n x x ---≤。
四． 设有方程组112233131232a x b a x b a x b ⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦⎣⎦当参数a 满足什么条件时，雅可比方法对任意的初始向量都收敛。
写出与雅可比方法对应的高斯赛德尔迭代公式。
（12分） 五．用欧拉预估校正法求解初值问题 '2 (00.2)(0)1x y y x y y ⎧=-≤≤⎪⎨⎪=⎩ 取h=0.1，小数点后保留5位。
西北工业大学数值分析(附答案)
西北工业大学数值分析习题集第一章 绪 论1. 设x >0,x 的相对误差为δ,求ln x 的误差.2. 设x 的相对误差为2％,求nx 的相对误差.3.下列各数都是经过四舍五入得到的近似数,即误差限不超过最后一位的半个单位,试指出它们是几位有效数字:*****123451.1021,0.031,385.6,56.430,7 1.0.x x x x x =====⨯4.利用公式求下列各近似值的误差限:********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 .若取(五位有效数字),试问计算100Y 将有多大误差7.求方程25610x x -+=的两个根,使它至少具有四位有效数字(≈.8.当N 充分大时,怎样求211Ndx x +∞+⎰9.正方形的边长大约为100㎝,应怎样测量才能使其面积误差不超过1㎝210. 设212S gt =假定g 是准确的,而对t 的测量有±秒的误差,证明当t增加时S 的绝对误差增加,而相对误差却减小.11. 序列{}n y 满足递推关系1101n n y y -=-(n=1,2,…),若0 1.41y ≈(三位有效数字),计算到10y 时误差有多大这个计算过程稳定吗12.计算61)f =,1.4≈,利用下列等式计算,哪一个得到的结果最好3--13.()ln(f x x =,求f (30)的值.若开平方用六位函数表,问求对数时误差有多大若改用另一等价公式ln(ln(x x =-计算,求对数时误差有多大14. 试用消元法解方程组{101012121010;2.x x x x +=+=假定只用三位数计算,问结果是否可靠15. 已知三角形面积1sin ,2s ab c =其中c 为弧度,02c π<<,且测量a ,b ,c的误差分别为,,.a b c ∆∆∆证明面积的误差s ∆满足.s a b cs a b c ∆∆∆∆≤++第二章 插值法1.根据定义的范德蒙行列式,令2000011211121()(,,,,)11n n n n n n n n n x x x V x V x x x x x x x xx x ----==证明()n V x 是n 次多项式,它的根是01,,n x x -,且101101()(,,,)()()n n n n V x V x x x x x x x ---=--.2.当x = 1 , -1 , 2 时, f (x)= 0 , -3 , 4 ,求f (x )的二次插值多项式.3.给出f (x )=ln x 的数值表用线性插值及二次插值计算ln 的近似值.4.给出cos x ,0°≤x ≤90°的函数表,步长h =1′=(1/60)°,若函数表具有5位有效数字,研究用线性插值求cos x 近似值时的总误差界.5. 设0k x x kh =+,k =0,1,2,3,求032max ()x x x l x ≤≤.6.设j x 为互异节点(j =0,1,…,n ),求证:i) 0()(0,1,,);nkkj jj x l x x k n =≡=∑ii) 0()()1,2,,).nk jj j xx l x k n =-≡0(=∑7.设[]2(),f x C a b ∈且()()0f a f b ==,求证21()()().8maxmax a x ba xb f x b a f x ≤≤≤≤≤-"8.在44x -≤≤上给出()xf x e =的等距节点函数表,若用二次插值求xe的近似值,要使截断误差不超过610-,问使用函数表的步长h 应取多少9. 若2n n y =,求4n y ∆及4n y δ.10. 如果()f x 是m 次多项式,记()()()f x f x h f x ∆=+-,证明()f x 的k 阶差分()(0)k f x k m ∆≤≤是m k -次多项式,并且()0(m lf x l +∆=为正整数).11. 证明1()k k k k k k f g f g g f +∆=∆+∆.12. 证明110010.n n kkn n k k k k f gf g f g g f --+==∆=--∆∑∑13. 证明1200.n j n j y y y -=∆=∆-∆∑14. 若1011()n n n n f x a a x a x a x --=++++有n 个不同实根12,,,n x x x ,证明{10,02;, 1.1()n k njk n a k n j jx f x -≤≤-=-=='∑15. 证明n 阶均差有下列性质:i)若()()F x cf x =,则[][]0101,,,,,,n n F x x x cf x x x =;ii) 若()()()F x f x g x =+,则[][][]010101,,,,,,,,,n n n F x x x f x x x g x x x =+.16. 74()31f x x x x =+++,求0172,2,,2f ⎡⎤⎣⎦及0182,2,,2f ⎡⎤⎣⎦.17. 证明两点三次埃尔米特插值余项是(4)22311()()()()/4!,(,)k k k k R x f x x x x x x ++=ξ--ξ∈并由此求出分段三次埃尔米特插值的误差限.18. 求一个次数不高于4次的多项式()P x ,使它满足(0)(1)P P k =-+并由此求出分段三次埃尔米特插值的误差限.19. 试求出一个最高次数不高于4次的函数多项式()P x ,以便使它能够满足以下边界条件(0)(0)0P P ='=,(1)(1)1P P ='=,(2)1P =.20. 设[](),f x C a b ∈,把[],a b 分为n 等分,试构造一个台阶形的零次分段插值函数()n x ϕ并证明当n →∞时,()n x ϕ在[],a b 上一致收敛到()f x .21. 设2()1/(1)f x x =+,在55x -≤≤上取10n =,按等距节点求分段线性插值函数()h I x ,计算各节点间中点处的()h I x 与()f x 的值,并估计误差.22. 求2()f x x=在[],a b 上的分段线性插值函数()h I x ,并估计误差.23. 求4()f x x =在[],a b 上的分段埃尔米特插值,并估计误差.24. 给定数据表如下:试求三次样条插值()S x 并满足条件i) (0.25) 1.0000,(0.53)0.6868;S S '='= ii)(0.25)(0.53)0.S S "="=25. 若[]2(),f x C a b ∈,()S x 是三次样条函数,证明 i)[][][][]222()()()()2()()()bbbbaaaaf x dx S x dx f x S x dx S x f x S x dx"-"="-"+""-"⎰⎰⎰⎰;ii) 若()()(0,1,,)i i f x S x i n ==,式中ix 为插值节点,且01n a x x x b=<<<=,则[][][]()()()()()()()()()baS x f x S x dx S b f b S b S a f a S a ""-"="'-'-"'-'⎰.26. 编出计算三次样条函数()S x 系数及其在插值节点中点的值的程序框图(()S x 可用式的表达式). 第三章 函数逼近与计算1.(a)利用区间变换推出区间为[],a b 的伯恩斯坦多项式.(b)对()sin f x x =在[]0,/2π上求1次和三次伯恩斯坦多项式并画出图形,并与相应的马克劳林级数部分和误差做比较.2.求证:(a)当()m f x M ≤≤时,(,)n m B f x M ≤≤. (b)当()f x x =时,(,)n B f x x =.3.在次数不超过6的多项式中,求()sin 4f x x =在[]0,2π的最佳一致逼近多项式.4.假设()f x 在[],a b 上连续,求()f x 的零次最佳一致逼近多项式. 5.选取常数a ,使301max x x ax≤≤-达到极小,又问这个解是否唯一6.求()sin f x x =在[]0,/2π上的最佳一次逼近多项式,并估计误差.7. 求()xf x e =在[]0,1上的最佳一次逼近多项式.8.如何选取r ,使2()p x x r =+在[]1,1-上与零偏差最小r 是否唯一9. 设43()31f x x x =+-,在[]0,1上求三次最佳逼近多项式. 10. 令[]()(21),0,1n n T x T x x =-∈,求***0123(),(),(),()T x T x T x T x .11. 试证{}*()nT x 是在[]0,1上带权ρ=.12. 在[]1,1-上利用插值极小化求11()f x tg x -=的三次近似最佳逼近多项式.13. 设()xf x e=在[]1,1-上的插值极小化近似最佳逼近多项式为()n L x ,若nf L ∞-有界,证明对任何1n ≥,存在常数n α、n β,使11()()()()(11).n n n n n T x f x L x T x x ++α≤-≤β-≤≤14. 设在[]1,1-上234511315165()128243843840x x x x x x ϕ=-----,试将()x ϕ降低到3次多项式并估计误差.15. 在[]1,1-上利用幂级数项数求()sin f x x =的3次逼近多项式,使误差不超过.16. ()f x 是[],a a -上的连续奇(偶)函数,证明不管n 是奇数或偶数,()f x 的最佳逼近多项式*()n n F x H ∈也是奇(偶)函数.17. 求a 、b 使[]22sin ax b x dx π+-⎰为最小.并与1题及6题的一次逼近多项式误差作比较.18. ()f x 、[]1(),g x C a b ∈,定义 ()(,)()();()(,)()()()();b baaa f g f x g x dxb f g f x g x dx f a g a =''=''+⎰⎰问它们是否构成内积19. 用许瓦兹不等式估计6101x dx x +⎰的上界,并用积分中值定理估计同一积分的上下界,并比较其结果.20. 选择a ,使下列积分取得最小值:1122211(),x ax dx x ax dx----⎰⎰.21. 设空间{}{}10010121,,,span x span x x 1ϕ=ϕ=,分别在1ϕ、2ϕ上求出一个元素,使得其为[]20,1x C ∈的最佳平方逼近,并比较其结果.22. ()f x x =在[]1,1-上,求在{}2411,,span x x ϕ=上的最佳平方逼近.23.sin (1)arccos ()n n x u x +=是第二类切比雪夫多项式,证明它有递推关系()()()112n n n u x xu x u x +-=-.24. 将1()sin 2f x x=在[]1,1-上按勒让德多项式及切比雪夫多项式展开,求三次最佳平方逼近多项式并画出误差图形,再计算均方误差.25. 把()arccos f x x =在[]1,1-上展成切比雪夫级数.26. 用最小二乘法求一个形如2y a bx=+的经验公式,使它与下列数据拟合,并求均方误差.27. 观测物体的直线运动,得出以下数据:求运动方程.28. 在某化学反应里,根据实验所得分解物的浓度与时间关系如下:用最小二乘拟合求()y f t =.29. 编出用正交多项式做最小二乘拟合的程序框图. 30. 编出改进FFT 算法的程序框图.31. 现给出一张记录{}{}4,3,2,1,0,1,2,3k x =,试用改进FFT 算法求出序列{}k x 的离散频谱{}k C (0,1,,7).k =第四章 数值积分与数值微分1.确定下列求积公式中的待定参数,使其代数精度尽量高,并指明所构造出的求积公式所具有的代数精度: (1)101()()(0)()hh f x dx A f h A f A f h --≈-++⎰; (2)21012()()(0)()hh f x dx A f h A f A f h --≈-++⎰;(3)[]1121()(1)2()3()/3f x dx f f x f x -≈-++⎰;(4)[][]20()(0)()/1(0)()h f x dx h f f h ah f f h ≈++'-'⎰.2.分别用梯形公式和辛普森公式计算下列积分:(1)120,84xdx n x =+⎰; (2)1210(1),10x e dx n x --=⎰;(3)1,4n =⎰; (4),6n =.3.直接验证柯特斯公式具有5次代数精度. 4. 用辛普森公式求积分1x e dx-⎰并计算误差.5.推导下列三种矩形求积公式: (1)2()()()()()2ba f f x dxb a f a b a 'η=-+-⎰; (2)2()()()()()2ba f f x dxb a f b b a 'η=---⎰;(3)3()()()()()224baa b f f x dx b a f b a +"η=-+-⎰.6.证明梯形公式和辛普森公式当n →∞时收敛到积分()baf x dx⎰.7.用复化梯形公式求积分()ba f x dx⎰,问要将积分区间[],a b 分成多少等分,才能保证误差不超过ε(设不计舍入误差)8. 1xedx-⎰,要求误差不超过510-.9.卫星轨道是一个椭圆,椭圆周长的计算公式是S a =θ,这里a 是椭圆的半长轴,c 是地球中心与轨道中心(椭圆中心)的距离,记h 为近地点距离,H 为远地点距离,6371R =公里为地球半径,则(2)/2,()/2a R H h c H h =++=-.我国第一颗人造卫星近地点距离439h =公里,远地点距离2384H =公里,试求卫星轨道的周长.10. 证明等式3524sin3!5!n nnnππππ=-+-试依据sin(/)(3,6,12)n n n π=的值,用外推算法求π的近似值.11. 用下列方法计算积分31dyy ⎰并比较结果.(1) 龙贝格方法;(2) 三点及五点高斯公式;(3) 将积分区间分为四等分,用复化两点高斯公式.12. 用三点公式和五点公式分别求21()(1)f x x =+在x =,和处的导数值,并估计误差.()f x 的值由下表给出:第五章 常微分方程数值解法1. 就初值问题0)0(,=+='y b ax y 分别导出尤拉方法和改进的尤拉方法的近似解的表达式，并与准确解bx ax y +=221相比较。
西北工业大学数值分析(附答案)
西北工业大学数值分析习题集第一章 绪 论1. 设x >0,x 的相对误差为δ,求ln x 的误差.2. 设x 的相对误差为2％,求nx 的相对误差.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 -=( n=1,2,…)计算到100Y .27.982(五位有效数字),试问计算100Y 将有多大误差?7. 求方程25610x x -+=的两个根,使它至少具有四位有效数字27.982).8. 当N 充分大时,怎样求211N 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. 计算61)f =, 1.4≈,利用下列等式计算,哪一个得到的结果最好?3--13.()ln(f x x =,求f (30)的值.若开平方用六位函数表,问求对数时误差有多大?若改用另一等价公式ln(ln(x x =-计算,求对数时误差有多大?14. 试用消元法解方程组{101012121010;2.x x x x +=+=假定只用三位数计算,问结果是否可靠?15. 已知三角形面积1sin ,2s ab c =其中c 为弧度,02c π<<,且测量a ,b ,c 的误差分别为,,.a b c ∆∆∆证明面积的误差s ∆满足.s a b cs a b c ∆∆∆∆≤++第二章 插值法1. 根据(2.2)定义的范德蒙行列式,令2000011211121()(,,,,)11n n n n n n n n n x x x V x V x x x x x x x xx x ----==L L L L L L L L L证明()n V x 是n 次多项式,它的根是01,,n x x -L ,且101101()(,,,)()()n n n n V x V x x x x x x x ---=--L L .2. 当x = 1 , -1 , 2 时, f (x)= 0 , -3 , 4 ,求f (x )的二次插值多项式.3.4. 给出cos x ,0°≤x ≤90°的函数表,步长h =1′=(1/60)°,若函数表具有5位有效数字,研究用线性插值求cos x 近似值时的总误差界.5. 设0k x x kh =+,k =0,1,2,3,求032max ()x x x l x ≤≤.6. 设jx 为互异节点(j =0,1,…,n ),求证:i) 0()(0,1,,);nk k j j j x l x x k n =≡=∑Lii) 0()()1,2,,).n k jj j xx l x k n =-≡0(=∑L7. 设[]2(),f x C a b ∈且()()0f a f b ==,求证21()()().8max max a x ba xb f x b a f x ≤≤≤≤≤-"8. 在44x -≤≤上给出()xf x e =的等距节点函数表,若用二次插值求xe 的近似值,要使截断误差不超过610-,问使用函数表的步长h 应取多少?9. 若2n n y =,求4n y ∆及4n y δ.10. 如果()f x 是m 次多项式,记()()()f x f x h f x ∆=+-,证明()f x 的k 阶差分()(0)k f x k m ∆≤≤是m k -次多项式,并且()0(m l f x l +∆=为正整数).11. 证明1()k k k k k k f g f g g f +∆=∆+∆.12. 证明110010.n n kkn n k k k k f gf g f g g f --+==∆=--∆∑∑13. 证明1200.n j n j y y y -=∆=∆-∆∑14. 若1011()n nn n f x a a x a x a x --=++++L 有n 个不同实根12,,,n x x x L ,证明{10,02;, 1.1()n k njk n a k n j jx f x -≤≤-=-=='∑15. 证明n 阶均差有下列性质: i)若()()F x cf x =,则[][]0101,,,,,,n n F x x x cf x x x =L L ;ii) 若()()()F x f x g x =+,则[][][]010101,,,,,,,,,n n n F x x x f x x x g x x x =+L L L .16. 74()31f x x x x =+++,求0172,2,,2f ⎡⎤⎣⎦L 及0182,2,,2f ⎡⎤⎣⎦L . 17. 证明两点三次埃尔米特插值余项是(4)22311()()()()/4!,(,)k k k k R x f x x x x x x ++=ξ--ξ∈并由此求出分段三次埃尔米特插值的误差限.18. 求一个次数不高于4次的多项式()P x ,使它满足(0)(1)P P k =-+并由此求出分段三次埃尔米特插值的误差限. 19. 试求出一个最高次数不高于4次的函数多项式()P x ,以便使它能够满足以下边界条件(0)(0)0P P ='=,(1)(1)1P P ='=,(2)1P =.20. 设[](),f x C a b ∈,把[],a b 分为n 等分,试构造一个台阶形的零次分段插值函数()n x ϕ并证明当n →∞时,()n x ϕ在[],a b 上一致收敛到()f x .21. 设2()1/(1)f x x =+,在55x -≤≤上取10n =,按等距节点求分段线性插值函数()h I x ,计算各节点间中点处的()h I x 与()f x 的值,并估计误差. 22. 求2()f x x =在[],a b 上的分段线性插值函数()h I x ,并估计误差. 23. 求4()f x x =在[],a b 上的分段埃尔米特插值,并估计误差.试求三次样条插值并满足条件 i) (0.25) 1.0000,(0.53)0.6868;S S '='= ii)(0.25)(0.53)0.S S "="=25. 若[]2(),f x C a b ∈,()S x 是三次样条函数,证明 i)[][][][]222()()()()2()()()bbbbaaaaf x dx S x dx f x S x dx S x f x S x dx"-"="-"+""-"⎰⎰⎰⎰;ii) 若()()(0,1,,)i i f x S x i n ==L ,式中i x 为插值节点,且01n a x x x b =<<<=L ,则[][][]()()()()()()()()()baS x f x S x dx S b f b S b S a f a S a ""-"="'-'-"'-'⎰.26. 编出计算三次样条函数()S x 系数及其在插值节点中点的值的程序框图(()S x 可用(8.7)式的表达式).第三章 函数逼近与计算1. (a)利用区间变换推出区间为[],a b 的伯恩斯坦多项式.(b)对()sin f x x =在[]0,/2π上求1次和三次伯恩斯坦多项式并画出图形,并与相应的马克劳林级数部分和误差做比较. 2. 求证:(a)当()m f x M ≤≤时,(,)n m B f x M ≤≤. (b)当()f x x =时,(,)n B f x x =.3. 在次数不超过6的多项式中,求()sin 4f x x =在[]0,2π的最佳一致逼近多项式.4. 假设()f x 在[],a b 上连续,求()f x 的零次最佳一致逼近多项式.5. 选取常数a ,使301max x x ax≤≤-达到极小,又问这个解是否唯一?6. 求()sin f x x =在[]0,/2π上的最佳一次逼近多项式,并估计误差.7. 求()xf x e =在[]0,1上的最佳一次逼近多项式.8. 如何选取r ,使2()p x x r =+在[]1,1-上与零偏差最小?r 是否唯一? 9. 设43()31f x x x =+-,在[]0,1上求三次最佳逼近多项式.10. 令[]()(21),0,1n n T x T x x =-∈,求***0123(),(),(),()T x T x T x T x .11. 试证{}*()nT x 是在[]0,1上带权ρ=的正交多项式.12. 在[]1,1-上利用插值极小化求11()f x tg x -=的三次近似最佳逼近多项式. 13. 设()xf x e =在[]1,1-上的插值极小化近似最佳逼近多项式为()n L x ,若n f L ∞-有界,证明对任何1n ≥,存在常数n α、n β,使11()()()()(11).n n n n n T x f x L x T x x ++α≤-≤β-≤≤14. 设在[]1,1-上234511315165()128243843840x x x x x x ϕ=-----,试将()x ϕ降低到3次多项式并估计误差. 15. 在[]1,1-上利用幂级数项数求()sin f x x =的3次逼近多项式,使误差不超过0.005.16. ()f x 是[],a a -上的连续奇(偶)函数,证明不管n 是奇数或偶数,()f x 的最佳逼近多项式*()n n F x H ∈也是奇(偶)函数.17. 求a 、b 使[]22sin ax b x dx π+-⎰为最小.并与1题及6题的一次逼近多项式误差作比较.18. ()f x 、[]1(),g x C a b ∈,定义()(,)()();()(,)()()()();bbaaa f g f x g x dxb f g f x g x dx f a g a =''=''+⎰⎰问它们是否构成内积?19. 用许瓦兹不等式(4.5)估计6101x dx x +⎰的上界,并用积分中值定理估计同一积分的上下界,并比较其结果.20. 选择a ,使下列积分取得最小值:1122211(),x ax dx x ax dx----⎰⎰.21. 设空间{}{}10010121,,,span x span x x 1ϕ=ϕ=,分别在1ϕ、2ϕ上求出一个元素,使得其为[]20,1x C ∈的最佳平方逼近,并比较其结果.22. ()f x x =在[]1,1-上,求在{}2411,,span x x ϕ=上的最佳平方逼近.23.sin (1)arccos ()n n x u x +=是第二类切比雪夫多项式,证明它有递推关系()()()112n n n u x xu x u x +-=-.24. 将1()sin 2f x x=在[]1,1-上按勒让德多项式及切比雪夫多项式展开,求三次最佳平方逼近多项式并画出误差图形,再计算均方误差.25. 把()arccos f x x =在[]1,1-上展成切比雪夫级数.26.2y a bx =+.27.用最小二乘拟合求.29. 编出用正交多项式做最小二乘拟合的程序框图. 30. 编出改进FFT 算法的程序框图. 31. 现给出一张记录{}{}4,3,2,1,0,1,2,3k x =,试用改进FFT 算法求出序列{}k x 的离散频谱{}k C (0,1,,7).k =L第四章 数值积分与数值微分1. 确定下列求积公式中的待定参数,使其代数精度尽量高,并指明所构造出的求积公式所具有的代数精度:(1)101()()(0)()hh f x dx A f h A f A f h --≈-++⎰; (2)21012()()(0)()hh f x dx A f h A f A f h --≈-++⎰;(3)[]1121()(1)2()3()/3f x dx f f x f x -≈-++⎰;(4)[][]20()(0)()/1(0)()hf x dx h f f h ah f f h ≈++'-'⎰.2. 分别用梯形公式和辛普森公式计算下列积分:(1)120,84xdx n x =+⎰; (2)1210(1),10x e dx n x --=⎰;(3)1,4n =⎰;(4),6n =.3. 直接验证柯特斯公式(2.4)具有5次代数精度.4. 用辛普森公式求积分10xedx-⎰并计算误差.5. 推导下列三种矩形求积公式:(1)2()()()()()2ba f f x dxb a f a b a 'η=-+-⎰; (2)2()()()()()2ba f f x dxb a f b b a 'η=---⎰;(3)3()()()()()224baa b f f x dx b a f b a +"η=-+-⎰. 6. 证明梯形公式(2.9)和辛普森公式(2.11)当n →∞时收敛到积分()baf x dx⎰.7. 用复化梯形公式求积分()baf x dx⎰,问要将积分区间[],a b 分成多少等分,才能保证误差不超过ε(设不计舍入误差)?8.1x e dx-,要求误差不超过510-.9. 卫星轨道是一个椭圆,椭圆周长的计算公式是S a =θ,这里a 是椭圆的半长轴,c 是地球中心与轨道中心(椭圆中心)的距离,记h 为近地点距离,H 为远地点距离,6371R =公里为地球半径,则(2)/2,()/2a R H h c H h =++=-.我国第一颗人造卫星近地点距离439h =公里,远地点距离2384H =公里,试求卫星轨道的周长.10. 证明等式3524sin3!5!n nn n ππππ=-+-L试依据sin(/)(3,6,12)n n n π=的值,用外推算法求π的近似值.11. 用下列方法计算积分31dyy ⎰并比较结果.(1) 龙贝格方法;(2) 三点及五点高斯公式;(3) 将积分区间分为四等分,用复化两点高斯公式.12. 用三点公式和五点公式分别求21()(1)f x x =+在x =1.0,1.1和1.2处的导数值,并估计误()f x第五章 常微分方程数值解法1. 就初值问题0)0(,=+='y b ax y 分别导出尤拉方法和改进的尤拉方法的近似解的表达式，并与准确解bx ax y +=221相比较。
西工大2020年4月《计算机操作系统》作业机考参考答案
西工大2020年4月《计算机操作系统》作业机考参考答案试卷总分:100 得分:96本作业每次试题内容相同，只是题目和选项顺序是随机调整的，大家可放心下载使用一、单选题(共50 道试题,共100 分)1. 一作业进入内存后，则所属该作业的进程初始时处于（）状态。
A.运行B.等待C.就绪D.收容正确答案:C2.用户在一次计算过程中，或者一次事物处理中，要求计算机完成所做的工作的集合，这是指（）。
A.进程B.程序C.作业D.系统调用正确答案:C3. 下列方法中哪一个破坏了“循环等待”条件？（）A.银行家算法B.一次性分配策略（即预分配策略）C.剥夺资源法D.资源有序分配正确答案:D4. 现代操作系统的基本特征是（）、资源共享和操作的异步性。
A.多道程序设计B.中断处理C.程序的并发执行D.实现分时与实时处理正确答案:C5. 操作系统中采用缓冲技术的目的是为了增强系统（）的能力。
A.串行操作B.控制操作C.重执操作D.并行操作正确答案:D6. 进程间的基本关系为（）。
A.相互独立与相互制约B.同步与互斥C.并行执行与资源共享D.信息传递与信息缓冲正确答案:C7. 碎片现象的存在使得（）。
A.内存空间利用率降低B.内存空间利用率提高C.内存空间利用率得以改善D.内存空间利用率不影响正确答案:A8. 进程和程序是两个既有联系又有区别的概念，下面描述中，（）是错误的。
A.进程是动态的，程序是静态的B.一个程序可对应多个进程C.进程有生命周期，而程序相对是永久的D.程序是可以并发执行的正确答案:D9. 外存（如磁盘）上存放的程序和数据（）。
A.可由CPU直接访问B.必须在CPU访问之前移入内存C.必须由文件系统管理D.必须由进程调度程序管理正确答案:B10. 只能在核心态下执行的指令是（）。
A.读时钟日期B.屏蔽所有中断C.改变文件内容D.调用库函数正确答案:B11. 下述几种引起进程永远等待的情况中，只有（）是属于死锁问题。
西北工业大学计算方法第五周作业答案第五章答案
J
( x1,
x2 )
取得极小值必然有
J
( x1, x1
x2 )
0,
J ( x1, x2 ) 0; 即 x2
J
( x1 , x1
x2 )
2( x1
x2
1)
2( x1
x2
2)
6(3x1
x2
1)
J
( x1 , x2
x2
)
2( x1
x2
1)
2( x1
x2
2)
2(3x1
x2
1)
化简得到正规法方程组:
a b
34.9 102.8
求解得到 a 0.95, b 3.11 ；进而拟合关系为 0.95 3.11 。
计算决定系数：
1
∑ ∑
̄
1
0.007 48.3675
0.99985
由此可知，本问题中拟合得到的线性关系显著（显著或不显著）。
基于线性模型得到的正规方程组为：
5 3
3 1.9
ln I0 a
0.028638 0.304447
上述正规方程组的解为： ln I0 1.7179; a 2.8726 ；
模型方程为： I 5.5728e2.8726t ；
利用上述模型方程计算不同 ti 时刻的发射强度预测值 并填入下表：
ti
或
x1 0.7500, x2 1.2500 。
解法 2：通过最小二乘法原理求解：
设上述三个式子在取定 x1, x2 时的误差为 e1, e2 , e3 ，上述误差的平方和
3
J (x1, x2 ) ei2 (x1 x2 1)2 (x1 x2 2)2 i 1
(3x1 x2 1)2
西北工业大学数值分析(附答案)
西北工业大学数值分析习题集第一章 绪 论1. 设x >0,x 的相对误差为δ,求ln x 的误差.2. 设x 的相对误差为2％,求nx 的相对误差.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 -=( n=1,2,…)计算到100Y .27.982(五位有效数字),试问计算100Y 将有多大误差?7. 求方程25610x x -+=的两个根,使它至少具有四位有效数字27.982).8. 当N 充分大时,怎样求211N 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. 计算61)f =, 1.4≈,利用下列等式计算,哪一个得到的结果最好?3--13.()ln(f x x =,求f (30)的值.若开平方用六位函数表,问求对数时误差有多大?若改用另一等价公式ln(ln(x x =-计算,求对数时误差有多大?14. 试用消元法解方程组{101012121010;2.x x x x +=+=假定只用三位数计算,问结果是否可靠?15. 已知三角形面积1sin ,2s ab c =其中c 为弧度,02c π<<,且测量a ,b ,c 的误差分别为,,.a b c ∆∆∆证明面积的误差s ∆满足.s a b cs a b c ∆∆∆∆≤++第二章 插值法1. 根据(2.2)定义的范德蒙行列式,令2000011211121()(,,,,)11n n n n n n n n n x x x V x V x x x x x x x xx x ----==证明()n V x 是n 次多项式,它的根是01,,n x x -,且 101101()(,,,)()()n n n n V x V x x x x x x x ---=--.2. 当x = 1 , -1 , 2 时, f (x)= 0 , -3 , 4 ,求f (x )的二次插值多项式.3.4. 给出cos x ,0°≤x ≤90°的函数表,步长h =1′=(1/60)°,若函数表具有5位有效数字,研究用线性插值求cos x 近似值时的总误差界.5. 设0k x x kh =+,k =0,1,2,3,求032max ()x x x l x ≤≤.6. 设jx 为互异节点(j =0,1,…,n ),求证:i) 0()(0,1,,);nk k j j j x l x x k n =≡=∑ii) 0()()1,2,,).n k jj j xx l x k n =-≡0(=∑7. 设[]2(),f x C a b ∈且()()0f a f b ==,求证21()()().8max max a x ba xb f x b a f x ≤≤≤≤≤-"8. 在44x -≤≤上给出()xf x e =的等距节点函数表,若用二次插值求xe 的近似值,要使截断误差不超过610-,问使用函数表的步长h 应取多少?9. 若2n n y =,求4n y ∆及4n y δ.10. 如果()f x 是m 次多项式,记()()()f x f x h f x ∆=+-,证明()f x 的k 阶差分()(0)k f x k m ∆≤≤是m k -次多项式,并且()0(m l f x l +∆=为正整数).11. 证明1()k k k k k k f g f g g f +∆=∆+∆.12. 证明110010.n n kkn n k k k k f gf g f g g f --+==∆=--∆∑∑13. 证明1200.n jn j yy y -=∆=∆-∆∑14. 若1011()n n n n f x a a x a x a x --=++++有n 个不同实根12,,,n x x x ,证明{10,02;, 1.1()n k njk n a k n j jx f x -≤≤-=-=='∑15. 证明n 阶均差有下列性质: i)若()()F x cf x =,则[][]0101,,,,,,n n F x x x cf x x x =;ii) 若()()()F x f x g x =+,则[][][]010101,,,,,,,,,n n n F x x x f x x x g x x x =+.16. 74()31f x x x x =+++,求0172,2,,2f ⎡⎤⎣⎦及0182,2,,2f ⎡⎤⎣⎦.17. 证明两点三次埃尔米特插值余项是(4)22311()()()()/4!,(,)k k k k R x f x x x x x x ++=ξ--ξ∈并由此求出分段三次埃尔米特插值的误差限.18. 求一个次数不高于4次的多项式()P x ,使它满足(0)(1)P P k =-+并由此求出分段三次埃尔米特插值的误差限. 19. 试求出一个最高次数不高于4次的函数多项式()P x ,以便使它能够满足以下边界条件(0)(0)0P P ='=,(1)(1)1P P ='=,(2)1P =.20. 设[](),f x C a b ∈,把[],a b 分为n 等分,试构造一个台阶形的零次分段插值函数()n x ϕ并证明当n →∞时,()n x ϕ在[],a b 上一致收敛到()f x .21. 设2()1/(1)f x x =+,在55x -≤≤上取10n =,按等距节点求分段线性插值函数()h I x ,计算各节点间中点处的()h I x 与()f x 的值,并估计误差. 22. 求2()f x x =在[],a b 上的分段线性插值函数()h I x ,并估计误差. 23. 求4()f x x =在[],a b 上的分段埃尔米特插值,并估计误差.试求三次样条插值并满足条件 i) (0.25) 1.0000,(0.53)0.6868;S S '='= ii)(0.25)(0.53)0.S S "="=25. 若[]2(),f x C a b ∈,()S x 是三次样条函数,证明 i)[][][][]222()()()()2()()()bbbbaaaaf x dx S x dx f x S x dx S x f x S x dx"-"="-"+""-"⎰⎰⎰⎰;若()()(0,1,,)i i f x S x i n ==,式中i x 为插值节点,且01n a x x x b =<<<=,则[][][]()()()()()()()()()baS x f x S x dx S b f b S b S a f a S a ""-"="'-'-"'-'⎰.26. 编出计算三次样条函数()S x 系数及其在插值节点中点的值的程序框图(()S x 可用(8.7)式的表达式).第三章 函数逼近与计算1. (a)利用区间变换推出区间为[],a b 的伯恩斯坦多项式.(b)对()sin f x x =在[]0,/2π上求1次和三次伯恩斯坦多项式并画出图形,并与相应的马克劳林级数部分和误差做比较. 2. 求证:(a)当()m f x M ≤≤时,(,)n m B f x M ≤≤. (b)当()f x x =时,(,)n B f x x =.3. 在次数不超过6的多项式中,求()sin 4f x x =在[]0,2π的最佳一致逼近多项式.4. 假设()f x 在[],a b 上连续,求()f x 的零次最佳一致逼近多项式.5. 选取常数a ,使301max x x ax≤≤-达到极小,又问这个解是否唯一?6. 求()sin f x x =在[]0,/2π上的最佳一次逼近多项式,并估计误差.7. 求()xf x e =在[]0,1上的最佳一次逼近多项式.8. 如何选取r ,使2()p x x r =+在[]1,1-上与零偏差最小?r 是否唯一?9. 设43()31f x x x =+-,在[]0,1上求三次最佳逼近多项式.10. 令[]()(21),0,1n n T x T x x =-∈,求***0123(),(),(),()T x T x T x T x .11. 试证{}*()nT x 是在[]0,1上带权ρ=的正交多项式.12. 在[]1,1-上利用插值极小化求11()f x tg x -=的三次近似最佳逼近多项式. 13. 设()xf x e =在[]1,1-上的插值极小化近似最佳逼近多项式为()n L x ,若n f L ∞-有界,证明对任何1n ≥,存在常数n α、n β,使11()()()()(11).n n n n n T x f x L x T x x ++α≤-≤β-≤≤14. 设在[]1,1-上234511315165()128243843840x x x x x x ϕ=-----,试将()x ϕ降低到3次多项式并估计误差. 15. 在[]1,1-上利用幂级数项数求()sin f x x =的3次逼近多项式,使误差不超过0.005.16. ()f x 是[],a a -上的连续奇(偶)函数,证明不管n 是奇数或偶数,()f x 的最佳逼近多项式*()n n F x H ∈也是奇(偶)函数.17. 求a 、b 使[]220sin ax b x dx π+-⎰为最小.并与1题及6题的一次逼近多项式误差作比较.1g x C a b∈(),f x、[],定义18.()()(,)()();()(,)()()()();bbaaa f g f x g x dxb f g f x g x dx f a g a =''=''+⎰⎰问它们是否构成内积?19. 用许瓦兹不等式(4.5)估计6101x dx x +⎰的上界,并用积分中值定理估计同一积分的上下界,并比较其结果.20. 选择a ,使下列积分取得最小值:1122211(),x ax dx x ax dx----⎰⎰.21. 设空间{}{}10010121,,,span x span x x 1ϕ=ϕ=,分别在1ϕ、2ϕ上求出一个元素,使得其为[]20,1x C ∈的最佳平方逼近,并比较其结果.22. ()f x x =在[]1,1-上,求在{}2411,,span x x ϕ=上的最佳平方逼近.23.sin (1)arccos ()n n x u x +=是第二类切比雪夫多项式,证明它有递推关系()()()112n n n u x xu x u x +-=-.24. 将1()sin 2f x x=在[]1,1-上按勒让德多项式及切比雪夫多项式展开,求三次最佳平方逼近多项式并画出误差图形,再计算均方误差.25. 把()arccos f x x =在[]1,1-上展成切比雪夫级数.26.2y a bx =+.27.用最小二乘拟合求.29. 编出用正交多项式做最小二乘拟合的程序框图. 30. 编出改进FFT 算法的程序框图. 31. 现给出一张记录{}{}4,3,2,1,0,1,2,3k x =,试用改进FFT 算法求出序列{}k x 的离散频谱{}k C (0,1,,7).k =第四章 数值积分与数值微分1. 确定下列求积公式中的待定参数,使其代数精度尽量高,并指明所构造出的求积公式所具有的代数精度:(1)101()()(0)()hh f x dx A f h A f A f h --≈-++⎰; (2)21012()()(0)()hh f x dx A f h A f A f h --≈-++⎰;(3)[]1121()(1)2()3()/3f x dx f f x f x -≈-++⎰;(4)[][]20()(0)()/1(0)()hf x dx h f f h ah f f h ≈++'-'⎰.2. 分别用梯形公式和辛普森公式计算下列积分:(1)120,84xdx n x =+⎰; (2)1210(1),10x e dx n x --=⎰;(3)1,4n =⎰;(4),6n =.3. 直接验证柯特斯公式(2.4)具有5次代数精度.4. 用辛普森公式求积分10xedx-⎰并计算误差.5. 推导下列三种矩形求积公式:(1)2()()()()()2ba f f x dxb a f a b a 'η=-+-⎰; (2)2()()()()()2ba f f x dxb a f b b a 'η=---⎰;(3)3()()()()()224baa b f f x dx b a f b a +"η=-+-⎰. 6. 证明梯形公式(2.9)和辛普森公式(2.11)当n →∞时收敛到积分()baf x dx⎰.7. 用复化梯形公式求积分()baf x dx⎰,问要将积分区间[],a b 分成多少等分,才能保证误差不超过ε(设不计舍入误差)?8.1x e dx-,要求误差不超过510-.9. 卫星轨道是一个椭圆,椭圆周长的计算公式是S a =θ,这里a 是椭圆的半长轴,c 是地球中心与轨道中心(椭圆中心)的距离,记h 为近地点距离,H 为远地点距离,6371R =公里为地球半径,则(2)/2,()/2a R H h c H h =++=-.我国第一颗人造卫星近地点距离439h =公里,远地点距离2384H =公里,试求卫星轨道的周长. 10. 证明等式3524sin3!5!n nn n ππππ=-+-试依据sin(/)(3,6,12)n n n π=的值,用外推算法求π的近似值.11. 用下列方法计算积分31dyy ⎰并比较结果.(1) 龙贝格方法;(2) 三点及五点高斯公式;(3) 将积分区间分为四等分,用复化两点高斯公式.12. 用三点公式和五点公式分别求21()(1)f x x =+在x =1.0,1.1和1.2处的导数值,并估计误()f x第五章 常微分方程数值解法1. 就初值问题0)0(,=+='y b ax y 分别导出尤拉方法和改进的尤拉方法的近似解的表达式，并与准确解bx ax y +=221相比较。
西北工业大学_计算方法作业_答案
西工大计算方法作业答案参考答案 第一章1 *1x =1.7； *2x =1.73； *3x =1.732 。
2．3． （1） ≤++)(*3*2*1x x x e r 0.00050； （注意：应该用相对误差的定义去求） （2） ≤)(*3*2*1x x x e r 0.50517； （3） ≤)/(*4*2x x e r 0.50002。
4．设6有n 位有效数字，由6≈2.4494……，知6的第一位有效数字1a =2。
令3)1()1(1*1021102211021)(-----⨯≤⨯⨯=⨯=n n r a x ε 可求得满足上述不等式的最小正整数n =4，即至少取四位有效数字，故满足精度要求可取6≈2.449。
5. 答：（1）*x (0>x )的相对误差约是*x 的相对误差的1/2倍； （2）n x )(* 的相对误差约是*x 的相对误差的n 倍。
6． 根据********************sin 21)(cos 21sin 21)(sin 21sin 21)(sin 21)(c b a c e c b a c b a b e c a c b a a e c b S e r ++≤ =******)()()(tgcc e b b e a a e ++ 注意当20*π<<c 时，0**>>c tgc ，即1*1*)()(--<c tgc 。
则有)()()()(****c e b e a e S e r r r r ++<7．设20=y ，41.1*0=y ，δ=⨯≤--2*001021y y由 δ1*001*111010--≤-=-y y y y ，δ2*111*221010--≤-=-y y y yδ10*991*10101010--≤-=-y y y y即当0y 有初始误差δ时，10y 的绝对误差的绝对值将减小1010-倍。
而11010<<-δ，故计算过程稳定。
西工大2021年4月机考《计算方法》作业参考答案非免费
西工大2021年4月机考《计算方法》作业试卷总分:100 得分:100答案网叫福到（这四个字的拼音）一、单选题(共50 道试题,共100 分)1. {A.0B.1C.2D.{<i正确答案:D2. {A.5次B.2次C.1次D.3次正确答案:3.{A.{<imgB.{<imgC.2D.1正确答案:4.3.141580是π的有( )位有效数字的近似值。
A.6B.5C.4D.7正确答案:5.舍入误差是(?? ?)产生的误差。
A.只取有限位数B.模型准确值与用数值方法求得的准确值C.观察与测量D.数学模型准确值与实际值正确答案:6.反幂法是用来求矩阵( )特征值及相应特征向量的一种向量迭代法。
A.按模最大B.按模最小C.全部D.任意一个正确答案:7. {A.{<img s>B.{<imgC.{<img sD.0正确答案:8.求解线性方程组Ax=b的LLT分解法中，A须满足的条件是( )。
A.对称阵B.正定矩阵C.任意阵D.各阶顺序主子式均不为零正确答案:9.幂法的收敛速度与特征值的分布（）。
A.有关B.不一定C.无关D.无法判断正确答案:10. {A.3B.2C.1D.4正确答案:11. {A.0B.{<img s>C.2D.1正确答案:12.{A.舍入B.观测C.模型D.截断正确答案:13.解线性方程组Ax=b的迭代格式x(k+1)=Mx(k)+f收敛的充要条件是( )。
A.{<img s>B.{<img s">C.{<imgD.{<img sr>正确答案:14. {A.2B.{<img >C.3D.1正确答案:15.求解线性方程组Ax=b的LLT分解法中，A须满足的条件是( )。
A.对称阵B.各阶顺序主子式均大于零C.任意阵D.各阶顺序主子式均不为零正确答案:16.三点的高斯求积公式的代数精度为( ).A.2B.5C.3D.4正确答案:17. 三点的高斯型求积公式的代数精度为(????)。
计算方法各习题及参考答案
第二章 数值分析2.1 已知多项式432()1p x x x x x =-+-+通过下列点：试构造一多项式()q x 通过下列点：答案：54313()()()3122q x p x r x x x x x =-=-++-+． 2.2 观测得到二次多项式2()p x 的值：表中2()p x 的某一个函数值有错误，试找出并校正它．答案：函数值表中2(1)p -错误，应有2(1)0p -=．2.3 利用差分的性质证明22212(1)(21)/6n n n n +++=++ ．2.4 当用等距节点的分段二次插值多项式在区间[1,1]-近似函数xe 时，使用多少个节点能够保证误差不超过61102-⨯． 答案：需要143个插值节点．2.5 设被插值函数4()[,]f x C a b ∈，()3()h H x 是()f x 关于等距节点01n a x x x b =<<<= 的分段三次艾尔米特插值多项式，步长b ah n-=．试估计()3||()()||h f x H x ∞-．答案：()443||()()||384h M f x H x h ∞-≤．第三章 函数逼近3.1 求()sin ,[0,0.1]f x x x =∈在空间2{1,,}span x x Φ=上最佳平方逼近多项式，并给出平方误差．答案：()sin f x x =的二次最佳平方逼近多项式为-522sin ()0.832 440 710 1.000 999 10.024 985 1x p x x x ≈=-⨯+-，二次最佳平方逼近的平方误差为0.122-1220(sin )())0.989 310 710x p x dx δ=-=⨯⎰．3.2 确定参数,a b c 和，使得积分2121(,,)[I a b c ax bx c -=++-⎰取最小值．答案：810, 0, 33a b c ππ=-== 3.3 求多项式432()251f x x x x =+++在[1,1]-上的３次最佳一致逼近多项式()p x ．答案：()f x 的最佳一致逼近多项式为323()74p x x x =++． 3.4 用幂级数缩合方法，求() (11)x f x e x =-≤≤上的３次近似多项式6,3()p x ，并估计6,3||()()||f x p x ∞-．答案：236,3()0.994 574 650.997 395 830.542 968 750.177 083 33p x x x x =+++， 6,3||()()||0.006 572 327 7f x p x ∞-≤3.5 求() (11)xf x e x =-≤≤上的关于权函数()x ρ=的三次最佳平方逼近多项式3()S x ，并估计误差32||()()||f x S x -和3||()()||f x S x ∞-．答案：233()0.994 5710.997 3080.542 9910.177 347S x x x x =+++，32||()()||0.006 894 83f x S x -=，3||()()||0.006 442 575f x S x ∞-≤．第四章 数值积分与数值微分4.1 用梯形公式、辛浦生公式和柯特斯公式分别计算积分1(1,2,3,4)n x dx n =⎰，并与精确值比较．答案：计算结果如下表所示4.2 确定下列求积公式中的待定参数，使得求积公式的代数精度尽量高，并指明所确定的求积公式具有的代数精度． （１）101()()(0)()hh f x dx A f h A f A f h --≈-++⎰（２）11211()[(1)2()3()]3f x dx f f x f x -≈-++⎰ （３）20()[(0)()][(0)()]2h h f x dx f f h h f f h α''≈++-⎰答案：（１）具有三次代数精确度（２）具有二次代数精确度（３）具有三次代数精确度．4.3 设10h x x =-，确定求积公式12300101()()[()()][()()][]x x x x f x dx h Af x Bf x h Cf x Df x R f ''-=++++⎰中的待定参数,,,A B C D ，使得该求积公式的代数精确度尽量高，并给出余项表达式．答案：3711,,,20203020A B C D ====-，(4)6()[]1440f R f h η=，其中01(,)x x η∈．4.4 设2()P x 是以0,,2h h 为插值点的()f x 的二次插值多项式，用2()P x 导出计算积分30()hI f x dx =⎰的数值积分公式h I ，并用台劳展开法证明：453(0)()8h I I h f O h '''-=+． 答案：3203()[(0)3(2)]4h h I p x dx h f f h ==+⎰．4.5 给定积分10sin xI dx x =⎰（１）运用复化梯形公式计算上述积分值，使其截断误差不超过31102-⨯． （２）取同样的求积节点，改用复化辛浦生公式计算时，截断误差是多少？（３）要求的截断误差不超过610-，若用复化辛浦生公式，应取多少个节点处的函数值？ 答案：（１）只需7.5n ≥，取９个节点，0.946I ≈（２）4(4)46111|[]||()|()0.271102880288045n b a R f h f η--=-≤=⨯ （３）取７个节点处的函数值．4.6 用变步长的复化梯形公式和变步长的复化辛浦生公式计算积分10sin xI dx x =⎰．要求用事后误差估计法时，截断误不超过31102-⨯和61102-⨯． 答案：使用复化梯形公式时，80.946I T ≈=满足精度要求；使用复化辛浦生公式时，40.946 083I s ≈=满足精度要求．4.7（１）利用埃尔米特插值公式推导带有导数值的求积公式2()()[()()][()()][]212ba b a b a f x dx f a f b f b f a R f --''=+--+⎰，其中余项为 5(4)()[](), (,)4!30b a R f f a b ηη-=∈． （２）利用上述公式推导带修正项的复化梯形求积公式020()[()()]12Nx N N x h f x dx T f x f x ''≈--⎰，其中 0121[()2()2()2()()]2N N N hT f x f x f x f x f x -=+++++ ，而 00, (0,1,2,,), i N x x ih i N Nh x x =+==- ．4.8 用龙贝格方法计算椭圆2214x y +=的周长，使结果具有五位有效数字． 答案：49.6884l I =≈．4.9确定高斯型求积公式0011()()()x dx A f x A f x ≈+⎰的节点0x ，1x 及系数0A ，1A ．答案：00.289 949x =，10.821 162x =，00.277 556A =，10.389 111A =．4.10 验证高斯型求积公式00110()()()x e f x dx A f x A f x +∞-≈+⎰的系数及节点分别为0001 2 2A A x x ===-=+第五章 解线性方程组的直接法5.1 用按列选主元的高斯-若当消去法求矩阵A 的逆矩阵，其中111210110A -⎛⎫ ⎪= ⎪ ⎪-⎝⎭． 答案： 1110331203321133A -⎛⎫ ⎪ ⎪⎪=- ⎪ ⎪⎪-- ⎪⎝⎭5.2 用矩阵的直接三角分解法解方程组1234102050101312431701037x x x x ⎛⎫⎛⎫⎛⎫⎪ ⎪ ⎪⎪ ⎪ ⎪= ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪⎝⎭⎝⎭⎝⎭ 答案： 42x =，32x =，21x =，11x =．5.3 用平方根法(Cholesky 分解法)求解方程组12341161 4.25 2.750.51 2.75 3.5 1.25x x x -⎛⎫⎛⎫⎛⎫ ⎪⎪ ⎪-=- ⎪⎪ ⎪ ⎪⎪ ⎪⎝⎭⎝⎭⎝⎭答案： 12x =，21x =，31x =-．5.4 用追赶法求解三对角方程组123421113121112210x x x x ⎛⎫⎛⎫⎛⎫ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪= ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪⎝⎭⎝⎭⎝⎭ 答案：42x =，31x =-，21x =，10x =．第六章 解线性代数方程组的迭代法６.１ 对方程1212123879897x x x x x x x -+=⎧⎪-+=⎨⎪--=⎩作简单调整，使得用高斯－赛得尔迭代法求解时对任意初始向量都收敛，并取初始向量(0)[0 0 0]T x =，用该方法求近似解(1)k x+，使(1)()3||||10k k x x +-∞-≤． 答案：近似解为(4)[1.0000 1.0000 1.0000]Tx =．６.２ 讨论松弛因子 1.25ω=时，用SOR 方法求解方程组121232343163420412x x x x x x x +=⎧⎪+-=⎨⎪-+=-⎩ 的收敛性．若收敛，则取(0)[0 0 0]T x=迭代求解，使(1)()41||||102k k x x +-∞-<⨯． 答案：方程组的近似解为*1 1.50001x =，*2 3.33333x =，*3 2.16667x =-．６.３ 给定线性方程组Ax b =，其中111221112211122A ⎛⎫ ⎪ ⎪⎪= ⎪ ⎪ ⎪ ⎪⎝⎭，证明用雅可比迭代法解此方程组发散，而高斯－赛得尔迭代法收敛．６.４ 设有方程组112233302021212x b x b x b -⎛⎫⎛⎫⎛⎫ ⎪⎪ ⎪= ⎪⎪ ⎪ ⎪⎪ ⎪-⎝⎭⎝⎭⎝⎭，讨论用雅可比方法和高斯－赛得尔方法解此方程组的收敛性．如果收敛，比较哪种方法收敛较快．答案：雅可比方法收敛，高斯－赛得尔方法收敛,且较快．６.５ 设矩阵A 非奇异．求证：方程组Ax b =的解总能通过高斯－赛得尔方法得到．６.６ 设()ij n nA a ⨯=为对称正定矩阵，对角阵1122(,,,)nn D diag a a a = ．求证：高斯－赛得尔方法求解方程组1122D AD x b --=时对任意初始向量都收敛．第七章 非线性方程求根例７.４ 对方程230xx e -=确定迭代函数()x ϕ及区间[,]a b ，使对0[,]x a b ∀∈，迭代过程1(), 0,1,2,k x x k ϕ+== 均收敛，并求解．要求51||10k k x x -+-<． 答案：若取2()x x ϕ=，则在[1,0]-中满足收敛性条件，因此迭代法121, 0,1,2,k x k x k +== 在(1,0)-中有惟一解．取00.5x =-，*70.458960903x x ≈=-．取2()x x ϕ=，在[0,1上满足收敛性条件，迭代序列121, 0,1,2,k x k x k +== 在[0,1]中有惟一解．取00.5x =，*140.910001967x x ≈=- 在[3,4]上，将原方程改写为23xe x =，取对数得2ln(3)()x x x ϕ==．满足收敛性条件，则迭代序列21ln(3), 0,1,2,k k x x k +== 在[3,4]中有惟一解．取0 3.5x =， *16 3.733067511x x ≈=．例７.６ 对于迭代函数2()(3)x x c x ϕ=+-，试讨论：（１）当c 为何值时，1()k k x x ϕ+=产生的序列{}k x（２）c 取何值时收敛最快？（３）取1,2c =-()x ϕ51||10k k x x -+-<．答案：（１）(c ∈时迭代收敛．（２）c =时收敛最快．（３）分别取1, 2c =--，并取0 1.5x =，计算结果如下表７.７所示表７.７例７.13 设不动点迭代1()k x x ϕ+=的迭代函数()x ϕ具有二阶连续导数，*x 是()x ϕ的不动点，且*()1x ϕ'≠，证明Steffensen 迭代式21(), (), 0,1,2,()2k k k k k k k k k k k y x z x k y x x x z y xϕϕ+==⎧⎪=-⎨=-⎪-+⎩二阶收敛于*x ．例7.15 设2()()()()()x x p x f x q x f x ϕ=--，试确定函数()p x 和()q x ，使求解()0f x =且以()x ϕ为迭代函数的迭代法至少三阶收敛．答案：1()()p x f x ='，31()()2[()]f x q x f x ''=' 例７.19 设()f x 在[,]a b 上有高阶导数，*(,)x a b ∈是()0f x =的(2)m m ≥重根，且牛顿法收敛，证明牛顿迭代序列{}k x 有下列极限关系：111lim2k kk k k k x x m x x x -→∞-+-=-+．第八章 矩阵特征值8.1 用乘幂法求矩阵A 的按模最大的特征值与对应的特征向量，已知5500 5.51031A -⎛⎫ ⎪=- ⎪ ⎪-⎝⎭，要求(1)()611||10k k λλ+--<，这里()1k λ表示1λ的第k 次近似值．答案：15λ≈，对应的特征向量为[5,0,0]T-；25λ≈-，对应的特征向量为[5,10,5]T --． 8.2 用反幂法求矩阵110242012A -⎛⎫⎪=-- ⎪ ⎪-⎝⎭的按模最小的特征值．知A 的按模较大的特征值的近似值为15λ=，用5p =的原点平移法计算1λ及其对应的特征向量．答案：(1) A 的按模最小的特征值为30.2384428λ≈(2) 1 5.1248854λ≈，对应的特征向量为(8)[0.242 4310, 1 ,0.320 011 7]T U =--．8.3 设方阵A 的特征值都是实数，且满足121, ||||n n λλλλλ>≥≥> ，为求1λ而作原点平移，试证：当平移量21()2n p λλ=+时，幂法收敛最快． 8.4 用二分法求三对角对称方阵1221221221A ⎛⎫ ⎪ ⎪= ⎪ ⎪ ⎪⎝⎭的最小特征值，使它至少具有２位有效数字．答案：取5 2.234375λ≈-即有２位有效数字．8.5 用平面旋转变换和反射变换将向量[2 3 0 5]T x =变为与1[1 0 0 0]Te =平行的向量．答案：203/2/00001010/0T ⎛⎫⎪- ⎪=⎪--⎝0.324 442 8400.486 664 26200.811 107 1040.486 664 2620.812 176 04800.298 039 92200100.811 107 1040.298 039 92200.530 266 798H --⎛⎫⎪--⎪= ⎪ ⎪⎪--⎝⎭8.6 若532644445A -⎛⎫⎪=- ⎪ ⎪-⎝⎭，试把A 化为相似的上Hessenberg 阵，然后用QR 方法求A 的全部特征值．第九章 微分方程初值问题的数值解法9.1 用反复迭代（反复校正）的欧拉预估－校正法求解初值问题0, 0<0.2(0)1y y x y '+=≤⎧⎨=⎩，要求取步长0.1h =，每步迭代误差不超过510-． 答案： [4]11(0.1)0.904 762y y y ≈==，[4]22(0.2)0.818 594y y y ≈==9.2 用二阶中点格式和二阶休恩格式求初值问题2, 0<0.4(0)1dy x y x dx y ⎧=+≤⎪⎨⎪=⎩的数值解(取步长0.2h =，运算过程中保留五位小数)．答案：用二阶中点格式，取初值01y =计算得0n =时，1211.000 00, 1.200 00, (0.2)=1.240 00K K y y ==≈ 1n =时，1221.737 60, 2.298 72, (0.4)=1.699 74K K y y ==≈用二阶休恩格式，取初值01y =计算得0n =时，1211.000 00, 1.266 67, (0.2)=1.240 00K K y y ==≈ 1n =时，1221.737 60, 2.499 18, (0.4)=1.701 76K K y y ==≈9.3 用如下四步四阶阿达姆斯显格式1123(5559379)/24n n n n n n y y h f f f f +---=+-+-求初值问题, (0)1y x y y '=+=在[0,0.5]上的数值解．取步长0.1h =，小数点后保留８位．答案：4(0.4)0.583 640 216y y ≈=，5(0.5) 1.797 421 984y y ≈=． 9.4 为使二阶中点公式1(,(,))22n n n n n n h hy y hf x y f x y +=+++，求解初值问题 , (0)y y y aλλ'=-⎧⎨=⎩为实常数绝对稳定，试求步长h 的大小应受到的限制条件． 答案：2h λ≤．9.5 用如下反复迭代的欧拉预估－校正格式(0)1(1)()111(,)[(,)(,)]2 0,1,2,; 0,1,2,nn n n k k n n n n n n y y hf x y h y y f x y f x y k n +++++⎧=+⎪⎪=++⎨⎪⎪==⎩，求解初值问题sin(), 01(0)1x y e xy x y '⎧=<≤⎨=⎩时，如何选择步长h ，使上述格式关于k 的迭代收敛． 答案：2h e<时上述格式关于k 的迭代是收敛的．9.6 求系数,,,a b c d ，使求解初值问题0(,), ()y f x y y x a '==的如下隐式二步法221()n n n n n y ay h bf cf df +++=+++的误差阶尽可能高，并指出其阶数．答案：系数为142,,33a b d c ====，此时方法的局部截断误差阶最高，为五阶5()O h ．9.7 试用欧拉预估－校正法求解初值问题, (0)=1, 0<0.2()/, (0)2dyxy z y dxx dz x y z z dx⎧=-⎪⎪≤⎨⎪=+=⎪⎩，取步长0.1h =，小数点后至少保留六位．答案：由初值00(0)1, (0)2y y z z ====可计算得110.800 000z 2.050 000y =⎧⎨=⎩ ， 11(0.1)0.801 500(0.1) 2.046 951y y z z ≈=⎧⎨≈=⎩ 220.604 820z 2.090 992y =⎧⎨=⎩ ， 22(0.2)0.604 659(0.2) 2.088 216y y z z ≈=⎧⎨≈=⎩。
西工大NOJ答案完全版
西工大NOJ答案完全版#include<stdio.h>int main(){int a,b,sum;scanf("%d%d",&a,&b);sum=a+b;printf("%d\n",sum);return 0;}#include<stdio.h>#define PI 3.int main(){double r,h,l,s,sq,vq,vz;scanf("%lf%lf",&r,&h);l=2*PI*r;s=PI*r*r;sq=4*PI*r*r;vz=PI*r*r*h;printf("%.2lf\n%.2lf\n%.2lf\n%.2lf\n%.2lf\n",l,s,sq,vq,vz);ret urn 0;}#include<stdio.h>int main(){double ma,eng,c,sum,ave;scanf("%lf%lf%lf",&ma,&eng,&c);sum=ma+eng+c;ave=sum/3;printf("%lf\n%lf\n",sum,ave);return 0;}#include<stdio.h>int main(){int a,b,c,m;scanf("%d%d%d",&a,&b,&c);if (a>b) m=a;else m=b;if (m<c) m=c;printf("%d",m);return 0;}#include<stdio.h>int main(){int n;scanf("%d",&n);if((1000<n<)&&(n/1000==n%10)&&(n/100%10==n/10%10)) printf("yes\n");else if((100<n<=1000)&&(n/100==n%10)) printf("yes\n");else if((10<n<=100)&&(n/10==n%10)) printf("yes\n");else if(0<n<=10) printf("yes\n");else printf("no\n");return 0;}#include<stdio.h>int main(){double l,bon;scanf("%lf",&l);if(l<=10) bon=l*0.1;else if(l<20) bon=1+(l-10)*0.075; else if(l<40) bon=1.75+(l-20)*0.05; else if(l<60) bon=2.75+(l-40)*0.03; else if(l<100) bon=3.35+(l-60)*0.015; else bon=3.95+(l-100)*0.01;printf("%lf\n",bon);return 0;}#include<stdio.h>int main(){double d,m;scanf("%lf",&d);if(d<=2) m=7;else if(d<=15){if(d-2==(int)(d-2)) m=7+(d-2)*1.5;else m=7+((int)(d-2)+1)*1.5;}else if(d-15==(int)(d-15)) m=26.5+(d-15)*2.1; else m=26.5+((int)(d-15)+1)*2.1;printf("%lf\n",m);return 0;}#include<stdio.h>int main(){int y,m,d,Days,sum;scanf("%d-%d-%d",&y,&m,&d);if((y%4==0&&y%100!=0)||(y%400==0)) Days=29; else Days=28;switch(m){case 1:sum=d;break;case 2:sum=31+d;break;case 3:sum=31+Days+d;break;case 4:sum=62+Days+d;break;case 5:sum=92+Days+d;break; case 6:sum=123+Days+d;break; case 7:sum=153+Days+d;break; case 8:sum=184+Days+d;break; case 9:sum=215+Days+d;break; case 10:sum=245+Days+d;break; case 11:sum=276+Days+d;break; case 12:sum=307+Days+d;break; }printf("%d\n",sum);return 0;}#include<stdio.h>int main(){int i;scanf("%d",&i);if(i>=90) printf("A\n");else if(i>=80) printf("B\n"); else if(i>=70) printf("C\n");else if(i>=60) printf("D\n");else printf("E\n");return 0;}#include<stdio.h>int main(){double x,y;scanf("%lf,%lf",&x,&y);if((x-2)*(x-2)+(y-2)*(y-2)<=1) printf("10");else if((x-2)*(x-2)+(y+2)*(y+2)<=1) printf("10"); else if((x+2)*(x+2)+(y-2)*(y-2)<=1) printf("10"); else if((x+2)*(x+2)+(y+2)*(y+2)<=1) printf("10"); else printf("0");return 0;}#include<stdio.h>int main(){double l,x,r;scanf("%lf %lf",&l,&r);while((2*l*l*l-4*l*l+3*l-6)!=0&&(2*r*r*r-4*r*r+3*r-6)!=0){x=(l+r)/2;if((2*l*l*l-4*l*l+3*l-6)*(2*x*x*x-4*x*x+3*x-6)<=0) r=x;else l=x;}if(2*l*l*l-4*l*l+3*l-6==0) printf("%.2lf",l);else printf("%.2lf",r);return 0;}#include<stdio.h>#include<math.h>int main(){int i=800,t=2,cnt=0,sum=0;double e=-1;while(i>=500){while(t<=i-1){if(i%t==0) break;t++;}if(t==i) e=pow(-1,cnt),sum=sum+e*i,cnt++; i--;t=2;}printf("%d %d",cnt,sum);return 0;}#include<stdio.h>#include<math.h>int main(){int a=1;double b=1,pi=0,c=1;while(fabs(c)>=1e-6)pi=pi+c,b=b+2,a=-a,c=a/b;pi=pi*4;printf("%lf\n",pi);return 0;}#include<stdio.h>int main(){int a1=1,a2=1,n=2,sum=2,t; while(sum<=100){t=a1;a1=a2;a2=t+2*a2;sum=sum+a2;n++;}printf("%d\n",n-1);while(sum<=1000){t=a1;a1=a2;a2=t+2*a2;sum=sum+a2;n++;}printf("%d\n",n-1); while(sum<=){t=a1;a1=a2;a2=t+2*a2;sum=sum+a2;n++;}printf("%d\n",n-1);}#include<stdio.h>int main(){int x,a,s,n=1;scanf("%d %d",&x,&a); s=x;if(a!=0){for(;n<a;n++){s=s*x;if(s>=1000) s=s/100%10*100+s/10%10*10+s%10; }}printf("%d\n",s);return 0;}#include<stdio.h>int main(){int m,n,s;scanf("%d",&n);s=n*n*n;printf("%d*%d*%d=%d=",n,n,n,s);for(m=1;s!=n*m;m++);if(n%2==1){for(s=-n/2;s<n/2;s++)printf("%d+",m+2*s);printf("%d",m+n/2*2);}else{for(s=-n/2;s<n/2-1;s++)printf("%d+",m+s*2+1); printf("%d",m+(n/2-1)*2+1);}return 0;}#include<stdio.h>int main(){char a,b,c,x,y,z;a='A',b='B',c='C',x='X',y='Y',z='Z';printf("%c=%c\n",a,z);printf("%c=%c\n",b,x);printf("%c=%c\n",c,y);return 0;}#include<stdio.h>int main(){int a,b,t;scanf("%d %d",&a,&b);if(a>b)t=a,a=b,b=t;for(;a<b;a++){for(t=2;t<a;t++)if(a%t==0) break; if(t==a)printf("%d ",a);}return 0;}#include<stdio.h>int main(){int n=1;double a1=1,a2=2,a3,sum=2; while(n<=19){a3=a1+a2;sum=sum+a3/a2;a1=a2;a2=a3;n++;}printf("%lf\n",sum);return 0;}#include<stdio.h>#include<math.h> int main(){double a;int n=0;scanf("%lf",&a);a=fabs(a);if(a<=1)printf("0\n"); else{while(a>1){a=a/10;n++;}printf("%d\n",n);}return 0;}#include<stdio.h> int main(){int a=1,b=0,t,m,n=0; scanf("%d",&t); while(n<t){m=b;b=3*a+2*b;a=m;n++;}printf("%d %d",a,b); return 0;}#include<stdio.h> #include<math.h> int main(){int n;int f(int n);scanf("%d",&n);printf("%d\n",f(n));}int f(int n){int a;if(n==1||n==2)a=0;else if(n==3)a=1;else if(n==4)a=3;elsea=f(n-1)*2+pow(2,n-4)-f(n-4); return a;}#include<stdio.h>int main(){int n,x=1234,y=1,a,b,c,d,e,f,g,h,i,j;scanf("%d",&n);for(;x<&&y<;x++){y=x*n;a=x/%10;b=x/1000%10;c=x/100%10;d=x/10%10;e=x%10;f=y/%10;g=y/1000%10;h=y/100%10;i=y/10%10;j=y%10;if(a==b||a==c||a==d||a==e||a==f||a==g||a==h||a==i||a==j) continue;if(b==c||b==d||b==e||b==f||b==g||b==h||b==i||b==j) continue;if(c==d||c==e||c==f||c==g||c==h||c==i||c==j) continue;if(d==e||d==f||d==g||d==h||d==i||d==j) continue;if(e==f||e==g||e==h||e==i||e==j) continue;if(f==g||f==h||f==i||f==j) continue;if(g==h||g==i||g==j) continue;if(h==i||h==j) continue;if(i==j) continue;printf("%05d/%05d=%d\n",y,x,n); }return 0;}#include<stdio.h>#include<stdlib.h>#include<math.h>int main(){int m,n,i;double x,s=0;scanf("%d%d",&n,&m);for(i=n;i<=m;i++){x=pow(i,2.0);s=s+1/x;}printf("%.5lf\n",s);return 0;}#include<stdio.h>int main(){int x,y,a,b,L;double t;scanf("%d%d%d%d%d",&x,&y,&a,&b,&L); if(a==b) printf("impossible\n");else if(x>y){if(a>b)y=L-x+y,t=(double)y/(a-b);else y=x-y,t=(double)y/(b-a);if((int)t==t)printf("%d\n",(int)t);else printf("%lf\n",t);}else {if(a>b)y=y-x,t=(double)y/(a-b);else y=L-y+x,t=(double)y/(b-a);if((int)t==t)printf("%d\n",(int)t);else printf("%lf\n",t);}return 0;}#include<stdio.h>int _max(int a,int b){return a>b?a:b;}int a[20];int f[20][20];int main(){int n,i,j,s=0;scanf("%d",&n);for(i=0;i<n;i++)scanf("%d",&a[i]);for(j=1;j<n;j++)f[0][0]=a[0],f[0][j]=f[0][j-1]*a[j]; for(i=1;i<n;i++){f[i][i-1]=1;for(j=i;j<n;j++)f[i][j]=f[i][j-1]*a[j];}for(i=0;i<n;i++)for(j=i;j<n;j++)s=_max(s,f[i][j]);if(s==0)printf("-1\n");else printf("%d\n",s);return 0;}#include<stdio.h>int main(){int x=192,y,z,a,b,c,d,e,f,g,h,i; for(;x<328;x++){y=2*x;z=3*x;a=x/100%10;b=x/10%10;c=x%10;d=y/100%10;e=y/10%10;f=y%10;g=z/100%10;h=z/10%10;i=z%10;if(a==b||a==c||a==d||a==e||a==f||a==g||a==h||a==i||a==0) continue;if(b==c||b==d||b==e||b==f||b==g||b==h||b==i||b==0) continue;if(c==d||c==e||c==f||c==g||c==h||c==i||c==0) continue;if(d==e||d==f||d==g||d==h||d==i||d==0) continue;if(e==f||e==g||e==h||e==i||e==0) continue;if(f==g||f==h||f==i||f==0) continue;if(g==h||g==i||g==0) continue;if(h==i||h==0) continue;printf("%d %d %d\n",x,y,z);}return 0;}#include<stdio.h>int main(){int a,b,c,sum=10;scanf("%d%d%d",&a,&b,&c);for(;sum<=100;sum++){if(sum%3==a&&sum%5==b&&sum%7==c){ printf("%d\n", sum);break;}}if(sum==101)printf("-1\n");return 0;}#include<stdio.h>#include<math.h>int main(){int ce,y,m,n,a=0;scanf("%d",&n);for(ce=0;;ce+=100){for(y=ce+1;y<ce+100;y+=2){for(m=3;m<sqrt(y);m+=2){if(y%m==0) break;}if(m>=sqrt(y)) break;}if(y==ce+101) a++;if(a==n) break;}printf("%d %d\n",ce,ce+99);return 0;}#include<stdio.h>int main(){int n,i;scanf("%d",&n);for(i=1;i<=n;i++){if(i%7==0) printf("%d ",i);else if(i/1000%10==7||i/100%10==7||i/10%10==7||i%10==7) printf("%d ",i);}return 0;}#include<stdio.h>double a[];int main(){int n,i;double ave,sum=0;scanf("%d",&n);for(i=0;i<n;i++){scanf("%lf",&a[i]);sum=sum+a[i];}ave=sum/n;for(i=0,sum=0;i<n;i++)sum=sum+(a[i]-ave)*(a[i]-ave); printf("%lf\n",sum);return 0;}#include<stdio.h>int f[];int main(){int n,i;scanf("%d",&n);for(i=0;i<n;i++){int a,b;scanf("%d%d",&a,&b);f[i]=a+b;}for(i=0;i<n;i++){if(f[i]>100)f[i]=f[i]/10%10*10+f[i]%10,printf("%d\n",f[i]);else printf("%d\n",f[i]);}return 0;}#include<stdio.h>#include<math.h>int getbit(int n,int k){n=n>>k-1;return n&1;}int main(){int n,k;scanf("%d%d",&n,&k);printf("%d",getbit(n,k));return 0;}#include<stdio.h>int _min(int n){int i=0;if(n>=100) i=n/100%10,n=n/10%10*10+n%10;if(n>=50) i+=1,n-=50;if(n>=10) i+=n/10%10,n%=10; if(n>=5) i+=1,n-=5;if(n>=2) i+=n/2,n%=2;if(n==1) i+=1;return i;}int f[100];int main(){int n,i,sum;scanf("%d",&n);for(i=0;i<n;i++)scanf("%d",&f[i]);for(sum=0,i=0;i<n;i++)sum=sum+_min(f[i]);printf("%d\n",sum);return 0;}#include<stdio.h>int fceil(double x){if(x>0){if((int)x==x) return x; else return (int)x+1;}else return (int)x;}int main(){double x;scanf("%lf",&x);printf("%d\n",fceil(x)); return 0;}#include<stdio.h>int _move(int value,int n) {int t;if(n>0)t=value>>n|value<<32-n;else {n=-n;t=value<<n|value>>32-n;}return t;}int main(){int value,n;scanf("%d%d",&value,&n);printf("%d\n",_move(value,n));return 0;}#include<stdio.h>void timesum(int AH,int AM,int AS,int BH,int BM,int BS){ int h=0,m=0,s=0;if(AS+BS>59) s=(AS+BS)-60,m++;else s=AS+BS;if(AM+BM>59) m+=((AM+BM)-60),h++;else m=AM+BM;h+=(AH+BH);printf("%d %d %d\n",h,m,s);}int main(){int AH,AM,AS,BH,BM,BS;scanf("%d%d%d%d%d%d",&AH,&AM,&AS,&BH,&BM,& BS);timesum(AH,AM,AS,BH,BM,BS);return 0;}#include<stdio.h>void f(long int n){int m;m=n%10;printf("%d",m);n/=10;if(n!=0) f(n);}int main(){long int n;scanf("%d",&n);f(n);return 0;}#include<stdio.h>int getfloor(double x) {if(x>0) return (int)x;else return (int)x-1;}int main(){double x;scanf("%lf",&x);printf("%d\n",getfloor(x)); return 0;}#include<stdio.h>inline int xchg(unsigned char n){n=n>>4|n<<4;return n;}int main(){unsigned char n;scanf("%d",&n);printf("%d\n",xchg(n));return 0;}#include<stdio.h>void QuickSort(int A[100],int s,int m) {int i,t;for(;s<=m;s++){for(i=s;i<=m;i++)if(A[s]<A[i]){t=A[s];A[s]=A[i];A[i]=t;}printf("%d ",A[s]);}}int main(){int s,m,n,i,A[100]; scanf("%d",&n);for(i=0;i<n;i++)scanf("%d",&A[i]); scanf("%d%d",&s,&m); for(i=0;i<s;i++)printf("%d ",A[i]); QuickSort(A,s,m);for(i=m+1;i<n;i++)printf("%d ",A[i]);return 0;}#include<stdio.h>double avg(int A[100],int s,int e) {double sum;for(sum=0;s<=e;s++)sum=sum+A[s];return sum;}int main(){int A[100],s,e,i,n;double a;scanf("%d",&n);for(i=0;i<n;i++)scanf("%d",&A[i]);scanf("%d%d",&s,&e);a=avg(A,s,e)/(e-s+1); printf("%lf\n",a); return 0;}#include<stdio.h>int main(){int n,a[100],i,j,s=0,t=0; scanf("%d",&n);for(i=0;i<n;i++) scanf("%d",&a[i]);for(i=0;i<n-1;i++)for(j=i+1;j<n;j++) {if(a[i]>a[j]) s++;else if(a[i]<a[j]) t++; }s=s>t?t:s;printf("%d\n",s); return 0;}#include<stdio.h>void SelectionSort(int A[],int s,int m) {int i,t,n;for(n=s;n<=s+m-1;n++){for(i=n;i<=s+m-1;i++)if(A[n]<A[i]){t=A[n];A[n]=A[i];A[i]=t;}printf("%d ",A[n]);}}int main(){int s,m,n,i,A[100];scanf("%d",&n);for(i=0;i<n;i++)scanf("%d",&A[i]);scanf("%d%d",&s,&m);for(i=0;i<s;i++)printf("%d ",A[i]);SelectionSort(A,s,m);for(i=m+1;i<n;i++)printf("%d ",A[i]);return 0;}#include<stdio.h>void SelectionSort(int A[],int s,int m) {int i,t,n;for(n=s;n<=s+m-1;n++){for(i=n;i<=s+m-1;i++)if(A[n]<A[i]){t=A[n];A[n]=A[i];A[i]=t;}printf("%d ",A[n]);}}int main(){int s,m,n,i,A[100]; scanf("%d",&n);for(i=0;i<n;i++)scanf("%d",&A[i]); scanf("%d%d",&s,&m); for(i=0;i<s;i++)printf("%d ",A[i]); SelectionSort(A,s,m); for(i=m+1;i<n;i++) printf("%d ",A[i]); return 0;}#include<stdio.h>int main(){int a[100][100],n,i,j; scanf("%d",&n);for(i=0;i<n;i++)for(j=0;j<n;j++)scanf("%d",&a[i][j]); for(i=0;i<n;i++)for(j=0;j<n;j++){if(j<n-i-1) printf(" "); else printf("%d ",a[i][j]); if(j==n-1) printf("\n"); }return 0;}#include<stdio.h>int main(){int a[100][100],n,i,j; scanf("%d",&n);for(i=0;i<n;i++)for(j=0;j<n;j++)scanf("%d",&a[i][j]); for(i=0;i<n;i++)for(j=0;j<n;j++){if(j<i) printf(" ");else printf("%d ",a[i][j]); if(j==n-1) printf("\n"); }return 0;}#include<stdio.h>int main(){int a[100][100],n,i,j; scanf("%d",&n);for(i=0;i<n;i++)for(j=0;j<n;j++)scanf("%d",&a[i][j]);for(i=0;i<n;i++)for(j=0;j<n;j++){if(j>n-i-1) printf(" ");else printf("%d ",a[i][j]);if(j==n-1) printf("\n");}return 0;}#include<stdio.h>int main(){int a[100][100],n,m,i,j,s1,s2; scanf("%d%d",&n,&m);for(i=0;i<n;i++)for(j=0;j<m;j++)scanf("%d",&a[i][j]);for(j=0,s1=0;j<m;j++) s1=s1+a[0][j]+a[n-1][j]; for(i=1,s2=0;i<n-1;i++) {s1=s1+a[i][0]+a[i][m-1]; for(j=1;j<m-1;j++)s2=s2+a[i][j];}printf("%d\n",s1-s2); return 0;}#include<stdio.h>int main(){int a[100],b[100],i,j,n; scanf("%d",&n);for(i=0;i<n;i++){scanf("%d",&a[i]);b[i]=a[i];}for(i=1;i<n;i+=2)printf("%d ",b[i]);return 0;}#include<stdio.h>#include<string.h>int main(){char a[10][10];char b[10];int i,j;for(i=0;i<10;i++)scanf("%s",&a[i]);for(i=0;i<9;i++)for(j=0;j<9-i;j++)if(strcmp(a[j],a[j+1])>0)strcpy(b,a[j]),strcpy(a[j],a[j+1]),strcpy(a[j+1],b);for(i=0;i<10; i++)printf("%s ",a[i]);return 0;}#include<stdio.h>int main(){char a[],b[1000][10];int i,j=0,k=0,m;gets(a);for(i=0;;i++){if(a[i]==' '||a[i]==','||a[i]=='.'||a[i]=='\0'){for(m=0;j<i;j++,m++) b[k][m]=a[j];k++;for(;;i++) if(a[i+1]!=' '&&a[i+1]!=','&&a[i+1]!='.'){ j=i+1;break;}}if(a[i]=='\0')break;}for(i=k-1;i>=0;i--) printf("%s ",b[i]);return 0;}#include<stdio.h>#include<string.h>int main(){charA[][5]={"I","II","III","IV","V","VI","VII","VIII","IX","X","XI"," XII"};char B[5];int i=0;while((B[i]=getchar())!='\n')i++;B[i]='\0'; //这一句不能落下，否则结果都是13，不过为啥？for(i=0;;i++){if(strcmp(A[i],B)==0){printf("%d\n",i+1);break;}}return 0;}#include<stdio.h>#include<string.h>void stringmerge(char s1[],char s2[]) {strcat(s1,s2);}int main(){char s1[],s2[];gets(s1);gets(s2);stringmerge(s1,s2);puts(s1);return 0;}#include<stdio.h>void deletechar(char s[],char c) {int i;for(i=0;;i++){if(s[i]==c)continue;if(putchar(s[i])=='\0')break;}}int main(){char s[],c;gets(s);scanf("%c",&c);deletechar(s,c);return 0;}#include<stdio.h>#include<string.h>int main(){char a[1000];int i,n;gets(a);n=strlen(a);for(i=0;i<n;){if(a[i]=='y'&&a[i+1]=='o'&&a[i+2]=='u') {printf("%s","we");i+=3;continue;}putchar(a[i]);i++;}return 0;}#include<stdio.h>。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
参考答案 第一章1 *1x =1.7； *2x =1.73； *3x =1.732 。
2．3． （1） ≤++)(*3*2*1x x x e r 0.00050； （注意：应该用相对误差的定义去求） （2） ≤)(*3*2*1x x x e r 0.50517； （3） ≤)/(*4*2x x e r 0.50002。
4．设6有n 位有效数字，由6≈2.4494……，知6的第一位有效数字1a =2。
令3)1()1(1*1021102211021)(-----⨯≤⨯⨯=⨯=n n r a x ε 可求得满足上述不等式的最小正整数n =4，即至少取四位有效数字，故满足精度要求可取6≈2.449。
5. 答：（1）*x (0>x )的相对误差约是*x 的相对误差的1/2倍；（2）nx )(* 的相对误差约是*x 的相对误差的n 倍。
6． 根据********************sin 21)(cos 21sin 21)(sin 21sin 21)(sin 21)(c b a c e c b a c b a b e c a c b a a e c b S e r ++≤=******)()()(tgcc e b b e a a e ++ 注意当20*π<<c 时，0**>>c tgc ，即1*1*)()(--<c tgc 。
则有)()()()(****c e b e a e S e r r r r ++<7．设20=y ，41.1*0=y ，δ=⨯≤--2*001021y y 由 δ1*001*111010--≤-=-y y y y ，δ2*111*221010--≤-=-y y y yMδ10*991*10101010--≤-=-y y y y即当0y 有初始误差δ时，10y 的绝对误差的绝对值将减小1010-倍。
而11010<<-δ，故计算过程稳定。
8. 变形后的表达式为：（1）)1ln(2--x x =)1ln(2-+-x x（2）arctgx x arctg -+)1(=)1(11++x x arctg（3）1ln )1ln()1(ln 1--++=⎰+N N N N dx x N N=ΛΛ+-+-+32413121)1ln(N N N N 1ln )11ln()1(-+++=N N N N =1)1ln()11ln(-+++N NN （4）xx sin cos 1-=x x cos 1sin +=2x tg第二章1．绝对误差限312110-⨯, 对分8次2． (1) 隔根区间[0, 0.8]；(2) 等价变形 )2ln(x x -=； 迭代公式Λ,2,1)2ln(1=-=-n x x n n 。
(3) 收敛性论证：用局部收敛性定理论证。
3． (1) 7210-=x x ；(2) 2/)7(lg +=x x ； (3) 31+=x x ；4． 143)(2++='x x x f牛顿迭代公式为：ΛΛ143122231++-++-='-=+n n n n n n n n n n x x x x x x )x (f )x (f x x 列表计算根的近似值为6．⎪⎪⎭⎫⎝⎛+=+==+223123132)(n n nn n n x a x x a x x x ϕ 证明：2()3,()6f x x f x x '''==Q当0x >时，()0,()0;f x f x '''>>当0x <时，()0,()0;f x f x '''<<因此，对于0>a，当0x ≥00()()0f x f x ''>，牛顿迭代法收敛，当0x ∈时，)23001022022)033x x a x x x x+==>1x ≥1x。
对于0a <，当00x <<时，00()()0f x f x ''>，牛顿迭代法收敛；当0x ∈时，)201022)03x x x x=<1x <1x当0a =时，迭代变为312233k k k k k x x x x x +=-=该迭代发对于任何0x R ∈均收敛。
第三章1． x 1=2，x 2=1，x 3=1/22． ⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡---=-3132132310313101A 3． L = ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-153012001 ， U = ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--2400410321 y 1 =14, y 2 = -10, y 3 = -72x 1 =1, x 2 =2, x 3 =34. x 1≈-4.00, x 2≈3.00, x 3≈2.005. B 的特征值为：0，0，0，ρ(B)=0<1(E －B 1)-1B 2的特征值为：0，2，2，ρ[(E －B 1)-1B 2]=2>1. 6. x (5)=(0.4999, 1.0004, -0.4997)T 7.∣a ∣>2第四章1.5615.3)7(1≈λ 相应近似特征向量为 c =( 2258 ， 1268 ， 2258 ）T ，（ 0≠c ）第五章1． 取0x =100、1x =121用线性插值时，115≈10.7143；取0x =100、1x =121、2x =144用二次插值时，115≈10.7228。
2．选取插值节点为：0x =1.4、1x =1.5、2x =1.6，)54.1(f ≈1.9447。
3．利用∑=+ω'=pj j p j p )x ()x (f ]x ,x ,x [f 0110Λ，并注意当n p ≤时，对p j ,,1,0Λ=，0)(=j x f ，故有 n p x x x f p ≤=0],,,[10Λ 而1+=n p 时，)()(11++'=n n x x f ω，故有11],,,[10+==n p x x x f p Λ，4. )(3x L =)(3x N =)926913(5123-+-x x x5. （1）用反插值法得根的近似值*α=0.3376； （2）用牛顿迭代法得根的近似值*α=0.337667。
6. 令311)3(10))()((!3)(max 11-+-≤≤≤---+-k k k x x x x x x x x x f k k ξ可求得h ≤0.2498（或h ≤0.2289）。
详解：由题义知，所采用的是三点等距插值，由误差公式：(3)2()()()()()3!1sin ()()()3!1()()()6k k k k k k k k k f R x x x h x x x x h x x h x x x x h x x h x x x x h ξξ=----+≤----+≤----+ 令 ()()()()k k k g x x x h x x x x h =----+由 ()0g x '=得：223()0k x x h --=得 ()g x 的驻点为：3k xx h =±% 故，{}11113max ()max (),(),(),()()k k k k k x x x g x g x g x g x g x g x -+-+≤≤===%%所以，3321()6R x ≤=令 3310-≤ 解得：13310)0.2498h -≤≈7. (1) 5982)(233+-+-=x x x x H22)4(3)2()1)((!41)(--=x x f x R ξ )2,1(∈ξ (2)61592)(233-+-=x x x x H)3()2)(1)((!41)(2)4(3---=x x x f x R ξ )3,1(∈ξ第六章1．正规方程组为⎪⎪⎭⎫ ⎝⎛493330⎪⎪⎭⎫ ⎝⎛21x x =⎪⎪⎭⎫⎝⎛2973 3888.21≈x ， 4456.02≈x2．正规方程组为⎪⎪⎭⎫⎝⎛7277699532753275⎪⎪⎭⎫ ⎝⎛b a =⎪⎪⎭⎫ ⎝⎛5.3693214.271 9726.0≈a ， 0500.0≈b20500.09726.0x y +=3． 取对数at I I -=0ln ln 相应的正规方程组为⎪⎪⎭⎫ ⎝⎛--03.25.35.37⎪⎪⎭⎫ ⎝⎛a I 0ln =⎪⎪⎭⎫ ⎝⎛-1858.09890.172825.1ln 0=I ， 8882.2≈a 6308.50≈Ite I 8882.26308.5-=4．正规方程组为⎪⎪⎭⎫⎝⎛6092.31781.31781.34⎪⎪⎭⎫ ⎝⎛b a =⎪⎪⎭⎫⎝⎛9607.124.14 4864.2≈a ， 4016.1≈b x y ln 4016.14864.2+=。