高等工程数学试题解答

合集下载

《高等工程数学》科学出版社吴孟达版习题答案(1-8章)

《高等工程数学》――科学出版社版习题答案：第一章习题（P26） 1．略2．在R 4中，求向量a ＝[1,2,1,1]T ,在基a 1 ＝ [1 , 1, 1, 1]T , a 2 ＝ [1 , 1, －1,－1]Ta 3 ＝ [1 , －1, 1, －1]T a 4 ＝ [1 , －1,－1, 1]T 下的坐标。

解：其坐标为：x ＝( 5/4, 1/4, －1/4，－1/4 )T 3．在R 2×2中，求矩阵12A=03⎡⎤⎢⎥⎣⎦，在基 111B =11⎡⎤⎢⎥⎣⎦，211B =10⎡⎤⎢⎥⎣⎦，311B =00⎡⎤⎢⎥⎣⎦，410B =00⎡⎤⎢⎥⎣⎦下的坐标。

解：其坐标为：x ＝( 3, －3, 2，－1 )T4．试证：在R 2×2中，矩阵111B =11⎡⎤⎢⎥⎣⎦，211B =01⎡⎤⎢⎥⎣⎦，311B =10⎡⎤⎢⎥⎣⎦，410B =11⎡⎤⎢⎥⎣⎦线性无关。

证明：设 k 1B 1+ k 2B 2+ k 3B 3+ k 4B 4=0000⎡⎤⎢⎥⎣⎦，只要证明k 1= k 2 = k 3= k 4 =0即可。

余略。

5．已知R 4中的两组基：T T T T 1234=[1,0,0,0],=[0,1,0,0],=[0,0,1,0],=[0,0,0,1]αααα和T T T T 1234=[2,1,1,1],=[0,3,1,0],=[5,3,2,1],=[6,6,1,3]ββββ－求由基1234{,,,}αααααB =到基1234{,,,}βββββB =的过渡矩阵，并求向量1234[,,,]x x x x ξ=在基1234{,,,}βββββB =的坐标。

解：基1234{,,,}αααααB =到基1234{,,,}βββββB =的过渡矩阵是：2056133611211013⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦－向量1234[,,,]x x x x ξ=在基1234{,,,}βββββB =的坐标是：11234205612927331336112923x 112190018101373926x x x x ⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦－－－－－1＝－－27－－6．设R[x]n 是所有次数小于n 的实系数多项式组成的线性空间，求多项式p(x) = 1+ 2x n －1在基{1，(x －1)，(x －1)2，(x －1)3，….，(x －1)n －1}的坐标。

中南大学高等工程数学试题及参考答案.docx

中南大学专业硕士“高等工程数学Ⅰ”考试试卷（开卷）考试日期： 2014 年月日时间 100 分钟注：解答全部写在答题纸上一、填空题 ( 本题 24 分，每小题 3 分 )1111324（1）如果Ax b, A 161，矩阵 A 1， A，利用 Gauss-Seidel 迭253113344代法求解此方程组是否收敛；答案：953，收敛，212解析： || A ||1为列范数，等于各列绝对值之和的最大值，||A ||为行范数，等于各行绝对值之和的最大值， A 为严格对角占优矩阵，根据课本P143定理 5.4.12 知， Jacobi 和 G-S 均收敛。

（ 2）利用迭代法求解非线性方程 f ( x) 2x e x0 的根，取初值 x0 0.5 。

给出一个根的存在区间，在该区间上收敛的迭代函数为；答案： [-1 ，0] ，g( x) 1 e x2解析：1 1 xf (1)20,f(0)10 ，故在[-10]g(x)e，根据课本P93定理 4.2.32e1可知迭代函数收敛的条件：（1）在[-1,0] 上一阶导数存在；（ 2）x [1,0] ，均有 | g(x) |[-1,0]；（3）| g' ( x) |max 1 ，2 1e x在[-1,0]上收敛。

