数值分析试题(A)06.5.27.
- 格式:doc
- 大小:192.50 KB
- 文档页数:5
一、单项选择题(每小题3分,共15分) 1、用Simpson 公式求积分1401x dx +⎰的近似值为 ( ).A.2924 B.2429C.65D. 562、已知(1)0.401f =,且用梯形公式计算积分2()f x dx ⎰的近似值10.864T =,若将区间[0,2]二等分,则用递推公式计算近似值2T 等于( ). A.0.824 B.0.401 C.0.864 D. 0.8333、设3()32=+f x x ,则差商0123[,,,]f x x x x 等于( ).A.0B.9C.3D. 64的近似值的绝对误差小于0.01%,要取多少位有效数字( ). A.3 B.4 C.5 D. 25、用二分法求方程()0=f x 在区间[1,2]上的一个实根,若要求准确到小数 点后第四位,则至少二分区间多少次( ).A.12B.13C.14D. 15二、填空题(每小题4分,共40分)1、对于迭代函数2()=(3)ϕ+-x x a x ,要使迭代公式1=()ϕ+k k x x则a 的取值范围为 .2、假设按四舍五入的近似值为2.312,则该近似值的绝对误差限为 .3、迭代公式212(3)=,03++>+k k k k x x a x a x a收敛于α= (0)α>. 4、解方程4()530f x x x =+-=的牛顿迭代公式为 . 5、设()f x 在[1,1]-上具有2阶连续导数,[1,1]x ∀∈-,有1()2f x ''≤,则()f x 在[1,1]-上的线性插值函数1()L x 在点0处的误差限1(0)R ≤______.6、求解微分方程初值问题2(0)1'=-⎧⎨=⎩y xy yy ,0x 1≤≤的向前Euler 格式为 .7、设310131013A -⎛⎫⎪=-- ⎪ ⎪-⎝⎭,则A ∞= .8、用梯形公式计算积分112-⎰dx x 的近似值为 . 9、设12A 21+⎡⎤=⎢⎥⎣⎦a 可作Cholesky 分解,则a 的取值范围为 . 10、设(0)1,(0.5) 1.5,(1)2,(1.5) 2.5,(2) 3.4f f f f f =====,若1=h ,则用三点公式计算(1)'≈f .三、解答题(共45分) 1、给定数据用复化Simpson 公式计算1.381.30()f x dx ⎰的近似值,并估计误差,小数点后保留3位. (8分)2、用直接三角分解法求线性代数方程组⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡432631531321321x x x 的解. (8分) 3、求()λx ,使得迭代公式1()()λ+=+k k k k f x x x x 求方程2()31=+-f x x x 的根的相应迭代序列{}k x 具有平方收敛. (5分)4、已知数据试对数据用最小二乘法求出形如=+y x b的拟合曲线. (8分) 5、已知(2)8f -=,(0)4f =,(2)8=f ,试求二次拉格朗日插值多项式. (8分) 6、设矩阵A 如下,根据谱半径判断用Jacobi 迭代法求解方程组Ax b =的敛散性.(8分)1102111221012A ⎡⎤-⎢⎥⎢⎥⎢⎥=--⎢⎥⎢⎥⎢⎥-⎢⎥⎣⎦一、单项选择题(每小题3分,合计15分) 1、A 2、D 3、C 4、C 5、D 二、填空题(每小题3分,合计30分) 1、0<<a ; 2、31102-⨯; 3;4、4135345++-=-+k k k k k x x x x x ; 5、14; 6、1(2)+=+-n n n n n y y h x y y ; 7、5;8、34-; 9、3>a ;10、1.2;三、计算题(合计55分) 1、给定数据用复化Simpson 公式计算 1.381.30()f x dx ⎰的近似值,并估计误差,小数点后保留3位. (8分)解: 401024S [()4()()]6-=++x x f x f x f x ………… 1分 1.38 1.30(3.624 4.20 5.19)6-=+⨯+ 0.341= ………… 2分20422012234S [()4()()][()4()()]66--=+++++x x x xf x f x f x f x f x f x =0.342 ………… 6分2211[]15-≈-I S S S =-⨯40.6710 ………… 8分 2、用直接三角分解法求线性代数方程组⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡432631531321321x x x 的解. (8分) 解:设111213212223313233u u u 123100135l 100u u 136l l 100u ⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥=*⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦………… 1分 111=u ,212=u ,313=u ,121=l ,131=l 122=u ,223=u ,132=l133=u ,133=l …………6分所以⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=111011001L ,⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=100210321U …………7分 由b Ly =得Ty )1,1,2(=;由y Ux =得Tx )1,1,1(-=. ………… 8分3、求()λx ,使得迭代公式1()()λ+=+k k k k f x x x x 求方程2()31=+-f x x x 的根的相应迭代序列{}k x 具有平方收敛.(6分)解:要使迭代序列具有平方收敛,则()0ϕ'*=x ………… 2分 而()()()ϕλ=+f x x x x ,即 ………… 3分 2()()()()10()λλλ''**-**+=*f x x x f x x …………4分 而()0*=f x 则有()1()λ'*=-*f x x ………… 5分所以()()23λ'=-=--x f x x ………… 6分4、已知数据试对数据用最小二乘法求出形如=+ay x b的拟合曲线. (8分) 解:因为11=+b x y a a ,令0111,,,====b a a y x x a a y……2分 则有法方程01461061410⎛⎫⎛⎫⎛⎫= ⎪⎪ ⎪⎝⎭⎝⎭⎝⎭a a ……5分解出014,1==-a a ,则1,4=-=-a b ……7分 所以1=4-y x……8分5、已知(2)8f -=,(0)4f =,(2)8=f ,试求二次拉格朗日插值多项式. (7分)解:01()(2)8l x x x =- …………2分 211()(4)4l x x =-- …………4分21()(2)8l x x x =+ …………6分 2012()()(2)()(0)()(2)L x l x f l x f l x f =-++24=+x …………7分6、设矩阵A 如下,根据谱半径判断用Jacobi 迭代法求解方程组Ax b =的敛散性.(8分)1102111221012A ⎡⎤-⎢⎥⎢⎥⎢⎥=--⎢⎥⎢⎥⎢⎥-⎢⎥⎣⎦解:100010001D ⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦,00010021002L ⎡⎤⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎣⎦,10021002000U ⎡⎤⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦…………3分1100211()0221002J B D L U -⎡⎤⎢⎥⎢⎥⎢⎥=+=⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦…………5分 2102111()0222102J E B λλλλλλ⎡⎤-⎢⎥⎢⎥⎢⎥-=--=-=⎢⎥⎢⎥⎢⎥-⎢⎥⎣⎦…………6分()2J B ρ=…………7分 所以用Jacobi 迭代法求解方程组Ax b =收敛 …………8分。
一. 填空题(本大题共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) 迭代法和直接法是数值计算中常用的两种方法。
迭代法通过反复迭代逼近解,直到满足所需精度为止;而直接法则通过一系列代数运算直接得到解。
数值分析试卷及答案数值分析试卷一、选择题(共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. 数值解线性方程组的常用迭代法有_________方法和_________方法。
数值分析试题院系: 专业: 分数:姓名: 学号: 日期:2006.5.27 一、 填空题(每空2分,共20分) 1.设1221A -⎡⎤=⎢⎥-⎣⎦,则A 的奇异值1_____.σ= 2. 已知2()P x 是用极小化插值法得到的sin x 在[0,3]上的二次插值多项式,则2()P x 的 截断误差上界为2()sin ()R x x P x =-≤_________. 3. 设42()231f x x x =++和节点,0,1,2,2k k x k ==则015[,,,]________f x x x = 和40()_________f x ∆=.4.如下两种计算1e -近似值的方法中哪种方法能够提供较好的近似。
_____方法1: 19101!n en --=⎛⎫≈ ⎪⎝⎭∑ 方法2:19101(9)!n e n --=⎛⎫≈ ⎪-⎝⎭∑5. 已知α是非线性方程f (x )=0的二重根,试构造至少二阶收敛的迭代格式__________________.6.给出求解线性方程组1231231238892688x x x x x x x x x -++=-⎧⎪-+=⎨⎪-+-=⎩ 的收敛的Jacobi 迭代格式(分量形式)______________________及相应的迭代矩阵______________________。
7. 解线性方程组Ax=b 的简单迭代格式(1)()k k xB xg +=+收敛的充要条件是__________.8. 下面Matlab 程序所解决的数学问题为____________________. function x=fun(A,b)n=length(b);x=zeros(n,1); x(n)=b(n)/A(n,n);for i =n-1:-1:1x(i )=(b(i )-A(i ,i +1:n)* x(i +1:n))/A(i ,i);end二、(15分) 已知方程组Ax=b ,即12121.000122x x x x +=⎧⎨+=⎩有解x =(2,0)T,(1) 求()cond A ∞;(2) 求右端项有小扰动的方程组12121.00012.00012x x x x +=⎧⎨+=⎩的解x x +∆;(3) 计算b b∞∞∆和x x∞∞∆,结果说明了什么问题。
一、填空题(每题3分,共30分)1. 用1415.3近似π,有效位数为 ① 。
2. 若干个浮点数做连加运算,按 ② 安排运算时,计算误差小。
3. 对称正定矩阵⎪⎪⎪⎭⎫ ⎝⎛----=1911215412416A 做Cholesky 分解,得⎪⎪⎪⎭⎫ ⎝⎛-=b a L 3214,那么,=a ③ ,=b ④ 。
4. 用部分选主元的Doolittle 分解法分解矩阵⎪⎪⎪⎪⎪⎭⎫⎝⎛2103673285213234,经过第一轮分解后得到⎪⎪⎪⎪⎪⎭⎫⎝⎛2104367321852413234,在第二轮分解时,应选择第 ⑤ 行作为主元行。
5.以这三点为节点的二次Newton 插值多项式为 ⑦ 。
6. Cotes 系数)(n k C 只跟将积分区间等分的份数有关,而跟 ⑧ ,和 ⑨ 都无关。
7. 用Jacobi 迭代求解线性方程组⎪⎩⎪⎨⎧=++-=-+=+-9353258462321321321x x x x x x x x x 取初始值T x)0,0,0()0(=,则=)1(x ⑩ 。
二、设序列{}n y 满足关系式11-=n n y ny ,假设在求0y 时的误差为ε,求计算10y 的误差,并讨论计算的稳定性?(8分)三、用紧凑格式的Doolittle 分解法求解线性方程组(10分)⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=⎪⎪⎪⎪⎪⎭⎫ ⎝⎛⎪⎪⎪⎪⎪⎭⎫ ⎝⎛---35303311111066133154602134321x x x x四、已知数据表①求最小二乘拟合函数2210)(x a x a a x P ++=和拟合误差。
(保留5位有效数字)(10分)②求二次Lagrange 插值多项式)(2x L (7分)③比较)(x P 和)(2x L ,并对结果做出说明(5分)五、证明⎩⎨⎧∈+-+∈+++=].2,1[458]1,0[223)(2323x x x x x x x x x f 是以)34,2(),8,1(),2,0(为节点的三次样条插值函数。
数值分析试题及答案一、选择题(每题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. 简述数值分析中病态问题的特点及其对算法的影响。
期末考试试卷(A 卷)2007学年第二学期 考试科目: 数值分析 考试时间:120 分钟学号 姓名 年级专业一、判断题(每小题2分,共10分)1. 用计算机求1000100011n n=∑时,应按照n 从小到大的顺序相加。
( )2. 为了减少误差,进行计算。
( )3. 用数值微分公式中求导数值时,步长越小计算就越精确。
( )4. 采用龙格-库塔法求解常微分方程的初值问题时,公式阶数越高,数值解越精确。
( )5. 用迭代法解线性方程组时,迭代能否收敛与初始向量的选择、系数矩阵及其演变方式有关,与常数项无关。
( )二、填空题(每空2分,共36分)1. 已知数a 的有效数为0.01,则它的绝对误差限为________,相对误差限为_________.2. 设1010021,5,1301A x -⎡⎤⎡⎤⎢⎥⎢⎥=-=-⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦则1A =_____,2x =______,Ax ∞=_____.3. 已知53()245,f x x x x =+-则[1,1,0]f -= ,[3,2,1,1,2,3]f ---= .4. 为使求积公式11231()((0)f x dx A f A f A f -≈++⎰的代数精度尽量高,应使1A = ,2A = ,3A = ,此时公式具有 次的代数精度。
5. n 阶方阵A 的谱半径()A ρ与它的任意一种范数A 的关系是 .6. 用迭代法解线性方程组AX B =时,使迭代公式(1)()(0,1,2,)k k XMX N k +=+=K 产生的向量序列{}()k X收敛的充分必要条件是 .7. 使用消元法解线性方程组AX B =时,系数矩阵A 可以分解为下三角矩阵L 和上三角矩阵U 的乘积,即.A LU = 若采用高斯消元法解AX B =,其中4221A -⎡⎤=⎢⎥⎣⎦,则L =_______________,U =______________;若使用克劳特消元法解AX B =,则11u =____;若使用平方根方法解AX B =,则11l 与11u 的大小关系为_____(选填:>,<,=,不一定)。
线封密三峡大学试卷班级姓名学号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分)学完《数值分析》这门课程后,请你简述一下“插值、逼近、拟合”三者的区别和联系.。
数值分析试卷及答案**注意:以下是一份数值分析试卷及答案,试卷和答案分别按照题目和解答的格式排版,以确保整洁美观,语句通顺。
**---数值分析试卷一、选择题(每题2分,共20分)1. 数值分析是研究如何用计算机处理数值计算问题的一门学科。
以下哪个选项不是数值分析的应用领域?A. 金融风险评估B. 天气预测C. 数据挖掘D. 图像处理2. 在数值计算中,稳定性是指算法对于输入数据的微小扰动具有较好的性质。
以下哪个算法是稳定的?A. 高斯消元法B. 牛顿迭代法C. 不动点迭代法D. 雅可比迭代法二、填空题(每题3分,共30分)1. 下面关于插值多项式的说法中,不正确的是:一般情况下,插值多项式的次数等于插值点的个数减1。
2. 线性方程组中,如果系数矩阵A是奇异的,则该方程组可能无解或有无穷多解。
......三、解答题(共50分)1. 请给出用割线法求解非线性方程 f(x) = 0 的迭代格式,并选择合适的初始值进行计算。
解:割线法的迭代公式为:x_(k+1) = x_k - f(x_k) * (x_k - x_(k-1)) / (f(x_k) - f(x_(k-1)))选择初始值 x0 = 1,x1 = 2 进行计算:迭代1次得到:x2 = x1 - f(x1) * (x1 - x0) / (f(x1) - f(x0))迭代2次得到:x3 = x2 - f(x2) * (x2 - x1) / (f(x2) - f(x1))继续迭代直至满足精度要求。
2. 对于一个给定的线性方程组,高斯消元法可以用来求解其解空间中的向量。
请简要描述高斯消元法的基本思想并给出求解步骤。
高斯消元法的基本思想是通过一系列的行变换将线性方程组化为上三角形式,然后再通过回代求解方程组的未知数。
求解步骤如下:步骤1:将方程组表示为增广矩阵形式,即将系数矩阵和常数向量连接在一起。
步骤2:从第一行开始,选取第一个非零元素作为主元,然后通过行变换将其它行的该列元素消去。
数值分析试题院系: 专业: 分数:姓名: 学号: 日期:2006.5.27 一、 填空题(每空2分,共20分) 1.设1221A -⎡⎤=⎢⎥-⎣⎦,则A 的奇异值1_____.σ= 2. 已知2()P x 是用极小化插值法得到的sin x 在[0,3]上的二次插值多项式,则2()P x 的 截断误差上界为2()sin ()R x x P x =-≤_________. 3. 设42()231f x x x =++和节点,0,1,2,2k k x k ==则015[,,,]________f x x x = 和40()_________f x ∆=.4.如下两种计算1e -近似值的方法中哪种方法能够提供较好的近似。
_____方法1: 19101!n en --=⎛⎫≈ ⎪⎝⎭∑ 方法2:19101(9)!n e n --=⎛⎫≈ ⎪-⎝⎭∑5. 已知α是非线性方程f (x )=0的二重根,试构造至少二阶收敛的迭代格式__________________.6.给出求解线性方程组1231231238892688x x x x x x x x x -++=-⎧⎪-+=⎨⎪-+-=⎩ 的收敛的Jacobi 迭代格式(分量形式)______________________及相应的迭代矩阵______________________。
7. 解线性方程组Ax=b 的简单迭代格式(1)()k k xB xg +=+收敛的充要条件是__________.8. 下面Matlab 程序所解决的数学问题为____________________. function x=fun(A,b)n=length(b);x=zeros(n,1); x(n)=b(n)/A(n,n);for i =n-1:-1:1x(i )=(b(i )-A(i ,i +1:n)* x(i +1:n))/A(i ,i);end二、(15分) 已知方程组Ax=b ,即12121.000122x x x x +=⎧⎨+=⎩有解x =(2,0)T,(1) 求()cond A ∞;(2) 求右端项有小扰动的方程组12121.00012.00012x x x x +=⎧⎨+=⎩的解x x +∆;(3) 计算b b∞∞∆和x x∞∞∆,结果说明了什么问题。
三、(15分) 已知函数值表 2101201210iix y --⎧⎨⎩在函数空间{}21,H span x =中求最佳平方逼近多项式,并估计误差22δ。
(注:取小数点后四位)四、(15分) 已知函数值表 1234() 1.12.62.81.6i i x f x =⎧⎨=⎩用二次多项式计算x =0.26时函数的较好近似值,并估计误差.五、(15分)(1)求[0,1]区间上关于权函数()ln x x ρ=-的首项系数为1的正交多项式01(),()x x ϕϕ。
(2)构造带权()ln x x ρ=-的高斯型求积公式1110()()()x f x dx A f x ρ≈⎰(3) 导出此高斯型求积公式的截断误差。
六、(10分)已知近似数x =10的绝对误差限为0.05,试求函数()f x =的相对误差限.七、(10分) 用Householder方法求矩阵1113001A ⎡⎢=-⎢⎥⎥⎦的正交分解,即A=QR 。
数值分析答案一、 填空题(每空2分,共20分)1. 3 .2. 90.1464= 3. 015[,,,]0f x x x = 和40()3f x ∆=4. 方法25. 12()'()k k k k f x x x f x +=-6。
Jacobi 迭代格式(1)()()123(1)()()213(1)()()312(62)/9(8)/8(8)/8k k k k k k k k k x x x x x x x x x +++⎧=+-⎪=++⎨⎪=-+-⎩迭代矩阵210991108811088J B -⎡⎤⎢⎥⎢⎥=⎢⎥⎢⎥-⎢⎥⎣⎦7. ()1B ρ< 8.解上三角形方程组Ax=b二、 (15分)(1)4414414410 1.0001101010() 2.0001(2.000110)410Acond A AA--∞∞∞⎡⎤-⨯=⎢⎥-⎣⎦==⨯⨯≈⨯(2) []11Tx x +∆=(3)[][][][][]22,0.00010()112011TTTTTb b x x x x =∆=∆=+∆-=-=-0.005%b b∞∞∆=和50%x x∞∞∆=虽然方程组右端项扰动的相对误差仅为0.005%,然而此小扰动引起解的相对误差却高达50%,这是由于”系数矩阵的条件数比较大,方程组是病态的”,从而导致上述结果. 三、(15分)212()1,()x x x ϕϕ==222122210(2)11(1),,,120111102Y ⎡⎤-⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥-⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥Φ=Φ==⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦2222158510435,1034237583() 1.65720.4286357(,)(,)(,)5836420.2286357a ab b s x x xY Y a Y b Y δ⎡⎤⎡⎤⎡⎤⎡⎤⎡⎤⎢⎥==⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦⎣⎦-⎣⎦=-=-=-Φ-Φ=-⨯+⨯=四、(15分)(1)建立如下差商表()1 1.12 2.6 1.53 2.80.20.6541.61.20.7i i x f x ---一阶差商二阶差商22() 1.1 1.5(1)0.65(1)(2)(0.26) 1.1 1.50.740.650.74 1.740.84694N x x x x N =+----=+⨯-⨯⨯=-22() 2.60.2(2)0.7(2)(3)(0.26) 2.60.2 1.740.7 1.74 2.74 1.08523N x x x x N =+----=+⨯-⨯⨯=-(2)22220.261(0.26)(0.26)((0.26)(0.26))140.261(0.84694 1.08532)0.058814R f N N N -=-≈---=-+=-22220.264(0.26)(0.26)((0.26)(0.26))140.264(0.84694 1.08532)0.297214R f N N N -=-≈---=-+=-五、 (15分) (1)由首1正交多项式的构造公式,可得0()1x ϕ=,01000(,())()()((),())x x x x x x x ϕϕϕϕϕ=-11121001100011ln 1,ln ,ln 49ln (,())1/411,()((),())144ln xdx x xdx x xdx x xdxx x x x x x xdxϕϕϕϕ=-=-=--====--⎰⎰⎰⎰⎰,(2) 11110011,ln ,()()()44x A xdx x f x dx f ρ==-≈⎰⎰(3) Gauss 型求积公式的截断误差为(2)(2)1122(2)(2)12(2)()1()1()(ln )()(ln )()2!42!4()11()1117(ln ln ln )()()2!2162!9816288ffR f x x dx x x dxffx x x x x dx f ξηηηη=--=--=--+=-+=⎰⎰⎰六、(10分)()0.05,'()e x f x ≤=n=20()(())'()()1(())()()()10.05(())()0.000252010r r x e f x f x e x e f x e x f x f x nxe f x e x n x=≈==≈≤=⨯七、(10分)(1,1,,(2,0,0),(3,1,112293125121312612122311212251326100233TTTTT x y u x y uu H I u u H A RQ H ==-=-=⎡---⎢⎡⎡⎤⎢⎢⎢⎥=-=-=--⎢⎢⎥⎢⎥⎢⎢⎥⎢⎥⎣⎦⎢⎥⎢⎥⎣⎦--⎢⎥⎣⎦⎡---⎡⎤-⎢⎢⎢⎢-=-===--⎢⎢⎥---⎢⎥⎣⎦⎣⎢⎢⎢⎥⎢⎥⎦。