第1章2范数

合集下载

第一章 误差与范数

第一章 误差与范数

1.1 误差的来源例1.1.1 用差商ha f h a f a f )()()(-+≈'求x x f ln )(=在3=x 处导数的近似值.取1.0=h ,1000.0=h ,h =0.000 000 000 000 001和h =0.000 000 000 000 000 1分别用MATLAB 软件计算,取十五位数字计算.解 在MATLAB 工作窗口输入下面程序>>a=3;h=0.1;y=log(a+h)-log(a);yx=y/h运行后得yx = 0.32789822822991 将此程序中h 改为0.000 1,运行后得yx = 0.33332777790385后者比前者好.再取h = 0.000 000 000 000 001,运行后得yx = 0.44408920985006不如前者好.取h = 0.000 000 000 000 000 1,运行后得yx = 0算出的结果反而毫无价值.例1.1.2 分别求方程组b AX =在下列情况时的解,其中A ⎪⎪⎭⎫⎝⎛=011111.. (1)⎪⎪⎭⎫ ⎝⎛=22b ;(2)⎪⎪⎭⎫⎝⎛=0122.b . 解 (1) 首先将方程组b AX =化为同解方程b A X 1-=,然后在MATLAB 工作窗口输入程序>> b=[2,2]';A=[1,1;1,1.01]; X=A\b运行后输出当⎪⎪⎭⎫ ⎝⎛=22b 时,b AX =的解为⎪⎪⎭⎫⎝⎛=02X ;(2)同理可得,当⎪⎪⎭⎫ ⎝⎛=0122.b 时,b AX =的解为⎪⎪⎭⎫⎝⎛=11X .例1.1.3 计算e 的近似值. 解 泰勒级数e +++++++=!!4!3!21432 n x x x x x n x)(∞<<-∞x , 取1=x ,得e +++++++=!1!41!31!2111  n . (1.2)这是一个无限过程,计算机无法求到精确值.只能在(1.2)取有限项时计算,再估计误差.如果取有限项!!!!)(  n s n 1413121111++++++=作为e 的值必然会有误差,根据泰勒余项定理可知其截断误差为e !)1()1( +=-n e s n θ)10(<<θ.如果取(1.2)的前九项,输入程序>> n =8;s=1;S =1; for k=1:ns=s*k;S=S+1/s, ends,S,R=3/(s*(n+1)) 或>>S1=1+1+1/2+1/(1*2*3)+1/(1*2*3*4)+1/(1*2*3*4*5)+1/(1*2*3*4*5*6)+1/(1*2*3*4*5*6*7)+1/(1*2*3*4*5*6*7*8),R1=3/(1*2*3*4*5*6*7*8*9)运行后结果S = R =2.71827876984127 8.267195767195768e-006 因为截断误差为e ),10(101968.267!93!)18()1(6-8<<⨯≈<+=-θθ  e s 所以e 的近似值e ≈≈++++++++=!81!71!61!51!41!31!2111)1(8 s 2.718 28.1.2 误差和有效数字例1.2.1 取282.718作为e 的四舍五入近似值时,求其绝对误差和相对误差. 解 在MATLAB 工作窗口输入程序>>juewu=exp(1)-2.71828运行后输出结果为juewu = 1.828 459 045 505 326e-006例1.2.2 计算⎰π20sin x xd x 的近似值,并确定其绝对误差和相对误差.解 因为被积函数xxsin 的原函数不是初等函数,故用泰勒级数求之.++-+-=!!!!sin 9 75 386x x x x x x 421 )(∞<<-∞x , (1.5) 这是一个无限过程,计算机无法求到精确值.可用(1.5)的前四项!!!75 36x x x -+-421代替被积函数xxsin ,得 ⎰π=20sin x x y d ⎰π≈20(x !!!14275 36x x x -+-)d x =!7)2(!5)2(!3)2(275375 3⋅π-⋅π+⋅π-π=y ˆ. 根据泰勒余项定理和交错级数收敛性的判别定理,得到绝对误差!99)2(ˆ9⋅<-=πyy R = WU , 在MA TLAB 命令窗口输入计算程序如下:>>syms xf=1-x^2/(1*2*3)+x^4/(1*2*3*4*5)-x^6/(1*2*3*4*5*6*7)y=int(f,x,0,pi/2),y1=double(y)y11=pi/2-(pi/2)^3/(3*3*2)+(pi/2)^5/(5*5*4*3*2)-(pi/2)^7/(7*7*6*5*4*3*2)inf=int(sin(x)/x,x,0,pi/2) ,infd=double(inf) WU =(pi/2)^9/(9*9*8*7*6*5*4*3*2), R =infd-y11因为运行后输出结果为: =y 1.370 762 168 154 49,yˆ=1.370 744 664 189 38,=R 1.750 396 510 491 47e-005, WU = 1.782 679 830 970 664e-005410-<.所以,yˆ的绝对误差为=ε410-,故⎰π=20sin x xy d 7 1.370≈x .yˆ的相对误差为 =r ε71.37010ˆ4-=y ε<0.007 3%.1.3 误差估计的基本方法例1.3.4 设计三种算法求方程01522=-+x x 在)3,2(的一个正根*x 的近似值,并研究每种算法的误差传播情况.解 为解已知方程,我们可以设计如下三种算法,然后将计算结果与此方程的精确解5.2*=x 比较,观察误差的传播.算法1 将已知方程化为同解方程=x 2215x -.取初值20=x ,按迭代公式21215k k x x -=+依次计算 ,,,,21n x x x ,结果列入表1–3中.算法2 将已知方程化为同解方程1215+=x x .取初值20=x ,按迭代公式 12151+=+k k x x依次计算 ,,,,21n x x x ,结果列入表1–3中.算法3 将已知方程化为同解方程141522+-+-=x x x x x .取初值20=x ,按迭代公式为1415221+-+-=+k k kk k x x x x x 依次计算 ,,,,21n x x x ,结果列入表1–3中.我们为这三种算法的计算编写两套MATLAB 程序如下: (1)MATLAB 主程序function [k,juecha,xiangcha,xk]= liti112(x0,x1,limax) % 输入的量--x0是初值, limax 是迭代次数和精确值x; % 输出的量--每次迭代次数k 和迭代值xk,% --每次迭代的绝对误差juecha 和相对误差xiangcha , x(1)=x0;for i=1:limaxx(i+1)=fl(x(i));%程序中调用的fl.m juecha = abs(x(i)-x1);xiangcha = juecha /( abs(x(i))+eps);xk=x(i);k=i-1;[(i-1),juecha,xiangcha,xk] end(2)MATLAB 调用函数程序及其计算结果①算法2的MATLAB 调用函数程序function y1=fl(x)y1=15/(2*x+1);② 将MATLAB 主程序和调用函数程序分别命名liti112.m 和fl.m ,分别保存为M 文件,然后在MATLAB 工作窗口输入命令>> [k,juecha,xiangcha,xk]= liti112(2,2.5,100) ③运行后输出计算结果列入表1–3和表 1-4中.④将算法2的MATLAB 调用函数程序的函数分别用y1=15-2*x^2和y1=x-(2*x^2+x-15)/(4*x+1)代替,得到算法1和算法3的调用函数程序,将其保存,运行后将三种算法的前8个迭代值821,,,x x x 列在一起(见表 1-3),进行比较.将三种算法的821,,,x x x 对应的绝对误差和相对误差的值列在一起(见表 1-4),进行比较.1.4 数值计算中应注意的问题例1.4.1 求数)181(71915-+⨯=-x 的近似值. 解 (1)直接用MATLAB 命令>> x=(7^15)*(sqrt(1+8^(-19))-1)运行后输出结果x = 0问题出现在两个相近的数1981-+与1相减时,计算机运行程序>>sqrt(1+8^(-19))-1运行后输出结果ans = 0由于计算机硬件只支持有限位机器数的运算,因此在计算中可能引入和传播舍入误差.因为有效数字的严重损失,导致输出18119-+-的结果为0,计算机不能再与数157继续进行真实的计算,所以,最后输出的结果与x 的精确值不符.(2)如果化为18187)181(71919151915++⨯=-+⨯=---x ,再用MATLAB 命令>> x=(7^15)*( (8^(-19))/(sqrt(1+8^(-19))+1))运行后输出结果x = 1.6471e-005 这是因为18119-+-化为18181919++--后,计算机运行程序>> x= (8^(-19))/(sqrt(1+8^(-19))+1)运行后的结果为x =3.4694e-018 由于有效数字的损失甚少,所以运算的结果-18103.4694⨯再与157继续计算,最后输出的结果与x 的精确值相差无几.例1.4.2 求数)13030ln(2--=y 的近似值. 解 (1)直接用MATLAB 程序>> x=30;x1= sqrt(x^2-1)运行后输出结果x1 = 29.9833 输入MATLAB 程序>> x=30; x1=29.9833;y=log(x-x1)运行后输出结果y = -4.0923(2)因为)13030ln(2--中的30=x 很大,如果采用倒数变换法111221-+=--=x x x x z ,即130301ln)13030ln(22-+=--)190030ln(-+-=.输入MATLAB 程序>> x=30;y=-log(x+sqrt(x^2-1))运行后输出结果y = -4.0941(3)输入MA TLAB 程序>> x=30; y=log(x-sqrt(x^2-1))运行后输出结果y = -4.0941 可见,(2)计算的近似值比(1)的误差小.参加计算的数,有时数量级相差很大.如果不注意采取相应的措施,在它们的加减法运算中,绝对值很小的那个数经常会被绝对值较大的那个数“吃掉”,不能发挥其作用,造成计算结果失真.例1.4.4 请在16位十进制数值精度计算机上利用软件MATLAB 计算下面的两个数0.30.1111111111111111*++=x 和0.30.11111111111111111*++=y将计算结果与准确值比较,解释计算结果.解 在MATLAB 工作窗口输入下面程序>> x=111111*********+0.1+0.3, y=1111111111111111+0.1+0.3运行后输出结果x = 1.111111*********e+014,y =1.111111*********e+015 从输出的结果可以看出,x *x =,而y *y ≠.为什么*y 仅仅比*x 多一位1,而y *y ≠呢?这是因为计算机进行运算时,首先要把参加运算的数写成绝对值小于1而“阶码”相同的数,这一过程称为数的“对阶”.在16位十进制数值精度计算机上利用软件MATLAB 计算这两个数,把运算的数*x 写成浮点规格化形式为,151515*103000**********.0001010000000000000.000100111111111111111.0⨯+⨯+⨯=x在16位十进制数值精度计算机上,三项的数都表示为小数点后面16位数字的数与1510之积,所以,计算机没有对数进行截断,而是按原来的三个数进行计算.因此,计算的结果x *x =.而161616*10030000000000000.00010010000000000000.000101111111111111111.0⨯+⨯+⨯=y三项的数都表示写成绝对值小于1而“阶码”都为1610的数以后,第一项的纯小数的小数点后面有16位数字.但是,后两项数的纯小数的小数点后面有17位数字,超过了16位十进制数值精度计算机的存储量,计算机对后两项的数都进行截断最后一位,即后两项的数都是16位机上的零,再进行计算,所以计算结果与实际不符.五、向量和矩阵的范数例1.5.1: 用matlab 求下列向量的2,1,5,,-∞∞范数。

向量范数3-1,3-2,3-3

向量范数3-1,3-2,3-3
X
A
X AX


X x1 , x2 , , xn R n
T
试证上述函数是向量范数,称为向量的加权范数或椭圆范数。 证明 因为A是正定对称矩阵,故存在可逆矩阵P,使得
P T AP I
从而
A P
X
A

1 2 A
T 1
P P
T T 1 2
1

1 T
1 2
P 1 B T B
证明 易验证条件(i)和(ii)成立,现验证条件(iii)也 成立。 下面用到了Chauchy-Schwarz不等式。
x y
2 2
x y , x y ( x, x ) ( x, y ) ( y , x ) ( y , y )
x
2 2
2 x
2y2源自 y2 2定理对 x ( x , x ,, x )T C n C n R 分别定义三个函数 1 2 n
x
x
1

x
i 1
n i 1
n
i
1 2
1-范数,
)
2
( xi
2
2-范数(或Euclid范数)
x

max xi
1 i n
∞-范数(或最大值范数)。
它们均构成范数。 说明:在同一个向量空间,可以定义多种向量范数,而对 于同一个向量,不同定义的范数,其大小可能不同。

AX

AX H A X H A X
即矩阵范数与向量范数相容
算子范数
定义 设

即由向量范数构造矩阵范数



分别是 C m 和 C n

第一章 误差与范数

第一章 误差与范数

习 题1.11. 分别求方程组b AX=在下列情况时的解,其中⎪⎪⎭⎫⎝⎛=01.2222A .(1)⎪⎪⎭⎫ ⎝⎛=22b ;(2)⎪⎪⎭⎫ ⎝⎛=01.22b . 2. 用差商h a f h a f a f )()()('-+≈,求x x f ln )(=在5=x 处导数)5('f 的近似值,取十五位数字计算,h 分别取0.000 1和0.000 000 01.习 题1.21. 为了使 4.420=的近似值的相对误差小于%1.0,试问取几位有效数字?2. 写出光速10*10)001000.0925997.2(⨯±=c s /cm 的绝对误差和相对误差. 3. 取655923.141作为π的四舍五入近似值时,求其有效数字、绝对误差和相对误差.4. 问722,141.3,142.3分别作为π的近似值各具有几位有效数字?习 题1.31.测得电压)2110(±=V V ,电流)5.020(±=I A ,则由欧姆定律得Ω==5.5I VR ,求R 的绝对误差和相对误差.2. 设nx y =,求y 的相对误差与x 的相对误差的关系式.3. 设计三种算法求方程0322=-+x x 在)2,1(的一个正根*x 的近似值,并研究每种算法的误差传播情况.4. 已测得某场地长*x 的值110=x m ,宽*y 的值80=y m ,已知x 的绝对误差是2.0m ,y 的绝对误差是1.0m ,求场地面积xy S =的绝对误差和相对误差.5. 正方形的边长约为200cm ,怎样测量才能使其面积误差不超过2cm 1?1. 求数)274(41915-+⨯=-x 的近似值.2. 利用四位数学用表求2cos 1-的近似值,要求至少有2位有效数字.3. 请在16位十进制数值精度计算机上利用软件MATLAB 计算下面的两个数0.30.1333333333333333*++=x 和0.30.133333333333333333*++=y将计算结果与准确值比较,解释计算结果.4. 利用恒等变换,使下列表达式的计算结果比较精确.(1)1,11>>--+x x x x x ; (2)e12-x ; (3)x xsin cos 1-. 5. 设x y ln =.当)0(>≈a a x 时,如果已知对数a ln 的绝对误差限为n-⨯1021,试估计真数a 的相对误差限及有效数字位数.6. 指出下列各题的合理计算途径(对给出具体数据的,请算出结果):(1)o1cos 1-(三角函数值取四位有效数字);(2))13030ln(2--;(3)x xsin cos 1- (其中x充分小); (4)127x .习 题 1.51.设⎪⎪⎪⎪⎪⎭⎫⎝⎛=1095791068565778710A ,当b T)31,33,23,32(=有微小误差)1.0,1.0,1.0,1.0(--=δb T 时,估计方程组b Ax =解的变化.2. 用MATLAB 求矩阵A 的2范数条件数、1范数条件数、∞范数条件数、Frobenius 范数条件数、条件数倒数的估计值,其中3. 用MATLAB 求第2题中矩阵A 的2范数、1范数、∞范数、Frobenius 范数、谱半径.4. 用MATLAB 求下列向量的2范数、1范数、∞范数、5范数、-∞范数:(1))5,2,1,1(-=X ;(2) )12,45,06.0,89,4.0(--=Y T. 5. 线性方程组求解与性态讨论.求b Ax =的解向量x ,其中然后把b 扰动为b ˆT )99.30,01.33,99.22,01.32(=,再求解b xA ˆˆ=. 计算xx x -ˆ(使用1范数或∞范数),讨论方程组性态..98755.03.02.064321543210⎪⎪⎪⎪⎪⎭⎫⎝⎛---=A .31332332,9105791068565778710⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=⎪⎪⎪⎪⎪⎭⎫⎝⎛=b A。

矩阵论简明教程第二版(张凯院著)课后答案下载

矩阵论简明教程第二版(张凯院著)课后答案下载

矩阵论简明教程第二版(张凯院著)课后答案下载《矩阵论简明教程》是xx年科学出版社出版的图书。

以下是要与大家分享的矩阵论简明教程第二版(张凯院著),供大家参考!点击此处下载???矩阵论简明教程第二版(张凯院著)课后答案???矩阵的相似变换,范数理论,矩阵分析,矩阵分解,特征值的估计与表示,广义逆矩阵,矩阵的直积以及线性空间与线性变换。

各章均配有习题,书末有习题解答与提示。

与传统矩阵论教材不同的是,《矩阵论简明教程》不是从较抽象的线性空间与线性变换开始,而是以较具体的矩阵相似变换理论作为基础来介绍矩阵理论的主要内容,以达到由浅入深的目的,并使读者在较短时间内掌握近现代矩阵理论相当广泛而又很基本的内容。

学习过工科线性代数课程的读者均可阅读《矩阵论简明教程》。

[1]第一章矩阵的相似变换1.1特征值与特征向量1.2相似对角化1.3Jordan标准形介绍1.4IHamilton-CayIey定理1.5向量的内积1.6酉相似下的标准形习题1第2章范数理论2.1向量范数2.2矩阵范数2.2.1方阵的范数2.2.2与向量范数的相容性2.2.3从属范数2.2.4长方阵的范数2.3范数应用举例2.3.1矩阵的谱半径2.3.2矩阵的条件数习题2第3章矩阵第4章矩阵分解第5章特征值的估计与表示第6章广义逆矩阵第7章矩阵的直积第8章线性空间与线性变换习题解答与提示参考文献1.实分析与复分析WalterRudin著课后习题答案机械工业出版社2.计算机专业英语教程第4版金志权课后习题答案电子工业出版社3.矩阵论简明教程第二版徐仲张凯院著课后答案科学出版社。

线性代数 第六版

线性代数 第六版

线性代数第六版简介线性代数是一门研究向量空间及其上的线性变换的学科。

在数学领域中,它是一门重要而基础的学科,被广泛应用于各个领域,如物理学、计算机科学、经济学等等。

本文将介绍《线性代数第六版》这本书的内容和特点。

作者本书的作者是Gilbert Strang。

Gilbert Strang是美国的一位著名数学家,现任麻省理工学院应用数学系教授。

他的主要研究领域是应用数学和数值分析,特别是在线性代数的教学和应用领域有着丰富的经验。

内容概述本书是一本线性代数的教材,共分为十二个章节。

以下是每个章节的简略概述:1.第一章介绍了向量和矩阵的基本概念,包括向量的几何解释、矩阵运算和矩阵的性质。

2.第二章讨论了线性方程组和矩阵的消元法,以及矩阵的秩和求解线性方程组的方法。

3.第三章介绍了矩阵的逆和逆矩阵的性质,以及逆矩阵的求解方法。

4.第四章讨论了线性变换和坐标变换,以及线性变换对于矩阵的表示。

5.第五章介绍了特征值和特征向量的概念,以及对角化和相似矩阵的性质。

6.第六章讨论了正交向量和正交矩阵,以及正交矩阵的特性和应用。

7.第七章介绍了复向量空间和复数域上的线性代数,包括复数的运算和复向量的性质。

8.第八章讨论了对称矩阵和二次型,以及对称矩阵的对角化和奇异值分解。

9.第九章介绍了线性相关性和线性无关性,并讨论了向量空间的基与维数。

10.第十章讨论了正交补空间和投影运算,以及最小二乘问题的求解方法。

11.第十一章介绍了复数域上的正交矩阵和正交变换,以及复数域上的最小二乘问题。

12.第十二章讨论了内积空间和希尔伯特空间,包括内积、范数和正交性的概念。

特点《线性代数第六版》有以下几个特点:•简洁明了的叙述风格,易于理解和学习。

•丰富的例子和练习,帮助读者掌握概念和方法。

•强调线性代数与实际问题的联系,注重应用层面的讲解。

•提供了大量的实际应用案例,帮助读者将理论知识应用到实际中。

•给出了详细的解题步骤和解答,方便读者自学和复习。

复向量二范数-概述说明以及解释

复向量二范数-概述说明以及解释

复向量二范数-概述说明以及解释1.引言1.1 概述复向量二范数是线性代数中的重要概念,对于复数域上的向量进行度量和定义具有重要意义。

在实际应用中,复向量二范数被广泛用于信号处理、图像处理、统计学等领域。

本文将从理解向量及范数开始,逐步介绍复向量的定义和复向量二范数的计算方法。

通过深入探讨,读者将能够全面了解复向量二范数的重要性及其在不同应用领域的实际意义。

同时,我们也将展望未来研究方向,希望能够激发更多学者对复向量二范数的研究和应用。

1.2 文章结构文章结构部分是关于整篇文章的框架和组织方式的介绍,它用来向读者展示本文的主要内容和逻辑结构。

在本篇文章中,文章结构部分应当包括以下内容:文章结构部分内容:本文主要分为引言、正文和结论三个部分。

1. 引言部分包括概述、文章结构和目的。

在这部分中,我们将介绍复向量二范数的背景和意义,以及本文的主要内容和目的。

2. 正文部分分为三个小节:理解向量及范数、复向量的定义和复向量的二范数。

在这部分中,我们将介绍向量、范数以及复向量的定义和性质,最终重点介绍复向量的二范数。

3. 结论部分包括总结复向量二范数的重要性、应用领域和展望未来研究。

在这部分中,我们将总结复向量二范数的意义,探讨其在实际应用中的作用,同时也展望未来对复向量二范数相关研究的发展方向。

通过以上文章结构的介绍,读者可以清晰地了解本文的结构和内容安排,帮助他们更好地理解和阅读整篇文章。

1.3 目的本文旨在深入探讨复向量二范数的概念和重要性,帮助读者更加全面地理解和应用复向量的数学理论。

通过对复向量的定义和二范数的推导,我们可以更好地理解复数空间中的向量运算和性质。

此外,通过这篇文章的阐述,我们也希望为读者提供一些关于复数空间在实际问题中的应用领域和未来研究方向的思考,以启发更多对复向量二范数的研究与探讨。

通过本文的阅读,读者可以更深入地了解复向量二范数的内涵和意义,为进一步深入学习和研究复数空间打下坚实的基础。

运筹与优化第1-2章

运筹与优化第1-2章

◆锥、凸锥 定义2 . 1 . 2 设有集合C ⊂ E n,若对C 中每 一点x,当λ取任何非负数时,都有 λx∈C ,则称C 为锥,又若C 为凸集,则 称C为凸锥。
lim x
k→∞
(k )
=x
序列若存在极限,则任何子序列有相 同的极限,极限唯一
四、聚点
定义1.3.5设{x(k)}是Rn中的一个向量序 (k j ) 列,如存在一个子序列{x } ,使
k j →∞
lim x
(k j )
ˆ = x
ˆ 则称 x 是序列{x(k)}的一个聚点 • 如无穷序列有界,即存在正数M,使得 对所有k均有|| x(k) ||≤ M,则该序列必 有聚点
八、梯度
• 函数f在x处的梯度为n维列向量:
⎡ ∂ f ( x ) ∂f ( x ) ∇f ( x ) = ⎢ , , ∂x 2 ⎣ ∂x1
∂f ( x ) ⎤ , ⎥ ∂x n ⎦
T
九、Hesse矩阵
∂ f ( x) [∇ f ( x )]ij = ,1 ≤ i , j ≤ n ∂x i ∂x j
n×k
第2章 凸集与凸函数
§2.1凸集(Convex Set)
一、凸集的概念

实空间上定义了+、数乘、内积
x(1)和 x(2)的凸组合(convex combination x(1) x(2) x(1) x(2)
●例2.1.1 集合H = { x | pTx = α}为凸集, 其中,p为n维列向量, α为实数。
• 设每人每天需要各种食品的数量分别为x1、 x2、…、xn • Min c1x1+ c2x2+ …+ cnxn a11x1 +a12x2+ …+ a1nxn≥b1 a21x1 +a22x2+ …+ a2nxn≥b2 s.t. … am1x1 +am2x2+ …+ amnxn≥bm x1 、x2、 … 、 xn≥ 0

高等工程数学知到章节答案智慧树2023年南京理工大学

高等工程数学知到章节答案智慧树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.有界区域上的弦振动方程定解问题可以用傅里叶积分变换法求解。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
Axx T xx T
Axx T xx T
T
消去非0数
xx
||xxT||,即 得证明。
24
T
Axx
A xx
T
方阵谱半径与范数关系
定理:对任意的正数ε>0,存在某个矩阵范数||A|| 使得
A ( A)
定理:对任何一种矩阵范数||A||都有
k 1k
lim A
1
绝对值不等式
根据向量范数定义容易导出类似于绝对值不等式:
a b a b a b x y x y x y
定义不同的向量范数就可以得到不同的不等式!
2
Minkowski不等式
向量范数定义为:
x
p
n p xi i 1
1 p
p 1
如果下式成立则向量x,y相互正交。 0向量与任 何向量与此 正交!
4
( x, y ) 0
一些常用的向量范数
在向量空间Rn可以定义很多向量范数,其中有一些常用的:
2范数: x
2
2 x1ቤተ መጻሕፍቲ ባይዱ x2
2 xn
xT x
x1 x2 x x T x x1 xn
m x

x

M x

11
范数的等价性
定理:同一个向量空间中任意两个不同范数 ||*||α、||*||β都相互等价。
定义向量范数的目的就是为了研究向量序列的 收敛性问题!
12
向量序列的收敛性定义
利用向量范数可以简化向量序列的收敛性定义,给 向量序列的研究带来方便,特别是讲到多元方程组的 迭代法收敛性时,常常要考虑向量序列。 定义:对向量空间Rn中的向量序列
A max aij 5
1i m j 1
19
求矩阵2范数需要计算ATA特征值:
5 0 0 AT A 0 9 6 , I AT A ( 3)( 5)( 15) 0 0 6 9
因此得:
T A2 ( A A) max 15 1 2
( A) max j
j
特征值可 能是复数, 故要取绝 对值!
o
23
方阵谱半径与范数关系
定理:方阵谱半径不大于任何一种矩阵范数。即
( A) A
证明:设A的绝对值最大的特征值为λ,对应的特征向 量为x,注意到特征向量是非0的!
( A) , Ax x, x 0
20
定理5.8 矩阵2-范数与F-范数具有正交不变性。 即若U、V 都是n阶正交方阵,则有:
UAV
2
A2
UAV
F
AF
这种情况表明在正 交变换下2、F范数 不改变原来矩阵的 范数!但对无穷范 数、1范数却不具有 这种性质!
21
矩阵序列的收敛性定义
利用矩阵范数可以简化矩阵序列的收敛性定义,给 矩阵序列的研究带来方便,特别是讲到多元方程组的 迭代法收敛性时,常常要考虑矩阵序列。 定义:对矩阵空间Rn×n中的矩阵序列
列和矩阵范数:它定义每列矩阵元素绝对值之和中 最大者!
A 1 max aij
1 j n i 1
15
n
几个常用的矩阵范数
无穷矩阵范数:它定义每行矩阵元素绝对值之和中 最大者!
A max aij
1i n j 1
n
谱矩阵范数:它定义方阵ATA特征值之中最大者的 算术平方根!
利用内积与2范数可以证明著名的不等式:
( x, y) x 2 y
n 2
2
n 2 n 2 xi yi xi yi i 1 i 1 i 1
7
例题1:设m×n阶矩阵A为列满秩的,并且向量空间
Rm中有一个向量范数||*||α,证明下式为向量空间Rn中的范数。
Ax A x
18
例题:求下面矩阵的各种范数。
2 1 0 A 0 2 3 1 2 0
解:这些范数包括:frobenius、2、1、∞
1 2
A 1 max aij 5
1 j n i 1
n
m
AF
2 aij 23 i, j
x y

k Ax

k x
Ay


A( x y ) Ax


x

y

根据向量范数的定义可知:f(x)是一个Rn中的范数。
9
范数的等价性
同一个向量空间中两个不同范数虽然表达式不同, 但它们之间仍然有关系---等价关系。 定义:同一个向量空间中两个不同范数||*||α、 ||*||β,如果存在两个正数m、M使得下式成立则 这两个范数相互等价。
17
矩阵算子范数
A 2 max
x0
Ax x
2
2
max Ax
x 2 1
2
A 1 max
x0
Ax 1 x1 Ax x
max Ax 1
x 1 1
可以证明如下 三个结果,因 此这三种范数 都为算子范数!
A

max
x0
max Ax
x 1


定理:如果一个矩阵范数为算子范数,则必有
A(1) , A(2) ,
, A( k ) ,
如果存在矩阵空间Rn×n中的矩阵A,使得对某 种范数成立
lim A( k ) A 0
k
则称这个矩阵序列收敛于矩阵A。
表明:随着k 的增大,矩阵 A(k)与矩阵A 越来越接近!
22
方阵谱半径
方阵谱半径是指能包括它的所有特征值的最小圆 盘的半径即:
向量范数
在计算方法中经常要进行理论分析,经常涉及到 线性代数中学习过的向量空间,包括Rn、Rm×n。 定义:向量范数定义 向量范数是从向量空间Rn到实数集R的一个满足如下 三个条件的函数||*||:Rn→R:
(1) x 0; x 0 x 0 (2) kx k x (3) x y x y
x y
p
x
p
y
p
3
向量内积
定义:向量内积定义 向量内积是从向量空间Rn到实数集R的一个满足如下 四个条件的二元函数(*,*):Rn→R:
(1) ( x, y) ( y, x ) (2) (kx, y ) k ( x, y ) (3) ( x y, z ) ( x, z ) ( y, z ) (4) ( x, x ) 0;( x, x ) 0 x 0
x (1) , x (2) ,
, x( k ) ,
如果存在向量空间Rn中的向量x,使得对某种 范数成立
lim x ( k ) x 0
k
则称这个向量序列收敛于向量x。
表明:随着k 的增大,向量 x(k)与向量x越 来越接近!
13
矩阵范数
在计算方法中经常要进行理论分析,还经常涉及到 矩阵的度量问题,我们只讨论方阵范数定义。 定义:矩阵范数定义 矩阵范数是从矩阵空间Rn×n到实数集R的一个满足 如下四个条件的一元函数||*||:Rn×n→R:
x2
x1 x2 xn xn
5
一些常用的向量范数
无穷范数:
x
1范数:

max x1 , x2 ,
, xn max xi
1i n
x 1 x1 x2
xn xi
i 1
n
6
Cauchy-Schwarz不等式
A 2 max ( AT A) max ( AAT )
16
由向量范数导出矩阵范数
定理:设||*||是向量空间Rn中的一个向量范数,则
A max
x 0
Ax x
max Ax
x 1
确定了矩阵空间Rn×n中的一个矩阵范数,称它为算子范数。 由向量的2范数、无穷范数、1范数都可以导出相 应的矩阵范数,称为矩阵的2范数或谱范数、无穷 范数、1范数或列范数,它们都为算子范数。
x

Ax

, xR
n
证明:根据向量范数的定义来证明。
f ( x) x

Ax

, x Rn
是定义在向量空间Rn中的一个函数;它满足:
(1)
f ( x) 0 , x Rn
(2)
f ( x) 0 Ax 0 Ax 0 x 0
8
(3)
(4)
f (kx ) kAx
(1) (2) (3) (4)
A 0; A 0 A 0 kA k A A B A B AB A B
14
几个常用的矩阵范数
Frobenius矩阵范数:它定义所有矩阵元素的平方 和的算术平方根!
A
F

aij
i 1 j 1
n
n
2
trace( AT A) trace( AAT )
k
( A)
定理:如果方阵A的谱半径小于1,则方阵A的幂序列总收 敛于0阵!即
lim Ak 0
k
25
相关文档
最新文档