故 g( x)2（3）设事件A发生的概率为p，在 n 次重复试验中事件m np近似服A 发生次数为m，当 n 充分大时，m )m(1n从的分布为；答案:N (0,1)解析：课本 P187 定理 7.2.4（4）设x1 , x2 , x3 , x4[ 1,1] ，若数值积分公式1 f (x)dx A1 f ( x1 ) A2 f ( x2 ) A3 f ( x3 )A4 f ( x4 ) 的代数精度大于11，则A1A2A3A4；答案： 21解析：令 f ( x) 1 ，可得1dx2A1 A2A3A4。

1（ 5）已知y f ( x) 通过点(x i, y i), i0,1,2,3 ，则其Lagrange插值基函数l2( x)；答案： l 2 ( x)(x x0 )( x x1)( x x3 ) ( x2 x0 )( x2x1 )( x2x3 )解析：课本 P20 拉格朗日插值基函数的定义（式 2.3.2）。

大学工程数学考试题及答案

大学工程数学考试题及答案一、选择题（每题2分，共10分）1. 下列哪个选项是微积分的基本定理？A. 积分中值定理B. 洛必达法则C. 牛顿-莱布尼茨公式D. 泰勒级数展开答案：C2. 在概率论中，随机变量X服从二项分布B(n, p)，其中n=10，p=0.3，那么E(X)等于多少？A. 2B. 3C. 4D. 5答案：A3. 线性代数中，一个矩阵A可逆的充分必要条件是什么？A. 行列式非零B. 秩等于A的阶数C. A的所有特征值非零D. 所有选项都是答案：D4. 在复数域中，下列哪个表达式表示复数的共轭？A. z + z*B. z - z*C. |z|^2D. z * z*答案：B5. 傅里叶级数在工程数学中的应用之一是？A. 信号处理B. 量子力学C. 统计物理D. 所有选项都是答案：A二、填空题（每题3分，共15分）6. 函数f(x) = sin(x)的一阶导数是_________。

答案：cos(x)7. 矩阵的特征值是_________。

答案：λ8. 拉普拉斯变换的逆变换通常使用_________。

答案：拉普拉斯逆变换9. 随机变量X和Y相互独立，且P(X=x) = 2x，P(Y=y) = 3y，则P(X+Y=4)等于_________。

答案：1/410. 曲线y = x^2在点(1,1)处的切线斜率是_________。

答案：2三、解答题（共75分）11. （15分）证明函数f(x) = e^x在实数域上是单调递增的。

答案：由于f'(x) = e^x > 0对于所有实数x，因此f(x)在实数域上是单调递增的。

12. （20分）解线性方程组：\[\begin{align*}x + 2y &= 5 \\3x - y &= 4\end{align*}\]答案：使用高斯消元法或克拉默法则，解得 \( x = 2, y = 1.5 \)。

13. （20分）计算下列定积分：\[\int_{0}^{1} x^2 dx\]答案：使用基本积分公式，得到 \( \frac{1}{3}x^3 \) 在0到1的积分为 \( \frac{1}{3} \)。

高等工程数学考试部分答案1

工程硕士学位课程考试
高等工程数学试题
注意：每位考生只要选做以下两部分试题，答案必须写在答题纸上
矩阵分析部分
一．（6分）设求值。

解：参考试题2第一题
二．（8分）已知函数矩阵：，求矩阵
解：参考试题2第二题
三．（10分）设向量
与，令，
（1）求的一组基和维数；（2）求维数。

解：参考试题2第三题
四．（10分）设，
1．求的Jordan标准形及最小多项式；
解: 矩阵的最小多项式为, Jordan标准形为
2。

求解初值问题
解：参考试题2第四题（2）小题
五．（8分）设与是线性空间的两个基，为从基到的过渡矩阵，为的一个线性变换，在基下的矩阵，求线性变换在基下的矩阵。

解: 由题意有
所以由第一式有
把第二式和第三式代入得到
把第一式代入左边得到
从而有, 所以
六．（8分）设且可逆，，求证：的特征值都是正数。

证明: 因为为正规矩阵, 所以酉矩阵与对角矩阵. 即存在酉矩阵, 使得, 其中为对角矩阵, 从而
所以的元素全为实数. 设为任意一个特征值, 是属于的特征向量, 则有
得证.。

高等工程数学第二章习题及答案

第2章线性代数方程组数值解法研究n 阶线性方程组Ax b =的数值解法.()ij A a =是n n⨯矩阵且非奇异，12(,,,)Tn x x x x = ,12(,,,)Tn b b b b =两类数值方法：(1) 直接法：通过有限次的算术运算，若计算过程中没有舍入误差，可以求出精确解的方法.Ax b Gx d == 等价变换G 通常是对角矩阵、三角矩阵或者是一些结构简单的矩阵的乘积.(2) 迭代法：用某种极限过程去逐次逼近方程组的解的方法.(1)()i i Ax b x Bx k x Bx k +==+−−−−−→=+ 等价变换建立迭代格式，0,1,i =一、向量范数与矩阵范数 1. 向量范数【定义】若对nK 上任一向量x ，对应一个非负实数x ，对任意,nx y R ∈及K α∈，满足如下条件(向量范数三公理) (1) 非负性：0x ≥，且0x =的充要条件是0x =；(2)齐次性：x xαα=；(3)三角不等式：x y x y+≤+.则称x为向量x的范数.常用的向量范数： (1) 1—范数11nii x x ==∑(2) 2—范数12221()ni i x x ==∑(3) ∞—范数1max ii nxx ∞≤≤=(4) 一般的p —范数11()pnpi pi xx ==∑2. 矩阵范数【定义】若n nK ⨯上任一矩阵()ij n n A a ⨯=，对应一个非负实数A ，对任意的,n nA B K ⨯∈和K α∈，满足如下条件(矩阵范数公理)：(1) 非负性：0A ≥，且0A =的充要条件是0A =；(2)齐次性：A Aαα=；(3)三角不等式：A B A B +≤+；(4)乘法不等式：AB A B≤.则称A为矩阵A的范数.矩阵范数与向量范数是相容的：Ax A x≤向量范数产生的从属范数或算子范数：10max maxx x AxA Ax x=≠==常见从属范数：(1) 1—范数111max ||nij j ni A a ≤≤==∑(2) ∞—范数11max ||nij i nj A a ∞≤≤==∑(3) 2—范数2A =谱半径1()max ||H i i n A A ρλ≤≤=，iλ为H A A 的特征值.H A 为A 的共轭转置. 注：矩阵A 的谱半径不超过A 的任一范数，即()A A ρ≤范数等价性定理：,s t x x为n R 上向量的任意两种范数，则存在常数12,0c c >，使得12,ns t s c x x c x x R ≤≤ ∀∈.注：矩阵范数有同样的结论. 【定理2.1】是任一向量范数，向量序列()k x 收敛于向量*x 的充要条件是()*0,k x x k -→ →∞二、 Gauss 消去法 1．顺序Gauss 消去法将方程Ax b =写成如下形式11112211,121122222,11122,1n n n n n n n n nn n n n a x a x a x a a x a x a x a a x a x a x a ++++++=⎧⎪+++=⎪⎨⎪⎪+++=⎩其中记,1,1,2,,.i n i a b i n +==消元过程：第一次消元：设110a ≠，由第2,3,,n 个方程减去第一个方程乘以1111/(2,3,,)i i m a a i n == ，则将方程组中第一个未知数1x消去，得到同解方程11112211,1(1)(1)(1)22222,1(1)(1)(1)22,1n n n n n n n nn n n n a x a x a x a a x a x a a x a x a ++++++=⎧⎪ ++=⎪⎨⎪⎪ ++=⎩其中， (1)11,2,3,,;2,3,,,1ijij i j a a m a i n j n n =-==+ . 1111/i i m a a =，2,3,,i n = .第二次消元：设(1)220a ≠，.由第2,3,,n 个方程减去方程组中的第2个方程乘以(1)(1)2222/(3,4,,)i i m a a i n == ，则将方程组第2个未知数2x 消去，得到同解方程11112213311,1(1)(1)(1)(1)2222322,1(2)(2)(2)33333,1(2)(2)(2)33,1n n n n n n n n n nnn n n n a x a x a x a x a a x a a x a a x a x a a x a x a ++++++++=⎧⎪ +++=⎪⎪ ++=⎨⎪⎪⎪ ++=⎩其中(2)(1)(1)22, 3,4,,; 3,4,,,1ij ij i j a a m a i n j n n =-==+ . (1)(1)2222/i i m a a =，3,4,,i n = .经过1n -次消元后，原方程组变成等价方程组11112213311,1(1)(1)(1)(1)2222322,1(2)(2)(2)33333,1(1)(1),1n n n n n n n n n n n nn n n n a x a x a x a x a a x a a x a a x a x a a x a +++--+++++=⎧⎪ +++=⎪⎪ ++=⎨⎪⎪⎪ =⎩其中()(1)(1), 1,2,,k k k ij ij ik ij a a m a i k k n --=-=++ ， 1,2,,,1j k k n n =+++ .(1)(1)/k k ik ik kkm a a --=，1,2,,i k k n =++ ；1,2,,1k n =- .回代过程：(1)(1),1(1)(1)(1),1,,1/[]/,1,2,,2,1.n n n n n m n i i i ii n i j j i j j i x a a x a a x a i n n --+---+=+⎧=⎪⎨=-=--⎪⎩∑计算量：按常规把乘除法的计算次数合在一起作为Gauss 消去法总的计算量，而略去加减法的计算次数. 在消去过程中，对固定的消去次数(1,2,,1)k k n =- ，有：除法(1)(1),,/,1,1,,k k ik i k k k m a a i k k n --= =++ 共计n k -次；乘法(1),,1,2,,;1,2,,,1k ik k j m a i k k n j k k n n - =++ =+++ 共计()(1)n k n k --+次.因此，消去过程总的计算量为1311[()(1)]3n k M n k n k n k n-==--++-≈∑ 回代过程的乘除法计算次数为21()2n n +.与消去法计算量相比可以略去不计.所以， Gauss 消去法总的计算量大约为313n .2. Gauss-Jordan 消去法Gauss-Jordan 消去法是Gauss 消去法的一种变形.此方法的第一次消元过程同Gauss 消去法一样，得到(1)(1)(1)(1)11112213311,1(1)(1)(1)(1)22223322,1(1)(1)(1)(1)32233333,1(1)(1)(1)(1)2233,1,,,,n n n n n n n n n nn nn n n n a x a x a x a x a a x a x a x a a x a x a x a a x a x a x a ++++⎧++++=⎪ +++=⎪ +++=⎨ +++= ⎪⎪⎪⎪⎩其中，(1)11,2,,,1jj a a j n n ==+ . 第二次消元：设(1)220a ≠，由第1,3,4,,n 个方程减去第2个方程乘以(1)(1)2222/(1,3,4,,)i i m a a i n == ,则得到同解方程组(2)(2)(2)11113311,1(1)(2)(2)(2)22223322,1(2)(2)(2)33333,1(2)(2)33,1,,,n n n n n n n n n nnn n n n a x a x a x a a x a x a x a a x a x a a x a x a +++++ +++= +++= ++= ++= (2),⎧⎪⎪⎪⎨⎪⎪⎪⎩继续类似的过程，在第k 次消元时，设(1)k kk a -，将第i 个方程减去第k 个方程乘以(1)(1)/k k ik ik kk m a a --=，这里1,3,4,1,1,,i k k n =-+ .经过1n -次消元，得到(2)1111,1(1)(2)2222,1(2)(2)33,1,,,n n n n n a x a a x a a x a +++⎧ =⎪ =⎪⎪ ⎨⎪⎪⎪ =⎩其中()(1)(1),1,2,,1,1,,k k k ij ij ik kj a a m a i k k n --=-=-+ ；1,2,,,1; 1,2,,1j n n k n =+=- .此时，求解回代过程为(1)(1),1/,1,2,,n i i i n iix a a i n --+= = 经统计，总的计算量约为312M n ≈次乘除法. 从表面上看Gauss-Jordan 消去法似乎比Gauss 消去法好，但从计算量上看Gauss －Jordan 消去法明显比Gauss消去法的计算量要大，这说明用Gauss-Jordan 消去法解线性方程组并不可取.但用此方法求矩阵的逆却很方便. 3．列选主元Gauss 消去法在介绍Gauss 消去法时，始终假设(1)0k kk a -≠，称(1)k kka -为主元.若(1)0k kka -=，显然消去过程无法进行.实际上，既使(1)0k kka -≠，但(1)k kka -很小时，用它作除数对实际计算结果也是很不利的.称这样的(1)k kka -为小主元.【例2.2】设计算机可保证10位有效数字，用消元法解方程1112120.3100.7,0.9,x x x x -⎧⨯+=⎪⎨ +=⎪⎩【解】经过第一次消元：第2个方程减去第1个方程乘以212111/m a a =得1112(1)(1)222230.3100.7x x a x a -⎧⨯+=⎪⎨ =⎪⎩其中(1)1222222111/0.333333333310a a a a =-=-⨯，(1)123323211113(/)0.233333333310a a a a a =-⋅=-⨯于是解得(1)(1)223221/0.7000000000,0.0000000000,x a a x ⎧==⎪⎨=⎪⎩而真解为120.2,0.7x x = =注：造成结果失真的主要因素是主元素11a太小，而且在消元过程中作了分母，为避免这个情况发生，应在消元之前，作行交换.【定义】若 (1)(1)||max ||k k k r k ik k i na a --≤≤=，则称(1)||k k r k a - 为列主元素. k r 行为主元素行，这时可将第 k r行与第k 行进行交换，使(1)||k k r k a - 位于交换后的等价方程组的 (1)k kk a - 位置，然后再施实消去法，这种方法称为列选主元Gauss 消去法或部分主元Gauss 消去法.【例2.3】应用列选主元Gauss 消去法解上述方程. 【解】因为2111a a >，所以先交换第1行与第2行，得1211120.9,0.3100.7,x x x x -⎧+=⎪⎨⨯+=⎪⎩ 然后再应用Gauss 消去法，得到消元后的方程组为1220.9,0.7.x x x ⎧+=⎨=⎩回代求解，可以得到正确的结果.即120.2,0.7x x = =.三、三角分解法设方程组Ax b =的系数矩阵A 的顺序主子式不为零.即1112121222110,1,2,,.kk k k k kka a a a a a k n a a a ∆=≠=在Gauss 消去法中，第一次消元时，相当于用单位下三角阵211131111010010n m L m m -⎡⎤⎢⎥- ⎢⎥⎢⎥=- ⎢⎥ ⎢⎥⎢⎥- ⎢⎥⎣⎦ ，左乘方程组Ax b =，得11A x b =，其中11121(1)(1)122211(1)200n n n nn a a a a a A L a a -(1)⎡⎤⎢⎥ ⎢⎥==⎢⎥ ⎢⎥⎢⎥ ⎣⎦ ，1(1)(1)111,11,1,1(,,,)Tn n n n b L b a a a -+++== .第二次消元时，相当于用单位下三角阵1232210101001n L m m - ⎡⎤⎢⎥ ⎢⎥⎢⎥= - ⎢⎥⎢⎥⎢⎥ - ⎢⎥⎣⎦0 ，左乘方程组11A x b =，得22A x b =其中11121(1)(1)22211(2)(2)221333(2)(2)300000n n n n nn a a a a a A L L A a a a a --⎡⎤ ⎢⎥ ⎢⎥⎢⎥== ⎢⎥⎢⎥ ⎢⎥ ⎢⎥⎣⎦ ，11(1)(2)(2)2211,12,13,1,1(,,,,).Tn n n n n b L L b a a a a --++++==经过1n -次消元，最后得到等价方程组11n n A x b --=其中11121(1)222111111221(1)n n n n n n nn a a a a a A L L L L A a (1)--------⎡⎤⎢⎥ ⎢⎥==⎢⎥⎢⎥⎢⎥ ⎣⎦1111(1)(1)112221,12,1,1(,,,)n Tn n n n n n n b L L L L b a a a --------+++==注意到1n A -是一个上三角阵，记111111221n n n U A L L L L A -------==则121()n A L L L U LU -==其中，121n L L L L -= . 不难验证21313212_1111n n nn m L m m m m m ⎡⎤⎢⎥ ⎢⎥⎢⎥= ⎢⎥ ⎢⎥⎢⎥ 1 ⎢⎥⎣⎦是单位下三角阵.于是解线性方程组Ax b =，就转化为解方程 LUx b =，若令Ux y =就得到一个与 Ax b =等价的方程组Ly b Ux y =⎧⎨=⎩【定理2.2】若 A 为 n 阶方阵，且 A 的所有顺序主子式0k ∆≠，1,2,,k n = .则存在唯一的一个单位下三角矩阵 L 和一个上三角矩阵 U ，使A LU =.在上述过程中，若不假设A 的顺序主子式都不为零，只假设A 非奇异，那么Gauss 消去法将不可避免要应用两行对换的初等变换.第一次消元，将第1行与第1r 行交换，相当于将方程组Ax b =左乘矩阵11r P ：1111r r P Ax P b=经第一次消元得11111111r r L P Ax L P b--=即系数矩阵为11111r A L P A-=，其中110111r P ⎡⎢ ⎢ 1= 1 0 1 ⎣0 0 ⎤⎥⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎦1 列 1r列类似地，经1n -次消元，有121111111,22,11n n n n n r n n r r A L P L P L P A----------= .如果预先知道每一个(1,2,,1)iir P i n =- ，则在消元之前就全部作交换，得 1211,2,1,n n n r n r r A P P P A PA----== ，其中，1211,2,1,n n n r n r r P P P P ----= .即原方程变为PAx Pb =然后再消元，相当于对PA 做三角分解PA LU =由以上讨论，可得结论【定理2.3】若A 非奇异，则一定存在排列矩阵 P ，使得 PA 被分解为一个单位下三角阵和一个上三角1 行1行r阵的乘积，即PA LU =成立.这时，原方程组Ax b = 等价于 PAx Pb =，即等价于求解LUx Pb =令Ux y =则Ly Pb =实际求解时，先解方程组Ly Pb =,再根据 y 求解 Ux y =,即得原方程组Ax b =的解. 这种求解方法称为三角分解法.常用三角分解方法有以下几种. 1．Doolittle 分解方法假设系数矩阵A 不需要进行行交换，且三角分解是唯一的. 记21121110n n l L l l ⎡⎤⎢⎥ ⎢⎥=⎢⎥ ⎢⎥ ⎢⎥⎣⎦ , 11121222n n nn u u u u u U u ⎡⎤⎢⎥ ⎢⎥=⎢⎥ ⎢⎥ 0 ⎣⎦ 于是有1112111121222212222112111110n n n n n n n n nn a a a u u u u u a a a l l l a a a ⎡⎤ ⎡⎤⎢⎥⎢⎥ ⎢⎥⎢⎥=⎢⎥⎢⎥ ⎢⎥⎢⎥ ⎢⎥⎢⎥ ⎣⎦⎣⎦ nn u ⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥0 ⎣⎦从前面讨论A 的LU 分解过程可看出,L 、U 的元素都是用有关的(1)k ij a -来表示的，而它们的计算较麻烦.现在给出直接从系数矩阵A ，通过比较等式的两边逐步把L 和U 构造出来的方法，而不必利用Gauss 消去法的中间结果(1)k ij a -.计算步骤： (1) 由L 阵的第1行分别乘U 阵的各列，先算出U 阵的第1行元素 11,1,2,,j j u a j n = = .然后，由L 阵的各行分别去乘U 阵的第1列，算出L 阵的第1列元素1111/,2,3,,i i l a a i n = = .(2)现假设已经算出U 阵的前1r -行元素，L 阵的前1r -列元素，下面来算U 阵的第r 行元素，L 阵的第r 列元素.由L 阵的第r 行分别乘U 阵的第j 列(,1,,)j r r n =+ ，得11r ij rk kj rjk a l u u -==+∑所以，得U 阵的第r 行元素11,,1,,r rj rj rk kj k u a l u j r r n-==- =+∑ .再由L 阵的第i 行(1,2,,)i r r n =++ 分别去乘U 阵的第r 列,得11r ir ik kr ir rrk a l u l u -==+∑,所以，得L 阵的第r 列元素11[]/,1,2,,.r ir ir ik kr rr k l a l u u i r r n -==- =++∑取1,2,,r n = 逐步计算，就可完成三角分解A LU =；(3)解与Ax b = 等价的方程组Ly b Ux y =⎧⎨=⎩逐次用向前代入过程先解Ly b = 得1111,2,3,,.i i i ij j j y b y b l y i n -==⎧⎪⎨=- =⎪⎩∑然后再用逐次向后回代过程解Ux y =得1/,()/,1,2,,2,1.n n nn n i i ij j ii j i x y u x y u x u i n n =+=⎧⎪⎨=- =--⎪⎩∑2．Crout 分解方法仍假设系数矩阵A 不需要进行行交换，且三角分解是唯一的.即ˆA L=ˆU .与Doolittle 分解方法的区别在111212122211n n n n nn a a a a a a a a a ⎡⎤ ⎢⎥ ⎢⎥=⎢⎥ ⎢⎥⎢⎥ ⎣⎦ 1122ˆˆl l ⎡⎤ 0⎢⎥ ⎢⎥⎢⎥ ⎢⎥⎢⎥⎣⎦ 122ˆ1ˆ10n u u ⎡⎤⎢⎥ ⎢⎥⎢⎥ ⎢⎥ 1 ⎣⎦ 比较两边，则可推导出与Doolittle 分解方法类似的公式，不过Crout 分解方法是先算ˆL 的第r 列，然后再算ˆU的第r 行.3．Cholesky 分解方法若 A 为对称正定矩阵，则有 ˆT U L =，即11()()TT T A LDL LD LD LL ===其中L 为下三角阵. 进一步展开为1121111211112122221222221212n n n n n n nn n n nn a a a l l l l a a a l l l l l l l a a a ⎡⎤⎡⎤ ⎢⎥⎢⎥ 0 ⎢⎥⎢⎥=⎢⎥⎢⎥ ⎢⎥⎢⎥ ⎢⎥ ⎢⎥⎣⎦⎣⎦ 0nn l ⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥ ⎣⎦ 比较两边对应元素，容易得到12121()r rr rr rk k l a l -==-∑ ，11()/r ir ir ik rk rrk l a l l l -==-∑ 1,2,,;1,2,,.r n i r r n ==++Cholesky 分解的优点：不用选主元. 由21rrr rk k a l ==∑ 可以看出||1,2,,.rk l k r ≤=这表明中间量rk l得以控制，因此不会产生由中间量放大使计算不稳定的现象. Cholesky 分解的缺点：需要作开方运算. 改进的Cholesky 分解：改为使用分解T A LDL =即11121121121221222121111n n n n n n n n nn a a a d l l l d a a a l l d a a a ⎡⎤ 1 ⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥ 1 1 ⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥ ⎢⎥⎢⎥⎢⎥ ⎢⎥ ⎣⎦⎣⎦⎣⎦ 2n l ⎡⎤⎢⎥ ⎢⎥⎢⎥ ⎢⎥ 1⎣⎦其中21ˆl 1ˆn l 2ˆn l ˆnn l 1ˆn u12111()/r r rr rk k k r ir ir ik k rk rk d a l d l a l d l d-=-=⎧=-⎪⎪⎨⎪=-⎪⎩∑∑，1,2,,;1,2,,.r n i r r n ==++Cholesky 分解方法或平方根法：应用Cholesky 分解可将Ax b =分解为两个三角形方程组T Ly b L x y ⎧= ⎪⎨= ⎪⎩分别可解得111111/,()/.i i i ik k ii k y b l y b l y l i n -=⎧=⎪⎨=-, =2,3,,⎪⎩∑和1/,()/1,.n n nn n i i ki k ii k i x y l x y l x l i n n =+⎧=⎪⎨=-, =--2,,2,1⎪⎩∑改进的Cholesky 分解方法或改进的平方根法：应用改进的Cholesky 分解，将方程组Ax b =分解为下面两个方程组1,,T Ly b L x D y -= ⎧⎨= ⎩同理可解得1111,,2,3,,.i i i ik k k y b y b l y i n ==⎧=⎪⎨=- =⎪⎩∑和1/,/,1,2,,2,1.n n n n i i i ki k k i x y d x y d l x i n n =+⎧=⎪⎨=- =--⎪⎩∑ 4．解三对角方程组的追赶法若()ij n n A a ⨯=满足1||||,1,2,,.nii ij j j ia a i n =≠> =∑则称A 为严格对角占优矩阵.若A 满足1||||,1,2,,.nii ij j j ia a i n =≠≥ =∑且其中至少有一个严格不等式成立，则称A 为弱对角占优矩阵.现在考虑Ax d = 的求解，即11112222211111n n n n n n n n n b c x d a b c x d a b c x d d a b x -----⎡⎤⎡⎤⎡⎤ ⎢⎥⎢⎥⎢⎥ ⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥ = ⎢⎥⎢⎥⎢⎥ ⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦ 系数矩阵A 满足条件11||||0,||||||,,0,2,3,, 1.||||0,i i i i i n n b c b a c a c i n b a ⎧>>⎪≥+ ≠=-⎨⎪>>⎩采用Crout 分解方法11112222221111n n n n n n n b c a b c a b c a b βαβγαγα---⎡⎤ ⎡⎤⎢⎥ 1 ⎢⎥⎢⎥ ⎢⎥⎢⎥ = ⎢⎥⎢⎥ ⎢⎥ ⎢⎥ ⎢⎥⎢⎥⎣⎦ ⎣⎦ 1n β-⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥1 ⎢⎥⎢⎥ 1 ⎣⎦其中，,,i i i αβγ为待定系数.比较上式两边可得到111111,;,,2,3,,;,2,3,, 1.i i i i i i i i i b c a b i n c i n ααβγγβααβ-= == =+ == =-进而可导出1111111,2,3,,.,/,,2,3,,./(),2,3,, 1.i i i i i i ii i i i a i n b c b b i n c b i n γαβααββαβ--⎧= =⎪= =⎪⎨=- =⎪⎪=- =-⎩由此可看出，真正需要计算的是(1,2,,1)i n β=- ，而i α可由,i i b a 和1i β-产生.因此，实现了A 的Crout 分解后，求解Ax d =就等价于解方程组Ly dUx y =⎧⎨=⎩从而得到解三对角方程组的追赶法公式： (1) 计算i β的递推公式：1111/,/(),2,3,, 1.i i i i i c b c b i n ββαβ-⎧=⎪⎨=- =-⎪⎩(2) 解方程组Ly d =：11111/()/(),2,3,,.i i i i i i i y d b y d a y b a i n β--⎧=⎪⎨=-- =⎪⎩(3) 解方程组Ux y =：1,1,2,,2,1.n n i i i i x y x y x i n n β+⎧=⎪⎨=- =--⎪⎩追赶法的乘除法次数是66n -次.将计算121n βββ-→→→ 及12n y y y →→→ 的过程称之为“追”的过程，将计算方程组Ax d =的解121n n x x x x -→→→→ 的过程称之为“赶”的过程.四、迭代法将Ax b =改写为一个等价的方程组 x Bx k =+建立迭代公式 (1)(),0,1,2,.i i x Bx k i +=+ =称矩阵B 为迭代矩阵.【定义】如果对固定的矩阵B及向量k，对任意初始猜值向量(0)x ，迭代公式(1)()i i +()i()*lim i i x x →+∞=成立，其中*x 是一确定的向量，它不依赖于(0)x 的选取.则称此迭代公式是收敛的，否则称为发散的.如果迭代收敛，则应有**,x Bx k =+1. 收敛性()()*,0,1,2,i i x x i ε=- =为第i步迭代的误差向量.则有(1)(1)*()*()(),0,1,2,.x x B x x B i εε++=-=-==所以，容易推出()(0),0,1,2,,i i B i εε= =其中，(0)(0)*xxε=-为初始猜值的误差向量.设n nB K ⨯∈，lim 0i i B →+∞=⇔ ()1B ρ<.迭代法收敛基本定理：下面三个命题是等价的 (1) 迭代法(1)()i i x Bx k +=+收敛；(2)()1B ρ<；(3) 至少存在一种矩阵的从属范数⋅，使1B <注：当条件()1B ρ<难以检验时，用1B 或B ∞等容易求出的范数，检验11B <或1B∞<来作为收敛的充分条件较为方便.常用迭代法如下. 2．Jacob 迭代考察线性方程组Ax b =，设A 为非奇异的n 阶方阵，且对角线元素0ii a ≠(1,2,,)i n = .此时，可将矩阵A 写成如下形式A D L U =++,1122(,,,)nn D diag a a a = ，21313212000n n a L a a a a ⎡⎤⎢⎥ ⎢⎥⎢⎥= ⎢⎥ ⎢⎥⎢⎥ 0 ⎢⎥⎣⎦ ,12131232000n n a a a a a U ⎡⎤ ⎢⎥ ⎢⎥⎢⎥= 0 ⎢⎥ ⎢⎥⎢⎥ ⎢⎥⎣⎦ ,建立Jacobi 迭代公式(1)1()1(),i i x D L U x D b +--=-++迭代矩阵11()J B D L U I D A --=-+=-J B 的具体元素为112111122122221200n n J n n nn nn a a a a a a B a a a a a a ⎡⎤ - -⎢⎥⎢⎥⎢⎥- - ⎢⎥=⎢⎥⎢⎥ ⎢⎥⎢⎥- - 0 ⎢⎥⎣⎦ Jacobi 迭代法的分量形式如下1(1)()()111(),j n i i i jj jm m jm m m m j jj xb a x a x a -+==+=--∑∑1,2,,;0,1,2,.j n i = =3．Gauss-Seidel 迭代容易看出，在Jacobi 迭代法中，每次迭代用的是前一次迭代的全部分量()(1,2,,)i jx j n = .实际上，在计算(1)i j x +时，最新的分量(1)(1)(1)121,,,i i i j x x x +++- 已经算出，但没有被利用.事实上，如果Jacobi 迭代收敛，最新算出的分量一般都比前一次旧的分量更加逼近精确解，因此，若在求(1)i j x+时，利用刚刚计算出的新分量(1)(1)(1)121,,,i i i j x x x+++- ，对Jacobi 迭代加以修改，可得迭代公式1(1)(1)()111(),j ni i i jj jm m jm m m m j jj xb a x a x a -++==+=--∑∑1,2,,;0,1,2,.j n i = =矩阵形式(1)1()1()(),0,1,2,.i i x D L Ux D L b i +--=-++-+=1()G B D L U -=--+注：（1）两种迭代法均收敛时，Gauss-Seidt 迭代收敛速度更快一些.（2）但也有这样的方程组，对Jacobi 迭代法收敛，而对Gauss-Seidel 迭代法却是发散的. 【例2.4】分别用Jacobi 迭代法和Gauss-Seidel 迭代法求解下面的方程组121232342,46,4 2.x x x x x x x ⎧- =⎪-+-=⎨⎪-+=⎩初始猜值取0(0,0,0)x =. 【解】 Jacobi 迭代公式为(1)()12(1)()()213(1)()321(2),41(6),0,1,2,41(2),4i i i i i i i x x x x x i x x +++⎧=+⎪⎪⎪=++=⎨⎪⎪=+⎪⎩迭代计算4次的结果如下 (1)(2)(3)(4)(0.5,1.5,0.5),(0.875,1.75,0.875),(0.938,1.938,0.938),(0.984,1.969,0.984).T T T T x x x x ====Gauss-Seidel 迭代公式为(1)()12(1)(1)()213(1)(1)321(2),41(6),0,1,2,41(2),4i i i i i i i x x x x x i x x +++++⎧=+⎪⎪⎪=++=⎨⎪⎪=+⎪⎩迭代计算4次的结果如下(1)(2)(3)(4)(0.5,1.625,0.9063),(0.9063,1.9532,0.9883),(0.9883,2.0,0.9985),(0.9985,1.999,0.9998).T T T T x x x x ====从这个例子可以看到，两种迭代法作出的向量序列(){}i x 逐步逼近方程组的精确解*(1,2,1)T x =,而且Gauss-Seidel 迭代法收敛速度较快.一般情况下，当这两种迭代法均收敛时，Gauss-Seidt 迭代收敛速度更3．超松弛迭代法为了加快迭代的收敛速度，可将Gauss-Seidel 迭代公式改写成1(1)()(1)()11(),j ni i i i jjj jm m jm m m m jjj xx b a x a x a -++===+--∑∑ 1,2,,;0,1,2,.j n i = =并记1(1)(1)()11(),j ni i i jj jm m jm m m m jjj rb a x a x a -++===--∑∑称 (1)i j r + 为 1i + 步迭代的第 j 个分量的误差向量.当迭代收敛时，显然有所有的误差向量(1)0(),1,2,,.i j r i j n +→→∞=为了获得更快的迭代公式，引入因子R ω∈，对误差向量 (1)i j r + 加以修正，得超松弛迭代法（简称SOR 方法）(1)()(1),0,1,2,.i i i j j j x x r i ω++=+ =即1(1)()(1)()1(),j ni i i i jjj jm mjm m m m jjjxx b a xa x a ω-++===+--∑∑1,2,,;0,1,2,.j n i = =适当选取因子ω，可望比Gauss-Seidel 迭代法收敛得更快.称ω为松弛因子.特别当1ω=时，SOR 方法就是Gauss-Seidel 迭代法.写成矩阵向量形式(1)1()1()[(1)](),j i x D L D U x D L b ωωωωω+--=+--++0,1,2,.i =迭代矩阵为1()[(1)].B D L D U ωωωω-=+--实际计算时，大部分是由计算经验或通过试算法来确定opt ω的近似值.所谓试算法就是从同一初始向量出发，取不同的松驰因子ω迭代相同次数(注意：迭代次数不应太少)，然后比较其相应的误差向量()()i i r b Ax =-(或()(1)i i x x --)，并取使其范数最小的松弛因子ω作为最佳松弛因子opt ω的近似值.实践证明，此方法虽然简单，但往往是行之有效的. 4．迭代收敛其它判别方法：用迭代法收敛基本定理来判断收敛性时，当n 较大时，迭代矩阵的谱半径计算比较困难，因此，人们试图建立直接利用矩阵元素的条件来判别迭代法的收敛定理. （1）若方程组Ax b =中的系数矩阵A 是对称正定阵，则 Gauss-Seidel 迭代法收敛. 对于SOR 方法，当02ω<< 时迭代收敛（2）若A 为严格对角占优阵，则解方程组 Ax b = 的Jacobi 迭代法，Gauss －Seidel 迭代法均收敛. 对于SOR 方法，当01ω<< 时迭代收敛.【例2.5】设线性方程组为121221,32,x x x x ⎧+=-⎪⎨+=⎪⎩建立收敛的Jacobi 迭代公式和Gauss －Seidel 迭代公式. 【解】对方程组直接建立迭代公式，其Jacobi 迭代矩阵为0230J B -⎡⎤=⎢⎥- ⎣⎦，显见谱半径()1J B ρ=>，故Jacobi 迭代公式发散.同理Gauss －Seidel 迭代矩阵为0206G B -⎡⎤=⎢⎥ ⎣⎦，谱半径()61G B ρ=>，故Gauss －Seidel 选代公式也发散. 若交换原方程组两个方程的次序，得一等价方程组121232,21,x x x x ⎧+=⎪⎨+=-⎪⎩其系数矩阵显然对角占优，故对这一等价方程组建立的Jacobi 迭代公式，Gauss －Seidel 迭代公式皆收敛. （3）SOR 方法收敛的必要条件是 02ω<<【定理2.5】如果A 是对称正定阵，且02ω<<，则解Ax b =的SOR 方法收敛.注：当(0,2)ω∈ 时，并不是对任意类型的矩阵A ，解线性方程组Ax b =的SOR 方法都是收敛的.当SOR 方法收敛时，通常希望选择一个最佳的值opt ω使SOR 方法的收敛速度最快.然而遗憾的是，目前尚无确定最佳超松弛因子opt ω的一般理论结果.实际计算时，大部分是由计算经验或通过试算法来确定opt ω的近似值.所谓试算法就是从同一初始向量出发，取不同的松驰因子ω迭代相同次数(注意：迭代次数不应太少)，然后比较其相应的误差向量()()i i r b Ax =-(或()(1)i i x x --)，并取使其范数最小的松弛因子ω作为最佳松弛因子opt ω的近似值.实践证明，此方法虽然简单，但往往是行之有效的.【例2.6】求解线性方程组Ax b =，其中10.3000900.308980.30009100.4669110.274710.30898A - -- -0.46691 0= - -- 00.274711(5.32088,6.07624,8.80455,2.67600).T b ⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥ - ⎣⎦ =-分别利用Jacobi 迭代法，Gauss －Seidel 迭代法，SOR 迭代法求解. 【解】其结果列入下表中，方程组精确解(五位有效数字)为*(8.4877,6.4275, 4.7028,4.0066).T x =-Jacobi 迭代法计算结果i()1i x()2i x ()3i x ()4i x ()2||||i r0 012.3095 1 5.3209 6.0762 －8.8046 2.6760 5.3609 27.97113.5621 －5.2324 1.90143.631820 8.4872 6.4263 －4.7035 4.0041 0.0041 218.48606.4271 －4.7050 4.0063 0.0028Gauss-Seidel 迭代法计算结果i()1i x()2i x()3i x()4i x()2||||i r0 012.3095 1 5.3209 7.6730 －5.2220 2.8855 3.6202 28.51506.1933 －5.1201 3.90040.49098 8.4832 6.4228 －4.7064 4.0043 0.0078 98.48556.4252－4.70554.00550.0038SOR 迭代法计算结果（1.16ω=）i()1i x()2i x()3i x()4i x()2||||i r0 012.3095 1 6.1722 9.1970 －5.2320 3.6492 3.6659 29.69416.1177 －4.8999 4.43351.33136 8.4842 6.4253 －4.7005 4.4047 0.0051 78.48686.4288－4.70314.00650.0016计算结果表明，若求出精确到小数点后两位的近似解，Jacobi 迭代法需要21次，Gauss －Seidel 迭代法需要9次，而SOR 迭代法(选松弛因子 1.16ω=)仅需要7次，起到加速作用.5．误差分析【定理2.6】设 *x 是方程 Ax b = 的惟一解，v ⋅ 是某一种向量范数，若对应的迭代矩阵其范数1v B <，则迭代法(1)(),0,1,2,.i i xBx k i +=+ = 收敛，且产生向量序列(){}i x 满足()*()(1)||||||||||||1||||i i i vv vvB x x x x B --≤--()*(1)(0)||||||||||||1||||i i vv vvB x x x x B -≤--【证明】由迭代收敛基本定理的(3)知，迭代法(1)(),0,1,2,.i i x Bx k i +=+ =收敛到方程的解*x .于是，由迭代公式立即得到(1)*()*(1)()()(1)(),().i i i i i i x x B x x x x B x x ++--=--=-为书写方便把v 范数中v 略去，有估计式(1)*()*||||||||||||,i i x x B x x +-≤⋅-(1)()()(1)||||||||||||.i i i i x x B x x +--≤⋅-再利用向量范数不等式||||||||||||x y x y -≥-于是得第一个不等式()(1)(1)()()*(1)*()*||||||||||||||||||||(1||||)||||,i i i i i i i B x x x x x x x x B x x -++ -≥-≥--- ≥--再反复递推即第二个不等式.注：（1）若事先给出误差精度ε，利用第二个不等式可得到迭代次数的估计(1)(0)(1||||)ln ln ||||||||v v v B i B x x ε⎡⎤->⎢⎥-⎣⎦ （2）在||||v B 不太接近1的情况下，由第一个不等式，可用()(1)||||i i v x x ε--<作为控制迭代终止的条件，并取 ()i x 作为方程组 Ax b = 的近似解.但是在||||v B 很接近1时，此方法并不可靠.一般可取1,2,v =∞或F .【例2.7】用Jacobi 迭代法解方程组123123123202324,812,231530.x x x x x x x x x ⎧++=⎪++=⎨⎪-+=⎩问Jacobi 迭代是否收敛?若收敛，取(0)(0,0,0)T x =，需要迭代多少次，才能保证各分量的误差绝对值小于610-?【解】 Jacobi 迭代的分量公式为(1)()()123(1)()()213(1)()()3121(2423)201(12),0,1,2,81(3022),15i i i i i i i i i x x x x x x i x x x +++⎧=--⎪⎪⎪=-- =⎨⎪⎪=-+⎪⎩Jacobi 迭代矩阵J B 为130102011088210155J B ⎡⎤ - -⎢⎥⎢⎥⎢⎥=- -⎢⎥⎢⎥⎢⎥- ⎢⎥⎣⎦，由5251||||max ,,1208153J B ∞⎧⎫==<⎨⎬⎩⎭知，Jacobi 迭代收敛. 因设(0)(0,0,0)Tx =，用迭代公式计算一次得(1)(1)(1)12363,, 2.52x x x = = =而(1)(0)|||| 2.x x ∞-=于是有6110(1)13ln ln 13.23i -⎡⎤⋅-⎢⎥>=⎢⎥⎢⎥⎣⎦所以，要保证各分量误差绝对值小于610-，需要迭代14次.【例2.8】用Gauss －Seidel 迭代法解例2.11中的方程组，问迭代是否收敛?若收敛，取(0)(0,0,0)Tx =，需要迭代多少次，才能保证各分量误差的绝对值小于610-?【解】 Gauss －Seidel 迭代矩阵G B 为102403601()03025524000G B D L U - - ⎡⎤⎢⎥=-+= -⎢⎥⎢⎥ 38 -3⎣⎦显然1||||14G B =<，所以迭代收敛. Gauss －Seidel 迭代分量公式为(1)()()123(1)(1)()213(1)(1)(1)3121(2423),201(12),0,1,2,81(3022),15i i i i i i i i i x x x x x x i x x x ++++++⎧=--⎪⎪⎪=-- =⎨⎪⎪=-+⎪⎩因取(0)(0,0,0)T x =，故迭代一次得(1)(1)(1)1231.2, 1.35, 2.11x x x = = =于是有(1)(0)|||| 2.11x x ∞-=，计算得6110(1)14ln ln 10.2.114i -⎡⎤⋅-⎢⎥>=⎢⎥⎢⎥⎣⎦所在，要保证各分量误差绝对值小于610-，需要迭代11次.。

高等工程数学试题及参考答案-工程硕士

中南大学工程硕士“高等工程数学”考试试卷考试日期：2010年 4 月日时间110分钟注：解答全部写在答题纸上一、填空题(本题24分，每小题3分) 1. 若方程0)(=x f 可表成)(x xϕ=，且在[,]a b 内有唯一根*x ，那么)(x ϕ满足，则由迭代公式)(1n n x x ϕ=+产生的序列{}n x 一定收敛于*x 。

（)(x ϕ满足：1()[,]x C a b ϕ∈,且[,]x a b ∀∈有()[,]x a b ϕ∈, '()1x L ϕ≤<；）2. 已知二元非线性函数221122120()24,(2,2)Tf x x x x x x x X =-++-=，该函数从X 0 出发的最速下降方向为（最速下降方向为：()4,2Tp =-）； 3．已知二元非线性函数221122120()24,(2,2)Tf x x x x x x x X =-++-=，该函数从X 0 出发的Newton 方向为（Newton 方向为： ()2,0Tp =-）； 4．已知)(x f y =在区间],[b a 上通过点(,),0,1,2,,i i x y i n =L ，则其三次样条插值函数)(x S 是满足（（1）在每个小区间是次数不超过3次的多项式，（2）在区间[,]a b 上二阶导数连续，（3）满足插值条件(),0,1,2,,i i S x y i n ==L ）；5．设某个假设检验问题的拒绝域为W ，且当原假设H 0成立时，样本值12(,,,)n X X X L 落入W 的概率为0.15，则犯第一类错误的概率为________（0.15）；6．在实际问题中求某参数的置信区间时，总是希望置信水平愈大愈好，而置信区间的长度愈短愈好。

但当增大置信水平时，则相应的置信区间长度总是变长； 7．取步长2.0=h ，解]1,0[,1)0(2'∈⎩⎨⎧=-=x y yx y 的Euler 法公式为：（1(2)0.60.2,0,1,2,,5n n n n n n y y h x y y x n +=+-=+=L ）；8．对实际问题进行建模求解时可能出现的误差有：（模型误差，观测误差，方法误差，舍入误差。

《高等工程数学》习题三参考答案

1i n
1 P{ X 1 x, X 2 x, , X n x} 1 (1 P{ X x}) n 1 (1 F ( x)) n ;
因为 X ( n ) max X i ，所以 FX ( n ) ( x) P{ X ( n ) x} P{ X 1 x} P{ X n x} F ( x ) 。
11. 解：因 X ~ N (80,20 2 ) ，样本容量为 100，所以 X ~ N (80,4) ，
3
P{ X 80 3} P{
X 3 2
3 3 } 2(1 ( )) 2 * (1 - normcdf(3/2)) 0.1336 。 2 2
3 )， 10
12. 解：设 X 1 , X 2 , , X 10 和 Y1 , Y2 , , Y15 为 N ( 20,3) 两独立样本，则 X ~ N (20,
2
2 ( n) ， X
X1 ~ t ( n) ， X2 / n
所以 X
2
X1 /1 ~ F (1, n) 。 X2 / n
9. 解：MATLAB 命令为（1）norminv(0.99)；（2）norminv(0.04)；（3）chi2inv(0.975,15);(4) chi2inv(0.025,15);(5) chi2inv(0.95,50);(6) chi2inv(0.95,100);(7) tinv(0.975,19);(8) tinv(0.975,99); (9) finv(0.95,2,6);(10) finv(0.05,3,40);(11) finv(0.05,2,6);(12) finv(0.01,3,40) 10.解：因 X ~ N (1,4) ，样本容量为 16，所以 X ~ N (1,

《高等工程数学》习题一参考答案

也是一组标准正交基。 19. 解：系数矩阵 A
2 1 1 1 3 1 0 0 1 4 ，可得基础解系为 1 1 1 0 1 0 1 1 1 5
f1 (0,1,1,0,0) , f 2 (1,1,0,1,0) ， f 3 (4,5,0,0,1) ，Schmidt 正交化得，
1
13.按 P21 欧氏空间定义 2.1，逐条验证， 1）不满足第（2 ）条，（4）条，故不是欧氏空间; 不满足第（4）条，故不是欧氏空间；3）都满足，故是欧氏空间。 14. 按 P21 欧氏空间定义 2.1，逐条验证，都满足，故是欧氏空间。 15. 设向量 ( x1 , x2 , x3 , x4 ) 与三个向量正交，则有
所以对两组基有相同坐标的非零向量可取为 (c, c, c,c)(c 0). 5. 由第 7 页子空间定义可得，1）向量满足加法和数乘封闭，是子空间；2）向量不满足加法或数乘封闭，故而不是子空间。注：从几何上看，子空间过原点，而不过原点的都不是。 6. 两个向量组生成相同子空间的充分必要条件是这两个向量组等价，即可以互相线性表示。解：因对应分量不成比例，故 α1 (1,1,0,0), α2 (1,0,1,1) ， β1 (1,1,0,0), β2 (1,0,1,1) 线性
2
T1T2 ( x1 , x2 ) T1[T2 ( x1 , x2 )] T1 ( x1 , x2 ) ( x2 , x1 ) T2T1 ( x1 , x2 ) T2 [T1 ( x1 , x2 )] T2 ( x2 , x1 ) ( x2 , x1 )
11.略。 12. 解：1）因为 T ( x1 , x2 , x3 ) ( 2 x1 x2 , x2 x3 , x1 ) ，按照 P18 (1.21)，可知

《高等工程数学》科学出版社版习题答案(第三章)

《高等工程数学》――科学出版社版习题答案（第三章）（此习题答案仅供学员作业时参考。

因时间匆忙，有错之处敬请指正，谢谢！）（联系地址：yangwq@ ）P501．自己验证范数的三个条件2．自己验证范数的三个条件3．（1）122222212111121()||||(||)||||||||||||||||(||1),||||||(11)||n n n n nk k i j k k k i j k n k k T nx x x x C x x x x x x x Cauchy Schwartz x x x I I x x I R x ==≠===∈==+∙≥=-=∙=<>≤∙==∈∑∑∑∑∑设，，...，，则有－－(*)另由不等式，有－－(**)其中，，...，1所以由(*)和(**)式有：212||||||||x x ≤≤(（2）121111211111()||||max ||||||||()||max ||||||||max ||||||||||||||||n n nk k k n k n i k k nn i k k n i x x x x C x x x x x x x x x x x x n x n x x x n x ∞≤≤=≤≤∞≤≤=∞∞=∈==≤=≤=≤=≤≤∑∑ 设，，...，，则有－－(*)另外对，，...，的任一分量有所以有：－－(**)所以由(*)和(**)式有：（3）12211212()||||max ||||||()||max ||||||||||||n n k k n n i k k nx x x x C x x x x x x x x x x x x ∞≤≤≤≤∞=∈==≤==≤=≤=设，，...，，则有－－(*)因对，，...，的任一分量有所以有：－－(**)所以由(*)和(**)式有：2||||||||x x ∞∞≤≤4．已知1321i A i -⎡⎤=⎢⎥+⎣⎦试求第12|||| ||||||||A A A A ρ∞，，，（）解：12222||||max{2||||max{3121365531215511655||176650(16)(1)5511||||413||(1)162421(H H A A i i i A A i i i i I A A i A iI A i λλλλλλλλλλλλλρ∞====+-+⎡⎤⎡⎤⎡⎤==⎢⎥⎢⎥⎢⎥-+-⎣⎦⎣⎦⎣⎦----==-+-=---+-=-+--==-+-=-----因所以)1A ==+5．证明：(1)211H U U U II U =因是酉矩阵，所以＝而单位矩阵的特征值为，所以(2) 222222)))))H H H H H H H H H H H H H H H H H H H H H H H H H U U U IUA UA A U UA A AUA A AU AU U A AU U A A UA A AU AU AU A U AU U AU U A UU AU U A A UA A U AU U AU U AUA =========因是酉矩阵，所以＝（）（所以（）（（）即矩阵与（）（相似，所以有相同的特征值即（）（（）即矩阵与（）（相似，所以有相同的特征值即6． ||||=1||||=1111-1||||=max||||=max||||=11||||=||||||||||||||||||||e e I Ie e I A A A A A A ---=≤∴≥7． (1)证明：假设I －A 不可逆，则|I-A|=0，即1是A 的特征值，所以 ()1()()1A A A A A A ρρρ≥≤<又因为对的任一范数，都有所以由题设知矛盾，所以I －A 可逆(2) 由||||||||||||||||1||||||||||||||||||||11||||1||||1||||I A I A I I A I A A II A I I A AI A I I A A I I A A I A A I A I A A A I A A --⇒--⇒-+-∴-=+-≤+-≤+-∴---≤<-≤- -1-1-1-1-1-1-1-1-1-1-1-1()()=()-()=()=()()()()()()()由得() 得证明：1||||0.9 ()<1lim k k A A A O ρ→∞=∴∴=（2）2||(2)() ()=2||1|lim 2k k c c I B cc c c B c c c c B O λλλλλρλ→∞---=--=-+∴--∴<= 当|时，9．(1) 解： 21334(4)(1)22()41A A λλλλλλρ--=--=-+--∴=>故发散(2) 因为收敛半径为：R=5，所以收敛10．解： 1210.10.80.60.30.90.534140.90771300.511773490111425371773311311702835377k k A A A A λλ∞=⎡⎤=⎢⎥⎣⎦==-⎡⎤⎢⎥-⎡⎤⎡⎤=⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦⎢⎥-⎢⎥⎣⎦⎡⎤⎡⎤⎡⎤⎢⎥-⎡⎤⎢⎥⎢⎥==⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎢⎥-⎣⎦⎣⎦⎢⎥⎣⎦∑设的特征值为，，所以(1) 222sin 2cos(2)sin 2sin()sin cos()sin cos 2sin(2)cos 2cos()cos sin()cos t tt At t t t e te e e e te e t t t t At t t t t t t t t At t t t t ⎡⎤⎢⎥⎢⎥=⎢⎥⎢⎥⎣⎦⎡⎤⎢⎥⎢⎥=⎢⎥⎢⎥⎣⎦-⎡⎤⎢⎥⎢⎥=⎢⎥-⎢⎥⎣⎦(2) 222sin 2cos(2)sin 2sin()sin cos()sin cos 2sin(2)cos 2cos()cos sin()cos tt t At t t t e te e e e te e t t t t At t t t t t t t t At t t t t ⎡⎤⎢⎥⎢⎥=⎢⎥⎢⎥⎣⎦⎡⎤⎢⎥⎢⎥=⎢⎥⎢⎥⎣⎦⎡⎤⎢⎥-⎢⎥=⎢⎥⎢⎥-⎣⎦12．解： A 的特征值为：－1，1，22221166110221102211sin 2(2sin 2sin )(sin 22sin )33sin()00sin 0sin 021cos 2(cos 2cos )(co 33cos()t t t t t t t At t t t t t t t t ee e e e e e e e e e e e e e e tt t t t At t t t t t At ------⎡⎤⎢⎥⎢⎥⎢⎥=+-⎢⎥⎢⎥⎢⎥-+⎢⎥⎣⎦⎡⎤--⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎣⎦-=（4-3-）（2-3+）（）（）（）（）s 2cos )0cos 000cos t t t t ⎡⎤-⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦13．解： A 的特征值为：1，1，42ln 1110240.5 1.50.50.50.5 2.5A -⎡⎤⎢⎥=-⎢⎥⎢⎥--⎣⎦-⎡⎤⎢⎥=--⎢⎥⎢⎥--⎣⎦-ln4-2ln4+6ln4-1-ln4+33。

