《高等工程数学》吴孟达版习题答案(第二章)
高等数学第二章答案
高等数学第二章答案【篇一:高等数学第二章复习题及答案】>第二章一、填空题f(a?x)?f(a?x)?x?0xf(3?h)?f(3)?2、设f?(3)?2,则lim。
h?0______________2h1、设f(x)在x?a可导,则lim。
3、设f(x)?e,则limh?0?1xf(2?h)?f(2)?。
_____________hcosx?,f?(x0)?2,(0?x0?),则f(x0)?。
_______________________1?sinx2dy?5、已知x2y?y2x?2?0,则当经x=1、y=1时,。
dx_______________4、已知f(x)?6、f(x)?xex,则f???(ln2)?_______________。
__________7、如果y?ax(a?0)是y?x2?1的切线,则a?。
8、若f(x)为奇函数,f?(x0)?1且,则f?(?x0)?9、f(x)?x(x?1)(x?2)?(x?n),则f?(0)?10、y?ln(1?3?x),则y??11、设f?(x0)??1,则limx?0______________________________________________________。
x。
?___________f(x0?2x)?f(x0?x)_________________________12、设x?y?tany,则dy?。
13、设y?y???(0)?。
_______________14、设函数y?f(x)由方程xy?2lnx?y4所确定,则曲线y?f(x)在点(1,1)处的切线方程是______________________。
1???xcos15、f(x)??x??0_______________________x?0x?0。
,其导数在x?0处连续,则?的取值范围是16、知曲线y?x3?3a2x?b与x轴相切,则b2可以通过a表示为二、选择题。
高等数学课后习题答案第二章
=
1 4
1 tan
x 2
sec 2
x 2
5、设、 y =
1 2π D 1 2π D
e
−
( x−a)2 2D
,其中 a, D 是常数,求出使导数 y ′( x ) = 0 的 x 值
( x −a ) 2 2D
解: y ′ =
e
−
( x − a )2 2D
3、证明: (1) 、可导的偶(奇)函数的导数是奇函数(偶) (2) 、可导的周期函数的导数是具有相同周期的函数 证明:设 f ( x ) 是偶函数,且可导 则
f ( x) = f ( − x ) f (− x + ∆x ) − f (− x ) f ( x − ∆x ) − f ( x ) = lim = − f ′( x ) ∆x → 0 ∆x ∆x
[1 − ( x + ∆x ) 2 ] − (1 − x 2 ) − 2 x∆x − (∆x) 2 = lim = −2 x ∆x → 0 ∆x → 0 ∆x ∆x −b ) 2a
:
3、 设函数 f ( x) = ax 2 + bx + c , 其中 a, b, c 是常数, 求 f ′( x) , f ′(0) , f ′( −1) , f ′( 解
f ′(− x ) = lim
∆x →0
表明 f ′( x) 是奇函数。 设 f ( x) = f ( x + T )
f ′( x + T ) = lim
∆x →0
f ( x + T + ∆x ) − f ( x + T ) f ( x + ∆x ) − f ( x ) = lim = f ′( x) ∆ x → 0 ∆x ∆x
高等数学二教材答案
高等数学二教材答案
在高等数学二这门课程中,教材中的习题和答案对学生的学习非常
重要。
通过解答教材中的问题,学生可以巩固知识、提高解题能力,
并更好地理解课程的重点。
下面是我为您整理的高等数学二教材答案,希望对您的学习有所帮助。
第一章:多元函数微分学
1.1 函数的极限与连续
1.2 偏导数与全微分
1.3 多元复合函数的导数
1.4 隐函数及其导数
1.5 微分中值定理及其应用
第二章:重积分
2.1 二重积分的概念与性质
2.2 二重积分的计算方法
2.3 三重积分的概念与性质
2.4 三重积分的计算方法
2.5 曲线曲面积分的计算方法
第三章:曲线积分与曲面积分
3.1 曲线积分的概念与性质
3.2 第一类曲线积分和第二类曲线积分的关系3.3 曲线积分的计算方法
3.4 曲面积分的概念与性质
3.5 曲面积分的计算方法
第四章:向量场的微积分
4.1 向量场及其积分
4.2 格林公式
4.3 散度与旋度
4.4 散度定理和斯托克斯定理
4.5 保守场与调和场
第五章:常微分方程
5.1 常微分方程的基本概念
5.2 一阶常微分方程的解法
5.3 高阶常微分方程的解法
5.4 线性方程组的解法
5.5 常微分方程的应用
以上是高等数学二教材中各章节知识点的答案和解析。
通过仔细阅读教材和解答,希望您能更好地理解和掌握数学知识,并能够顺利完成习题。
在解题过程中,建议您思路清晰,注意计算细节,同时注重理论与实际应用的结合。
希望这些答案对您的学习有所帮助。
祝您学习进步!。
高等工程数学知到章节答案智慧树2023年南京理工大学
高等工程数学知到章节测试答案智慧树2023年最新南京理工大学第一章测试1.有限维线性空间上范数1,范数2之间的关系是参考答案:等价2.赋范线性空间成为Banach空间,需要范数足?参考答案:完备性3.标准正交系是一个完全正交系的充要条件是满足Parseval等式参考答案:对4.在内积空间中,可以从一组线性无关向量得到一列标准正交系参考答案:对5.矩阵的F范数不满足酉不变性参考答案:错6.与任何向量范数相容的矩阵范数是?参考答案:算子范数7.正规矩阵的谱半径与矩阵何种范数一致参考答案:矩阵2范数8.矩阵收敛,则该矩阵的谱半径参考答案:小于19.矩阵幂级数收敛,则该矩阵的谱半径参考答案:小于110.正规矩阵的条件数等于其最大特征值的模与最小特征值的模之商参考答案:对第二章测试1.l矩阵不变因子的个数等于( )参考答案:矩阵的秩2.Jordan标准形中Jordan块的个数等于( )参考答案:初等因子的个数3.Jordan块的对角元等于其( )参考答案:初等因子的零点4.n阶矩阵A的特征多项式等于( )参考答案:A的n个不变因子的乘积;A的n阶行列式因子5.下述条件中,幂迭代法能够成功处理的有( )参考答案:主特征值是实r重的;主特征值有两个,是一对共轭的复特征值;主特征值有两个,是一对相反的实数;主特征值只有一个6.n阶矩阵A的特征值在( )参考答案:A的n个行盖尔圆构成的并集与n个列盖尔圆构成的并集的交集中;A的n 个列盖尔圆构成的并集中;A的n个行盖尔圆构成的并集中7.不变因子是首项系数为1的多项式参考答案:对8.任意具有互异特征值的矩阵,其盖尔圆均能分隔开参考答案:错9.特征值在两个或两个以上的盖尔圆构成的连通部分中分布是平均的参考答案:错10.规范化幂迭代法中,向量序列uk不收敛参考答案:错第三章测试1.二阶方阵可作Doolittle分解参考答案:错2.若矩阵A可作满秩分解A=FG,则F的列数为A的()参考答案:秩3.矩阵的满秩分解不唯一.参考答案:对4.酉等价矩阵有相同的奇异值.参考答案:对5.求矩阵A的加号逆的方法有()参考答案:满秩分解;Greville递推法;奇异值分解;矩阵迭代法6.若A为可逆方阵,则参考答案:对7.用A的加号逆可以判断线性方程组Ax=b是否有解?参考答案:对8.A的加号逆的秩与A的秩相等参考答案:对9.若方阵A是Hermite正定矩阵,则A的Cholesky分解存在且唯一.参考答案:对10.是Hermite标准形.参考答案:错第四章测试1.()是利用Gauss消去法求解线性方程组的条件.参考答案:系数矩阵的顺序主子式均不为0;所有主元均不为02.关于求解线性方程组的迭代解法, 下面说法正确的是().参考答案:J法和GS法的敛散性无相关性;若系数矩阵A对称正定, 则GS迭代法收敛3.如果不考虑舍入误差, ()最多经n步可迭代得到线性方程组的解.参考答案:共轭梯度法4.关于共轭梯度法, 下面说法正确的是()参考答案:B和C都对5.下面哪些是求解线性方程组的迭代解法().参考答案:共轭梯度法;最速下降法6.若系数矩阵A对称正定, 则()参考答案:可用Cholesky法求解线性方程组7.任意线性方程组都可以通过三角分解法求解.参考答案:错8.最速下降法和共轭梯度法的区别在于选取的搜索方向不同.参考答案:对9.广义逆矩阵法可用于任意线性方程组的求解.参考答案:对10.Gauss消去法和列主元素法的数值稳定性相当.参考答案:错第五章测试1.对于凸规划,如果x为问题的KKT点,则其为原问题的全局极小点参考答案:对2.对于无约束规划问题,如果海塞阵非正定,我们可采用哪种改进牛顿法求解原问题?参考答案:构造一对称正定矩阵来取代当前海塞阵,并一该矩阵的逆乘以当前梯度的负值作为方向3.共轭梯度法中,为参考答案:FR公式4.内点罚函数法中常用的障碍函数有参考答案:倒数障碍函数;对数障碍函数5.广义乘子罚函数的优点是在罚因子适当大的情形下,通过修正拉格朗日乘子就可逐步逼近原问题的最优解?参考答案:对6.分子停留在最低能量状态的概率随温度降低趋于( ).参考答案:17.模拟退火算法内循环终止准则可采用的方法.参考答案:固定步数;由接受和拒绝的比率控制迭代步8.背包问题是组合优化问题吗?参考答案:对9.单纯形算法是求解线性规划问题的多项式时间算法.参考答案:错10.对于难以确定初始基本可行解的线性规划问题,我们引入人工变量后,可采用哪些方法求解原问题?参考答案:两阶段法;大M法第六章测试1.如果不限定插值多项式的次数,满足插值条件的插值多项式也是唯一的()参考答案:错2.改变节点的排列顺序,差商的值不变()参考答案:对3.Hermite插值只能用插值基函数的方法求解()参考答案:错4.在最小二乘问题中,权系数越大表明相应的数据越重要()参考答案:对5.加窗傅里叶变换时频窗的长宽比是信号自适应的()参考答案:错6.傅里叶变换域的点和时间域上的点是一一对应的()参考答案:错7.若f(t)的傅里叶变换为,则 f(2t)的傅里叶变换为 ( )参考答案:8.小波函数对应了()参考答案:高通滤波器第七章测试1.有界区域上的弦振动方程定解问题可以用傅里叶积分变换法求解。
高等工程数学第二章习题及答案
第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次.。
高等工程数学课后习题(吴孟达)(2021年整理)
高等工程数学课后习题(吴孟达)(word版可编辑修改)
编辑整理:
尊敬的读者朋友们:
这里是精品文档编辑中心,本文档内容是由我和我的同事精心编辑整理后发布的,发布之前我们对文中内容进行仔细校对,但是难免会有疏漏的地方,但是任然希望(高等工程数学课后习题(吴孟达)(word版可编辑修改))的内容能够给您的工作和学习带来便利。
同时也真诚的希望收到您的建议和反馈,这将是我们进步的源泉,前进的动力。
本文可编辑可修改,如果觉得对您有帮助请收藏以便随时查阅,最后祝您生活愉快业绩进步,以下为高等工程数学课后习题(吴孟达)(word版可编辑修改)的全部内容。
9、证明比较简单:按照F分布定义即可证明.
第七章
第八章
正交试验设计的目的是如何科学、合理地安排试验,使之能在很多的试验条件中选出代表性强的少数几个试验条件,并通过较少次数的试验就能取得最好的结果。
《高等工程数学》科学出版社 吴孟达版习题答案(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]T a 3 = [1 , -1, 1, -1]T a 4 = [1 , -1,-1, 1]T 下的坐标。
解:其坐标为:x =( 5/4, 1/4, -1/4,-1/4 )T 3.在R2×2中,求矩阵12A=03⎡⎤⎢⎥⎣⎦,在基 111B =11⎡⎤⎢⎥⎣⎦,211B =10⎡⎤⎢⎥⎣⎦,311B =00⎡⎤⎢⎥⎣⎦,410B =00⎡⎤⎢⎥⎣⎦下的坐标。
解:其坐标为:x =( 3, -3, 2,-1 )T 4.试证:在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=[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 =的坐标是:6.设R[x]n 是所有次数小于n 的实系数多项式组成的线性空间,求多项式p(x) = 1+ 2x n -1在基{1,(x -1),(x -1)2,(x -1)3,….,(x -1)n -1}的坐标。
高等数学 线性代数 习题答案第二章
第二章习题2-11. 证明:若lim n →∞x n =a ,则对任何自然数k ,有lim n →∞x n +k =a .证:由lim n n x a →∞=,知0ε∀>,1N ∃,当1n N >时,有n x a ε-<取1N N k =-,有0ε∀>,N ∃,设n N >时(此时1n k N +>)有n k x a ε+-<由数列极限的定义得 lim n k x x a +→∞=.2. 证明:若lim n →∞x n =a ,则lim n →∞∣x n ∣=|a|.考察数列x n =(-1)n ,说明上述结论反之不成立.证:lim 0,,.使当时,有n x n x aN n N x a εε→∞=∴∀>∃>-<而 n n x a x a -≤- 于是0ε∀>,,使当时,有N n N ∃>n n x a x a ε-≤-< 即 n x a ε-<由数列极限的定义得 lim n n x a →∞=考察数列 (1)nn x =-,知lim n n x →∞不存在,而1n x =,lim 1n n x →∞=,所以前面所证结论反之不成立。
3. 证明:lim n →∞x n =0的充要条件是lim n →∞∣x n ∣=0.证:必要性由2题已证,下面证明充分性。
即证若lim 0n n x →∞=,则lim 0n n x →∞=,由lim 0n n x →∞=知,0ε∀>,N ∃,设当n N >时,有0 0n n n x x x εεε-<<-<即即由数列极限的定义可得 lim 0n n x →∞=4. 利用夹逼定理证明:(1) lim n →∞222111(1)(2)n n n ⎛⎫+++ ⎪+⎝⎭ =0; (2) lim n →∞2!n =0. 证:(1)因为222222111112(1)(2)n n n n n n n n n n++≤+++≤≤=+而且 21lim0n n →∞=,2lim 0n n→∞=, 所以由夹逼定理,得222111lim 0(1)(2)n n n n →∞⎛⎫+++= ⎪+⎝⎭ . (2)因为22222240!1231n n n n n<=<- ,而且4lim 0n n →∞=, 所以,由夹逼定理得2lim 0!nn n →∞= 5. 利用单调有界数列收敛准则证明下列数列的极限存在. (1) x 1>0,x n +1=13()2n nx x +,n =1,2,…; (2) x 1x n +1,n =1,2,…;(3) 设x n 单调递增,y n 单调递减,且lim n →∞(x n -y n )=0,证明x n 和y n 的极限均存在.证:(1)由10x >及13()2n n nx x x =+知,有0n x >(1,2,n = )即数列{}n x 有下界。
《高等工程数学》科学出版社版习题答案(第三章)
《高等工程数学》――科学出版社版习题答案(第三章)(此习题答案仅供学员作业时参考。
因时间匆忙,有错之处敬请指正,谢谢!) (联系地址: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。
高等工程数学习题答案
高等工程数学习题答案【篇一:高等工程数学考试题及参考解答(仅供参考)】xt>一、填空题(每小题3分,共15分)2x12???x101,设总体x服从正态分布n(0,4),而(x1,x2?,x15)是来自x的样本,则u?222(x11???x15)服从的分布是_______ .解:f(10,5).?是总体未知参数?的相合估计量的一个充分条件是2,?n?)??, limvar(??)?0.解:lime(?nnn??n??3,分布拟合检验方法有_______ 与____ ___. 解:?检验、柯尔莫哥洛夫检验. 4,方差分析的目的是_______ .解:推断各因素对试验结果影响是否显著.22?1二、单项选择题(每小题3分,共15分)1,设总体x~n(1,9),(x1,x2,?,x9)是x的样本,则(a)x?1x?1~n(0,1);(b)~n(0,1); 31x?1~n(0,1). ~n(0,1);(d92(c)2,若总体x?n(?,?),其中?已知,当样本容量n保持不变时,如果置信度1??减小,则?的2置信区间____b___ .(a)长度变大;(b)长度变小;(c)长度不变;(d)前述都有可能.3,在假设检验中,就检验结果而言,以下说法正确的是____b___ . (a)拒绝和接受原假设的理由都是充分的;(b)拒绝原假设的理由是充分的,接受原假设的理由是不充分的;(c)拒绝原假设的理由是不充分的,接受原假设的理由是充分的;(d)拒绝和接受原假设的理由都是不充分的.4,对于单因素试验方差分析的数学模型,设st为总离差平方和,se为误差平方和,sa为效应平方和,则总有___a___ .(a)st?se?sa;(b)sa?2??2(r?1);(c)sa/(r?1)?f(r?1,n?r);(d)sa与se相互独立.se/(n?r)?)=?[in?x(xx???0n;(b)cov(??)x?];(a)???2?1(c)?????n?p?1是?2的无偏估计;(d)(a)、(b)、(c)都对.22三、(本题10分)设总体x?n(?1,?)、y?n(?2,?),(x1,x2,?,xn1)和(y1,y2,?,yn2)分别是来自x和y的样本,且两个样本相互独立,和sx、sy分别是它们的样本均值和样本方差,证明2222(n1?1)sx?(n2?1)sy其中s??.n1?n2?22?t(n1?n2?2),证明:易知??n(?1??2,?2n1??2n2),u??n(0,1).由定理可知2(n1?1)sx?2由独立性和?分布的可加性可得2??(n1?1),22(n2?1)sy?2??2(n2?1).v?2(n1?1)sx?2?2(n2?1)sy?2??2(n1?n2?2).由u与v得独立性和t分布的定义可得??t(n1?n2?2).?1?2?, 0?x??,??1,??x?1,其中参数?(0???1) 四、(本题10分)设总体x的概率密度为f(x;?)??2(1??)??0, 其他,???;?,xn)是来自总体的一个样本,是样本均值,未知,(x1,x2,(1)求参数?的矩估计量?(2)证明4不是2?2的无偏估计量.解:(1)e(x)??????xf(x,?)dx???01xx1?dx??dx??,?2(1??)2?42??2?令?e(x),代入上式得到?的矩估计量为?(2)1. 2111?1?4e(42)?42?4[?()2]?4?dx?(??)2??dx?????,424?n?n因为d(x)?0,??0,所以 e(4)??.故42不是?的无偏估计量.五、(本题10分)设总体x服从[0,?](??0)上的均匀分布,(x1,x2,?xn)是来自总体x的一个样本,试求参数?的极大似然估计.解:x的密度函数为,0?x??;??f(x,?)??0,其他,?222似然函数为???n,0?xi??,i?1,2,?,n,l(?)??其它??0,??max?x,x,?,x?是?的显然??0时,l(?)是单调减函数,而??max?x1,x2,?,xn?,所以?12n极大似然估计.六、(本题10分)设总体x服从b(1,p)分布,(x1,x2,?xn)为总体的样本,证明是参数p的一个umvue.证明:x的分布律为f(x;p)?px(1?p)1?x,x?0,1.容易验证f(x;p)满足正则条件,于是???1i(p)?e?lnf(x;p)??.?pp(1?p)??另一方面2var()?1p(1?p)1, var(x)??nnni(p)即得方差达到c-r下界的无偏估计量,故是p的一个umvue.七、(本题10分)某异常区的磁场强度服从正态分布n(?0,?),由以前的观测可知?0?56.现有一台新仪器, 用它对该区进行磁测, 抽测了16个点, 得?61, s?400, 问此仪器测出的结果与以往相2解:设h0:???0?56.构造检验统计量22t???0~t(15), n确定拒绝域的形式?t?t??.由??0.05,定出临界值t?/2?t0.025?2.1315,从而求出拒绝域t?2.1315.?????2而n?16,?60,从而 |t|???0.8?2.1315,接受假设h0,即认为此仪器测222出的结果与以往相比无明显的差异.2八、(本题10分)已知两个总体x与y独立,x~(?1,?1),y~(?2,?2),?1, ?2, ?1, ?2未知,?12(x1,x2,?,xn)和(y1,y2,?,yn)分别是来自x和y的样本,求2的置信度为1??的置信区间.?2122分别表示总体x,y的样本方差,由抽样分布定理知解:设s12,s2p?f?/2(n1?1,n2?1)?f?f1??/2(n1?1,n2?1)??1??,则22??s12/s2?12s12/s2p??2???1??, ?f1??/2(n1?1,n2?1)?2f?/2(n 1?1,n2?1)?22??s12/s2s12/s2?12,所求2的置信度为1??的置信区间为 ??.?2f(n?1,n?1)f(n?1,n?1)2?/212?1??/21?九、(本题10分)试简要论述线性回归分析包括哪些内容或步骤.答:建立模型、参数估计、回归方程检验、回归系数检验、变量剔除、预测【篇二:高等工程数学试题答案】>一、设总体x具有分布律其中?(0???1)为未知参数,已知取得了样本值x1?1,x2?2,x3?1,求?的矩估计和最大似然估计.解:(1)矩估计:ex??2?2?2?(1??)?3(1??)2??2??314?(1?2?1)?33??5. 令ex?,得?6(2)最大似然估计:l(?)?????2?(1??)?2??2?2256dln(?)?10?4?12?5?0 d???5得?6二、(本题14分)某工厂正常生产时,排出的污水中动植物油的浓度x~n(10,1),今阶段性抽取10个水样,测得平均浓度为10.8(mg/l),标准差为1.2(mg/l),问该工厂生产是22否正常?(??0.05,t0.025(9)?2.2622,?0.025(9)?19.023,?0.975(9)?2.700)解:(1)检验假设h0:?=1,h1:?≠1;取统计量:??222(n?1)s2?20;拒绝域为:?2≤?21?2222?=2.70或≥(n?1)??(9)?(n?1)???0.975?0.025=19.023, 22经计算:??2(n?1)s22?09?1.22??12.96,由于?2?12.96?(2.700,19.023)2,1故接受h0,即可以认为排出的污水中动植物油浓度的方差为?2=1。
高等数学第二章课后习题答案
⾼等数学第⼆章课后习题答案第⼆章导数与微分1. ()().1,102-'=f x x f 试按定义求设200200(1)(1)10(1)10'(1)lim lim1020lim lim (1020)20x x x x f x f x f x xx x x x→?→?→?→-+?--?---==-?==?-=-?2. 下列各题中均假定()0x f '存在,按导数定义观察下列极限,指出此极限表⽰什么, 并将答案填在括号内。
⑴ ()()=?-?-→?xx f x x f x 000lim(0'()f x -);⑵ ()=→?xx f x 0lim ('(0)f ),其中()()存在;且0,00f f '= ⑶ ()()=--+→hh x f h x f h 000lim(02'()f x ).3. 求下列函数的导数:⑴ ='=y x y ,4则34x ⑵ ='=y x y ,32则1323x -⑶ ='=y xy ,1则3212x -- ⑷ ='=y x x y ,53则115165x 4.求曲线. 21,3 cos 程处的切线⽅程和法线⽅上点??=πx y'sin ,'()3y x y π=-==-2(1)0y +-=法线⽅程为1)23y x π-=-化简得3)0x π+-+= 5. 讨论函数=≠=0001sin 2x x xx y 在0=x 处的连续性和可导性. 20(0)01lim sin 0(0)()x f x f x→===因为有界量乘以⽆穷⼩所以函数在0x =处连续因为 20001sin(0)(0)1lim limlim sin 0x x x x f x f x x xx x→?→?→?+?-==?=所以函数在0x =处可导.6. 已知()()()()是否存在?⼜及求 0 ,0 0 ,0 2f f f x x x x x f '''<-≥=-+ 2'00(0)(0)(0)lim lim 0h h f h f h f h h+→+→++-==='0lim 1h h f h f hf h h-→-→++--===- ''(0)(0)f f +-≠Q '(0)f ∴不存在7. ()(). , 0sin x f x x x x x f '??≥<=求已知当0x <时, '()(sin )'cos f x x x ==; 当0x >时, '()()'1f x x ==;当0x =时'00(0)(0)(0)limlim 1h h f h f hf hh +→→+-===++ '00(0)(0)sin (0)limlim 1h h f h f h f h h-→-→+-===- '(0)1f ∴=综上,cos ,0'()1,0x x f x x8. 求下列函数的导数:(1);54323-+-=x x x y (2);1227445+-+=x xx y 2222222232242222csc cot (1)2csc 2'(1)2(1)csc cot 4csc (1)23(3)(3ln )(2ln )(2)'(3ln )(94)ln 32(3ln )x x x x xy x x x x x x x x x x x x x xx y x x x x x x x x x x -+-=+-+-=+++-++=+-+-+=+g g 2'364652'20282y x x x ---=--+(3);3253xxe x y +-= (4);1sec tan 2-+=x x y2'152ln 23x x y x e =-+ 2'2sec sec tan y x x x =+(5);log 3lg 2ln 2x x x y +-= (6)()();7432x x y -+= 123'ln10ln 2y x x x =-+'422y x =--(7);ln x xy =(8);cos ln 2x x x y = 21ln 'x xx y x-= 221'2ln cos cos ln sin y x x x x x x x x x =+- 21ln x x-= 22ln cos cos ln sin x x x x x x x x =+- (9);1csc 22 xxy +=2222csc cot (1)2csc 2'(1)x x x x x y x -+-=+g g 2222(1)csc cot 4csc (1)x x x x x x -+-=+ (10).ln 3ln 223xx x x y ++= 2232223(3)(3ln )(2ln )(2)x x x x x x x x y x x ++-++=+ 4222(94)ln 32(3ln )x x x x x x x x -+-+=+9. 已知. ,cos 21sin 4πρρ=+=d d 求因为1sin cos sin 2d d ρ=+-所以412422284d d πρπ?==+-=+10. .1轴交点处的切线⽅程与写出曲线x xx y -= 令0y =,得11x x ==-或因为2'1y x -=+,所以 11'2,'2x x y y ==-==曲线在(1,0)处的切线⽅程为2(1)y x =-,即220x y --=;曲线在(1,0)-处的切线⽅程为2(1)y x =+,即220x y -+=。
高等数学上册习题答案吴赣昌人民大学出版社高数理工类
第4章不定积分习题4-11.求下列不定积分:知识点:直接积分法的练习——求不定积分的基本方法。
思路分析:利用不定积分的运算性质和基本积分公式,直接求出不定积分!★(1)⎰思路: 被积函数52x-=,由积分表中的公式(2)可解。
解:532223x dx x C --==-+⎰★(2)dx-⎰思路:根据不定积分的线性性质,将被积函数分为两项,分别积分。
解:1141113332223()24dx x x dx x dx x dx x x C --=-=-=-+⎰⎰⎰⎰★(3)22x x dx +⎰() 思路:根据不定积分的线性性质,将被积函数分为两项,分别积分。
解:2232122ln 23x xxx dx dx x dx x C +=+=++⎰⎰⎰()★(4)3)x dx -思路:根据不定积分的线性性质,将被积函数分为两项,分别积分。
解:3153222223)325x dx x dx x dx x x C -=-=-+⎰⎰★★(5)4223311x x dx x +++⎰思路:观察到422223311311x x x x x ++=+++后,根据不定积分的线性性质,将被积函数分项,分别积分。
解:42232233113arctan 11x x dx x dx dx x x C x x++=+=++++⎰⎰⎰ ★★(6)221x dx x +⎰思路:注意到222221111111x x x x x +-==-+++,根据不定积分的线性性质,将被积函数分项,分别积分。
解:2221arctan .11x dx dx dx x x C x x =-=-+++⎰⎰⎰ 注:容易看出(5)(6)两题的解题思路是一致的。
一般地,如果被积函数为一个有理的假分式,通常先将其分解为一个整式加上或减去一个真分式的形式,再分项积分。
★(7)x dx x x x ⎰34134(-+-)2 思路:分项积分。
解:3411342x dx xdx dx x dx x dx x x x x --=-+-⎰⎰⎰⎰⎰34134(-+-)2 223134ln ||.423x x x x C --=--++ ★(8)23(1dx x -+⎰思路:分项积分。
高等工程数学科学出版社吴孟达版习题答案18章
《高等工程数学》――科学出版社版习题答案: 第一章习题(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}的坐标。
高等工程数学课后答案
第六章7、设X 1,X 2,…X n 为总体X~N (μ,σ2)的样本,求E[21)(x x ni i-∑=],D[21)(∑=-ni ix x ]。
解:E[21)(x x ni i -∑=]=(n-1)E[11-n 21)(x x ni i-∑=]=(n-1)σ2因为)1(~)(2212--∑=n X x xni iσ所以 D[21)(∑=-ni ix x ]= ])([212σ∑=-ni ix xD =σ22(n-1)8、设X 1,X 2,…X 5为总体X~N (0,1)的样本,(1)试确定常数c 1、d 1,使得)(~)()(2254312211n x x x d x x c χ++++并求出n ;(2)试确定常数c 2、d 2,使得),(~)()(2543222212n m F x x x d x x c +++。
解:(1)212)(1x x n S n i i -=∑=且总体为X~N (0,1),所以c 1=21,d 1=31因为2χ分布具有可加性,即若X i ~2χ(i=1,……k ),且各样本相互独立,则)(~121∑∑==ki i ki in xχ,所以n=2。
(2)因为)2,0(~21N x x +,)3,0(~)(543N x x x ++,)1,0(~221N x x +, )1,0(~3543N x x x ++且相互独立, 所以221]2[xx ++2543]3[x x x ++)2(~2χ 因)2(~22221χx x +,)1(~3)(22543χx x x ++,所)1,2(~)(2)(325432221F x x x x x +++,所以)1,2(,2322F d c =10、设X 1,X 2,…X n ,X n+1为总体X~N (μ,σ2)的样本的容量为n+1的样本,)(11~,1221x x n s x n x i n i i --==∑=试证:(1))1(~~1ˆ1---=+n t sx x n n T n (2))1,0(~21σn n N x x n +-+(3))1,0(~21σn n N x x -- 证明:(1)因为),(~),1(~~)1(),,(~212222σμχσσμN x n s n n N x n +-- 所以)1,0(~1),1,0(~121N nn xx n n N x x n n +-+-++σσ 所以)1(~)1(~)1(1221---+-+n t n sn n n x x n σσ,即)1(~~1ˆ1---=+n t s x x n n T n (2)因为),(~),,(~212σμσμN x nN x n + 所以)1,0(~21σnn N x x n +-+ (3)因为∑∑==--=-=-ni i n i i x n x n n x n x x x 21111111,011)(1)(1)11(22121=--=--=--∑∑∑===ni n i i n i i n n n x E n x E n n x n x n n E μμ2222221121)1()11(σσσnn nn n x n x n n D ni n i i -=+-=--∑∑== 所以)1,0(~21σnn N x x --15、设X 1,X 2,…X n ,1为总体X 的样本,如果X 具有下列密度函数(其中参数均未知)试分别求这些参数的矩估计量与极大似然估计量。
高等工程数学 试题 答案
《高等工程数学》试题一、 设总体X 具有分布律其中(01)θθ<<为未知参数,已知取得了样本值1231,2,1x x x ===,求θ的矩估计和最大似然估计.解:(1)矩估计:2222(1)3(1)23EX θθθθθ=+⨯-+-=-+14(121)33X =++=令EX X =,得5ˆ6θ=. (2)最大似然估计:2256()2(1)22L θθθθθθθ=⋅⋅-=-45ln()10120d d θθθθ=-= 得5ˆ6θ= 二、(本题14分)某工厂正常生产时,排出的污水中动植物油的浓度)1,10(~N X ,今阶段性抽取10个水样,测得平均浓度为10.8(mg/L ),标准差为1.2(mg/L ),问该工厂生产是否正常?(220.0250.0250.9750.05,(9) 2.2622,(9)19.023,(9) 2.700t αχχ====)解:(1)检验假设H 0:σ2=1,H 1:σ2≠1; 取统计量:2022)1(σχs n -=;拒绝域为:χ2≤)9()1(2975.0221χχα=--n =2.70或χ2≥2025.022)1(χχα=-n =19.023, 经计算:96.1212.19)1(2222=⨯=-=σχs n ,由于)023.19,700.2(96.122∈=χ2,故接受H 0,即可以认为排出的污水中动植物油浓度的方差为σ2=1。
(2)检验假设101010≠'='μμ:,:H H ; 取统计量:10/10S X t -=~ )9(2αt ;拒绝域为2622.2)9(025.0=≥t t ;1028.210/2.1108.10=-=t <2.2622 ,所以接受0H ', 即可以认为排出的污水中动植物油的平均浓度是10(mg/L )。
综上,认为工厂生产正常。
三、 在单因素方差分析中,因素A 有3个水平,每个水平各做4次重复实验,完成下列方差分析表,在显著水平0.05α=下对因素A 是否显著做检验。
高等数学专科教材课后答案
高等数学专科教材课后答案第一章:函数与极限1. a) 当x趋于正无穷时,f(x)趋近于0。
b) 当x趋于负无穷时,f(x)趋近于无穷大。
2. a) 函数f(x)在x = 2处有一个间断点。
b) 函数f(x)在x = 2处无定义。
3. a) 函数f(x)在x = 2处无极限。
b) 函数f(x)在x = 2处有一个有限极限。
4. a) 函数f(x)在x = 3处无极限。
b) 函数f(x)在x = 5处有一个有限极限。
5. a) 函数f(x)在x = 2处无极限。
b) 函数f(x)在x = 2处有一个无穷极限。
6. a) 函数f(x)在x = 0处无极限。
b) 函数f(x)在x = 0处有一个无穷极限。
7. a) 函数f(x)在x = 1处无极限。
b) 函数f(x)在x = 1处有一个有限极限。
第二章:导数与微分8. a) 函数f(x)在x = 2处不可导。
b) 函数f(x)在x = 2处可导。
9. a) 函数f(x)在x = 3处可导。
b) 函数f(x)在x = 3处不可导。
10. a) 函数f(x)在x = 0处可导。
b) 函数f(x)在x = 0处不可导。
11. a) 函数f(x)在x = 1处不可导。
b) 函数f(x)在x = 1处可导。
12. a) 函数f(x)在x = 2处可导。
b) 函数f(x)在x = 2处不可导。
13. a) 函数f(x)在x = 0处不可导。
b) 函数f(x)在x = 0处可导。
14. a) 函数f(x)在x = 1处可导。
b) 函数f(x)在x = 1处不可导。
第三章:微分中值定理与泰勒公式15. a) 函数f(x)在区间[a, b]上满足Rolle定理条件。
b) 函数f(x)在区间[a, b]上不满足Rolle定理条件。
16. a) 函数f(x)在区间[a, b]上满足Lagrange中值定理条件。
b) 函数f(x)在区间[a, b]上不满足Lagrange中值定理条件。
最新《高等工程数学》吴孟达版习题答案(第二章)
最新《高等工程数学》吴孟达版习题答案(第二章)《高等工程数学》――科学出版社版习题答案(第二章)(此习题答案仅供学员作业时参考。
因时间匆忙,有错之处敬请指正,谢谢!)(联系地址:yangwq@/doc/3318299940.html, ) P50 1.求下列矩阵的特征值、代数重数核几何重数,并判断矩阵是否可对角化(1)110020112-(2)011121213--(3)411030102-解:(1)特征值:1231(1)()λλλ=代数重数和几何重数均为,==2代数重数和几何重数均为2可对角化。
(2)特征值:1231(1)()λλλ=代数重数和几何重数均为,==2代数重数为2和几何重数为1不可对角化。
(3)特征值:123(1)λλλ===3代数重数为3、几何重数均为不可对角化。
2.求下列矩阵的不变因子、初等因子和Jordan 标准形(1)3732524103-----(2)413002 10-1 (3)1234012300120001(4)3000013000001100002000112-解:(1)不变因子是:123d d d i λλλ+=1,=1,=(-1)(-i)()初等因子是:i λλλ+(-1),(-i),()Jordan 标准形是:1000000i i ??-??(2)不变因子是:123d d d λ3=1,=1,=(-3)初等因子是:λ3(-3)Jordan 标准形是:310031003(3)不变因子是:1234d d d d λ4=1,=1,=1,=(-1)初等因子是:λ4(-1)Jordan 标准形是:1100011000110001(4)不变因子是:12345d d d d d λλλλλ=1,=1,=1,=(-2)(-3),=(-1)(-2)(-3)初等因子是:λλλλλ(-2),(-3),(-1),(-2),(-3)Jordan 标准形是:10000020000020*******0003??3.设(1)110A 0012-=22(2)33A 613--1=-7-11-(3)010A 111011=--求可逆矩阵P ,使得P -1AP 是Jordan 标准形解:(1)A 的特征值为1231λλλ=,==2 对应的特征向量是:121,ααTT=(,0,-1)=(0,0,1)二级根向量是:(2)2αT=(-1,1,0)(2)122101(,,0110002102P P AP ααα--??=??=??1)=0-1100(2)A 的特征值为123λλλ===2 对应的特征向量是:11αT=(,2,1)二级根向量和三级根向量是:(2)(3)11,ααT T=(1,3,3)=(0,2,2)(2)(3)111110(,,3232102102P P AP ααα-??=??=??1)=21200(3)此题数据不便于求解特征值,A 的特征多项式是:3210()|A|11121011f I λλλλλλλλ-=---=-??-??=-+4.试求第2题最小多项式。
高等数字习题答案
高等数字习题答案在高等数学的习题中,通常会涉及到微积分、线性代数、常微分方程、概率论与数理统计等主题。
由于题目的多样性和复杂性,这里我将提供一些典型的习题类型和解题方法,以及一些可能的解题步骤。
# 微积分习题类型:求导、求积分、微分方程求解、极值问题、泰勒级数展开等。
解题方法:- 求导:使用基本导数公式,如 \(\frac{d}{dx}x^n = nx^{n-1}\),或链式法则、乘积法则、商法则等。
- 求积分:利用基本积分公式,如 \(\int x^n dx =\frac{1}{n+1}x^{n+1} + C\),或者使用换元积分法、分部积分法等。
- 微分方程:根据方程类型(如一阶线性微分方程、二阶常系数微分方程)选择合适的解法。
- 极值问题:通过求导数的零点和二阶导数的符号来确定极值点。
示例:- 求函数 \(f(x) = 3x^2 + 2x - 5\) 的导数。
解:\(f'(x) = 6x + 2\)。
# 线性代数习题类型:矩阵运算、向量空间、线性变换、特征值与特征向量、二次型等。
- 矩阵运算:加减法、乘法、求逆矩阵、行列式计算等。
- 向量空间:基、维度、线性组合、线性相关性等。
- 特征值与特征向量:通过解特征方程 \(Av = \lambda v\) 来找到特征值和特征向量。
示例:- 求矩阵 \(A = \begin{bmatrix} 1 & 2 \\ 3 & 4 \end{bmatrix}\) 的逆矩阵。
解:\(A^{-1} = \frac{1}{-2} \begin{bmatrix} 4 & -2 \\ -3 & 1 \end{bmatrix}\)。
# 常微分方程习题类型:一阶微分方程、二阶微分方程、高阶微分方程、偏微分方程等。
解题方法:- 一阶微分方程:分离变量法、积分因子法、常数变易法等。
- 二阶微分方程:特征方程法、降阶法、幂级数法等。
示例:- 求解微分方程 \(y'' - 3y' + 2y = 0\)。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
《高等工程数学》――科学出版社版习题答案(第二章)
(此习题答案仅供学员作业时参考。
因时间匆忙,有错之处敬请指正,谢谢!) (联系地址:***************.cn ) P50
1. 求下列矩阵的特征值、代数重数核几何重数,并判断矩阵是否可对角化
(1)110020112⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦- (2)011121213⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦-- (3)411030102⎡⎤⎢⎥⎢⎥
⎢⎥⎣⎦
-
解:(1)特征值:
1231(1)()λλλ=代数重数和几何重数均为,==2代数重数和几何重数均为2
可对角化。
(2)特征值:
1231(1)()λλλ=代数重数和几何重数均为,==2代数重数为2和几何重数为1
不可对角化。
(3)特征值:
123(1)λλλ===3代数重数为3、几何重数均为
不可对角化。
2. 求下列矩阵的不变因子、初等因子和Jordan 标准形
(1)3732524103⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦
-----(2)413002⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦10-1 (3)1
2340
1230
0120
00
1⎡⎤⎢⎥⎢⎥
⎢⎥⎢⎥⎣⎦
(4)3
000
0130
0000110000200
01
12⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦
-
解:(1)不变因子是:123d d d i λλλ+=1,=1,=(-1)(-i)()
初等因子是:i λλλ+(-1),(-i),()
Jordan 标准形是:1000000i i ⎡⎤
⎢⎥⎢⎥⎢⎥-⎣⎦
(2)不变因子是:123d d d λ3
=1,=1,=(-3)
初等因子是:λ3
(-3)
Jordan 标准形是:310031003⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦
(3)不变因子是:1234d d d d λ4
=1,=1,=1,=(-1)
初等因子是:λ4
(-1)
Jordan 标准形是:11000
11000110001⎡⎤⎢⎥⎢
⎥⎢⎥
⎢⎥
⎣⎦
(4)不变因子是:12345d d d d d λλλλλ=1,=1,=1,=(-2)(-3),=(-1)(-2)(-3)
初等因子是:λλλλλ(-2),(-3),(-1),(-2),(-3)
Jordan 标准形是:1
0000020000
020*******
0003⎡⎤
⎢⎥⎢
⎥⎢⎥⎢⎥
⎢⎥⎢⎥⎣⎦
3. 设(1)110A 0012⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦-=22(2)33A 613⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦
--1=-7-11-(3)01
0A 111011⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦=-- 求可逆矩阵P ,使得P -
1AP 是Jordan 标准形
解:(1)A 的特征值为1231λλλ=
,==2 对应的特征向量是:121,ααT
T
=(,0,-1)=(0,0,1)
二级根向量是:(2)2αT
=(-1,1,0)
(2)
122101(,,0110002102P P AP ααα--⎡⎤
⎢⎥=⎢⎥
⎢⎥⎣⎦
⎡⎤
⎢⎥=⎢⎥
⎢⎥⎣⎦
1)=0-1100
(2)A 的特征值为123λλλ===2 对应的特征向量是:11αT
=(,2,1)
二级根向量和三级根向量是:(2)(3)11,ααT T
=(1,3,3)=(0,2,2)
(2)(3)111110(,,3232102102P P AP ααα-⎡⎤
⎢⎥=⎢⎥
⎢⎥⎣⎦
⎡⎤
⎢⎥=⎢⎥
⎢⎥⎣⎦
1)=21200
(3)此题数据不便于求解特征值,A 的特征多项式是:
3210()|A|11121011f I λλλλλλλλ-⎡⎤⎢⎥=---=-⎢⎥
⎢⎥-⎣⎦
=-+
4. 试求第2题 最小多项式。
解:(1)最小多项式是:A m ()i λλλλ=+(-1)(-i)() (2)最小多项式是:A m ()λλ=3
(-3) (3)最小多项式是:A m ()λλ=4(-1)
(4)最小多项式是:A m ()λλλλ=(-1)(-2)(-3)
5. 设10A 10⎡⎤
⎢⎥-⎢⎥⎢⎥⎣⎦
2=0101,计算方阵多项式8542
()34g A A A A A I -++-=2
解:因为:
854253232
()34
(245914)(21)(243710)
g λλλλλλλλλλλλλ-++-=+-+--++-+=2
而3
()(21)f λλλ=-+是A 的特征多项式 ,所以f (A )=0
故有2
34826()437100
956106134g A A A I --⎛⎫
⎪-+=- ⎪ ⎪-⎝⎭
=2
6. 设A 是可逆方阵,证明A -
1可表示为A 的方阵多项式。
证明:设A 是n 阶方阵,其特征多项式是:
1011()...n n n n f a a a a λλλλ--=++++
因A 可逆,所以0n a ≠(为什么?自己证明)
由1
011()...0n n n n f A a A a A a A a I --=++++= 得
112011(...)/n n n n A a A a A a I a ----=+++
所以A -1
可表示为A 的多项式。
7. 设0A ≠,0(2)k
A k =≥,证明A 不能与对角矩阵相似。
证明:由题设知,A 的最小多项式是:2
()A m λλ=,有重根,所以不能相似对角化。
8. 已知()p
A I p =为正整数,证明A 与对角矩阵相似。
证明:由题设知, ()1p
g λλ=-是A 的零化多项式,而多项式()1p
g λλ=-没有重根(为什么?自己证!!),所以A 的最小多项式没有重根,故与对角矩阵相似
9. 设2
A A =,试证A 的Jordan 标准形是diag{1,1,...,1,0, 0
证明:因为2
()g λλλ=-是A 的零化多项式,且是最小多项式,所以A 的特征值只能是0和1,且可对角化,所以A 的Jordan 标准形是diag{1,1,…,1,0,…,0} 10.
设方阵A 的特征多项式()f λ和最小多项式()m λ分别为:
(1)4
2
2
2
()(2)(3),()(2)(3)f m λλλλλλ=--=-- (2)3
3
2
()(3)(5),()(3)(5)f m λλλλλλ=--=-- 试确定A 的所有可能的Jordan 标准形 解:(1)A 的可能Jordan 标准形为
22212313⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦ 或21221
2313⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥
⎢⎥⎢⎥
⎣⎦
(2)A 的可能Jordan 标准形为
2212555⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥
⎢⎥⎢⎥
⎣⎦。