数值分析习题
- 格式:doc
- 大小:643.00 KB
- 文档页数:22
..数值分析试题集(试卷一)一( 10 分)已知 x 1* 1.3409 ,x 2* 1.0125 都是由四舍五入产生的近似值, 判断 x 1*x 2* 及 x 1* x 2*有几位有效数字。
二( 10 分)由下表求插值多项式x 01 2 y2 34 y1- 1三( 15 分)设 f ( x)C 4 [a,b] , H ( x )是满足下列条件的三次多项式H (a) f (a) , H (b) f (b) , H (c)f (c) , H (c) f (c)( a c b )求 f (x)H ( x) ,并证明之。
12四( 15 分)计算13 dx ,10 2。
x五( 15 分)在 [0,2]上取 x 0 0 , x 1 1 , x 22 ,用二种方法构造求积公式,并给出其公式的代数精度。
六( 10 分)证明改进的尢拉法的精度是 2 阶的。
七( 10 分)对模型 yy , 0 ,讨论改进的尢拉法的稳定性。
八( 15分)求方程 x 34x 2 7x 1 0 在 -1.2 附近的近似值,10 3。
-----------------------------------------------------------------------------------------------------------------------------(试卷二)一填空( 4*2 分)1 {k ( x) } k 0 是区间 [0, 1]上的权函数为( x) x 2 的最高项系数为 1 的正交多项式族,其中10 (x)1,则x0 ( x) dx ------------------- , 1 ( x) ------------------。
2 12 A,则 A1 4----------- ,( A) ----------------- 。
a 1 2 时, A 可作 LU 分解。
3 设 A,当 a 满足条件 ---------------- 14..4 设非线性方程 f ( x) (x33x23x1)( x 3) 0 ,其根 x1* 3 , x2*1,则求 x1* 的近似值时,二阶局部收敛的牛顿迭代公式是--------------------------- 。
第一章 绪论(12) 第二章 插值法(40-42)2、当2,1,1-=x 时,4,3,0)(-=x f ,求)(x f 的二次插值多项式。
[解]372365)1(34)23(21)12)(12()1)(1(4)21)(11()2)(1()3()21)(11()2)(1(0))(())(())(())(())(())(()(2221202102210120120102102-+=-++--=+-+-⨯+------⨯-+-+-+⨯=----+----+----=x x x x x x x x x x x x x x x x x x x y x x x x x x x x y x x x x x x x x y x L 。
3、给出x x f ln )(=的数值表用线性插值及二次插值计算54.0ln 的近似值。
X 0.4 0.5 0.6 0.7 0.8 x ln -0.916291 -0.693147 -0.510826 -0.357765 -0.223144[解]若取5.00=x ,6.01=x ,则693147.0)5.0()(00-===f x f y ,510826.0)6.0()(11-===f x f y ,则604752.182321.1)5.0(10826.5)6.0(93147.65.06.05.0510826.06.05.06.0693147.0)(010110101-=---=--⨯---⨯-=--+--=x x x x x x x x x y x x x x y x L ,从而6202186.0604752.19845334.0604752.154.082321.1)54.0(1-=-=-⨯=L 。
若取4.00=x ,5.01=x ,6.02=x ,则916291.0)4.0()(00-===f x f y ,693147.0)5.0()(11-===f x f y ,510826.0)6.0()(22-===f x f y ,则 217097.2068475.404115.2)2.09.0(5413.25)24.0(3147.69)3.01.1(81455.45)5.06.0)(4.06.0()5.0)(4.0()510826.0()6.05.0)(4.05.0()6.0)(4.0()693147.0()6.04.0)(5.04.0()6.0)(5.0(916291.0))(())(())(())(())(())(()(22221202102210120120102102-+-=+--+-⨯++-⨯-=----⨯-+----⨯-+----⨯-=----+----+----=x x x x x x x x x x x x x x x x x x x x x x y x x x x x x x x y x x x x x x x x y x L ,从而61531984.0217097.21969765.259519934.0217097.254.0068475.454.004115.2)54.0(22-=-+-=-⨯+⨯-=L补充题:1、令00=x ,11=x ,写出x e x y -=)(的一次插值多项式)(1x L ,并估计插值余项。
数值分析试题及答案汇总TTA standardization office【TTA 5AB- TTAK 08- TTA 2C】数值分析试题一、 填空题(2 0×2′) 1.⎥⎦⎤⎢⎣⎡-=⎥⎦⎤⎢⎣⎡-=32,1223X A 设x =是精确值x *=的近似值,则x 有 2 位有效数字。
2. 若f (x )=x 7-x 3+1,则f [20,21,22,23,24,25,26,27]= 1 , f [20,21,22,23,24,25,26,27,28]=0 。
3. 设,‖A ‖∞=___5 ____,‖X ‖∞=__ 3_____,‖AX ‖∞≤_15_ __。
4. 非线性方程f (x )=0的迭代函数x =(x )在有解区间满足 |’(x )| <1 ,则使用该迭代函数的迭代解法一定是局部收敛的。
5. 区间[a ,b ]上的三次样条插值函数S (x )在[a ,b ]上具有直到 2 阶的连续导数。
6. 当插值节点为等距分布时,若所求节点靠近首节点,应该选用等距节点下牛顿差商公式的前插公式 ,若所求节点靠近尾节点,应该选用等距节点下牛顿差商公式的 后插公式 ;如果要估计结果的舍入误差,应该选用插值公式中的 拉格朗日插值公式 。
7. 拉格朗日插值公式中f (x i )的系数a i (x )的特点是:=∑=ni i x a 0)( 1 ;所以当系数a i (x )满足 a i (x )>1 ,计算时不会放大f (x i )的误差。
8. 要使20的近似值的相对误差小于%,至少要取 4 位有效数字。
9. 对任意初始向量X (0)及任意向量g ,线性方程组的迭代公式x (k +1)=Bx (k )+g (k =0,1,…)收敛于方程组的精确解x *的充分必要条件是 (B)<1 。
10. 由下列数据所确定的插值多项式的次数最高是 5 。
11. 牛顿下山法的下山条件为 |f(xn+1)|<|f(xn)| 。
数值分析习题集及答案数值分析习题集及答案篇一:数值分析习题与答案第一章绪论习题一1.设x>0,x*的相对误差为δ,求f(x)=ln x的误差限。
解:求lnx 的误差极限就是求f(x)=lnx的误差限,由公式()有已知x*的相对误差,故即2.下列各数都是经过四舍五入得到的近似值,试指出它们有几位有效数字,并给出其误差限与相对误差限。
解:直接根据定义和式()()则得有5位有效数字,其误差限有2位有效数字,有5位有效数字,3.下列公式如何才比较准确?(1)(2),相对误差限满足,而解:要使计算较准确,主要是避免两相近数相减,故应变换所给公式。
(1)(2)4.近似数x*=,是 3位有数数字。
5.计算四个选项:取,利用:式计算误差最小。
第二、三章插值与函数逼近习题二、三 1.给定的数值表用线性插值与二次插值计算的近似值并估计误差限. 解:仍可使用n=1及n=2的Lagrange插值或Newton插值,并应用误差估计()。
线性插值时,用及两点,用Newton插值误差限,因,故二次插值时,用,,三点,作二次Newton 插值误差限,故2. 在-4≤x≤4上给出的等距节点函数表,若用二次,函数表的步长h插值法求的近似值,要使误差不超过应取多少? 解:用误差估计式(),令因得3. 若,求和.解:由均差与导数关系于是4. 若的值,这里p≤n+1.解:可知当而当P=n+1时于是得有互异,求,由均差对称性5.求证.解:解:只要按差分定义直接展开得6.已知的函数表求出三次Newton均差插值多项式,计算f()的近似值并用均差的余项表达式估计误差. 解:根据给定函数表构造均差表由式()当n=3时得Newton均差插值多项式N3(x)=+()+()() 由此可得f() N3()= 由余项表达式()可得由于7. 给定f(x)=cosx的函数表篇二:数值分析试题1参考答案参考答案 1 一、1.22.xn?1?xn?3.1, 0 4.7,f(xn)(n?0,1,?) ?f(xn)25 7(k?1)15(k)x2x11336. ? ,1(k?1)x2??x1(k?1)1220?2003??10?2?4二、(1) L?0?13?00?1??(2)1?0?120??,U??0100?5??4000?23100??0?? 3??4?1??l65?a65?(l61u15?l62u25?l63u35?l64u45);u55u56?a55?(l51u16?l52u26?l53u356?l54u46)三、先造差分表如下:(1)选x1?,x2?,x3?,x4?为节点,构造三次向前Newton插值多项式2y1?3y1N(x?th)?y1??y1?t(t?1)?t(t?1)(t?2) 31 2!3!将x1和h代入上式,则有N3(?)?25?2t?1/2*t(t?1)?5/6*t(t?1)(?2)由??解得t?,所以f()?N()?(2) 选x3?,x4?,x5?为节点,构造二次向前Newton插值式N2(x3?th)?y3??y3t?t(t?1)2!将x3和h代入上式,则有N2(?)?20?t?t(t?1) 由+=解得t=,所以 f()?N2()?(3)由f(?)3ht(t?1)(t?2)3!(,0?t?2)R2(x0?th)?f(?)3600有R(2(xi?)?(t?1)(t?2)?**maxt(t?1)(t?2)0?t?23!3!可知f(x)有两位整数,故能保证有两位有效数字。
数值分析期末考试一、 设80~=x ,若要确保其近似数的相对误差限为0.1%,则它的近似数x 至少取几位有效数字?(4分)解:设x 有n 位有效数字。
因为98180648=<<=,所以可得x 的第一位有效数字为8(1分) 又因为21101011000110821--⨯=<⨯⨯≤n ε,令321=⇒-=-n n ,可知x 至少具有3位有效数字(3分)。
二、求矩阵A 的条件数1)(A Cond (4分)。
其中⎥⎦⎤⎢⎣⎡=4231A 解:⎥⎦⎤⎢⎣⎡--=-5.05.1121A (1分) 1A =7(1分) 2711=-A (1分)249)(1=A Cond (1分)三、用列主元Gauss 消元法法求解以下方程组(6分)942822032321321321=++-=++--=+-x x x x x x x x x解:→⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---→⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----→⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----5.245.2405.35.230914220321821191429142821120321 ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---→⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---8175835005,245.24091425.33.2305.245.2409142(4分) 等价三角方程组为:⎪⎪⎩⎪⎪⎨⎧-=-=+-=++,8175835,5.245.24,942332321x x x x x x (1分)回代得1,3,5123==-=x x x (1分)四、设.0,2,3,1,103)(3210234=-===-+-=x x x x x x x x f 1)求以3210,,,x x x x 为节3次Lagrange 多项式;(6分) 2)求以3210,,,x x x x 为节3次Newton 多项式;(6分)3)给出以上插值多项式的插值余项的表达式(3分)解:由0,2,3,13210=-===x x x x 可得10)(,34)(,1)(,11)(3210-==-=-=x f x f x f x f即得: +------+------=))()(())()(()())()(())()(()()(312101320130201032103x x x x x x x x x x x x x f x x x x x x x x x x x x x f x L=------+------))()(())()(()())()(())()(()(23130321033212023102x x x x x x x x x x x x x f x x x x x x x x x x x x x f+-+--+-⨯-+-+--+-⨯-)03)(23)(13()0)(2)(1()1()01)(21)(31()0)(2)(3(11x x x x x x326610.)20)(30)(10()2)(3)(1()10()02)(32)(12()0)(3)(1(34x x x x x x x x x -+--=+--+--⨯-+---------⨯2)计算差商表如下:i x )(i x f 一阶差商 二阶差商 三阶差商1 -11 3 -1 5 -2 34 -7 4 0-10-225-1则=+-----+-+-=)2)(3)(1()3)(1(4)1(511)(3x x x x x x x N326610x x x -+--3))2)(3)(1())()()((!4)()(3210)4(3+--=----=x x x x x x x x x x x x f x R ξ五、给定方程组b Ax =,其中⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=100131w w w w A 。
]第一章 绪论姓名 学号 班级习题主要考察点:有效数字的计算、计算方法的比较选择、误差和误差限的计算。
1若误差限为5105.0-⨯,那么近似数有几位有效数字(有效数字的计算) 解:2*103400.0-⨯=x ,325*10211021---⨯=⨯≤-x x 故具有3位有效数字。
2 14159.3=π具有4位有效数字的近似值是多少(有效数字的计算) 解:10314159.0⨯= π,欲使其近似值*π具有4位有效数字,必需!41*1021-⨯≤-ππ,3*310211021--⨯+≤≤⨯-πππ,即14209.314109.3*≤≤π即取( , )之间的任意数,都具有4位有效数字。
3已知2031.1=a ,978.0=b 是经过四舍五入后得到的近似值,问b a +,b a ⨯有几位有效数字(有效数字的计算)解:3*1021-⨯≤-aa ,2*1021-⨯≤-b b ,而1811.2=+b a ,1766.1=⨯b a 2123****102110211021)()(---⨯≤⨯+⨯≤-+-≤+-+b b a a b a b a故b a +至少具有2位有效数字。
2123*****10210065.01022031.1102978.0)()(---⨯≤=⨯+⨯≤-+-≤-b b a a a b b a ab 故b a ⨯至少具有2位有效数字。
4设0>x ,x 的相对误差为δ,求x ln 的误差和相对误差(误差的计算)~解:已知δ=-**xx x ,则误差为 δ=-=-***ln ln xx x x x则相对误差为******ln ln 1ln ln ln xxx x xxx x δ=-=-5测得某圆柱体高度h 的值为cm h 20*=,底面半径r 的值为cm r 5*=,已知cm h h 2.0||*≤-,cm r r 1.0||*≤-,求圆柱体体积h r v2π=的绝对误差限与相对误差限。
(误差限的计算)解:*2******2),(),(h h r r r h r r h v r h v -+-≤-ππ绝对误差限为πππ252.051.02052)5,20(),(2=⨯⋅+⨯⋅⋅⋅≤-v r h v相对误差限为%420120525)5,20()5,20(),(2==⋅⋅≤-ππv v r h v 6设x 的相对误差为%a ,求nx y =的相对误差。
数值分析期末考卷一、选择题(每题4分,共40分)A. 插值法B. 拟合法C. 微分法D. 积分法A. 高斯消元法B. 高斯赛德尔迭代法C. 共轭梯度法D.SOR方法3. 下列哪个算法不是求解非线性方程的方法?A. 二分法B. 牛顿法C. 割线法D. 高斯消元法A. 梯形法B. 辛普森法C. 高斯积分法D. 复化求积法A. 欧拉法B. 龙格库塔法C.亚当斯法D. 高斯消元法A. 幂法B. 反幂法C. 逆迭代法D. QR算法A. 梯度下降法B. 牛顿法C. 共轭梯度法D. 高斯消元法A. 拉格朗日插值法B. 牛顿插值法C. 埃尔米特插值法D. 分段插值法A. 前向差分法B. 后向差分法C. 中心差分法D. 拉格朗日插值法A. 牛顿法B. 割线法C. 雅可比迭代法D. 高斯消元法二、填空题(每题4分,共40分)1. 数值分析的主要任务包括数值逼近、数值微积分、数值线性代数和______。
2. 在求解线性方程组时,迭代法的收敛速度与______密切相关。
3. 牛顿法的迭代公式为:x_{k+1} = x_k f(x_k)/______。
4. 在数值积分中,复化梯形公式的误差为______。
5. 求解常微分方程初值问题,龙格库塔法的阶数取决于______。
6. 矩阵特征值的雅可比方法是一种______方法。
7. 梯度下降法在求解无约束优化问题时,每次迭代的方向为______。
8. 拉格朗日插值多项式的基函数为______。
9. 数值微分中的中心差分公式具有______阶精度。
10. 在求解非线性方程组时,牛顿法的迭代公式为:x_{k+1} =x_k J(x_k)^{1}______。
三、计算题(每题10分,共60分)1. 给定数据点(1,2),(2,3),(3,5),(4,7),求经过这四个数据点的拉格朗日插值多项式。
2. 用牛顿迭代法求解方程x^3 2x 5 = 0,初始近似值为x0 = 2,计算前三次迭代结果。
《数值分析》练习题及答案解析一、单选题1. 以下误差公式不正确的是( D )A .()1212x x x x ∆-≈∆-∆B .()1212x x x x ∆+≈∆+∆C .()122112x x x x x x ∆≈∆+∆D .1122()x x x x ∆≈∆-∆ 2. 已知等距节点的插值型求积公式()()352kkk f x dx A f x =≈∑⎰,那么3kk A==∑( C )A .1 B. 2 C.3 D. 4 3.辛卜生公式的余项为( c )A .()()32880b a f η-''-B .()()312b a f η-''-C .()()()542880b a f η--D .()()()452880b a f η--4. 用紧凑格式对矩阵4222222312A -⎡⎤⎢⎥=-⎢⎥⎢⎥--⎣⎦进行的三角分解,则22r =( A ) A .1 B .12C .–1D .–25. 通过点()()0011,,,x y x y 的拉格朗日插值基函数()()01,l x l x 满足( D ) A .()00l x =0,()110l x = B . ()00l x =0,()111l x = C .()00l x =1,()111l x = D . ()00l x =1,()111l x =6. 用二分法求方程()0f x =在区间[],a b 上的根,若给定误差限ε,则计算二分次数的公式是n ≥( D )A .ln()ln 1ln 2b a ε-++ B. ln()ln 1ln 2b a ε-+-C.ln()ln 1ln 2b a ε--+ D. ln()ln 1ln 2b a ε--- 7.若用列主元消去法求解下列线性方程组,其主元必定在系数矩阵主对角线上的方程组是( B )A .123123123104025261x x x x x x x x x -+=⎧⎪-+=⎨⎪-+=-⎩ B 。
第一章 绪论习题主要考察点:有效数字的计算、计算方法的比较选择、误差和误差限的计算。
1 若误差限为5105.0-⨯,那么近似数0.003400有几位有效数字?(有效数字的计算) 2 14159.3=π具有4位有效数字的近似值是多少?(有效数字的计算)3 已知2031.1=a ,978.0=b 是经过四舍五入后得到的近似值,问b a +,b a ⨯有几位有效数字?(有效数字的计算)4 设0>x ,x 的相对误差为δ,求x ln 的误差和相对误差?(误差的计算)5测得某圆柱体高度h 的值为cm h 20*=,底面半径r 的值为cm r 5*=,已知cm h h 2.0||*≤-,cm r r 1.0||*≤-,求圆柱体体积h r v2π=的绝对误差限与相对误差限。
(误差限的计算)6 设x 的相对误差为%a ,求nx y =的相对误差。
(函数误差的计算)7计算球的体积,为了使体积的相对误差限为%1,问度量半径r 时允许的相对误差限为多大?(函数误差的计算)8 设⎰-=11dx e x eI x n n ,求证: (1))2,1,0(11 =-=-n nI I n n(2)利用(1)中的公式正向递推计算时误差逐步增大;反向递推计算时误差逐步减小。
(计算方法的比较选择)第二章 插值法习题主要考察点:拉格朗日插值法的构造,均差的计算,牛顿插值和埃尔米特插值构造,插值余项的计算和应用。
1 已知1)2(,1)1(,2)1(===-f f f ,求)(x f 的拉氏插值多项式。
(拉格朗日插值)2 已知9,4,10===x x x y ,用线性插值求7的近似值。
(拉格朗日线性插值)3 若),...1,0(n j x j =为互异节点,且有)())(())(()())(())(()(11101110n j j j j j j j n j j j x x x x x x x x x x x x x x x x x x x x x l ----------=+-+-试证明),...1,0()(0n k x x l xnj k jk j =≡∑=。
1.1 求下列各数的具有四位有效数字的近似值, 并指出其绝对误差限和相对误差限)1.0ln(,121,1011,1014321====x x x x1.2 下列各数都是对准确值进行四舍五入得到的近似值, 指出它们的绝对误差限、相对误差限和有效数字的位数。
3*5*4*3*2*1100.5,5000,50.31,3015.0,0315.0⨯=====x x x x x1.3 为了使31的近似值的相对误差不超过0.1%, 问应取几位有效数字?1.4 怎样计算下列各题才能使得结果比较精确?(1) x x sin )sin(-+ε,其中ε充分小 (2) ⎰++121N Nx dx,其中N 是充分大的正数(3)xxsin cos 1-,其中x 充分小(4) o 1cos 1- (5) 1001.0-e(6) )11010ln(84--1.5 求方程01562=+-x x 的两个根, 使至少具有四位有效数字。
2.1 证明方程043=-+x x 在区间[1,2]内有且仅有一个根。
如果用二分法求它具有五位有效数字的根,试问需对分多少次?(不必求根)2.2 用二分法求方程0134=+-x x 在[0.3, 0.4]内的一个根, 精度要求21021-⨯=ε。
2.3 找出下列方程的有根区间,选择适当的初始点用二分法求方程的根,精度要求210-=ε。
(1) 02=--x x ;(2) 06cos 2=-++-x e x x ; (3) 01tan =--x x ; (4) 0sin 2=--x e x 。
2.4 考虑方程032=-x e x ,将其改写为3xex ±=,取00=x ,用两种迭代公式迭代,分别收敛到1.0和-0.5附近的两个根(取精度要求310-=ε)。
2.5 为求方程0123=--x x 在5.1=x 附近的一个根,建立下列形式的迭代公式:(1) 2121111kk x x xx +=⇒+=+,;(2) 3212311k k x x x x +=⇒+=+,;(3) 111112-=⇒-=+k k x x x x ,。
试分析每一种迭代公式的收敛性。
2.6 考虑用迭代法求解下列方程:(1) )2(312x ex x+-=-;(2) xx -=5;(3) 27475.1--+=x x x 。
按所给的形式建立迭代公式,试确定区间[a, b ], 使迭代公式收敛, 并求出满足精度要求310-=ε的解。
2.7 用迭代法的思想,给出求22222+++++的迭代公式,并证明:222222lim=+++++∞→nn 。
2.8 能否用简单迭代法求解下列方程,如果不能,请给出收敛的迭代公式。
(1) )sin (cos 41x x x +=;(2) x x 24-=。
2.9 已知)(x x ϕ=在区间[a,b ]内有一个根,且当a<x<b 时,1)(>≥'k x ϕ。
试问如何将)(x x ϕ=化为收敛的迭代公式。
2.10 用Steffensen 加速迭代法求方程13-=x x 在[1,1.5]内的根2.11 试用Newton 法求方程032=-x e x 的根, 分别取初始点0.4,0.1,5.00-=x , 精度要求为310-=ε。
2.12 选择适当的初始点, 试用Newton 法求出满足精度要求为310-=ε的解(1) )2(312x ex x+-=-;(2) 0cos 102=+x x 。
2.13 导出计算)0(1>a a的Newton 迭代公式,使公式中即无开方又无除法运算。
2.14 设⎩⎨⎧-∞∈--∞∈=)0,(,),0[,)(1x x x x x f ,⎪⎩⎪⎨⎧-∞∈-∞∈=)0,(,),0[,)(32322x x x x x f ,函数)(1x f 和)(2x f 均有零点x =0,分别讨论用Newton 法解0)(1=x f 和0)(2=x f 是否收敛?收敛的阶是多少?2.15 用Newton 法设计一种不用除法的迭代公式,求正数c 的倒数,并证明:当初值0x 满足cx 200<<时,该迭代法收敛。
2.16 试用Newton 法解方程03=-a x , 导出求立方根3a 的迭代公式,讨论取什么初值可使迭代收敛。
2.17 为了简化,在Newton 迭代公式中用)(0x f '代替)(k x f ',试问这种迭代是几阶的? 2.18 设23)()(a x x f -=, 写出解0)(=x f 的Newton 迭代公式, 并证明迭代公式是线性收敛的。
2.19 设非线性方程0)3)(133()(23=+-+-=x x x x x f , 其根1,3*2*1=-=x x . 写出求*1x 的近似值时,二阶局部收敛的Newton 迭代公式和求*2x 的近似值时,二阶局部收敛的Newton 迭代公式。
2.20 设0)(=x f 有根,且)()(0∞<<-∞≤'≤<x M x f m 。
试证明:由)(1k k k x f x x λ-=+产生的序列k x 对任意的0x 和M20<<λ均收敛于根。
2.21 为用迭代法求方程0)(=x f 的根α,若将方程改写成)()(x g x Cf x x =+=, 其中C 为待定常数。
设)(x f '连续且0)(≠'αf , 试确定C , 使序列)(1k k x g x =+收敛于α,且尽可能收敛得快。
习题三3.1 考虑线性方程组6545953232121-=+-=-+-=-x x x x x x x (1) 用顺序Gauss 消去法求解该方程组; (2) 用LU 分解算法求解该方程组。
3.2 考虑线性方程组114231124342321321321=+-=-+=--x x x x x x x x x (1) 用顺序Gauss 消去法求解该方程组; (2) 用列主元Gauss 消去法求解该方程组。
3.3 考虑线性方程组53.048.104.368.203.130.1874.048.105.053.448.151.2321321321-=-+=-+=++x x x x x x x x x (1) 用顺序Gauss 消去法求解该方程组; (2) 用列主元Gauss 消去法求解该方程组;(3) 试检验⑴和⑵所得的两个解中, 哪个更接近准确解?(计算过程保留三位有效数字. 方程的准确解为: (1.35533, -1.29208, -0.252451)。
3.4 试用列主元Gauss-Jordan 消去法求下列矩阵的逆矩阵(用分数运算)⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--121011322)2(,011012111)1( 3.5 对下列矩阵作LU 分解⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--3010211123.6 设n n ij a A ⨯=)(, 试导出A=LU 分解(Crout 分解)的计算公式. 其中L 为下三角矩阵, U 为单位上三角矩阵, 即⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=111,211221222111nn nnn n u u u U l l l l l l L 并用此公式得到求解线性方程组Ax=b 的计算公式。
3.7 用追赶法解方程组453423124343232121=++=+++=++=+x x x x x x x x x x3.8 导出Crout 的形式⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡----1111111221112211n n nn n nn n u u l a l l a l b a c b b a c b追赶法的计算公式。
3.9 设0),(11≠=a a A ij , 经过一步Gauss 消去法得到⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎦⎤⎢⎣⎡=)2()2(2)2(2)2(222211,0nn n n Ta a a a A A a A其中α 试证明:(1) 若A 对称,则2A 对称;(2) 若A 对称正定,则2A 对称正定;3.10 设A 对称正定,[]ij a A =, 经1步Gauss 消去后约化为[])2()2(ija A=,试证:(1) n i a ii ,,2,1,0 =>, 且A 绝对值最大的元素在对角线上;(2) n i a a ii ii ,,2,1,)2( =≤;(3) ij nj i ijnj i a a max max ,1)2(,1≤≤≤≤≤。
3.11 试证明:单位下三角阵的逆、积还是单位下三角阵。
3.12 试证明:如果A 是对称正定阵,则1-A 也是对称正定矩阵且A 可唯一地写成形式L L A T=,其中L 是具有对角元的下三角阵。
3.13 设⎥⎦⎤⎢⎣⎡-=211a A , 试分析当a 为何值时可作L L T分解, 其中L 是对角线元素为正的下三角阵,并求矩阵L 。
3.14 设⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=22112aa A , 为使A 可分解为L L A T =, 其中L 为对角线元素为正的下三角阵, a 的取值范围是多少?若取a =1, 求矩阵L 。
3.15 已知Ax=f , 其中⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---=010,212112f ba b A(1) 试问参数a 和b 满足什么条件时,可选用平方根法求解该方程组? (2) 取b=0, a = -1, 试用追赶法求解该方程组。
3.16 设n R x ∈,试证明:∞∞≤≤xn xx 1∞∞≤≤xn xx 21211x xxn≤≤3.17 设n n R A ⨯∈,试证明:∑=≤≤=ni ij nj a A111max3.18 设n n R A ⨯∈,试证明:FF AAAn ≤≤213.19 设矩阵A 非奇异,求证AA11≥-3.20 设矩阵A 非奇异,求证nA λλ1)(Cond ≥其中n λλ,1分别是矩阵A 的最大最小特征值,且当A 为对称矩阵时,上式等号成立。
3.21 设矩阵A 非奇异,试证明:若11<-A A δ, 则)(A A δ+非奇异,且满足AA A A δδ-≤+-1)(13.22 方程组Ax=b , 设A 非奇异阵,b 有扰动b δ, 从而引起方程组解x 有扰动x δ, 试证明:bbA xxδδ)(Cond ≤3.23 求下面两个方程的解,并利用矩阵的条件数对解进行分析。
⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡++⎥⎦⎤⎢⎣⎡--⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡⎥⎦⎤⎢⎣⎡--43240179319240432405.1795.319240221121x x x x x x δδ, 3.24 已知⎥⎦⎤⎢⎣⎡=989999100A , 求pA 和∞=,2,1,)(Cond p A p 。