高等工程数学练习题及答案

高等工程数学练习题（2012年12月16）1. n 位男士和n 位女士排成一行，要求男女相间，求有多少种不同的排法？把n 个男、n 个女分别进行全排列，然后按乘法法则放到一起，而男女分别在前面,应该再乘2,即方案数为2·(n!) 个2. n 个人围圆圈坐下做游戏，求不同的坐法数？若某两人不愿坐在一起，有多少种不同的坐法? 若有3人总是坐在一起，又有多少种不同的坐法？ A: Q(n ,n)=(n-1)!;B: Q(n ,n)- Q(n-1,n-1) *2！=(n-1)!- (n-2)! *2！ C: Q(n-2,n-2) *3！= (n-3)! *3！3. 书架上有一部24卷的百科全书，现要从中取出5本，使得没有两本书是连续的，问有多少种不同的取法？ C(24-5+1,5)=C(20,5)4. 设{1,2,3,,(1)}.A n =+ （1）证明最大元素恰为j 的子集的个数是12j -；（2）证明：2112222 1.m m +++++=-A 、最大元素恰为j 的子集的个数，相当于前j-1个元素，每个元素出现或不出现的情况构成的所有子集的数量，每个元素出现或不出现2种可能，因此j-1个2相乘即为所有的情况，即12j -。

