文档库 最新最全的文档下载
当前位置:文档库 › 西安交通大学2004年计算方法A试题

西安交通大学2004年计算方法A试题

西安交通大学2004年计算方法A试题
西安交通大学2004年计算方法A试题

统计西安交大期末考试试题(含答案)

西安交大统计学考试试卷 一、单项选择题(每小题2 分,共20 分) 1.在企业统计中,下列统计标志中属于数量标志的是(C) A、文化程度 B、职业 C、月工资 D、行业 2.下列属于相对数的综合指标有(B ) A、国民收入 B、人均国民收入 C、国内生产净值 D、设备台数 3.有三个企业的年利润额分别是5000 万元、8000 万元和3900 万元,则这句话中有(B)个变量? A、0 个 B、两个 C、1 个 D、3 个 4.下列变量中属于连续型变量的是(A ) A、身高 B、产品件数 C、企业人数 D、产品品种 5.下列各项中,属于时点指标的有(A ) A、库存额 B、总收入 C、平均收入 D、人均收入 6.典型调查是(B )确定调查单位的 A、随机 B、主观 C、随意 D 盲目 7.总体标准差未知时总体均值的假设检验要用到(A ): A、Z 统计量 B、t 统计量 C、统计量 D、X 统计量 8.把样本总体中全部单位数的集合称为(A ) A、样本 B、小总体 C、样本容量 D、总体容量 9.概率的取值范围是p(D ) A、大于1 B、大于-1 C、小于1 D、在0 与1 之间 10.算术平均数的离差之和等于(A ) A、零 B、1 C、-1 D、2 二、多项选择题(每小题2 分,共10 分。每题全部答对才给分,否则不计分) 1.数据的计量尺度包括(ABCD ): A、定类尺度 B、定序尺度 C、定距尺度 D、定比尺度 E、测量尺度 2.下列属于连续型变量的有(BE ): A、工人人数 B、商品销售额 C、商品库存额 D、商品库存量 E、总产值 3.测量变量离中趋势的指标有(ABE ) A、极差 B、平均差 C、几何平均数 D、众数 E、标准差 4.在工业企业的设备调查中(BDE ) A、工业企业是调查对象 B、工业企业的所有设备是调查对象 C、每台设备是 填报单位D、每台设备是调查单位E、每个工业企业是填报单位 5.下列平均数中,容易受数列中极端值影响的平均数有(ABC ) A、算术平均数 B、调和平均数 C、几何平均数 D、中位数 E、众数 三、判断题(在正确答案后写“对”,在错误答案后写“错”。每小题1 分,共10 分) 1、“性别”是品质标志。(对)

《数值计算方法》试题集及答案(1-6) 2

