数值分析试卷及答案
- 格式:docx
- 大小:220.59 KB
- 文档页数:12
数值分析试题及答案一、选择题1. 下列哪个方法不适合用于求解非线性方程的根?A. 二分法B. 牛顿法C. 弦截法D. 正割法2. 当使用二分法求解非线性方程的根时,需要满足的条件是:A. 函数f(x)在区间[a, b]上连续B. 函数f(x)在区间[a, b]上单调递增C. 函数f(x)在区间[a, b]上存在根D. 函数f(x)在区间[a, b]上可导3. 数值积分是通过将定积分转化为求和的方法来近似计算积分值的过程。
下列哪个方法是常用的数值积分方法?A. 矩形法则B. 辛普森规则C. 梯形规则D. 高斯-勒让德法则4. 龙格-库塔法是常用于求解常微分方程的数值解法。
以下哪个选项是描述龙格-库塔法的特点?A. 该方法是一种多步法B. 该方法是一种多项式插值法C. 该方法是一种单步法D. 该方法是一种数值积分法5. 用有限差分法求解偏微分方程时,通常需要进行网格剖分。
以下哪个选项是常用的网格剖分方法?A. 多边形剖分法B. 三角剖分法C. 矩形剖分法D. 圆形剖分法二、解答题1. 将函数f(x) = e^x 在区间[0, 1]上用复化梯形规则进行数值积分,分为6个子区间,求得的近似积分值为多少?解:将区间[0, 1]等分为6个子区间,每个子区间的长度为h = (1-0)/6 = 1/6。
根据复化梯形规则的公式,近似积分值为:I ≈ (1/2) * h * [f(0) + 2f(1/6) + 2f(2/6) + 2f(3/6) + 2f(4/6) + 2f(5/6) +f(1)]≈ (1/2) * (1/6) * [e^0 + 2e^(1/6) + 2e^(2/6) + 2e^(3/6) + 2e^(4/6) +2e^(5/6) + e^1]2. 使用二分法求解方程 x^3 - 3x + 1 = 0 在区间[1, 2]上的根。
要求精确到小数点后三位。
解:首先需要判断方程在区间[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 =的牛顿-科特斯求积公式计算定积分1011I dx x=+⎰。
(10分) 五.用Newton 法求()cos 0f x x x =-=的近似解。
(10分) 六.试用Doolittle 分解法求解方程组:12325610413191963630x x x -⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥-=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥----⎣⎦⎣⎦⎣⎦ (10分)七.请写出雅可比迭代法求解线性方程组123123123202324812231530x x x x x x x x x ++=⎧⎪++=⎨⎪-+=⎩ 的迭代格式,并判断其是否收敛?(10分)八.就初值问题0(0)y yy y λ'=⎧⎨=⎩考察欧拉显式格式的收敛性。
(10分)《数值分析》(A )卷标准答案(2009-2010-1)一. 填空题(每小题3分,共12分) 1. ()1200102()()()()x x x x l x x x x x --=--; 2.7;3. 3,8;4. 2n+1。
数值分析期末试题及答案试题一:1. 简答题(共10分)a) 什么是数值分析?它的主要应用领域是什么?b) 请简要解释迭代法和直接法在数值计算中的区别。
2. 填空题(共10分)a) 欧拉方法是一种______型的数值解法。
b) 二分法是一种______法则。
c) 梯形法则是一种______型的数值积分方法。
3. 计算题(共80分)将以下函数进行数值求解:a) 通过使用二分法求解方程 f(x) = x^3 - 4x - 9 = 0 的近似解。
b) 利用欧拉方法求解微分方程 dy/dx = x^2 + 2x + 1, y(0) = 1 在 x = 1 处的解。
c) 使用梯形法则计算积分∫[0, π/4] sin(x) dx 的近似值。
试题二:1. 简答题(共10分)a) 请解释什么是舍入误差,并描述它在数值计算中的影响。
b) 请解释牛顿插值多项式的概念及其应用。
2. 填空题(共10分)a) 数值稳定性通过______号检查。
b) 龙格-库塔法是一种______计算方法。
c) 零点的迭代法在本质上是将方程______转化为______方程。
3. 计算题(共80分)使用牛顿插值多项式进行以下计算:a) 已知插值节点 (-2, 1), (-1, 1), (0, 2), (1, 4),求在 x = 0.5 处的插值多项式值。
b) 已知插值节点 (0, 1), (1, 2), (3, 7),求插值多项式,并计算在 x = 2 处的值。
c) 使用 4 阶龙格-库塔法求解微分方程 dy/dx = x^2 + 1, y(0) = 1。
答案:试题一:1. a) 数值分析是研究使用数值方法解决数学问题的一门学科。
它的主要应用领域包括数值微积分、数值代数、插值和逼近、求解非线性方程、数值积分和数值解微分方程等。
b) 迭代法和直接法是数值计算中常用的两种方法。
迭代法通过反复迭代逼近解,直到满足所需精度为止;而直接法则通过一系列代数运算直接得到解。
数值分析期末试题及答案一、选择题(每题5分,共20分)1. 在数值分析中,下列哪个算法不是用于求解线性方程组的?A. 高斯消元法B. 牛顿法C. 雅可比法D. 追赶法答案:B2. 插值法中,拉格朗日插值法属于:A. 多项式插值B. 样条插值C. 线性插值D. 非线性插值答案:A3. 以下哪个选项不是数值分析中的误差来源?A. 截断误差B. 舍入误差C. 计算误差D. 测量误差答案:C4. 在数值积分中,梯形法则的误差项是:A. O(h^2)B. O(h^3)C. O(h)D. O(1)答案:A二、填空题(每题5分,共20分)1. 牛顿插值法中,插值多项式的一般形式为:______。
答案:f(x) = a_0 + a_1(x-x_0) + a_2(x-x_0)(x-x_1) + ...2. 牛顿迭代法求解方程的根时,迭代公式为:x_{n+1} = x_n -f(x_n) / __________。
答案:f'(x_n)3. 在数值分析中,______ 用于衡量函数在区间上的近似积分值与真实积分值之间的差异。
答案:误差4. 线性方程组的解法中,______ 法是利用矩阵的LU分解来求解。
答案:克兰特三、解答题(每题10分,共60分)1. 给定函数f(x) = e^(-x),使用拉格朗日插值法,求x = 0.5时的插值值。
解答:首先选取插值节点x_0 = 0, x_1 = 0.5, x_2 = 1,对应的函数值分别为f(0) = 1, f(0.5) = e^(-0.5), f(1) = e^(-1)。
拉格朗日插值多项式为:L(x) = f(0) * (x-0.5)(x-1) / (0-0.5)(0-1) + f(0.5) * (x-0)(x-1) / (0.5-0)(0.5-1) + f(1) * (x-0)(x-0.5) / (1-0)(1-0.5)将x = 0.5代入得:L(0.5) = 1 * (0.5-0.5)(0.5-1) / (0-0.5)(0-1) + e^(-0.5) * (0.5-0)(0.5-1) / (0.5-0)(0.5-1) + e^(-1) * (0.5-0)(0.5-0.5) / (1-0)(1-0.5)计算得L(0.5) = e^(-0.5)。
数值分析试卷及答案数值分析试卷一、选择题(共10题,每题2分,共计20分)1. 数值分析的研究内容主要包括以下哪几个方面?A. 数值计算方法B. 数值误差C. 数值软件D. 数学分析答:A、B、C2. 下列哪种方法不属于数值积分的基本方法?A. 插值法B. 微积分基本公式C. 数值微积分D. 数值积分公式答:A3. 数值积分的目的是求解什么?A. 函数的导数B. 函数的原函数C. 函数的极值D. 函数的积分答:D4. 数值微分的目的是求解什么?A. 函数的导数B. 函数的原函数C. 函数的极值D. 函数的积分答:A5. 数值微分的基本方法有哪几种?A. 前向差分B. 后向差分C. 中心差分D. 插值法答:A、B、C6. 用数值方法求解方程的基本方法有哪几种?A. 迭代法B. 曲线拟合法C. 插值法D. 数值积分法答:A、B、C7. 用迭代法求方程的根时,当迭代结果满足何条件时可停止迭代?A. 当迭代结果开始发散B. 当迭代结果接近真实解C. 当迭代次数超过一定阈值D. 当迭代结果在一定范围内波动答:B8. 下列哪种插值方法能够确保经过所有给定数据点?A. 拉格朗日插值B. 牛顿插值C. 三次样条插值D. 二次插值答:A、B、C9. 数值解线性方程组的基本方法有哪几种?A. 直接法B. 迭代法C. 插值法D. 拟合法答:A、B10. 下列哪种方程求解方法适用于非线性方程?A. 直接法B. 迭代法C. 插值法D. 曲线拟合法答:B二、填空题(共5题,每题4分,共计20分)1. 数值积分的基本公式是_________。
答:牛顿-科特斯公式2. 数值微分的基本公式是_________。
答:中心差分公式3. 数值积分的误差分为_________误差和_________误差。
答:截断、舍入4. 用插值法求解函数值时,通常采用_________插值。
答:拉格朗日5. 数值解线性方程组的常用迭代法有_________方法和_________方法。
例1、 已知函数表求()f x 的Lagrange 二次插值多项式和Newton 二次插值多项式. 解:(1)插值基函数分别为()()()()()()()()()()1200102121()1211126x 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()1121213x x x x x x l x x x x x x x --+-===-+--+-故所求二次拉格朗日插值多项式为()()()()()()()()()()()2202()11131201241162314121123537623k k k L x y l x x x x x x x x x x x x x ==⎡⎤=-⨯--+⨯-+-+⨯+-⎢⎥⎣⎦=---++-=+-∑(2)一阶均差、二阶均差分别为[]()()[]()()[][][]010*********011201202303,11204,41234,,52,,126f 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 ---===-----===----===---故所求Newton 二次插值多项式为()()[]()[]()()()()()20010012012,,,35311126537623P 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 xx =++,[0,1]x ∈,试求()f x 在[0, 1]上关于()1x ρ=,{}span 1,x Φ=的最佳平方逼近多项式.解:若{}span 1,x Φ=,则0()1x ϕ=,1()x x ϕ=,且()1x ρ=,这样,有()()()()()()()()1120011011201100012101,11,,3123,,,,32269,324dx x dx xdx f x x dx f x x x dx ϕϕϕϕϕϕϕϕϕϕ========++==++=⎰⎰⎰⎰⎰ 所以,法方程为01123126119234a a ⎡⎤⎡⎤⎢⎥⎢⎥⎡⎤=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎢⎥⎢⎥⎣⎦⎣⎦,经过消元得01231162110123a a ⎡⎤⎡⎤⎢⎥⎢⎥⎡⎤=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎢⎥⎢⎥⎣⎦⎣⎦再回代解该方程,得到14a =,0116a =故,所求最佳平方逼近多项式为*111()46S x x =+ 例3、 设()xf x e =,[0,1]x ∈,试求()f x 在[0, 1]上关于()1x ρ=,{}span 1,x Φ=的最佳平方逼近多项式. 解:若{}span 1,x Φ=,则0()1x ϕ=,1()x x ϕ=,这样,有()()()()()()100012110101100100110,111,31,,2, 1.7183,1x x dx x dx xdx f e dx f xe dx ϕϕϕϕϕϕϕϕϕϕ===========⎰⎰⎰⎰⎰所以,法方程为0111 1.7183211123a a ⎡⎤⎢⎥⎡⎤⎡⎤=⎢⎥⎢⎥⎢⎥⎣⎦⎢⎥⎣⎦⎢⎥⎣⎦解法方程,得到00.8732a =,1 1.6902a =, 故,所求最佳平方逼近多项式为*1()0.8732 1.6902S x x =+例4、 用4n =的复合梯形和复合辛普森公式计算积分1⎰。
数值分析期末考试题及答案一、选择题(每题2分,共20分)1. 在数值分析中,下列哪个算法用于求解线性方程组?A. 牛顿法B. 高斯消元法C. 插值法D. 傅里叶变换答案:B2. 以下哪个选项不是数值分析中的误差类型?A. 舍入误差B. 截断误差C. 测量误差D. 累积误差答案:C3. 多项式插值中,拉格朗日插值法的特点是:A. 插值点必须等距分布B. 插值多项式的次数与插值点的个数相同C. 插值多项式是唯一的D. 插值多项式在插值点处的值都为1答案:B4. 在数值分析中,下列哪个方法用于求解非线性方程?A. 辛普森法则B. 牛顿迭代法C. 欧拉法D. 龙格-库塔法答案:B5. 以下哪个是数值稳定性的指标?A. 收敛性B. 收敛速度C. 条件数D. 误差传播答案:C二、简答题(每题10分,共20分)1. 简述高斯消元法求解线性方程组的基本原理。
答案:高斯消元法是一种直接解法,通过行变换将增广矩阵转换为上三角形式,然后通过回代求解线性方程组。
它包括三个基本操作:行交换、行乘以非零常数、行相加。
2. 解释什么是数值稳定性,并举例说明。
答案:数值稳定性是指数值解对输入数据小的扰动不敏感的性质。
例如,某些数值方法在计算过程中可能会放大舍入误差,导致结果不可靠,这样的方法就被认为是数值不稳定的。
三、计算题(每题15分,共30分)1. 给定线性方程组:\[\begin{align*}x + 2y - z &= 4 \\3x - y + 2z &= 1 \\-x + y + z &= 2\end{align*}\]使用高斯消元法求解该方程组,并给出解。
答案:首先将增广矩阵转换为上三角形式,然后回代求解,得到\( x = 1, y = 2, z = 1 \)。
2. 给定函数 \( f(x) = x^2 - 3x + 2 \),使用拉格朗日插值法在\( x = 0, 1, 2 \) 处插值,并求出插值多项式。
数值分析试题及答案一、选择题(每题2分,共20分)1. 以下哪个算法是数值分析中用于求解线性方程组的直接方法?A. 牛顿法B. 高斯消元法C. 梯度下降法D. 蒙特卡洛方法答案:B2. 插值法中,拉格朗日插值法和牛顿插值法的共同点是:A. 都是多项式插值B. 都使用差商C. 都只适用于等距节点D. 都需要预先知道所有数据点答案:A3. 在数值积分中,辛普森(Simpson)公式比梯形公式的误差:A. 更大B. 更小C. 相同D. 无法比较答案:B4. 以下哪个是数值稳定性分析中常用的方法?A. 条件数B. 收敛性C. 收敛速度D. 误差分析答案:A5. 在求解常微分方程的数值解时,欧拉方法属于:A. 单步法B. 多步法C. 隐式方法D. 显式方法答案:A6. 以下哪个是数值分析中求解非线性方程的迭代方法?A. 高斯-约当消元法B. 牛顿-拉弗森方法C. 雅可比迭代法D. 高斯-赛德尔迭代法答案:B7. 线性插值公式中,如果给定两个点\( (x_0, y_0) \)和\( (x_1, y_1) \),插值多项式是:A. \( y = y_0 + \frac{y_1 - y_0}{x_1 - x_0}(x - x_0) \)B. \( y = y_0 + \frac{y_1 - y_0}{x_0 - x_1}(x - x_0) \)C. \( y = y_0 + \frac{x - x_0}{x_1 - x_0}(y_1 - y_0) \)D. \( y = y_1 + \frac{x_1 - x}{x_1 - x_0}(y_0 - y_1) \)答案:C8. 以下哪个是数值分析中用于求解特征值问题的算法?A. 幂法B. 共轭梯度法C. 牛顿法D. 欧拉法答案:A9. 在数值微分中,使用有限差分法来近似导数时,中心差分法的误差:A. 与步长成正比B. 与步长的平方成正比C. 与步长的立方成正比D. 与步长的四次方成正比答案:B10. 以下哪个是数值分析中用于求解线性最小二乘问题的算法?A. 梯度下降法B. 牛顿法C. 奇异值分解法D. 共轭梯度法答案:C二、简答题(每题10分,共30分)1. 简述数值分析中病态问题的特点及其对算法的影响。
线封密三峡大学试卷班级姓名学号2011年春季学期《数值分析》课程考试试卷( A 卷)答案及评分标准注意:1、本试卷共3页;2、考试时间:120 分钟;3、姓名、学号必须写在指定地方;一、(16分)填空题1. 已知1125A ⎡⎤=⎢⎥⎣⎦,则1A 6= (1分),∞A 7= . (1分)2.迭代过程),1,0)((1 ==+n x x n n ϕ收敛的一个充分条件是迭代函数)(x ϕ满足1|)(|<'x ϕ. (2分)3. 设),,2,1,0(,,53)(2==+=k kh x x x f k 则差商0],,,[321=+++n n n n x x x x f .(2分)4. 设)(x f 可微,求方程)(x f x =根的牛顿迭代格式是.2,1,0,)(1)(1='---=+k x f x f x x x k k k k k (2分)5. 用二分法求方程01)(3=-+=x x x f 在区间]1,0[内的根,迭代进行二步后根所在区间为]75.0,5.0[.(2分)6.为尽量避免有效数字的严重损失,当1>>x 时,应将表达式x x -+1改写为xx ++11以保证计算结果比较精确.(2分)7. 将2111A ⎛⎫= ⎪⎝⎭作Doolittle 分解(即LU 分解),则100.51L ⎛⎫= ⎪⎝⎭(2分),2100.5U ⎛⎫= ⎪⎝⎭(2分)二、(10分)用最小二乘法解下列超定线性方程组:⎪⎩⎪⎨⎧=-=+=+2724212121x x x x x x 解:23222121,e e e x x ++=)(ϕ221221221)2()72()4(--+-++-+=x x x x x x由 ⎪⎪⎩⎪⎪⎨⎧=-+=∂∂=-+=∂∂0)1662(20)1323(2212211x x x x x x ϕϕ(8分)得法方程组 ⎩⎨⎧=+=+166213232121x x x x 7231=⇒x , 7112=x所以最小二乘解为: 7231=x 7112=x . (10分)三、(10分)已知)(x f 的函数值如下表25.15.001)(15.005.01---x f x用复合梯形公式和复合Simpson 公式求dx x f ⎰-11)(的近似值.解 用复合梯形公式,小区间数4=n ,步长5.0)]1(1[41=--⨯=h )]1())5.0()0()5.0((2)1([24f f f f f hT +++-+-=.线封密三峡大学试卷班级姓名学号25.1]2)5.15.00(21[25.0=++++-=(5分) 用复合Simpson. 小区间数2=n ,步长1)]1(1[21=--⨯=h)]1())5.0()5.0((4)0(2)1([62f f f f f hS ++-+⨯+-=33.168]2)5.10(45.021[61≈=+++⨯+-= (10分)四、(12分)初值问题 ⎩⎨⎧=>+='0)0(0,y x b ax y有精确解 bx ax x y +=221)(, 试证明: 用Euler 法以h 为步长所得近似解n y 的整体截断误差为n n n n ahx y x y 21)(=-=ε证: Euler 公式为:),(111---+=n n n n y x hf y y代入b ax y x f +=),(得:)(11b ax h y y n n n ++=-- 由0)0(0==y y 得:bh b ax h y y =++=)(001; 11122)(ahx bh b ax h y y +=++= )(3)(21223x x ah bh b ax h y y ++=++=……)()(12111---++++=++=n n n n x x x ah nbh b ax h y y (10分)因nh x n =,于是 )]1(21[2-++++=n ah bx y n n 2)1(2nn ah bx n -+==n n n bx x x a+-12∴n n n y x y -=)(ε)2(2112n n n n n bx x x abx ax +-+=-=n n n x x x a )(21--=n hx a 2 =221anh (12分)五、(10分) 取节点1,010==x x ,写出x e x y -=)(的一次插值多项式),(1x L 并估计插值误差.解: 建立Lagrange 公式为()=x L 110100101y x x x x y x x x x --+--=10101101-⨯--+⨯--=e x x x e x 11-+-=.(8分)())1)(0(!2)()()(11--''=-=x x y x L x y x R ξ )10(<<ξ ()811)0(max 2110≤--≤≤≤x x x(10分)六、(10分) 在区间]3,2[上利用压缩映像原理验证迭代格式,1,0,4ln 1==+k x x k k 的敛散性.解 : 在]3,2[上, 由迭代格式 ,1,0,4ln 1==+k x x k k , 知=)(x ϕx 4ln .因∈x ]3,2[时,]3,2[]12ln ,8[ln )]3(),2([)(⊂=∈ϕϕϕx (5分) 又1|1||)(|<='xx ϕ,故由压缩映像原理知对任意]3,2[0∈x 有收敛的迭代公式),1,0(,4ln 1 ==+k x x k k (10分)线封密三峡大学试卷班级姓名学号七、(10分)试构造方程组⎩⎨⎧=+=+423322121x x x x 收敛的Jacobi 迭代格式和Seidel Gauss -迭代格式,并说明其收敛的理由. 解:将原方程组调整次序如下:⎩⎨⎧=+=+324232121x x x x 调整次序后的方程组为主对角线严格占优方程组,故可保证建立的J 迭代格式和GS 迭代格式一定收敛.收敛的J 迭代格式为:⎪⎪⎩⎪⎪⎨⎧-=-=++)3(21)24(31)(1)1(2)(2)1(1k k k k x x x x .,1,0 =k (5分)收敛的GS 迭代格式为:⎪⎪⎩⎪⎪⎨⎧-=-=+++)3(21)24(31)1(1)1(2)(2)1(1k k k k x x x x .,1,0 =k (10分)八、(12分)已知43,21,41210===x x x 1)推导以这3个点作为求积节点在[0,1]上的插值型求积公式;2)指明求积公式所具有的代数精度.解:1)过这3个点的插值多项式)())(())(()())(())(()(121012002010212x f x x x x x x x x x f x x x x x x x x x p ----+----=+)())(())((2021201x f x x x x x x x x ----⎰⎰=∑=≈∴)()()(221010k k k x f A dx x p dx x f ,其中: ⎰⎰=----=----=32)4341)(2141()43)(21())(())((10201021100dx x x dx x x x x x x x x A ⎰⎰-=----=----=31)4321)(4121()43)(41())(())((10210120101dx x x dx x x x x x x x x A ⎰⎰=----=----=322143)(4143()21)(41())(())((10120210102dx x x dx x x x x x x x x A ∴所求的插值型求积公式为:⎰+-≈)]43(2)21()41(2[31)(10f f f dx x f (10分) 2)上述求积公式是由二次插值函数积分而来的,故至少具有2次代数精度,再将43,)(x x x f =代入上述求积公式,有:⎰+-==]43(2)21()41(2[3141333310dx x ⎰+-≠=])43(2)21(41(2[3151444410dx x 故上述求积公式具有3次代数精度. (12分)九、(10分)学完《数值分析》这门课程后,请你简述一下“插值、逼近、拟合”三者的区别和联系.。
二
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对称正定,,故也对称正定。
设迭代矩阵的特征值为,为相应的特征向量,则
与做内积,有
因正定,故,从而,格式收敛。
(2)设收敛到,则即,
即收敛到的解。
三
1 设且.求证:
证明以和为插值节点建立的不超过一次的插值多项式
应用插值余项公式有
2求一个次数不高于4次的多项式,使它满足.
解法一(待定参数法)满足的Hermite插值多项式为
设,令得
于是
解法二(带重节点的Newton插值法)建立如下差商表
这样可以写出Newton插值公式
3设,在上取,按等距节点求分段线性插值函数,计算各节点间中点处与的值,并估计误差.
解步长,.在区间上的线性插值函数
分段线性插值函数定义如下
,
各区间中点的函数值及插值函数值如表所示
估计误差:在区间上
而
令得的驻点,于是
故有结论
,
右端与无关,于是有
,
四
1确定参数和,使得积分取得最小值,并计算该最小值.
解本题实质上是求,关于权函数的二次最佳平方逼近多项式.选切比雪夫多项式为基函数进行计算:
于是得的二次最佳平方逼近多项式
进而有参数.
最小值就是平方误差:
2对彗星1968Tentax的移动在某个极坐标系下有如表所示的观察数据.
假设忽略来自行星的干扰,坐标应满足
其中为参数,为离心率,试用最小二乘法拟合和,并给出平方误差.
解由于关于参数和是非线性的,变形为,这样有下表的数据.
记,得拟合模型.
求解法方程组
得
进而有,拟合方程为
平方误差为
3 求函数在指定区间上关于的最佳平方逼近多项式.
解对做线性变换,即
利用勒让德正交多项式为基建立的一次最佳平方逼近多项式
的最佳平方逼近为
五
1确定中的待定参数,使其代数精确度尽量高,并指明求积公式所具有的代数精确度。
解令,代入公式两端并令其相等,得
解得
令,得
令,得故求积公式具有3次代数精确度。
2计算积分,若复化梯形公式,问区间应分多少等份才能使截断误差不超过?若改用复化辛普森公式,要达到同样精确度,区间应分多少等份?
解由于,故对复化梯形公式,要求
即。
取,即将区间分为213等份时,用复化梯形公式计算,截断误差
不超过。
用复化辛普森公式,要求
即。
取,即将区间等分为8等份时,复化辛普森公式可达精度。
3确定求积公式
中的系数,使代数精确度尽量高,并给出的表达式。
公式中。
解这是一个带权的且带导数值的求积公式。
为了积分方便,设该求积公式对
准确成立,得
化简得
解得
又因为
故求积公式
具有3次代数精确度。
下面估计求积公式的余项。
设在上三次插值多项式为,即满足。
因前述求积公式具有3次代数精确度,故它对于是准确成立的,且
因此有
注意到在上不变号,故余项
4已知。
(1)推导以这3个点作为求积节点在上的插值型求积公式;
(2)指明求积公式所具有的代数精确度;
(3)用所求公式计算。
解(1)过这3个点的插值多项式
故
其中
故所求的插值型求积公式为
(2)上述求积公式是由二次插值函数积分而来,故至少具有2次代数精确度。
再将代入上述求积公式,有
故上述求积公式具有3次代数精确度。
(3)
由于该求积公式具有3次代数精确度,从而为的精确度。
5设。
求证:(1)
(2)
(提示:直接使用泰勒展开即可得证)
七
1 对于迭代函数,试讨论:
(1)当为何值时,产生的序列收敛于;
(2)取何值时收敛最快?
(3)分别取计算的不动点,要求
解(1),根据定理7.3,当,亦即
时迭代收敛。
(2)由定理7.4知,当,即时迭代至少是二阶收敛的,收敛最快。
(3)分别取,并取,迭代计算结果如表7-4所示。
此时都达到。
事实上,
2(牛顿迭代法收敛性定理)设在上具有二阶连续导数,且满足条件
(1);
(2)在上;
(3)满足。
则由牛顿迭代法产生的序列单调收敛于在内的唯一实根,并且是平方收敛的。
证明因在上连续,由条件(1)知,方程在内有根。
又由条件(2)知在上恒正或恒负,所以在上严格单调,因而是在内的唯一实根。
条件(1)(2)共有四种情形:
(1)
(2)
(3)
(4)
仅就(1)进行定理证明,其余三种情况证明方法类似。
由,可知,再由知单增且。
又由牛顿迭代法知
由台劳展开的
其中介于,之间。
利用得
由以及前面证明的有
一般地,设,则必有且
再由台劳
及,得
根据归纳法原理数列单调下降有下界,因此有极限。
设,对迭代式两端取的极限,并利用,的连续性知即。
由上述证明知,有关系式,即对于单根,牛顿迭代法是平方收敛的。
3给定函数,对于一切,存在且,证明对于范围内的任意定数,迭代过程均收敛于的根.
证明由于,为单调增函数,故方程的根是唯一的(假定方程有根)。
迭代函数,。
由及得
,故。
由此可得
即
4设,试确定函数和,使求解且以为迭代函数的迭代法至少三阶收敛。
解要求三阶收敛到的根,根据定理7.4,应有于是由
得
故取
即迭代至少三阶收敛。