B:等比数列a1=2,q=2右侧为1+(2*（1-2^m ）/1-2)=2^(m+1)-2+1=2^(m+1)-1=左侧5. 证明等式: 22222012n n n n n n n ⎛⎫⎛⎫⎛⎫⎛⎫⎛⎫++++= ⎪ ⎪ ⎪ ⎪ ⎪⎝⎭⎝⎭⎝⎭⎝⎭⎝⎭C(n 0)*C(n 0)=C(n 0)*C(n n); C(n 1)*C(n 1)=C(n 1)*C(n n-1); ……C(n k)*C(n k)=C(n k)*C(n n-k) ,K=0~nC(n k)相当于（0 0）到直线{（n 0）(0 n)}上的某点（n-k,k ）的路径 C(n n-k) 相当于直线{（n 0）(0 n)}上的某点（n-k,k ）到（n n ）的路径根据乘法原理 C(n k)*C(n n-k)相当于（0 0）点通过直线{（n 0）(0 n)}上的某点（n-k,k ）到（n n ）的路径左侧为（0 0）点通过直线{（n 0）(0 n)}上所有点到（n n ）的路径相加由于（0 0）点到达（n n ）的所有路径均通过直线{（n 0）(0 n)}，所以根据加法原理左侧为（0 0）点到（n n ）的所有路径即等于（2n n ） 6. 证明恒等式： 11221011220n r n r n r n m r mn r m m m m m -+-+-+++⎛⎫⎛⎫⎛⎫⎛⎫⎛⎫⎛⎫⎛⎫⎛⎫⎛⎫++++=⎪ ⎪ ⎪ ⎪⎪ ⎪ ⎪ ⎪ ⎪--⎝⎭⎝⎭⎝⎭⎝⎭⎝⎭⎝⎭⎝⎭⎝⎭⎝⎭（-r-1,0）到(-1,i)路径为c(r+i,i)(-1,i)到(0,i)路径为1(0,i)到(n-m,m)路径为c(n-i.m-i)根据乘法原理，c(r+i,i) c(n-i.m-i)为（-r-1,0）经过(-1,i)到(0,i)再到(n-m,m)点的路径左侧为i 取0至m ，（-r-1,0）经过(-1,i)到(0,i)再到(n-m,m)点的路径之和，右侧围（-r-1,0）到(n-m,m)点的路径左右相等7. 求不定方程123(,)n x x x x r n r Z +++++=∈的非负整数解的个数；设r n ≥，求不定方程的正整数解的个数. C （n+r-1,r ）C(n+r-n-1,r-n)=c(r-1,r-n),相当于每盒先放一个球，球数量变成r-n ，再求解。