《计算方法》期中复习试题 一、填空题: 1、已知3.1)3(,2.1)2(,0.1)1(===f f f ,则用辛普生(辛卜生)公式计算求得 ?≈3 1 _________ )(dx x f ,用三点式求得≈')1(f 。 答案:2.367,0.25 2、1)3(,2)2(,1)1(==-=f f f ,则过这三点的二次插值多项式中2 x 的系数为 , 拉格朗日插值多项式为 。 答案:-1, )2)(1(21 )3)(1(2)3)(2(21)(2--------= x x x x x x x L 3、近似值*0.231x =关于真值229.0=x 有( 2 )位有效数字; 4、设)(x f 可微,求方程)(x f x =的牛顿迭代格式是( ); 答案 )(1)(1n n n n n x f x f x x x '--- =+ 5、对1)(3 ++=x x x f ,差商=]3,2,1,0[f ( 1 ),=]4,3,2,1,0[f ( 0 ); 6、计算方法主要研究( 截断 )误差和( 舍入 )误差; 7、用二分法求非线性方程 f (x )=0在区间(a ,b )内的根时,二分n 次后的误差限为 ( 1 2+-n a b ); 8、已知f (1)=2,f (2)=3,f (4)=5.9,则二次Newton 插值多项式中x 2系数为( 0.15 ); 11、 两点式高斯型求积公式?1 d )(x x f ≈( ?++-≈1 )] 321 3()3213([21d )(f f x x f ),代数精 度为( 5 ); 12、 为了使计算 32)1(6 )1(41310-- -+-+ =x x x y 的乘除法次数尽量地少,应将该表 达式改写为 11 ,))64(3(10-= -++=x t t t t y ,为了减少舍入误差,应将表达式 19992001-改写为 199920012 + 。

《数值计算方法》试题集及答案

《数值计算方法》复习试题 一、填空题: 1、????? ?????----=410141014A ,则A 的LU 分解为 A ??? ?????????=? ?????????? ?。 答案: ?? ????????--??????????--=1556141501 4115401411A 2、已知3.1)3(,2.1)2(,0.1)1(===f f f ,则用辛普生(辛卜生)公式计算求得 ?≈3 1 _________ )(dx x f ,用三点式求得≈')1(f 。 答案:, 3、1)3(,2)2(,1)1(==-=f f f ,则过这三点的二次插值多项式中2 x 的系数为 , 拉格朗日插值多项式为 。 答案:-1, )2)(1(21 )3)(1(2)3)(2(21)(2--------= x x x x x x x L 4、近似值*0.231x =关于真值229.0=x 有( 2 )位有效数字; 5、设)(x f 可微,求方程)(x f x =的牛顿迭代格式是( ); ( 答案 )(1)(1n n n n n x f x f x x x '--- =+ 6、对1)(3 ++=x x x f ,差商=]3,2,1,0[f ( 1 ),=]4,3,2,1,0[f ( 0 ); 7、计算方法主要研究( 截断 )误差和( 舍入 )误差; 8、用二分法求非线性方程 f (x )=0在区间(a ,b )内的根时,二分n 次后的误差限为 ( 1 2+-n a b ); 9、求解一阶常微分方程初值问题y '= f (x ,y ),y (x 0)=y 0的改进的欧拉公式为

( )] ,(),([2111+++++=n n n n n n y x f y x f h y y ); 10、已知f (1)=2,f (2)=3,f (4)=,则二次Newton 插值多项式中x 2系数为( ); 11、 两点式高斯型求积公式?1 d )(x x f ≈( ?++-≈1 )] 321 3()3213([21d )(f f x x f ),代数精 度为( 5 ); 12、 解线性方程组A x =b 的高斯顺序消元法满足的充要条件为(A 的各阶顺序主子式均 不为零)。 13、 为了使计算 32)1(6 )1(41310-- -+-+ =x x x y 的乘除法次数尽量地少,应将该表 达式改写为 11 ,))64(3(10-= -++=x t t t t y ,为了减少舍入误差,应将表达式 19992001-改写为 199920012 + 。 14、 用二分法求方程01)(3 =-+=x x x f 在区间[0,1]内的根,进行一步后根的所在区间 为 ,1 ,进行两步后根的所在区间为 , 。 15、 、 16、 计算积分?1 5 .0d x x ,取4位有效数字。用梯形公式计算求得的近似值为 ,用辛卜 生公式计算求得的近似值为 ,梯形公式的代数精度为 1 ,辛卜生公式的代数精度为 3 。 17、 求解方程组?? ?=+=+042.01532121x x x x 的高斯—塞德尔迭代格式为 ?????-=-=+++20/3/)51()1(1)1(2)(2)1(1 k k k k x x x x ,该迭 代格式的迭代矩阵的谱半径)(M ρ= 121 。 18、 设46)2(,16)1(,0)0(===f f f ,则=)(1x l )2()(1--=x x x l ,)(x f 的二次牛顿 插值多项式为 )1(716)(2-+=x x x x N 。 19、 求积公式 ?∑=≈b a k n k k x f A x x f )(d )(0 的代数精度以( 高斯型 )求积公式为最高,具 有( 12+n )次代数精度。

西安交通大学计算方法B上机试题

1.计算以下和式:01421181 84858616n n S n n n n ∞ =?? =--- ?++++??∑ ,要求: (1)若保留11个有效数字,给出计算结果,并评价计算的算法; (2)若要保留30个有效数字,则又将如何进行计算。 (1)题目分析 该题是对无穷级数求和,因此在使用matlab 进行累加时需要一个累加的终止条件。这里令?? ? ??+-+-+-+= 681581482184161n n n n a n n ,则 ()()1.016 1 6855844864816114851384128698161 681581482184161148113811282984161111<< ? ??? ????? ??++++++???? ????? ??++++++=??? ????? ??+-+-+-+??? ????? ??+-+-+-+=+++n n n n n n n n n n n n n n n n a a n n n n n n 故近似取其误差为1+≈k a ε,并且有m -1m -111021 21 ?=?=≈+βεk a , (2)算法依据 使用matlab 编程时用digits 函数和vpa 函数来控制位数。 (3)Matlab 运行程序 %%保留11位有效数字 k1=11; s1=0;%用于存储这一步计算值 for n=0:50 a=(1/16^n)*(4/(8*n+1)-2/(8*n+4)-1/(8*n+5)-1/(8*n+6)); n1=n-1; if a<=0.5*10^(1-k1) break end end; for i=0:1:n1 t=(1/16^i)*(4/(8*i+1)-2/(8*i+4)-1/(8*i+5)-1/(8*i+6)); s1=s1+t; end s11=vpa(s1,k1); disp('保留11位有效数字的结果为:');disp(s11); disp('此时n 值为:');disp(n1); %%保留30位有效数字 clear all; k2=30;

(完整word版)西工大计算方法试题参考(完整版).docx

2002-2003 第一学期 一.计算及推导( 5*8) 1.已知 x* 3.141, x ,试确定 x * 近似 x 的有效数字位数。 * * * 0.100 * * * 2.有效数 x 1 3.105, x 2 0.001, x 3 1 x 2 3 ,试确定 x x 的相对误差限。 3.已知 f ( x) 0.5 x 3 0.1x 2 ,试计算差商 f 0,1,2,3 4.给出拟合三点 A (0,1), B (1,0) 和 C (1,1) 的直线方程。 5.推导中矩形求积公式 b (b a) f ( a b ) 1 f '' ( )(b a)3 f (x)dx a 2 24 b n f (x)dx A i f ( x i ) a 6.试证明插值型求积公式 i 0 的代数精确度至少是 n 次。 7.已知非线性方程 x f (x) 在区间 a, b 内有一实根,试写出该实根的牛顿迭代 公式。 8.用三角分解法求解线性方程组 1 2 1 x 1 0 2 2 3 x 2 3 1 3 0 x 3 2 二.给出下列函数值表 0.4 0.5 0.6 0.7 0.8 x i 0.38942 0.47943 0.56464 0.64422 0.71736 f ( x i ) 要用二次插值多项式计算 f (0.63891) 的近似值,试选择合适的插值节点进行计 算,并说明所选用节点依据。 (保留 5 位有效数字)(12 分) 三. 已知方程 x ln x 0 在 (0,1) 内有一实根 ( 1)给出求该实根的一个迭代公式,试之对任意的初始近似 x 0 (0,1) 迭代法都收 敛,并证明其收敛性。 ( 2) x 0 0.5 试用构造的迭代公式计算 的近似值 x n ,要求 x n x n 1 10 3 。 四. 设有方程组

《计算方法》期末考试试题

《计算方法》期末考试试题 一 选 择(每题3分,合计42分) 1. x* = 1.732050808,取x =1.7320,则x 具有 位有效数字。 A 、3 B 、4 C 、5 D 、6 2. 取7 3.13≈(三位有效数字),则 ≤-73.13 。 A 、30.510-? B 、20.510-? C 、10.510-? D 、0.5 3. 下面_ _不是数值计算应注意的问题。 A 、注意简化计算步骤,减少运算次数 B 、要避免相近两数相减 C 、要防止大数吃掉小数 D 、要尽量消灭误差 4. 对任意初始向量)0(x 及常向量g ,迭代过程g x B x k k +=+)() 1(收敛的充分必要条件是_ _。 A 、11< B B 、1<∞ B C 、1)(

《数值计算方法》试题集及答案

《数值计算方法》复习试题 一、填空题: 1、????? ?????----=410141014A ,则A 的LU 分解为 A ??? ?????????=? ?????????? ?。 答案: ?? ????????--??????????--=1556141501 4115401411A 3、1)3(,2)2(,1)1(==-=f f f ,则过这三点的二次插值多项式中2 x 的系数为 ,拉 格朗日插值多项式为 。 答案:-1, )2)(1(21 )3)(1(2)3)(2(21)(2--------= x x x x x x x L 4、近似值*0.231x =关于真值229.0=x 有( 2 )位有效数字; 5、设)(x f 可微,求方程)(x f x =的牛顿迭代格式就是( ); 答案 )(1)(1n n n n n x f x f x x x '--- =+ 6、对1)(3 ++=x x x f ,差商=]3,2,1,0[f ( 1 ),=]4,3,2,1,0[f ( 0 ); 7、计算方法主要研究( 截断 )误差与( 舍入 )误差; 8、用二分法求非线性方程 f (x )=0在区间(a ,b )内的根时,二分n 次后的误差限为 ( 1 2+-n a b ); 10、已知f (1)=2,f (2)=3,f (4)=5、9,则二次Newton 插值多项式中x 2系数为( 0、15 ); 11、 解线性方程组A x =b 的高斯顺序消元法满足的充要条件为(A 的各阶顺序主子式均 不为零)。 12、 为了使计算 32)1(6 )1(41310-- -+-+ =x x x y 的乘除法次数尽量地少,应将该表 达式改写为 11 ,))64(3(10-= -++=x t t t t y ,为了减少舍入误差,应将表达式

计算方法试题

计算方法考试题(一) 满分70分 一、选择题:(共3道小题,第1小题4分,第2、3小题3分,共10分) 1、将A 分解为U L D A --=,其中),,(2211nn a a a diag D =,若对角阵D 非奇异(即),1,0n i a ii =≠,则b Ax =化为b D x U L D x 1 1)(--++=(1) 若记b D f U L D B 111 1),(--=+= (2) 则方程组(1)的迭代形式可写作 ) 2,1,0(1 )(1)1( =+=+k f x B x k k (3) 则(2)、(3)称 【 】 (A)、雅可比迭代。(B)、高斯—塞德尔迭代 (C)、LU 分解 (D)、Cholesky 分解。 2、记*x x e k k -=,若0lim 1≠=+∞→c e e p k k k (其中p 为一正数)称序列}{k x 是 【 】 (A)、p 阶收敛; (B)、1阶收敛; (C)、矩阵的算子范数; (D)、p 阶条件数。 3、牛顿切线法的迭代公式为 【 】 (A)、 ) () (1k x f x f x x k k k '- =+ (B)、 )()())((111--+--- =k k k k k k k x f x f x x x f x x 1 )() ()1()()()(x x f x f x f k i k i k i ??+=+ (D)、 )() ()()1(k k k x f x x -=+ 二、填空题:(共2道小题,每个空格2分,共10分) 1、设0)0(f =,16)1(f =,46)2(f =,则一阶差商 ,二阶差商=]1,2,0[f ,)x (f 的二次牛顿 插值多项式为 2、 用二分法求方程 01x x )x (f 3 =-+=在区间]1,0[内的根,进行第一步后根所在的区间为 ,进行第二步后根所在的区间 为 。 三、计算题:(共7道小题,第1小题8分,其余每小题7分,共50分) 1、表中各*x 都是对准确值x 进行四舍五入得到的近似值。试分别指出试用抛物插值计算115的近似值,并估计截断误差。 3、确定系数101,,A A A -,使求积公式 ) ()0()()(101h f A f A h f A dx x f h h ++-≈? -- (1) 具有尽可能高的代数精度,并指出所得求积公式的代数精度。

计算方法模拟试题及答案

计算方法模拟试题 一、 单项选择题(每小题3分,共15分) 1.近似值210450.0?的误差限为( )。 A . 0.5 B. 0.05 C . 0.005 D. 0.0005. 2. 求积公式)2(3 1 )1(34)0(31)(2 0f f f dx x f ++≈ ?的代数精确度为( )。 A. 1 B. 2 C. 3 D. 4 3. 若实方阵A 满足( )时,则存在唯一单位下三角阵L 和上三角阵R ,使LR A =。 A. 0det ≠A B. 某个0 det ≠k A C. )1,1(0det -=≠n k A k D. ),,1(0det n k A k =≠ 4.已知?? ?? ? ?????=531221112A ,则=∞A ( )。 A. 4 B. 5 C. 6 D 9 5.当实方阵A 满足)2(,221>>-=i i λλλλ,则乘幂法计算公式1e =( )。 A. 1+k x B. k k x x 11λ++ C. k x D. k k x x 11λ-+ 二、填空题(每小题3分,共15分) 1. 14159.3=π,具有4位有效数字的近似值为 。 2. 已知近似值21,x x ,则=-?)(21x x 。 3.已知1)(2-=x x f ,则差商=]3,2,1[f 。 4.雅可比法是求实对称阵 的一种变换方法。

5.改进欧拉法的公式为 。 三、计算题(每小题12分 ,共60分) 1. 求矛盾方程组; ??? ??=-=+=+2 42321 2121x x x x x x 的最小二乘解。 2.用列主元法解方程组 ??? ??=++=++=++4 26453426352321 321321x x x x x x x x x 3.已知方程组 ???? ? ?????=????????????????????----131********x x x a a a a (1) 写出雅可比法迭代公式; (2) 证明2

2016华工计算机计算方法(数值分析)考试试卷_共4页

考完试了,顺便把记得的题目背下来,应该都齐全了。我印象中也就只有这些题,题 目中的数字应该是对的,我也验证过,不过也不一定保证是对的,也有可能我也算错了。 还有就是试卷上面的题目可能没有我说的这么短,但是我也不能全把文字背下来,大概意 思就是这样吧。每个部分的题目的顺序可能不是这样,但总体就是这四大块。至于每道题 目的分值,我记得的就写出来了,有些题目没注意。我题目后面写的结果都是我考试时算 出来的,考完了也懒得验证了,可能不一定对,自己把握吧,仅供参考。 华南理工大学2016计算机计算方法(数值分析)考试试卷 一填空题(16分) 1.(6分)X* = 3.14,准确值x = 3.141592,求绝对误差e(x*) = ,相对误差e r(x*) = ,有效数位是。 2.(4分)当插值函数的n越大时,会出现龙格现象,为解决这个问题,分段函数不一个 不错的办法,请写出分段线性插值、分段三次Hermite插值和三次样条插值各自的特点。 3.(3分)已知x和y相近,将lgx – lgy变换成可以使其计算结果更准确。 4.(3分)已知2x3 – 3x2 +2 = 0,求牛顿迭代法的迭代式子。 解题思路:1. 这里的绝对误差和相对误差是没有加绝对值的,而且要注意是用哪个数减去哪个数得到的值,正负号会不一样;2. 可以从它们函数的连续性方面来说明;3. 只要满足课本所说的那几个要求就可以;这个记得迭代公式就可以直接写,记不住可以自己推导, 就是用泰勒展开式来近似求值得到的迭代公式。 我最终的结果是: 1.-0.001592 -0.000507 3 2.分段线性插值保证了插值函数的连续性,但是插值函数的一次导数不一定连续; 分段三次Hermite既保证了插值函数的连续性,也保证了其一次导数的连续性; 三次样条插值保证了插值函数及其一次导数和二次导数的连续性 3.lg(x/y) 4.x k+1 = x k – (2x3 – 3x2 +2)/(6x2 -6x) 二计算题(64分) 1.已知f(x) = x3 –x -1,用对分法求其在[0 , 2]区间内的根,误差要满小于0.2,需要对分多 少次?请写出最后的根结果。 解题思路:每次求区间的中值并计算其对应的函数值,然后再计算下一个区间中值及函数值,一直到两次区间中值的绝对值小于0.2为止。 我最终算得的对分次数是4,根的结果为11/8. 2.根据以下数据回答相应问题: x-2045 y51-31 (1)请根据以上数据构造Lagrange三次插值函数; (2)请列出差商表并写出Newton三次插值函数。 解题思路:(1) 直接按照书本的定义把公式列出来就可以了,这个要把公式记住了才行,不然也写不了;(2)差商表就是计算Newton三次插值函数过程中计算到的中间值及结

西安交通大学计算方法B大作业

计算方法上机报告 姓名: 学号: 班级:

目录 题目一------------------------------------------------------------------------------------------ - 4 - 1.1题目内容 ---------------------------------------------------------------------------- - 4 - 1.2算法思想 ---------------------------------------------------------------------------- - 4 - 1.3Matlab源程序----------------------------------------------------------------------- - 5 - 1.4计算结果及总结 ------------------------------------------------------------------- - 5 - 题目二------------------------------------------------------------------------------------------ - 7 - 2.1题目内容 ---------------------------------------------------------------------------- - 7 - 2.2算法思想 ---------------------------------------------------------------------------- - 7 - 2.3 Matlab源程序---------------------------------------------------------------------- - 8 - 2.4计算结果及总结 ------------------------------------------------------------------- - 9 - 题目三----------------------------------------------------------------------------------------- - 11 - 3.1题目内容 --------------------------------------------------------------------------- - 11 - 3.2算法思想 --------------------------------------------------------------------------- - 11 - 3.3Matlab源程序---------------------------------------------------------------------- - 13 - 3.4计算结果及总结 ------------------------------------------------------------------ - 14 - 题目四----------------------------------------------------------------------------------------- - 15 - 4.1题目内容 --------------------------------------------------------------------------- - 15 - 4.2算法思想 --------------------------------------------------------------------------- - 15 - 4.3Matlab源程序---------------------------------------------------------------------- - 15 - 4.4计算结果及总结 ------------------------------------------------------------------ - 16 - 题目五----------------------------------------------------------------------------------------- - 18 -

西交大计算方法上机报告

计算方法(B)实验报告 姓名: 学号: 学院: 专业:

实验一 三对角方程组Tx f =的求解 一、 实验目的 掌握三对角方程组Tx f =求解的方法。 二、 实验内容 求三对角方程组Tx f =的解,其中: 4 -1 -1 4 -1 -1 4 1 -1 4T ????????=?? ?? ???? , 3223f ?? ? ? ?= ? ? ??? 三、 算法组织 设系数矩阵为三对角矩阵 11222333111 b c a b c a b c a b c b n n n n T ---???????? =?????? ?????? 则方程组Tx f =称为三对角方程组。 设矩阵T 非奇异,T 可分解为T=LU ,其中L 为下三角矩阵,U 为单位上三角矩阵,记 1 1 212 313 1 1 1111 ,11n n n n n r l r l r L U l r l μμμμμ---???? ? ? ? ? ? ?== ? ? ? ? ? ? ? ? ? ?? ? ? ? 可先依次求出,L U 中的元素后,令Ux y =,先求解下三角方程组Ly f =得出 y ,再求解上三角方程组Ux y =。 追赶法的算法组织如下: 1.输入三对角矩阵T 和右端向量f ;

2.将Tx f =压缩为四个一维数组{}{}{}{}i i i i a b c d 、、、,{}{}{}i i i a b c 、、是T 的三对角线性方程组的三个对角,{}i d 是右端向量。将分解矩阵压缩为三个一维数组 {}{}{}i i i l r μ、、。 3.对T 做Crout 分解(也可以用Doolittle 分解)导出追赶法的计算步骤如下: 1111,b r c μ== for 2i n = 111, , ,i i i i i i i i i i i i i l a b a r r c y d l y μμ---==-==- end 4.回代求解x /n n n x y μ= for 11i n =- 1()/i i i i i x y c x μ+=- end 5. 停止,输出结果。 四、 MATLAB 程序 MATLAB 程序见附件1. 五、 结果及分析 实验结果为: (1.0000 1.0000 1.0000 1.0000)T x =

数值分析计算方法试题集及答案

数值分析复习试题 第一章 绪论 一. 填空题 1.* x 为精确值 x 的近似值;() **x f y =为一元函数 ()x f y =1的近似值; ()**,*y x f y =为二元函数()y x f y ,2=的近似值,请写出下面的公式:**e x x =-: *** r x x e x -= ()()()*'1**y f x x εε≈? ()() () ()'***1**r r x f x y x f x εε≈ ? ()()()() ()* *,**,*2**f x y f x y y x y x y εεε??≈?+??? ()()()()() ** * *,***,**222r f x y e x f x y e y y x y y y ε??≈ ?+??? 2、 计算方法实际计算时,对数据只能取有限位表示,这时所产生的误差叫 舍入误 差 。 3、 分别用2.718281,2.718282作数e 的近似值,则其有效数字分别有 6 位和 7 位;又取 1.73≈-21 1.73 10 2 ≤?。 4、 设121.216, 3.654x x ==均具有3位有效数字,则12x x 的相对误差限为 0.0055 。 5、 设121.216, 3.654x x ==均具有3位有效数字,则12x x +的误差限为 0.01 。 6、 已知近似值 2.4560A x =是由真值T x 经四舍五入得 到,则相对误差限为 0.0000204 . 7、 递推公式,??? ? ?0n n-1y =y =10y -1,n =1,2, 如果取0 1.41y ≈作计算,则计算到10y 时,误 差为 81 10 2 ?;这个计算公式数值稳定不稳定 不稳定 . 8、 精确值 14159265.3* =π,则近似值141.3*1=π和1415.3*2=π分别有 3

西交计算方法A上机大作业

计算方法A 上机大作业 1. 共轭梯度法求解线性方程组 算法原理:由定理3.4.1可知系数矩阵A 是对称正定矩阵的线性方程组Ax=b 的解与求解二次函数1()2 T T f x x Ax b x =-极小点具有等价性,所以可以利用共轭梯度法求解1()2 T T f x x Ax b x = -的极小点来达到求解Ax=b 的目的。 共轭梯度法在形式上具有迭代法的特征,在给定初始值情况下,根据迭代公式: (1)()()k k k k x x d α+=+ 产生的迭代序列(1)(2)(3)x x x ,,,... 在无舍入误差假定下,最多经过n 次迭代,就可求得()f x 的最小值,也就是方程Ax=b 的解。 首先导出最佳步长k α的计算式。 假设迭代点()k x 和搜索方向()k d 已经给定,便可以通过()()()() k k f x d φαα=+的极小化 ()()min ()()k k f x d φαα=+ 来求得,根据多元复合函数的求导法则得: ()()()'()()k k T k f x d d φαα=?+ 令'()0φα=,得到: ()() ()()k T k k k T k r d d Ad α=,其中()()k k r b Ax =- 然后确定搜索方向()k d 。给定初始向量(0)x 后,由于负梯度方向是函数下降最快的方向,故第一次迭代取搜索方向(0) (0)(0)(0)()d r f x b Ax ==-?=-。令 (1)(0)00x x d α=+ 其中(0)(0)0(0)(0) T T r d d Ad α=。第二次迭代时,从(1) x 出发的搜索方向不再取(1)r ,而是选取(1) (1)(0)0d r d β=+,使得(1)d 与(0)d 是关于矩阵A 的共轭向量,由此可 求得参数0β:

《计算方法》期末考试试题

一 选 择(每题3分,合计42分) 1. x* = ,取x =,则x 具有 位有效数字。 A 、3 B 、4 C 、5 D 、6 2. 取7 3.13≈(三位有效数字),则 ≤-73.13 。 A 、30.510-? B 、20.510-? C 、10.510-? D 、 3. 下面_ _不是数值计算应注意的问题。 A 、注意简化计算步骤,减少运算次数 B 、要避免相近两数相减 C 、要防止大数吃掉小数 D 、要尽量消灭误差 4. 对任意初始向量)0(x 及常向量g ,迭代过程g x B x k k +=+)() 1(收敛的充分必要条件是_ _。 A 、11< B B 、1<∞ B C 、1)(0 C 、f (a )f (b )<0 D 、f (a )f (b )>0 14. 由4个互异的数据点所构造的插值多项式的次数至多是____。

西安交通大学计算方法A实验报告

实验一 矩阵的分解 一、实验目的 掌握矩阵的分解原理和一般方法,学会利用矩阵分解直接求解线性方程组。 二、实验内容 求矩阵() 2020 =ij A α?的T LDL 分解与Cholesky 分解,其中 ,min(,),ij i i j i j i j α=?=? ≠? 。 三、问题分析 1. Cholesky 分解 Cholesky 分解是针对被分解矩阵为对称正定的情况给出的。 分解步骤如下: 11g =1111/y b g =,1111i i g g α= 2i n = ; DO 2j n = jj g = IF 0jj g < STOP ,JUMP TO (5) DO 1i j n =+ 1 1j ij ik kj k ij jj g g g g α-=??- ? ? ?=∑ ji ij g g = 1 1j i ik k k i jj b g y y g -=??- ? ? ?=∑ END DO END DO

2. T LDL 分解 T LDL 分解是针对Cholesky 分解中的开平方运算进行的改进。 分解步骤如下: 11i i r α=,1111/i i r r r =,11y b = 1i n = DO 2i n = DO j i n = 1 1i ij ij ik kj k r l r α-=??=- ??? ∑ /ji ij ii l r r = 1 1i i i ik k k y b l b -=??=- ??? ∑ END DO END DO 四、matlab 求解 分别写出T LDL 分解和Cholesky 分解的函数程序gaijinsqrt.m 和.cholesky m ,调用格 式如下: 1. [index,x,r]=gaijinsqrt(A,b) 参数说明: A 和b 分别是线性代数方程组Ax =b 的系数矩阵和右端向量;输出x 为解向量。 [index,x,g]=Cholesky(A,b) 参数说明: A 和b 分别是线性代数方程组Ax =b 的系数矩阵和右端向量;输出x 为解向量。 然后写出主程序2homework .m 如下: %生成矩阵A A=zeros(20,20); for i=1:20 for j=1:20 if i~=j if i>j A(i,j)=j; else A(i,j)=i; end

计算方法各章习题及答案

第二章 数值分析 2.1 已知多项式432()1p x x x x x =-+-+通过下列点: 试构造一多项式()q x 通过下列点: 答案:54313 ()()()3122 q 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]-近似函数x e 时,使用多少个节点能够保证误差不超过 61 102 -?. 答案:需要143个插值节点. 2.5 设被插值函数4()[,]f x C a b ∈,() 3()h H x 是()f x 关于等距节点 01n a x x x b =<<<=的分段三次艾尔米特插值多项式,步长b a h n -= .试估计() 3||()()||h f x H x ∞-. 答案:() 4 43||()()||384 h 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.1 22-1220 (sin )())0.989 310 710x p x dx δ=-=??. 3.2 确定参数,a b c 和,使得积分 2 1 2 1 (,,)[I a b c ax bx c -=++-?取最小值. 答案:810, 0, 33a b c ππ =- == 3.3 求多项式432()251f x x x x =+++在[1,1]-上的3次最佳一致逼近多项式 ()p x . 答案:()f x 的最佳一致逼近多项式为3 2 3 ()74 p x x x =++ . 3.4 用幂级数缩合方法,求() (11)x f x e x =-≤≤上的3次近似多项式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)x f 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 =? ,并与 精确值比较. 答案:计算结果如下表所示

相关文档 最新文档