工程数学试题及答案高级专

工程数学试题及答案高级专工程数学试题及答案（高级专科）一、单项选择题（每题3分，共30分）1. 极限的定义中，当x趋近于a时，f(x)的极限为L，意味着（）。

A. f(x) = LB. |f(x) - L| < ε，对任意的ε > 0，存在δ > 0，使得0 < |x - a| < δ时成立C. |f(x) - L| = 0D. f(x) ≠ L答案：B2. 函数f(x) = x^2在x=0处的导数为（）。

A. 0B. 1C. 2D. -1答案：B3. 以下哪个函数是奇函数？（）A. f(x) = x^2B. f(x) = x^3C. f(x) = x^4D. f(x) = x答案：B4. 以下哪个函数是周期函数？（）A. f(x) = e^xB. f(x) = sin(x)C. f(x) = ln(x)D. f(x) = x^2答案：B5. 以下哪个积分是发散的？（）A. ∫(0, +∞) e^(-x) dxB. ∫(0, +∞) x^2 dxC. ∫(0, +∞) e^x dxD. ∫(0, +∞) 1/x dx答案：D6. 以下哪个是二阶常系数线性微分方程？（）A. y'' + 2y' + y = 0B. y'' + 2y' + 3y = 0C. y'' + y' + y = 0D. y'' + y' = 0答案：A7. 以下哪个是二阶偏导数？（）A. ∂^2f/∂x∂yB. ∂^2f/∂x^2C. ∂^2f/∂y^2D. ∂^2f/∂x∂y^2答案：A8. 以下哪个是线性方程组的解？（）A. {x=1, y=2}B. {x=0, y=0}C. {x=1, y=1}D. {x=2, y=3}答案：C9. 以下哪个是矩阵的特征值？（）A. λ = 1B. λ = 2C. λ = 3D. λ = 4答案：A10. 以下哪个是傅里叶级数的系数？（）A. a_nB. b_nC. c_nD. d_n答案：A二、填空题（每题4分，共20分）11. 函数f(x) = sin(x)在x=π/2处的导数为______。

大学工程数学试题及答案

大学工程数学试题及答案一、选择题（每题3分，共30分）1. 以下哪个选项是微分方程的解？A. \( y = e^x \)B. \( y = e^{-x} \)C. \( y = x^2 \)D. \( y = \ln(x) \)答案：A2. 矩阵的行列式值表示了什么？A. 矩阵的面积B. 矩阵的体积C. 矩阵的旋转角度D. 矩阵的缩放因子答案：D3. 以下哪个是线性代数中的基本概念？A. 微分B. 积分C. 向量空间D. 极限答案：C4. 傅里叶变换用于解决什么问题？A. 微分方程B. 积分方程C. 信号处理D. 线性代数答案：C5. 欧拉公式 \( e^{ix} = \cos(x) + i\sin(x) \) 中，\( i \) 代表什么？A. 虚数单位B. 矩阵C. 行列式D. 向量答案：A6. 以下哪一项是拉普拉斯变换的基本性质？A. 线性性质B. 微分性质C. 积分性质D. 反演性质答案：A7. 泰勒级数展开是用于什么目的？A. 近似计算B. 精确计算C. 矩阵计算D. 向量计算答案：A8. 以下哪个函数是周期函数？A. \( y = x^2 \)B. \( y = e^x \)C. \( y = \sin(x) \)D. \( y = \ln(x) \)答案：C9. 以下哪一项是偏微分方程的解？A. \( u(x, y) = x^2 + y^2 \)B. \( u(x, y) = e^{x+y} \)C. \( u(x, y) = \ln(x+y) \)D. \( u(x, y) = \sin(x)\cos(y) \)答案：D10. 以下哪个选项是复数的性质？A. 可加性B. 可乘性C. 可除性D. 所有选项答案：D二、填空题（每题4分，共20分）1. 如果 \( f(x) = x^3 - 6x^2 + 11x - 6 \)，则 \( f'(x) \) 等于 _______。

答案：\( 3x^2 - 12x + 11 \)2. 矩阵 \( A = \begin{bmatrix} 1 & 2 \\ 3 & 4 \end{bmatrix} \) 的行列式 \( \det(A) \) 等于 _______。

工程数学试题及答案升专科

工程数学试题及答案升专科工程数学试题及答案（升专科）一、单项选择题（每题2分，共20分）1. 极限的定义中，当x趋近于a时，f(x)的极限为L，那么对于任意的ε>0，存在一个δ>0，使得当0<|x-a|<δ时，|f(x)-L|<ε。

以下哪个选项正确描述了这个定义？（）A. 当x趋近于a时，f(x)趋近于LB. 当x趋近于a时，f(x)可能趋近于LC. 当x趋近于a时，f(x)一定趋近于LD. 当x趋近于a时，f(x)不一定趋近于L答案：C2. 以下哪个函数是偶函数？（）A. f(x) = x^2 + xB. f(x) = x^3 - xC. f(x) = x^2 - 2x + 1D. f(x) = x^3 + 3x^2答案：C3. 以下哪个积分是定积分？（）A. ∫(1/x)dxB. ∫(1/x)dx from 1 to 2C. ∫(1/x)dx from 0 to 1D. ∫(1/x)dx from -1 to 1答案：B4. 以下哪个级数是收敛的？（）A. 1 + 1/2 + 1/3 + 1/4 + ...B. 1 - 1/2 + 1/3 - 1/4 + ...C. 1 + 1/2^2 + 1/3^2 + 1/4^2 + ...D. 1 - 1/2^2 + 1/3^2 - 1/4^2 + ...答案：C5. 以下哪个矩阵是可逆的？（）A. [1 2; 3 4]B. [1 0; 0 1]C. [1 1; 1 1]D. [2 3; 4 6]答案：B6. 以下哪个方程是二阶线性微分方程？（）A. y'' + 3y' + 2y = 0B. y''' + 3y'' + 2y' = 0C. y' + 3y = 0D. y'' + 3y = 0答案：A7. 以下哪个函数的导数是y' = 2x？（）A. y = x^2 + 1B. y = x^3 + 1C. y = x^2D. y = x^3答案：C8. 以下哪个函数是周期函数？（）A. f(x) = x^2B. f(x) = sin(x)C. f(x) = e^xD. f(x) = ln(x)答案：B9. 以下哪个是二重积分？（）A. ∫∫(x^2 + y^2)dxdyB. ∫(x^2 + y^2)dxC. ∫∫(x^2 + y^2)dAD. ∫(x^2 + y^2)dy答案：C10. 以下哪个是偏导数？（）A. ∂f/∂xB. df/dxC. ∂f/∂yD. dy/dx答案：A二、填空题（每题2分，共20分）11. 函数f(x) = x^3 - 3x在x=0处的导数是________。

大学工程数学试题及答案

大学工程数学试题及答案一、选择题（每题2分，共10分）1. 函数f(x)=\sin(x)在区间[0, π]上是：A. 增函数B. 减函数C. 先增后减D. 先减后增答案：C2. 以下哪个选项是二阶导数的几何意义？A. 切线的斜率B. 函数的增减性C. 函数的凹凸性D. 函数的极值点答案：C3. 复数z=3+4i的模是：A. 5B. 7C. √7D. √5答案：A4. 矩阵A=[1 2; 3 4]的行列式是：A. -2B. 2C. 0D. 5答案：B5. 以下哪个选项是泰勒级数展开的公式？A. f(x) = f(a) + f'(a)(x-a) + f''(a)(x-a)^2/2!B. f(x) = f(a) + f'(a)(x-a) + f''(a)(x-a)^2/3!C. f(x) = f(a) + f'(a)(x-a) + f''(a)(x-a)^2/4!D. f(x) = f(a) + f'(a)(x-a) + f''(a)(x-a)^2/1!答案：A二、填空题（每题3分，共15分）6. 函数f(x)=x^3-3x+2的导数是______。

答案：3x^2-37. 曲线y=x^2在点(1,1)处的切线方程是______。

答案：y-1=2(x-1)8. 向量(2, -3)和(1, 2)的点积是______。

答案：-49. 矩阵A=[1 0; 0 2]的逆矩阵是______。

答案：[1 0; 0 1/2]10. 函数f(x)=e^x的不定积分是______。

答案：e^x + C三、解答题（每题10分，共30分）11. 求函数f(x)=x^2-4x+3在区间[1, 3]上的定积分。

答案：∫(x^2-4x+3)dx从1到3 = (1/3x^3 - 2x^2 + 3x) | 从1到3 = 012. 证明函数f(x)=x^3在R上是单调递增的。

高等工程数学I 试题(A)与答案(2019.11.23)

2019年高等工程数学试题答案一、（15分）设210120003⎛⎫⎪= ⎪ ⎪⎝⎭A ，计算()ρA 、225max =x Ax 及()2cond A 。

解：12321012001;3003λλλλλλλ---=--=⇒===-I A ()3ρ=A 2||||()3是正规矩阵ρ∴== A A A 2222515max 5max 5155==∞===x x xAx AA ()2331是正规矩阵∴== A cond A 二、（10分）讲述一下求解矩阵A 的最靠近*λ的特征值的思路、步骤。

答：**对使用逆幂法，求出其按模最小的特征值再加上。

λλ-A I 000u v =≠任取*11()max()k k k k k u A I v u v u λ--⎧=-⎪⎨=⎪⎩*1()max()k k kk k A I u v u v u λ-⎧-=⎪⎨=⎪⎩即**()A I P A I LUλλ--=对进行选列主元的三角分解有1max()k k k kk k k Ly PvUu y u v u -⎧⎪=⎪⎪∴=⎨⎪⎪=⎪⎩1max()max()k ik i u x v x λλ*⎧→⎪⎪⎨⎪→⎪⎩－有三、（18分）已知矩阵200226044-⎛⎫⎪= ⎪ ⎪⎝⎭A ，求P 使得1-P AP 为A 的Jordan 标准型，同时需要求出A 的Jordan 标准型。

解：200226044λλλλ+-=-----I A ()()23+28λλ=-D 211D D ==()()23+28λλ=-d 211d d ==初等因子：()()2+2 8，λλ-Jordan 标准形：2128-⎛⎫ ⎪=- ⎪ ⎪⎝⎭J 1123212,[]8令--⎛⎫⎪==-= ⎪ ⎪⎝⎭P AP J P p p p 11121223332[032]512[0]228[011]∴=-∴=-=-=-==TT TAp p p Ap p p p Ap p p 15002131,2201使得-⎛⎫⎪ ⎪⎪∴=-= ⎪ ⎪- ⎪⎝⎭P P AP J四、(20分)已知241111212,212211⎛⎫⎛⎫⎪ ⎪=-= ⎪ ⎪ ⎪ ⎪---⎝⎭⎝⎭A b ，（1）求A 的满秩分解；（2）求A +；（3）判断Ax b =是否有解，有解时求极小范数解，无解时求极小范数最小二乘解。

高等工程数学考试题及参考解答(仅供参考)

考试题及参考解答(参考)一、填空题（每小题3分，共15分） 1，设总体X 服从正态分布(0,4)N ，而1215(,,)X X X 是来自X 的样本，则221102211152()X X U X X ++=++服从的分布是_______ .解：(10,5)F ．2，ˆnθ是总体未知参数θ的相合估计量的一个充分条件是_______ . 解：ˆˆlim (), lim Var()0n nn n E θθθ→∞→∞==． 3，分布拟合检验方法有_______ 与____ ___. 解：2χ检验、柯尔莫哥洛夫检验． 4，方差分析的目的是_______ .解：推断各因素对试验结果影响是否显著．5，多元线性回归模型=+Y βX ε中，β的最小二乘估计ˆβ的协方差矩阵ˆβCov()=_______ . 解：1ˆσ-'2Cov(β)=()X X ．二、单项选择题（每小题3分，共15分）1，设总体~(1,9)X N ，129(,,,)X X X 是X 的样本，则___B___ .（A ）1~(0,1)3X N -；（B ）1~(0,1)1X N -；（C ）1~(0,1)9X N -；（D ~(0,1)N ． 2，若总体2(,)XN μσ，其中2σ已知，当样本容量n 保持不变时，如果置信度1α-减小，则μ的置信区间____B___ .（A ）长度变大；（B ）长度变小；（C ）长度不变；（D ）前述都有可能.3，在假设检验中，就检验结果而言，以下说法正确的是____B___ . （A ）拒绝和接受原假设的理由都是充分的；（B ）拒绝原假设的理由是充分的，接受原假设的理由是不充分的；（C ）拒绝原假设的理由是不充分的，接受原假设的理由是充分的；（D ）拒绝和接受原假设的理由都是不充分的.4，对于单因素试验方差分析的数学模型，设T S 为总离差平方和，e S 为误差平方和，A S 为效应平方和，则总有___A___ .（A ）T e A S S S =+；（B ）22(1)AS r χσ-；（C ）/(1)(1,)/()A e S r F r n r S n r ----；（D ）A S 与e S 相互独立.5，在多元线性回归分析中，设ˆβ是β的最小二乘估计，ˆˆ=-εY βX 是残差向量，则___B____ . （A ）ˆn E ()=0ε；（B ）1ˆ]σ-''-εX X 2n Cov()=[()I X X ；（C ）ˆˆ1n p '--εε是2σ的无偏估计；（D ）（A ）、（B ）、（C ）都对.三、（本题10分）设总体21(,)XN μσ、22(,)Y N μσ，112(,,,)n X X X 和212(,,,)n Y Y Y 分别是来自X 和Y 的样本，且两个样本相互独立，X Y 、和22X Y S S 、分别是它们的样本均值和样本方差，证明12(2)X Y t n n +-，其中2221212(1)(1)2X Yn S n S S n n ω-+-=+-.证明：易知221212(,)X YN n n σσμμ--+，(0,1)X Y U N =．由定理可知22112(1)(1)Xn S n χσ--，22222(1)(1)Yn S n χσ--．由独立性和2χ分布的可加性可得222121222(1)(1)(2)XYn S n S V n n χσσ--=++-．由U 与V 得独立性和t 分布的定义可得12(2)X Y t n n =+-．四、（本题10分）设总体X 的概率密度为1, 0,21(;), 1,2(1)0, x f x x θθθθθ⎧<<⎪⎪⎪=≤<⎨-⎪⎪⎪⎩其他，其中参数01)θθ<<（未知，12()n X X X ，，，是来自总体的一个样本，X 是样本均值，（1）求参数；的矩估计量θθˆ（2）证明24X 不是2θ的无偏估计量．解：（1）101()(,)22(1)42x x E X xf x dx dx dx θθθθθθ+∞-∞==+=+-⎰⎰⎰，令()X E X =，代入上式得到θ的矩估计量为1ˆ22X θ=-．（2）222211141 (4)44[()]4()424E X EX DX EX DX DX n nθθθ⎡⎤==+=++=+++⎢⎥⎣⎦，因为()00D X θ≥>，，所以22(4)E X θ>．故24X 不是2θ的无偏估计量．五、（本题10分）设总体X 服从[0,](0)θθ>上的均匀分布，12(,,)n X X X 是来自总体X 的一个样本，试求参数θ的极大似然估计．解：X 的密度函数为1,0;(,)0,x f x θθθ≤≤⎧=⎨⎩其他, 似然函数为1,0,1,2,,,()0,n i x i n L θθθ<<=⎧⎪=⎨⎪⎩其它显然0θ>时，()L θ是单调减函数，而{}12max ,,,n x x x θ≥，所以{}12ˆmax ,,,n X X X θ=是θ的极大似然估计．六、（本题10分）设总体X 服从(1,)B p 分布，12(,,)n X X X 为总体的样本，证明X 是参数p 的一个UMVUE ．证明：X 的分布律为1(;)(1),0,1x x f x p p p x -=-=．容易验证(;)f x p 满足正则条件，于是21()ln (;)(1)I p E f x p p p p ⎡⎤∂==⎢⎥∂-⎣⎦．另一方面1(1)1Var()Var()()p p X X n n nI p -===，即X 得方差达到C-R 下界的无偏估计量，故X 是p 的一个UMVUE ．七、（本题10分）某异常区的磁场强度服从正态分布20(,)N μσ，由以前的观测可知056μ=．现有一台新仪器, 用它对该区进行磁测, 抽测了16个点, 得261, 400x s ==, 问此仪器测出的结果与以往相比是否有明显的差异(α=0.05)．附表如下：t 分布表 χ2分布表解：设0H ：560==μμ．构造检验统计量)15(~0t ns X t μ-=，确定拒绝域的形式2t t α⎧⎫>⎨⎬⎩⎭．由05.0=α，定出临界值1315.2025.02/==t t α，从而求出拒绝域{}1315.2>t ．而60,16==x n ，从而 ||0.8 2.1315t ===<，接受假设0H ，即认为此仪器测出的结果与以往相比无明显的差异．八、（本题10分）已知两个总体X 与Y 独立，211~(,)X μσ，222~(,)Y μσ，221212, , , μμσσ未知，112(,,,)n X X X 和212(,,,)n Y Y Y 分别是来自X 和Y 的样本，求2122σσ的置信度为1α-的置信区间.解：设布定理知的样本方差，由抽样分，分别表示总体Y X S S 2221 , []/2121/212(1,1)(1,1)1P F n n F F n n ααα---<<--=-，则222221211221/2122/212//1(1,1)(1,1)S S S S P F n n F n n αασασ-⎛⎫<<=- ⎪----⎝⎭，所求2221σσ的置信度为α-1的置信区间为 222212121/212/212//, (1,1)(1,1)S S S S F n n F n n αα-⎛⎫ ⎪----⎝⎭．九、(本题10分)试简要论述线性回归分析包括哪些内容或步骤．答：建立模型、参数估计、回归方程检验、回归系数检验、变量剔除、预测。

《高等工程数学》试题+解答

《高等工程数学》试题解答（工程硕士及进修生用 2003.1）考生注意：1、可不抄题,答案必须写在统一配发的专用答题纸上； 2、本试题可能用到的常数：5752961 64199509750950 . ,. ,....===u u u . 一、填空题 (每空3分,共30分)。

(1) ⎥⎥⎦⎤⎢⎢⎣⎡=010100001H ；(2) 1)(Cond 2=U ；(3) 7 3 , ；(4) )1 1 (~)()(221221,F X X X X -+；(5) X 2ˆ=θ； (6) 664≥n ；(7) e A SS SS SS +=.二、(10分)[解] 记)(21A A diag A ,=，则21A A ,的特征多项式为2)1()()(21-==λλλA A f f ， ∵ O I A ≠21 -，O I A ≠22 -，∴ 2)1()()(21-==λλλA A m m ，取)( )(21λλA A m m ,的最小公倍式，得 2)1()(-=λλA m ，故A 的Jordan 标准形为⎭⎬⎫⎩⎨⎧⎥⎦⎤⎢⎣⎡⎥⎦⎤⎢⎣⎡ 111111 , diag . 三、(10分)[解一] 记⎥⎥⎦⎤⎢⎢⎣⎡--=πππ021 A ,其特征值为πλ-=1 (二重根),记则令 ⎩⎨⎧=-=⇒⎩⎨⎧=-=-⇒⎩⎨⎧'='=t t a t t t a t t a t a a g f g f 1 0 1 101 1 1 1 c o s sin cos cos sin )()()()(πππππππλλλλ ∴ . ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡==t t t t t A f g A g g A g At sin 00cos sin 000sin )()()()()(sin 2 11πππππππ[解二] ∵ J A 2 2001200022ππ∆⎥⎥⎦⎤⎢⎢⎣⎡--= ∴ . ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡---==t t t t t t t t t t J t At sin 00cos 2sin 000sin )2(2sin 00)2(2cos 2)2(2sin 00022sin )2sin(sin πππππππππππ 四、(10分)[解] 对A 进行行初等变换故⎥⎥⎦⎤⎢⎢⎣⎡-=⎥⎥⎦⎤⎢⎢⎣⎡-==--21 12 121 1 211 121)(11 R L L , 从而A 有Doolittle 分解：五、(10分)[证] 将ω扩充为nV 的一个标准正交基 B } {n ααω,,,2 =则∴ T B =-==} {} {n n T T T ααααωω,,,,,,22 B P 其中} 1 1 1 {,,, -=diagP 为对称和正交矩阵，故T 是对称变换和正交变换。

工科高等数学试题及答案

工科高等数学试题及答案一、单项选择题（每题2分，共10分）1. 函数f(x)=x^2+3x+2的导数为：A. 2x+3B. 2x+6C. x^2+3D. 3x+2答案：A2. 以下哪个选项是不定积分的结果？A. ∫x dx = x^2 + CB. ∫x^2 dx = x + CC. ∫e^x dx = e^x + CD. ∫ln(x) dx = x + C答案：C3. 函数y=sin(x)的二阶导数为：A. -cos(x)B. cos(x)C. -sin(x)D. sin(x)答案：A4. 极限lim(x→0) (1+x)^(1/x)的值为：A. 0B. 1C. eD. -e5. 以下哪个选项是二阶偏导数？A. ∂^2f/∂x∂yB. ∂^2f/∂x^2C. ∂^2f/∂y^2D. ∂^2f/∂x答案：A二、填空题（每题3分，共15分）6. 函数f(x)=x^3-3x+2的极值点为______。

答案：x=17. 函数y=x^2+3x+2的拐点为______。

答案：x=-3/28. 函数y=e^x的不定积分为______。

答案：e^x + C9. 函数y=ln(x)的不定积分为______。

答案：xln(x) - x + C10. 函数y=cos(x)的二阶导数为______。

答案：-cos(x)三、解答题（每题10分，共40分）11. 求函数f(x)=x^3-6x^2+11x-6在区间[1,3]上的最大值和最小值。

答案：最大值在x=2处取得，为f(2)=2；最小值在x=3处取得，为f(3)=0。

12. 求极限lim(x→∞) (x^2-3x+2)/(x^3+2x^2-5x)。

13. 求函数y=x^2-4x+3的极值点。

答案：极小值点为x=2，极大值点不存在。

14. 求函数y=x^3-3x+1的拐点。

答案：拐点为x=1。

四、证明题（每题15分，共30分）15. 证明函数f(x)=x^3在R上是增函数。

高等工程数学试题解答

《高等工程数学》科学出版社 吴孟达版习题答案(1-8章)

中南大学高等工程数学试题及参考答案.docx

大学工程数学考试题及答案

高等工程数学考试部分答案1

高等工程数学第二章习题及答案

高等工程数学试题及参考答案-工程硕士

最新最全中南大学高等工程数学附答案

《高等工程数学》习题三参考答案

《高等工程数学》习题一参考答案

《高等工程数学》科学出版社版习题答案(第三章)

高等工程数学练习题及答案

工程数学试题及答案高级专

大学工程数学试题及答案

工程数学试题及答案升专科

大学工程数学试题及答案

高等工程数学I 试题(A)与答案(2019.11.23)

高等工程数学考试题及参考解答(仅供参考)

《高等工程数学》试题+解答

工科高等数学试题及答案

《高等工程数学》科学出版社吴孟达版习题答案(1-8章)