矩阵论文献翻译--5000字

合集下载

矩阵特征值、特征向量的研究【文献综述】

矩阵特征值、特征向量的研究【文献综述】

毕业论文文献综述数学与应用数学矩阵特征值、特征向量的研究一、前言部分数学作为一种研究问题的工具,大部分同学并未真正感受到它的实用价值,往往低估了数学对于学习知识及其解决问题的重要作用,或不会灵活运用数学这一工具去理解、解决问题.许多理论、规律、计算等若能灵活而有效地借助数学方法去剖析、推演,往往会有意外的收获[]1。

矩阵就是数学中的一小部分,英文名Matrix(SAMND矩阵)本意是子宫、控制中心的母体、孕育生命的地方,同时,在数学名词中,矩阵用来表示统计数据等方面的各种有关联的数据。

这个定义很好地解释了Matrix代码制造世界的数学逻辑基础。

在科学技术和工程应用中,矩阵理论的重要性和应用的广泛性是众所周知的,尤其是有了矩阵特征值、特征向量的各种求解及计算机的广泛使用和MATLAB等数学计算软件的迅猛普及为矩阵提供了更为广阔的发展和应用前景。

矩阵特征值、特征向量运用非常的广泛,在很多方面都有涉及。

本文将先从各种矩阵的特征值、特征向量求解方法和矩阵历史入手,从几个方面综述矩阵特征值、特征向量的应用[]2。

那什么是矩阵特征值、特征向量呢?定义:设A是N阶矩阵,如果数X和N维非零列向量x,使关系式Ax=Xx成立,那么,这样的数X就称为方阵A的特征值,非零向量x称为A的对应于特征值X的特征向量。

求特征值描述正方形矩阵的特征值的重要工具是特征多项式:说λ是A的特征值等价于λ) v = 0 (其中I是恒等矩阵)有非零解 (一个特征向量),因说线性系统 (A –iλ)=0。

此等价于行列式 det(A –i第一:运用MATLAB求解矩阵特征值、特征向量。

首先,我用下面的例子,来引导我们认识MATLAB在求解矩阵特征值、特征向量上的运用。

例1:对亏损矩阵进行 Jordan 分解[]5。

A=gallery(5) %MATLAB 设置的特殊矩阵,它具有五重特征值。

[VJ,DJ]=jordan(A); % 求出准确的特征值,使 A*VJ=VJ*D 成立。

矩阵研究毕业论文

矩阵研究毕业论文

0 引言为了利用矩阵研究线性变换, 希望能找到线性空间的基使线性变换在该基下的矩阵具有最简单的形式, 因此我们引进了特征值与特征向量. 特征值与特征向量在线性变换中起着举足轻重的作用, 充分利用特征值与特征向量的命题与性质对我们解题带来极大的帮助, 能使复杂的问题变的简单, 化简为易, 化繁为简. 本文就矩阵的特征值与特征向量在一些解题中的应用作了初步的探讨.1. 关于矩阵的特征值与特征向量的一般理论我们知道, 在有限维线性空间中, 取了一组基之后, 线性变换就可以用矩阵来表示. 为了利用矩阵来研究线性变换, 对于每个给定的线性变换, 我们希望能找到一组基使得它的矩阵具有最简单的形式. 从现在开始, 我们主要的来讨论, 在适当的选择基之后, 一个线性变换的矩阵可以化成什么样的简单形式. 为了这个目的, 先介绍特征值和特征向量的概念, 它们对于线性变换的研究具有基本的重要性.定义 1.1 设A 是数域P 上的一个n 阶方阵,若存在一个数P λ∈以及一个非零n 维列向量n x P ∈,使得Ax x λ=则称λ是矩阵A 的一个特征值,向量x 称为矩阵A 关于特征值λ的特征向量. 定义1.2 设A 是数域P 上一n 级矩阵, λ是一个文字. 矩阵A E -λ的行列式nnn n n n a a a a a a a a a A E ---------=-λλλλ212222111211,称为A 的特征多项式, 这是数域P 上的一个次多项式.设T 是n 维线性空间V 上的一个线性变换,求解T 的特征值与特征向量的方法可以分成一下三几步:1) 在线性空间V 中取一组基12,,,nξξξ, 写出/A 在这组基下的矩阵A ;2) 求出A 的特征多项式E Aλ-在数域P 中全部的根, 它们也就是线性变换/A 的全部特征值;3) 对于A 的每个特征值,j λ求其次线性方程组()0jI A X λ-=的一组基础解系:12,,,.t ηηη于是A 的属于jλ的全部特征值组成的集合是}{1122,0,1,2,,t t i i k k k k K k i t ηηη+++∈≠=例1 设V 是数域K 上3维线性空间,T 是V 上的一个线性变换,它在在V 的一个基1α,2α,3α下的矩阵A 是222214241A -⎛⎫ ⎪=-- ⎪⎪-⎝⎭,求A 的全部特征值与特征向量. 解: 因为特征多项式为2222214(3)(6)241I A λλλλλλ--⎛⎫ ⎪-=+-=-+ ⎪⎪+⎝⎭所以A 的全部特征值3(二重),-6.对于特征值3,解齐次线性方程组(3)0I A X -=,12312312322024402440x x x x x x x x x +-=⎧⎪+-=⎨⎪++=⎩得到一个基础解系:210-⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦, 201⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦因此,A 的属于3的两个线性无关的特征向量就是1122ζαα=-+,2132ζαα=+ 而A 的属于3的全部特征向量就是 .{}11221212,,,0k k k k K k k ζζ+∈且不全为对于特征值-6代入, 求出(6)0I A X --=的一个基础解系:122⎡⎤⎢⎥⎢⎥⎢⎥-⎣⎦.因此, A 的属于特征值-6的一个线性无关的特征向量就是312322ζααα=+-,而A 的属于特征值-6的全部特征向量是{}3,0k k K k ζ∈≠且.例2 设T 是复数域上n 维线性空间V 上的一个线性变换()1n >,它在V 的一个基12,,,nααα下的矩阵A 是1012210000010000001n n A ααααα--⎛⎫ ⎪ ⎪⎪= ⎪ ⎪ ⎪-----⎝⎭,求T 的特征值和特征向量. 解 :1012201221100001000100001000010000100001n n n n I A λλλλααααλαλαλαλαλαλα-------=-+--=--+令01221000100001000001001n n n D λαλαλαλαλα----=--+下面用数学归纳法求解()2n D n ≥当2n =时,22101.1D λαλαλαλα==++-+假设对于上述形式的1n -阶行列式,有012-132000100001000001001n n n D λαλαλαλαλα----=--+n-1n-2n-210=+++λαλαλα,对于n 阶行列式,把它第1行展开,得12102112111210121210000100010010(1)001000100101()(1)(1).n n n n n n n n n n n n D xλαλαλλαλαλαλλλλαλαλααλαλαλαλα+----+----=---+----+-=+++++--=++++根据数学归纳法原理,此命题对一切自然数2n ≥都成立. 故121210.n n n I A λλαλαλαλα---=++++即为T 的特征多项式.设12,,n λλλ 是I A λ-的全部复根. 对于1i n ≤≤,有111122201111,n n n n i i i i i i i ii i i n i A λλλλλλλλλλααλαλ-----⎛⎫⎛⎫⎛⎫⎪ ⎪ ⎪ ⎪ ⎪ ⎪⎪ ⎪ ⎪== ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪----⎝⎭⎝⎭⎝⎭ 因此12'(1,,,,)n i i i λλλ-(1i n ≤≤)是A 的属于特征值i λ的一个特征向量. 由于()()11,2,,110,2,3,,n i n I A n λ--⎛⎫-=-≠⎪⎝⎭而i I A λ-=,因此()1i rank I A n λ-=-. 从而齐次线性方程组()0i I A X λ-=的解空间的维数为(1)1n n --=. 于是A 的属于特征值i λ的所有特征向量组成的集合是{}21'(1,,,,)|,0.n i i i k k C k λλλ-∈≠从而T 的属于特征值i λ的全部特征向量是{}21'123()|,0.n i i i k k C k αλαλαλ-++++∈≠(1i n ≤≤)例2 在空间[]nP x (n>1)中(P 为实数域), 求微分运算D'()()f x f x ∂= 的 特征多项式,并证明:D 在任何一组基下的矩阵不可能是对角矩阵. 证:在[]nP x 中取一组基()211,,,,2!1!n x x x n --微分运算D 在此基下的矩阵为.0000100001000010⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡=DD 的特征多项式是.01000010001n D E λλλλλ=---=-从而D 的特征多项式为nλ. 因此D 的特征值为210n λλλ====.又D 的对应特征值0的奇次线性方程组()0A X -=的系数矩阵的秩为n-1,从而基础解系只含一个向量.它小于[]nP x 的维数n(n>1),故D 不可能同任何对角矩阵相似.所以微分运算D 在任何基下的矩阵都不可能是对角形. 2矩阵特征值与特征向量的五个应用2.1特征值与特征向量判断线性变换可对角化的应用定义2.1.1如果V 中存在一个基,使得线性变换A 在这个基下的的矩阵是对角矩阵,那么A 可对角化.由于线性变换A 在V 的不同基下的矩阵是相似的,因此线性变换A 可对角化当且仅当A 在V 的基下的矩阵A 可对角.定理2.1.1域F 上n 维线性空间V 上线性变换A 可对角化当且仅当A 有n 个线性无关的特征向量12,,,nξξξ,此时A 在基12,,,nξξξ下的矩阵A 为1000,00n λλ⎛⎫ ⎪⎪ ⎪⎝⎭其中i λ是i ξ所属的特征值(即i i i A ξλξ=),1,2,,.i n = 矩阵A 称为线性变换A 的标准形,除了主对角线上元素的排列次序外,A 的标准形是有A 唯一决定的.推论2.1.1 域F 上n 维线性空间V 上线性变换A 可对角化当且仅当V 中存在由A的特征向量组成的一个基.定义2.1.2设A 是域F 上线性空间V 上的一个线性变换,0λ是A 的一个特征值,令 {}00|,defV A V λααλαα==∈ .易验证V λ 是V 的一个子空间,称0V λ是A 的属于特征值0λ的特征子空间. 0V λ中全部非零向量就是A 的属于特征值0λ的全部特征向量. 由于()00000().V A I A Ker I A λααλαλααλ∈⇔=⇔-=⇔∈-因此 00().V Ker I A λλ=-即线性变换A 的属于特征值0λ的特征子空间等于线性变换0I A λ- 的核.设V 是域F 上n 维线性空间,V 上线性变换A 在V 的一个基12,,,nααα下的矩阵为A,λ是A 的一个特征值. 设σ是V 到nF 的一个同构映射,它把V 中向量对应于它在基12,,,nααα下的坐标,则()0V λσ等于n 元齐次线性方程组()00I A X λ-=的解空间,即矩阵A 的属于特征值0λ的特征子空间. 于是()()00dim V n rank I A λλ=-- .定理2.1.2设A 是域F 上n 维线性空间V 上的一个线性变换,则A 可对角化⇔A 有n 个线性无关的特征向量⇔V 中存在由A 的特征向量组成的一个基⇔A 的属于不同特征值的特征子空间的维数之和等于n 12,s V V V V λλλ⇔=⊕⊕⊕其中12,,,sλλλ 是A 的所有不同的特征值.例 3 设T 是复数域上n 维线性空间V 上的一个线性变换()1n >,它在V 的一个基12,,,nααα下的矩阵A 是1012210000010000001n n A ααααα--⎛⎫ ⎪ ⎪⎪= ⎪ ⎪ ⎪-----⎝⎭,称它是Frobennis 矩阵. 求T 的特征多项式和属于特征值i λ的全部特征向量(1,2,3,,)i n =;T 是否可对角化?令122221211112111n n n n n n P λλλλλλλλλ---⎛⎫⎪ ⎪ ⎪= ⎪⎪ ⎪⎝⎭情形112,,n λλλ两两不等. 此时0.p ≠从而P 的列向量组线性无关. 于是A 有n 个线性无关的特征向量,因此A 可对角化.此时{}112,,n p AP diag λλλ-=从而T 可对角化.情形 212,,n λλλ中有相等的. 此时0.p = 从而P 线性相关. 这时A 没有n 个线性无关的特征向量,因此A 不可对角化, 从而T 不可对角化.例4 设T 是数域K 上n 维线性空间V 上的对合变换(即T 满足2T I =),(1)证明T 有特征值,且它的特征值是1或-1.(2)判断T 是否可对角化;若可以对角化,请写出它的标准形. 解:设T 在V 的一个基12,,,nααα下的矩阵是A ,由2T I =,可得2A I =. 即A 是数域K 上的对合矩阵,设0λ是对合矩阵A 的一个特征值,则有0,α≠使0.A αλα=从而2200.A A αλαλα== 由于2A I =,因此20αλα=,即20(1)0.λα-=由于0,α≠因此2010.λ-= 即01.λ=± 当A I =时,1是A 的特征值,-1不是;当A I =-时,-1是A 的特征值,1不是; 当A I ≠±时,0.I A ±≠由于()()rank I A rank I A n -++=因此 ()().rank I A n rank I A n -=-+< 从而0.I A -=从而1是A 的一个特征值.同理可证,-1是A 的一个特征值.(1)从而,T 有特征值,且它的特征值是1或-1.(2)设().rank I A r +=由于()()rank I A rank I A n -++=,因此().rank I A n r -=- 属于特征值1的特征子空间1W 的维数为1dim ()();W n rank I A n n r r =--=--=属于特征值-1的特征子空间1W -的维数为1dim ()();W n rank I A n rank I A n r -=---=-+=-由于11dim dim (),W W r n r n -+=+-=因此A 可对角化.A 的相似标准形为{},.r n r diag I I --从而T 可对角化,且它的相似标准形为0,0rn r I I -⎛⎫ ⎪-⎝⎭其中().r rank I A =+2.2 特征值与特征向量在确定可对角化矩阵的应用当矩阵A 可对角化时,可根据A 的特征值和特征向量来确定它的元素.例 5 设3阶方阵A 的特征值1231,0,1,λλλ===-对应的特征向量分别是1231222,2,1.211ξξξ-⎛⎫⎛⎫⎛⎫ ⎪ ⎪ ⎪==-= ⎪ ⎪ ⎪ ⎪ ⎪ ⎪⎝⎭⎝⎭⎝⎭求A .分析:此题给了3阶矩阵A 的3个不相同的特征值及其对应的特征向量,那么矩阵A 可对角化,显然可用A 的特征值和特征向量来确定它的元素.解:由i ξ是方阵A 对应于特征值i λ 的特征向量,于是i i i A ξλξ=()1,2,3.i =令()123122221212P ξξξ-⎛⎫⎪==-- ⎪⎪⎝⎭,则112212219212P -⎛⎫ ⎪=- ⎪⎪--⎝⎭, ,PA PD =其中100000,001D ⎛⎫ ⎪= ⎪ ⎪-⎝⎭ 由上式可得:11021012,3220A PDP --⎛⎫ ⎪== ⎪⎪⎝⎭ 即为所求.2.3特征值与特征向量在n 阶矩阵的高次幂的求解中的应用当n 阶矩阵A 可对角化时,即矩阵A 可与对角阵相似时,可应用矩阵的特征值与特征向量计算其高次幂()k A k N *∈,且比较简单.当n 阶矩阵A 满足下面的四个条件之一时,即可对角化,即1.A PDP -=n 阶矩阵A 有n 个线性无关的特征向量. n 阶矩阵A 有n 个互不相等的特征值.n 阶矩阵A 的每个特征值的几何重数等于其代数重数. A 为是对称矩阵. 对于(){}11212,,,,,,,,n n A PDP P D diag ξξξλλλ-===其中12,,,nλλλ是A 的n 个互不相等的特征值,i ξ是A 的属于特征值i λ的特征向量()1,2,,.i n =例6 已知矩阵122212221A ⎛⎫ ⎪= ⎪⎪⎝⎭,求k A (其中k N *∈). 分析:矩阵的高次幂的求解一般是有技巧的,这里因为矩阵A 为实对称矩阵,故可对角化. 可按上面讨论的方法求之.解 因为,T A A =所以矩阵A 为实对称矩阵,故A 可对角化为D .()()212221251221I A λλλλλλ----=---=-+---故A 的特征值为1231,5,λλλ==-=当1λ=-时,解齐次线性方程()0,I A X --=求出一个基础解系:12111,001ηη--⎛⎫⎛⎫⎪ ⎪== ⎪ ⎪⎪ ⎪⎝⎭⎝⎭ 当5λ=时,可求()50A X λ-=的一个基础解系:311,1η⎛⎫ ⎪= ⎪ ⎪⎝⎭ 令111101,011P --⎛⎫ ⎪= ⎪ ⎪⎝⎭()1001,1,5010005D diag -⎛⎫ ⎪=--=- ⎪⎪⎝⎭ 则()11,1,5P AP D diag -==--则1A PDP -=于是()()()()()()()()1111111111111()()1001112111101010121301100511121515151152153k kkkkk k k k k k k k k k k A PP APP PP APP PP APP P P AP P AP PAP P -------------==⎛⎫----⎛⎫⎛⎫⎪ ⎪ ⎪ ⎪=--- ⎪ ⎪ ⎪ ⎪ ⎪⎪⎝⎭⎝⎭⎝⎭-+-+-+=-+-+-()()()()111151515215k kk k k k k k---⎛⎫⎪ ⎪+ ⎪ ⎪-+-+-+⎝⎭2.4 特征值与特征向量在求一些特殊数列通项公式的应用由一些特殊数列的递推公式,构造关系矩阵A ,并列出递推关系,当关系矩阵A 可对角化时,可利用A 的特征值与特征向量求解这些数列的通项公式.例7 斐波那契(Fibonacci )数列是0,1,1,2,3,5,8,13,它满足下列递推公式:21,n n n ααα++=+ 0,1,2,n=以及初始条件010, 1.αα== 求Fibonacci 数列的通项公式,并且求1lim.nn n αα→∞+解 由2111,,n n n n n ααααα++++=+⎧⎨=⎩ 可得21111.,10n n n n αααα+++⎛⎫⎛⎫⎛⎫= ⎪ ⎪ ⎪⎝⎭⎝⎭⎝⎭令11,10A ⎛⎫= ⎪⎝⎭ 1,0,1,2,n n n D n αα+⎛⎫== ⎪⎝⎭上式可写成1,n n D AD +=又由1001,0D αα⎛⎫⎛⎫== ⎪ ⎪⎝⎭⎝⎭ 所以0,.n n D A D n N *=∈于是求Fibonacci 数列的通项公式就只要去计算nA .可利用A 的相似标准形来求简化nA 的计算.211111122I A λλλλλλλ⎛⎫⎛---==--=-- ⎪ ⎪-⎝⎭⎝⎭于是A的特征值为12λλ==从而A 可对角化.对于特征值1λ,解奇次线性方程组()10,I A X λ-=求出一个基础解系:11,1λη⎛⎫= ⎪⎝⎭对于特征值2λ,可求出()20I A X λ-=的一个基础解系:22,1λη⎛⎫= ⎪⎝⎭ 令12,11P λλ⎛⎫= ⎪⎝⎭ 则1120,0P AP λλ-⎛⎫= ⎪⎝⎭从而12121121211212112010011101.1nn nn n n n n A P P λλλλλλλλλλλλλλ-++-⎛⎫⎛⎫⎛⎫⎛⎫== ⎪⎪ ⎪-⎝⎭⎝⎭⎭⎝-⎛⎫⎛⎫=⎪ ⎪-⎝⎭⎭由于110n n n A αα+⎛⎫⎛⎫= ⎪ ⎪⎝⎭⎝⎭因此))2121211110.n nn n n n nλαλλλλλ-⎛⎫⎛⎫==- ⎪⎪-⎝⎭⎝⎭⎡⎤⎢⎥=-⎢⎥⎝⎭⎝⎭⎣⎦即为Fibonacci 数列的通项公式. 于是211211112212111lim lim lim112nn nnnn nn n nnλλαλλαλλλλλλλ++→∞→∞→∞+⎛⎫- ⎪-⎝⎭==-⎛⎫- ⎪⎝⎭==例8已知()11,1,2i ii i ib cc b c--=⎧⎪⎨=+⎪⎩其中2,3,.i =设11,b c已知,求,.n nb c解由题可得1101,2,3,1122i ii ib bic c--⎛⎫⎛⎫⎛⎫⎪==⎪ ⎪⎪⎝⎭⎝⎭⎝⎭令01,1122B⎛⎫⎪=⎪⎝⎭则111,n nnb bBc c-⎛⎫⎛⎫=⎪ ⎪⎝⎭⎝⎭下面求1n B-.()111.11222I Bλλλλλ-⎛⎫-==-+⎪--⎝⎭因此B的全部特征值是11,.2-从而B可对角化.对于特征值1,解奇次线性方程组()0,I B X-=得到它的一个基础解系:11,1ξ⎛⎫= ⎪⎝⎭对于特征值1,2-解齐次线性方程组10,2I B X ⎛⎫--= ⎪⎝⎭得到它的一个基础解系:22.1ξ-⎛⎫= ⎪⎝⎭ 令12,11P -⎛⎫= ⎪⎝⎭ 则110.102P BP -⎛⎫⎪= ⎪-⎝⎭ 从而1111122111010210121211111130211122213111222n n n n n n n n B P P ---------⎛⎫⎪= ⎪-⎝⎭⎛⎫-⎛⎫⎛⎫ ⎪= ⎪ ⎪ ⎪⎛⎫--⎝⎭⎝⎭ ⎪ ⎪⎝⎭⎝⎭⎛⎫⎛⎫⎛⎫--+- ⎪ ⎪⎪⎝⎭⎝⎭ ⎪= ⎪⎛⎫⎛⎫ ⎪--+- ⎪ ⎪ ⎪⎝⎭⎝⎭⎝⎭因此22111111111112,3232111112.3232n n n n n n b b c c b c ----⎧⎡⎤⎡⎤⎛⎫⎛⎫=--++-⎪⎢⎥⎢⎥ ⎪ ⎪⎝⎭⎝⎭⎢⎥⎢⎥⎪⎣⎦⎣⎦⎨⎡⎤⎡⎤⎪⎛⎫⎛⎫=--++-⎢⎥⎢⎥ ⎪ ⎪⎪⎝⎭⎝⎭⎢⎥⎢⎥⎣⎦⎣⎦⎩2.5特征值与特征向量行列式计算中的应用用矩阵的特征值和特征向量计算三对角形的方法如下:设00000000000n a b c a b c a D a b ca =按第一行展开,得:12,n n n D aD cbD --=- 3,4,n =上式可写成21,n n n D aD cbD ++=- n N +∈由于2111,,n n n n n D aD cbD D D ++++=-⎧⎨=⎩ 令2111,,,10n n n n n n D D a cb d d A n N D D +++++-⎛⎫⎛⎫⎛⎫===∈ ⎪ ⎪ ⎪⎝⎭⎝⎭⎝⎭因此111,n nn n d Ad d A d +-==()2,3,n =其中2211D a cb d D a ⎛⎫-⎛⎫== ⎪⎪⎝⎭⎝⎭ 这样求nD 的问题就转化为nd 的问题,因而转化为求1,n A -即存在可逆矩阵P 使得 1P AP D -=(对角形),就可以算出1.n A -由201a cbI A a cb λλλλλ--==-+=-得A 的特征值12λλ==1) 若24a cb ≠① 若240,a cb -<则A 有两个不相等的复特征值12,,λλ在复数域上对应于12,λλ的特征向量分别为12,.ξξ取()12,P ξξ=则P 可逆 于是就有11111200n n n AP P λλ----⎛⎫=⎪⎝⎭所以111n n n n D d A d D+-⎛⎫== ⎪⎝⎭从而可求出nD .如果A 限制在实数域上,A 有复特征值,这时A 不可对角化.② 若240,a cb ->则A 有两个不同的特征值,则A 可对角化,按在复数域上的情况可求出nD2) 若24,a cb =这时A 有重根.若A 有两个线性无关的特征向量,则A 可对角化;若A 只有一个特征向量,这时可利用相似变换,把A 化若当标准形1100λλ⎛⎫ ⎪⎝⎭,可以算出1n A -,即可求出n D .例9 计算n 阶行列式:950004950004900.9500049n D =解:按第一行展开,得:12920,n n n D D D --=-()3,4,n =上式可写成21920,n n n D D D ++=-()n N +∈ 由2111920,,n n n n n D D D D D ++++=-⎧⎨=⎩ 令2111920,,,10n n n n n n D D d d A n N D D +++++-⎛⎫⎛⎫⎛⎫===∈ ⎪ ⎪ ⎪⎝⎭⎝⎭⎝⎭得111,n nn n d Ad d A d +-==()2,3,n =其中211619D d D ⎛⎫⎛⎫== ⎪ ⎪⎝⎭⎝⎭ 由于()()2920920451I A λλλλλλλ--==-+=---因此A 的特征值是124, 5.λλ==对于特征值14,λ=解其次线性方程组()40,I A X -=求出一个基础解系:14,1η⎛⎫= ⎪⎝⎭对于特征值25,λ=解其次线性方程组()50,I A X -=求出一个基础解系:25,1η⎛⎫= ⎪⎝⎭令45,11P ⎛⎫= ⎪⎝⎭ 则140,05P AP -⎛⎫= ⎪⎝⎭ 从而14005A P P-⎛⎫= ⎪⎝⎭111111111400545154011140554 5.4 4.554 5.4 4.5n n n n n n n n n n n n A P P---------⎛⎫= ⎪⎝⎭-⎛⎫⎛⎫⎛⎫= ⎪ ⎪ ⎪-⎝⎭⎝⎭⎝⎭⎛⎫--= ⎪--⎝⎭由于11619n n n D A D +-⎛⎫⎛⎫= ⎪ ⎪⎝⎭⎝⎭因此()11111161545.44.5549n n n n n n n D ----++⎛⎫=--=- ⎪⎝⎭例10 计算n 阶行列式:2120000121200012120000000210022n D ------=.解:将nD 按第一列展开得:1231232(2)22,n n n n n n n D D D D D D D ------=--+=+- ()4,5,6,n =上式可写成32122,n n n n D D D D +++=+-()n N *∈ 根据321221122,,,n n n n n n n n D D D D D D D D +++++++=+-⎧⎪=⎨⎪=⎩ 令323121*********,,100,5,0102n n n n n n n n D D D D D A D D D D ααα++++++-⎛⎫⎛⎫⎛⎫⎛⎫⎛⎫ ⎪ ⎪ ⎪ ⎪ ⎪===== ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪⎝⎭⎝⎭⎝⎭⎝⎭⎝⎭ 可得1,n n A αα+=11,n n A αα-=由于()()()2121011201I A λλλλλλλ---=-=-+-- 因此A 的特征值是1231,1, 2.λλλ==-= 对于特征值11,λ= 解其次线性方程组()0,I A X -=得到一个基础解系;111,1η⎛⎫ ⎪= ⎪ ⎪⎝⎭ 同理,分别可求231, 2.λλ=-=的一个特征向量23141,2,11ηη⎛⎫⎛⎫ ⎪ ⎪=-= ⎪ ⎪ ⎪ ⎪⎝⎭⎝⎭ 令114112,111P ⎛⎫ ⎪=- ⎪ ⎪⎝⎭ 则1100010002P AP -⎛⎫ ⎪=- ⎪ ⎪⎝⎭ 于是1100010002A P P -⎛⎫ ⎪=- ⎪ ⎪⎝⎭从而()()()()11111111111000100021001143361112010132611100220211233611121326202112n n n n n n n n n n n A P P -------+--⎛⎫ ⎪=- ⎪ ⎪⎝⎭⎛⎫-⎛⎫⎛⎫ ⎪ ⎪ ⎪=--- ⎪ ⎪ ⎪ ⎪ ⎪ ⎪- ⎪⎝⎭⎝⎭⎝⎭⎛⎫--⎛⎫ ⎪ ⎪ ⎪=-- ⎪ ⎪ ⎪- ⎪⎝⎭-⎝⎭于是()()()1121111123361011121325,62022112n n n n n n n n n D D D -+++--⎛⎫--⎛⎫⎛⎫⎛⎫ ⎪ ⎪ ⎪⎪ ⎪=-- ⎪ ⎪⎪ ⎪ ⎪ ⎪⎪- ⎪⎝⎭⎝⎭⎝⎭-⎝⎭从而()()()()()121013123 3.16 2.12562112263n n n nn n n n D -+⎛⎫ ⎪=-+-++-+-- ⎪ ⎪⎝⎭-=-++3.小结本文利用特征值与特征向量的一些命题和性质来探讨特征值与特征向量在一些解题计算中的应用,充分应用命题和性质给我们的解题带来很大的方便.参考文献[1] 大学数学系几何与代数教研室前代数小组.高等代数(第三版)[M].北京:高等教育出版社,2003.[2] 同济大学应用数学系. 工程数学- 线性代数(第4版) [M] . 北京:高等教育出版社,2003.[3] 奚传志. 矩阵的特征值与特征向量在行列式计算中的应用枣庄师专学报,1992年2期[4] 李淑花. 关于一类线性代数习题的快速解法[J]. 高等数学研究.[5] 谢国瑞. 线性代数及应用[M]. 北京:高等教育出版社,1999.[6] 戴华. 矩阵特征值反问题的若干进展[J]. 南京航空航天大学学报,1995.[7] 钱吉林.高等代数题解精粹[M].北京:中央民族大学出版社.[8]邵丽丽.矩阵的特征值和特征向量的研究.菏泽学院.计算机与信息工程系.山东菏泽(274015)[9] 朱凤娟.特征值与特征向量逆问题的研究[J].滨州学院学报2007.6 .[10] [英]S.巴比特. 科技工作者用矩阵方法[M] .北京:化学工业出版社.1984.126-137.[11]丘维声,高等代数(第二版)下册.北京:高等教育出版社[12] tephen H.Friedbeng等.Linear Algebra(4th Edition) [M].Prentice Hall/Pearson,1998.[13] Verler.W.J.Vectors Structures and Solutions of linear Matrix Equation, linear Algebra Appl;1975.[14]丘维声,高等代数(第二版)上册.北京:高等教育出版社[15] 熊全淹,线性代数[M].北京;高等教育出版社,1987.4.[16]丘维声,高等代数学习指导(下册).北京:清华大学出版社,2009[17]杨子胥,高等代数习题解(下册).济南:科学技术出版社,2009[18]丘维声,高等代数学习指导(上册).北京:清华大学出版社,2009致谢本学位论文是在我的指导老师张宝环老师的亲切关怀与细心指导下完成的.由于经验的匮乏,难免有许多考虑不周到的地方,从论文的选题、资料的搜集到论文的撰写编排整个过程中,张老师始终都给予了悉心的指导和不懈的支持,并为我指点迷津,帮助我开拓思路,精心点拨,热忱鼓励.张老师的一丝不苟的作风,严谨求实的态度,踏踏实实的精神,不仅授我以文,而且教我做人,给我以终生受益无穷之道.感谢老师们对我的教育培养.他们细心指导我的学习与研究.在此,我要向诸位老师深深地鞠上一躬.同时我要感谢同组的同学们,是我们相互的鼓励和支持才使得做论文的过程充满着快乐和感动.在此,我对所有帮助我的老师和同学们表达我衷心的感谢!。

矩阵分解的研究文献综述

矩阵分解的研究文献综述

矩阵分解的研究文献综述毕业论文文献综述数学与应用数学矩阵分解的研究一、前言部分(说明写作的目的,介绍有关概念,综述范围,扼要说明有关主题争论焦点)在近代数学、工程技术、经济理论管理科学中,大量涉及到矩阵理论的知识。

因此,矩阵理论自然就是学习和研究上述学科必不可少的基础之一。

矩阵理论发展到今天,已经形成了一整套的理论和方法,内容非常丰富。

矩阵分解对矩阵理论及近代计算数学的发展起了关键的作用。

寻求矩阵在各种意义下的分解形式,是对与矩阵有关的数值计算和理论都有着极为重要的意义。

因为这些分解式的特殊形式,一是能明显的反映出原矩阵的某些特征;二是分解的方法与过程提供了某些有效的数值计算方法和理论分析根据。

这些分解在数值代数和最优化问题的解决中都有着十分重要的角色以及在其他领域方面也起着必不可少的作用。

本文结合矩阵的基本知识原理,对矩阵分解的各种常用形式进行梳理、归纳,并举例进行说明。

矩阵的定义:由m n ?个数(1,2,,,1,2,,)ij a K i m j n ∈==K K 排成的m 行、n 列的长方形表111212122212n n m m mn a a a a a a a a a ?? ? ? ? ???K K M M O M K (1) 称为数域K 上的一个m n ?矩阵。

其中的ij a 称为这个矩阵的元。

两个矩阵相等就是它们对应位置的元全相等[1]。

矩阵通常用一个大写拉丁字母表示。

如(1)的矩阵可以被记为A .如果矩阵的行数m 与列数n 相等,则称它为n 阶方阵。

数域K 上所有m n ?矩阵的集合记为(),m n M K ,所有n 阶方阵的集合记为()n M K ,元全为0的矩阵称为零矩阵,记为0.矩阵A 的位于第i 行、第j 列的元简称为A 的(),i j 元,记为(),A i j 。

如果矩阵A 的(),i j 元是(1,2,,,1,2,,)ij a i m j n ==K K ,则可以写成()ij A a =。

矩阵理论应用论文

矩阵理论应用论文

高维随机矩阵理论在数组信号检测与估计中的应用摘要本文中,我们展示了高维随机矩阵理论在频谱中的要素、相关源的检测并解决了在大数组中的估计问题。

这些结果适用于样本空间的协方差矩阵R̂中所感测的数据。

可以看出,可以实现的检测样品尺寸大小小于传统方法所要求的。

如果确定了预定的方向,可以通过给R̂设置限制条件,包括从高维随机矩阵理论中提出的,可以得到更加准确的估计。

一组理论用来解决可行性问题。

讨论 了一些没有解决的问题。

问题声明我们认为,当p 很大时,检测映射在数列p (q<p )的传感器上的q 的数量以及他们的到达方向是个问题。

该模型的成像机制如下。

在每个时间t 的第j 个信号出现在场景中时,第i 个传感器的加性噪声和在第i 个传感器接收到的数据可以分别用平方可积的复数值随机变量序列S j (t)、N i (t)和X i (t)表示。

随机向量(S(t)=[S 1(t )……S q ])T,t ∈[0,+∞],ES (0)=0和奇异空间的协方差矩阵R S =ES(0)S(0)∗。

此外,假设随机变量序列(N i (t )|1≤i ≤p ,t ∈[0,+∞]),EN 1(0)=0和E |N 1(0)|2=σ2,σ2未知,与随机变量序列(S j (t )|1≤j ≤p ,t ∈[0,+∞])独立。

让N (t )=σW (t )=σ[W 1(t )……W p (t)]T (W i (t )被标准化)和X (t )=[X 1(t )……X p (t)]T 。

这些由阵列传感器收集的数据被建模成为随机向量的观测值X (t )=AS (t )+N(t)t ∈[0,+∞],A 是根据阵列的几何尺寸和信号参数的p*q 的矩阵,假设秩为q 。

在数据处理中的检测问题是从观测到的n 个快照(X (t i ))1≤i≤n 中估计q 。

根据上述假设,随机向量(X (t ))t∈[0,+∞]由空间的协方差R =EX (0)X (0)∗=AR S A ∗+σ2I p 决定,I p 表示p*p 的单位矩阵。

矩阵理论的应用

矩阵理论的应用

矩阵理论的应用摘要:矩阵是数学的基本概念之一。

作为线性代数的核心内容,矩阵广泛运用于各个领域,如数学建模、密码学、化学、通信和计算机科学等,解决了大量的实际问题。

关键词:矩阵;密码学;化学;数学建模;应用Abstract:Matrix is one of the fundamental conception in mathematics.As the core content in the linear algebra,It is used in various domains like mathematical modeling,cryptology,chemistry,communication&computer science,etc.and also solve a large amount of practical problems.Keyword:matrix,cryptology,chemistry,mathematical modeling,application. 一、引言矩阵理论在现代统计学的许多分支有着广泛的应用,成为统计学中不可缺少的工具,而且,随着研究的深入和应用的发展,矩阵与统计学之间的关系会越来越深刻。

一方面,统计学对矩阵研究提出了许多新的研究课题,刺激了有关矩阵理论研究的发展;另一方面,矩阵理论中的结果被越来越多地应用于统计学的理论研究及其应用中。

近三十年,许多统计学家致力于这方面的研究,并撰写了很多这方面的论文和著作,其中很多结论在统计学的研究中发挥着很大的作用。

矩阵理论在数值计算、线性规划、数据分析、科学试验、信号传输等重大领域也有着极其广泛的应用。

随着科技日新月异地进步,人类社会开始步入信息化、数字化时代,矩阵在生产实践中的应用越来越广泛,故矩阵理论的研究也就越来越重要。

二、矩阵理论在实际中的应用矩阵理论的应用是十分有必要,也是十分简便的。

它帮助我们解决了大量的实际问题,具体应用有如下几个方面:(1)在密码学中的应用古罗马时期,凯撒大帝为了避免信使在途中背杀以至于情报被敌军劫走,发明了一种方法,即,把明文中的每一个字母转化为英文字母表中的第四个字母。

矩阵理论论文

矩阵理论论文

矩阵分解在信号和图像处理方面的应用矩阵理论是一门发展完善、理论严谨、方法独特的理论基础课程,它对培养学生的逻辑能力、推理能力具有重要作用,但它又能广泛应用于各个领域。

矩阵理论主要内容包括线性空间、线性变换、范数理论;矩阵分析;矩阵分解;广义逆矩阵;特征值的估计以及广义特征值等。

用矩阵的理论和方法来处理现代工程技术中的各种问题已经越来越普遍。

下面简单介绍一下矩阵的奇异值分解在信号和图像处理方面的简单应用。

此方法近年来在数据降维和压缩,滤波器设网络节点估计、小波变换结果的后续处理等很多领域都获得了重要的应用。

在滤波器设计方面,VOZALIS等将SVD 用于协同滤波,他们的研究结果表明,SVD提高了协同滤波过程中预测的质量和精度。

而在消噪方面,LEHTOLA等利用SVD和数学形态学相结合,对心电信号(Electrocardiogram,ECG)进行处理,消除了噪声的影响,提高了心电图诊断的准确性。

同时奇异值分解已用于从孕妇皮肤测量信号中提取胎儿心电信号。

在另一些研究中SVD则被利用来实现特征提取和弱信号分离,如LIU等利用SVD从背景噪声强烈的振动信号中提取周期性冲击信息。

SVD在神经网络中也获得了应用,如TEOH等利用SVD实现了对隐层空间中模式的线性独立性分析,进而决定了隐层神经元节点的数目。

SVD的正交化特性在对小波和小波包变换结果的后续处理中也得到了有效的应用,如XIE等利用SVD对小波包分解后的肌电信号进行正交化处理,以获得代表肢体运动模式的最优特征,进而对肌电信号进行分类,用于对假肢的控制。

小波多分辨分析的本质就是把信号在一系列不同层次的空间上进行分解,获得相应的近似和细节信号,从而以不同的层次显示信号的各种概貌和细节特征[9],这种多分辨思想使得小波分析在很多领域获得了极为广泛的应用。

基于这种多分辨分析思想的思考,赵学智在SVD中提出了一种矩阵二分递推构造方法,根据该方法得到的SVD分解结果将分属于不同层次的空间,而且下一层次空间的基矢量是利用上一层次的近似基矢量而获得的,实现了利用SVD以不同的层次来展现信号的概貌和细部特征。

数学 外文翻译 外文文献 英文文献 矩阵

数学 外文翻译 外文文献 英文文献 矩阵

Assume that you have a guess U(n) of the solution. If U(n) is close enough to the exact solution, an improved approximation U(n + 1) is obtained by solving the linearized problemwhere have asolution.has. In this case, the Gauss-Newton iteration tends to be the minimizer of the residual, i.e., the solution of minUIt is well known that for sufficiently smallAndis called a descent direction for , where | is the l2-norm. The iteration iswhere is chosen as large as possible such that the step has a reasonable descent.The Gauss-Newton method is local, and convergence is assured only when U(0)is close enough to the solution. In general, the first guess may be outside thergion of convergence. To improve convergence from bad initial guesses, a damping strategy is implemented for choosing , the Armijo-Goldstein line search. It chooses the largestinequality holds:|which guarantees a reduction of the residual norm by at least Note that each step of the line-search algorithm requires an evaluation of the residualAn important point of this strategy is that when U(n) approaches the solution, then and thus the convergence rate increases. If there is a solution to the scheme ultimately recovers the quadratic convergence rate of the standard Newton iteration. Closely related to the above problem is the choice of the initial guess U(0). By default, the solver sets U(0) and then assembles the FEM matrices K and F and computesThe damped Gauss-Newton iteration is then started with U(1), which should be a better guess than U(0). If the boundary conditions do not depend on the solution u, then U(1) satisfies them even if U(0) does not. Furthermore, if the equation is linear, then U(1) is the exact FEM solution and the solver does not enter the Gauss-Newton loop.There are situations where U(0) = 0 makes no sense or convergence is impossible.In some situations you may already have a good approximation and the nonlinear solver can be started with it, avoiding the slow convergence regime.This idea is used in the adaptive mesh generator. It computes a solution on a mesh, evaluates the error, and may refine certain triangles. The interpolant of is a very good starting guess for the solution on the refined mesh.In general the exact Jacobianis not available. Approximation of Jn by finite differences in the following way is expensive but feasible. The ith column of Jn can be approximated bywhich implies the assembling of the FEM matrices for the triangles containing grid point i. A very simple approximation to Jn, which gives a fixed point iteration, is also possible as follows. Essentially, for a given U(n), compute the FEM matrices K and F and setNonlinear EquationsThis is equivalent to approximating the Jacobian with the stiffness matrix. Indeed, since putting Jn = K yields In many cases the convergence rate is slow, but the cost of each iteration is cheap.The nonlinear solver implemented in the PDE Toolbox also provides for a compromise between the two extremes. To compute the derivative of the mapping , proceed as follows. The a term has been omitted for clarity, but appears again in the final result below.The first integral term is nothing more than Ki,j.The second term is “lumped,” i.e., replaced by a diagonal matrix that contains the row j j = 1, the second term is approximated bywhich is the ith component of K(c')U, where K(c') is the stiffness matrixassociated with the coefficient rather than c. The same reasoning can beapplied to the derivative of the mapping . Finally note that thederivative of the mapping is exactlywhich is the mass matrix associated with the coefficient . Thus the Jacobian ofU) is approximated bywhere the differentiation is with respect to u. K and M designate stiffness and mass matrices and their indices designate the coefficients with respect to which they are assembled. At each Gauss-Newton iteration, the nonlinear solver assembles the matrices corresponding to the equationsand then produces the approximate Jacobian. The differentiations of the coefficients are done numerically.In the general setting of elliptic systems, the boundary conditions are appended to the stiffness matrix to form the full linear system: where the coefficients of and may depend on the solution . The “lumped”approach approximates the derivative mapping of the residual by The nonlinearities of the boundary conditions and the dependencies of the coefficients on the derivatives of are not properly linearized by this scheme. When such nonlinearities are strong, the scheme reduces to the fix-pointiter ation and may converge slowly or not at all. When the boundary condition sare linear, they do not affect the convergence properties of the iteration schemes. In the Neumann case they are invisible (H is an empty matrix) and in the Dirichlet case they merely state that the residual is zero on the corresponding boundary points.Adaptive Mesh RefinementThe toolbox has a function for global, uniform mesh refinement. It divides each triangle into four similar triangles by creating new corners at the midsides, adjusting for curved boundaries. You can assess the accuracy of the numerical solution by comparing results from a sequence of successively refined meshes.If the solution is smooth enough, more accurate results may be obtained by extra polation. The solutions of the toolbox equation often have geometric features like localized strong gradients. An example of engineering importance in elasticity is the stress concentration occurring at reentrant corners such as the MATLAB favorite, the L-shaped membrane. Then it is more economical to refine the mesh selectively, i.e., only where it is needed. When the selection is based ones timates of errors in the computed solutions, a posteriori estimates, we speak of adaptive mesh refinement. Seeadapt mesh for an example of the computational savings where global refinement needs more than 6000elements to compete with an adaptively refined mesh of 500 elements.The adaptive refinement generates a sequence of solutions on successively finer meshes, at each stage selecting and refining those elements that are judged to contribute most to the error. The process is terminated when the maximum number of elements is exceeded or when each triangle contributes less than a preset tolerance. You need to provide an initial mesh, and choose selection and termination criteria parameters. The initial mesh can be produced by the init mesh function. The three components of the algorithm are the error indicator function, which computes an estimate of the element error contribution, the mesh refiner, which selects and subdivides elements, and the termination criteria.The Error Indicator FunctionThe adaption is a feedback process. As such, it is easily applied to a lar gerrange of problems than those for which its design was tailored. You wantes timates, selection criteria, etc., to be optimal in the sense of giving the mostaccurate solution at fixed cost or lowest computational effort for a given accuracy. Such results have been proved only for model problems, butgenerally, the equid is tribution heuristic has been found near optimal. Element sizes should be chosen such that each element contributes the same to the error. The theory of adaptive schemes makes use of a priori bounds forsolutions in terms of the source function f. For none lli ptic problems such abound may not exist, while the refinement scheme is still well defined and has been found to work well.The error indicator function used in the toolbox is an element-wise estimate of the contribution, based on the work of C. Johnson et al. For Poisson'sequation –f -solution uh holds in the L2-normwhere h = h(x) is the local mesh size, andThe braced quantity is the jump in normal derivative of v hr is theEi, the set of all interior edges of thetrain gulation. This bound is turned into an element-wise error indicator function E(K) for element K by summing the contributions from its edges. The final form for the toolbox equation Becomeswhere n is the unit normal of edge and the braced term is the jump in flux across the element edge. The L2 norm is computed over the element K. This error indicator is computed by the pdejmps function.The Mesh RefinerThe PDE Toolbox is geared to elliptic problems. For reasons of accuracy and ill-conditioning, they require the elements not to deviate too much from beingequilateral. Thus, even at essentially one-dimensional solution features, such as boundary layers, the refinement technique must guarantee reasonably shaped triangles.When an element is refined, new nodes appear on its mid sides, and if the neighbor triangle is not refined in a similar way, it is said to have hanging nodes. The final triangulation must have no hanging nodes, and they are removed by splitting neighbor triangles. To avoid further deterioration oftriangle quality in successive generations, the “longest edge bisection” scheme Rosenberg-Stenger [8] is used, in which the longest side of a triangle is always split, whenever any of the sides have hanging nodes. This guarantees that no angle is ever smaller than half the smallest angle of the original triangulation. Two selection criteria can be used. One, pdead worst, refines all elements with value of the error indicator larger than half the worst of any element. The other, pdeadgsc, refines all elements with an indicator value exceeding a user-defined dimensionless tolerance. The comparison with the tolerance is properly scaled with respect to domain and solution size, etc.The Termination CriteriaFor smooth solutions, error equi distribution can be achieved by the pde adgsc selection if the maximum number of elements is large enough. The pdead worst adaption only terminates when the maximum number of elements has been exceeded. This mode is natural when the solution exhibits singularities. The error indicator of the elements next to the singularity may never vanish, regardless of element size.外文翻译假定估计值,如果是最接近的准确的求解,通过解决线性问题得到更精确的值当为正数时,( 有一个解,即使也有一个解都是不需要的。

数学专业英语第二版的课文翻译

数学专业英语第二版的课文翻译

2-A Why study geometry Many leading institutions of higher learning have recognized that positive benefits can be gained by all who study this branch of mathematics. This is evident from the fact that they require study of geometry as a prerequisite to matriculation in those schools. 许多居于领导地位的学术机构承认,所有学习这个数学分支的人都将得到确实的受益,许多学校把几何的学习作为入学考试的先决条件,从这一点上可以证明。

Geometry had its origin long ago in the measurement by the Babylonians and Egyptians of their lands inundated by the floods of the Nile River. The greek word geometry is derived from geo, meaning “earth” and metron, meaning “measure” . As early as 2000 . we find the land surveyors of these people re-establishing vanishing landmarks and boundaries by utilizing the truths of geometry . 几何学起源于很久以前巴比伦人和埃及人测量他们被尼罗河洪水淹没的土地,希腊语几何来源于geo ,意思是”土地“,和metron 意思是”测量“。

数学与应用数学专业英语翻译

数学与应用数学专业英语翻译

3-c (P38/P33)集合论在讨论任何数学的分支,如分析、代数、几何,最好使用符号和集合理论的专业术语。

在19世纪后期被布勒和康托尔发展的这个学科,已对20世纪数学发展有了深远的影响。

它具有统一的许多看似已断开连接的想法,并有助用优雅有系统的方式减少许多逻辑基础上的数学概念。

集合论彻底处理需要漫长的讨论,我们认为超出了这本书的范围。

幸运的是,基础的概念很少涉及数字,通过非正式的讨论他可能已经发展成为集合论的方法和想法的可操作的知识。

实际上,我们将讨论不是像一个有关精确术语的协议一样的一个新的理论,我们将或多或少的应用到熟悉的想法。

在数学中,'"集'"一词用于表示作为单个实体的对象的聚集,集合也被称为'"群'"部落,'一伙人'"'",团队',和‘选举团’这是所有集合的示例。

集合中的单个对象称为元素或集合的成员,他们都称为属于或包含于集合中。

反过来,集合包含或由其元素构成。

我们应该主要对数学对象的集合感兴趣:数集、曲线集、几何图形集等等。

在许多应用程序,通过对集合中的各个对象的性质进行没什么特别性质的假定去处理集合是方便快捷的。

这些称为抽象集。

抽象集理论已经发展到处理任意对象的集合,这种集合,并从这种普遍性理论派生其权力。

9-c(P80/P66)积分的基本属性由积分的定义,它可以推导出以下属性,证明有在第1.27 部分中。

定理1,线性性质对被积函数,如果f 和g 都在[a、b] 上可积,那么对于每一对常量c1 和c2有---。

而且有----------(积分等式成立)注解:利用数学归纳法,线性属性可推广,如下所示:如果--在【a,b】上可积那么对于所有实数-----有-----且---------(积分等式成立)定理2.区间可加性,如果以下三个积分的两个存在,那么第三个也存在,并且我们有...注意:在特别是,如果 f 是在[a、c]上是单调的和在[c、b]上也是单调的,且这两个积分存在,则【a,b】区间上的积分也存在且与其他两个积分的和相等定理3.平移不变性,如果f 在区间[a、b]上可积,那么对于任意实数c,都有… …。

【文献综述】应用矩阵的性质求解行列式

【文献综述】应用矩阵的性质求解行列式

【文献综述】应用矩阵的性质求解行列式文献综述数学与应用数学应用矩阵的性质求解行列式1.本课题的研究意义《高等代数》历来作为数学系各个专业的重要基础课,它在线性规划、离散数学、管理科学、计算机以及物理、化学等学科中也有极为广泛的应用;同时它也是学习相关专业课程的重要语言和工具。

矩阵理论是高等代数中的重要内容之一,而在矩阵理论中,方阵是最为重要的研究对象之一,方阵的可逆性在高等代数的许多领域有着举足轻重的作用。

在线性方程组的求解,线性空间结构问题,二次型的研究以及欧氏空间等等方面都可见其身影。

矩阵的可逆性研究离不开行列式的计算。

在行列式的计算中,当行列式转换成上三角行列式或者是下三角行列式,对角行列式和特殊行列式时计算就会会相对来说简单,于是在计算行列式时就尽量将其转化成三角型行列式,行列式的计算还有其他很多算法(降阶法,加边法,数学归纳法,按一行或一列展开法)还有些特殊的行列式还可以通过范德蒙公式来计算。

在行列式的计算中,运用了大量的矩阵的性质(矩阵的加法,减法,乘法,数乘,还用到了矩阵的分块,矩阵的秩)将行列式转换成三角型行列式。

所以本文从行列式和矩阵的相关性来阐述,运用矩阵的相关性质来求解行列式,以达到简化行列式,是复杂问题简单化,从而解出行列式。

2.目前国内、外的研究现状行列式的概念最早是由十七世纪日本数学家关孝和提出来的,他在1683年写了一部叫做《解伏题之法》的著作,标题的意思是“解行列式问题的方法”,书里对行列式的概念和它的展开已经有了清楚的叙述。

欧洲第一个提出行列式概念的是德国的数学家莱布尼茨。

德国数学家雅可比于1841年总结并提出了行列式的系统理论。

早在17世纪和18世纪初行列式就在解线性方程组中出现,1772年法国数学家范德蒙首先把行列式作为专门理论独立于线性方程之外进行研究,到了19世纪,是行列式理论形成和发展的重要时期,19世纪中叶出现了行列式的大量定理。

矩阵最早来于方程组的系数及常数所构成的方阵,这一概念由19世纪英国数学家凯利首先提出,林谨瑜运用分块矩阵的若干性质来解决行列式的计算。

矩阵分析小论文

矩阵分析小论文

浅谈正交矩阵与酉矩阵矩阵是数学中重要的基本概念,是高等代数的重要研究对象之一,也是数学与其它领域研究与应用的一个重要工具.矩阵是线性代数中的核心内容 ,而正交矩阵是一种较常用的矩阵 ,正交矩阵在矩阵论中占有重要地位,有着广泛的应用.对其本身的研究来说是富有创造性的领域.正交矩阵不仅在线性代数中,而且在理工各学科领域的数学方法中,如优化理论、计算方法、信息分析中都有着举足轻重的位置。

对矩阵性质的概括、改进和推广,以及对正交矩阵在数值分析中、矩阵分解中和对方程求解、数理统计中的应用的研究,对矩阵的理论研究有重要意义。

本文列举了正交矩阵与酉矩阵的一些常见的性质与定理,并对其应用进行了一些列举。

首先认识什么是正交矩阵,什么是酉矩阵。

酉矩阵的定义:n 阶复方阵U 的n 个列向量是U 空间的一个标准正交基,则U 是酉矩阵(Unitary Matrix)。

即若n 阶复矩阵A 满足条件:E A A AA H H ==(E 为单位矩阵,H A 表示“矩阵A 的共轭转置矩阵,即TH A A =”),则此时矩阵A 称为酉矩阵。

此时,容易验证,当矩阵A 、B 为酉矩阵时,则有如下的结论成立:(1)H A A =-1也为酉矩阵(2)1det =A(3)n n T U A ⨯∈,即T A 为酉矩阵(4)AB,BA 也均为酉矩阵正交矩阵的定义:正交矩阵是实数特殊化的酉矩阵。

如果实数矩阵A 满足E A A AA T T ==(E 为单位矩阵,T A 表示“矩阵A 的转置矩阵”),则n 阶实矩阵 A 称为正交矩阵。

此时,容易验证,当A 、B 为正交矩阵时,则有如下结论成立:(1)n n T E A A ⨯-∈=1,即1-A 、T A 均为正交矩阵(2)1det ±=A(3)AB,BA 也均为正交矩阵正交变换的定义:设A 是欧氏空间V 的一个线性变换,若A 保持向量的内积不变,即对于任意的α,β∈V 都有(A α,A β) = (α,β),则称A 为V 的正交变换。

矩阵论汉语

矩阵论汉语

PP波矩阵弦理论摘要:在平坦的背景上矩阵弦理论的简要回顾后,我们mulate矩阵字符串模型,在不同的PP波背景。

这将是完成两个常量和变量RR背景流量的情况下,为某些确切的字符串几何结构尝试。

我们表现出的字符串相互作用的非微扰的代表性,并展示如何特征值的隧道驱动器的WKB扩展到通常的微扰字符串也超对称的PP波背景情况下的互动。

诉讼RTN-研讨会“时空的量子结构和地理度量基本相互作用的性质“,鲁汶大学,2002年9月连接每矩阵弦理论矩阵弦理论[1](MST)概述turbative弦理论和中号(atrix)理论[2]解释如何理论字符串包括后者在第二个量化计划。

众所周知,11个维M理论compactified上一个圆S1(11)是给上升到ⅡA型字符串其中D -颗粒自称是自由的自然度的强耦合存储器在理论描述。

要在微扰政权能够达到一个弦理论,也许可以使用一个强/弱的两重性,即类型IIB级的S - selfduality。

因此,为了执行这样一个方案,一个是被迫进一步compactify理论另一个圆S1(9)能够达到一个T(9)对偶第IIB型弦理论。

Insodoing,在D -粒子映射到D型的字符串,其缠绕在第九方向计数其数量。

现在可以执行类型IIB级弦理论的S - selfduality在同一绕组在基本类型IIB级字符串映射的D -串第九届方向和弱耦合。

这些字符串后一个最后的t(9)对偶映射到弱耦合与固定光锥势头ⅡA型字符串。

请注意,我们才开始其中R与M理论的图片(11)= GSLS比R(9)是一个更大的圈-ⅡA型理论在强耦合机制- 当我们结束在M理论图片中,住宅(11)= GSLS是远远高于R(9)。

其实,链二元我们的表现可能已经得到一个翻转之间的M理论中的作用根本两个圆圈。

上述二重性链的结果是,微扰第IIA型字符串承认一个双描述包装类型IIB级ð字符串一旦光锥的势头字符串相关联的D弦清盘。

让我们通知从这个角度来看,整体的横向空间上八个立体几何中起着被动的角色,它可以作为一个任意的字符串背景下离开。

矩阵论论文

矩阵论论文

西安理工大学研究生课程论文课程名称:矩阵论任课教师:XXX论文/研究报告题目:线性变换在电路方程中的应用完成日期:2014年11月5日学科:Xxxx学号:XXXXXXX姓名:XXX成绩:线性变换在电路方程中的应用摘要:电路分析中的坐标变换和复杂绕组变压器分析中所用的变压器变换都是电路方程的线性变换。

根据矩阵理论,对坐标变换和变压器变换进行了统一阐释。

坐标变换本质是一个方阵和对角阵的相似变换,变压器变换的本质是新变量对旧变量的表示,当变换矩阵的逆阵等于它的转置(共轭转置)阵时,坐标变换和变压器变换数学表示是相同的。

通过对电路方程系数矩阵和三角阵的相似变换,同时得到了三相 abc 坐标系和任意速度旋转两相 dq0 坐标系、瞬时值复数分量 120 坐标系、前进 - 后退 FB0 坐标系之间的变换矩阵。

这有助于在更加基础的理论层面上揭示和理解电路方程线性变换的本质,也为提出电路方程线性变换的新类型提供了思路。

关键词:电路方程;线性变换;坐标变换;变压器变换引言在交流电机等电路分析中,常用的坐标变换是指三相静止 abc 坐标系任意速度旋转两相 dq坐标系、瞬时值复数分量 120 坐标系、前进 - 后退 FB坐标系,以及它们对应的特殊坐标系的变量之间的相互转换。

电路方程坐标变换的主要目的是使电压、电流、磁链方程系数矩阵对角化和非时变化,从而简化数学模型,使分析和控制变得简单、准确、易行。

还有一类电路方程变换,其目的是用旧变量表示出新变量,例如变压器中由原边变量利用变比变换而来的副边变量,把这类电路方程变换称为变压器变换。

坐标变换已有很多文献进行了阐述,但这些阐述大都是基于物理概念的。

变压器变换在复杂绕组变压器的分析中得到了应用,但只是针对具体问题对其方法的具体应用,没有明确提出变压器变换的概念。

这些文献对坐标变换和变压器变换都缺乏在数学层面上予以统一论述。

这种特殊和具体的阐述,不便于将之作普遍化和一般化的理解,也就妨碍了对它的推广和发展。

matrix theory(矩阵理论)

matrix theory(矩阵理论)

两个波导之间的耦合通过耦合矩阵可以表示为11nn n n b a t b a t κκ**++''⎛⎫⎛⎫⎛⎫= ⎪ ⎪ ⎪-⎝⎭⎝⎭⎝⎭, 其中t 为κ为, β为传播系数,2/n βπλ=,其中n 为折射率,λ为光波长。

这里假设有效折射率近似等于折射率[3]。

1nn n b ta a κ+''=+, 11n n n b a t a κ**++'=-+, 由上面的两个方程可得11()n nn a ta b κ+''=-+,2211[()]n nn b t a t b κκ*+''=-++, 上式可以写成矩阵的形式,表示为122111()n n nn n na a a tb b b t t κκ+*+''-⎛⎫⎛⎫⎛⎫⎛⎫=≡ ⎪ ⎪ ⎪ ⎪''-+⎝⎭⎝⎭⎝⎭⎝⎭P 。

矩阵P 描述了耦合去的耦合情况,可表示为2211()t t t κκ*-⎛⎫=⎪-+⎝⎭P 。

假设221t κ+=光信号在环中传输时,会经历相位的变化,而且还存在功率损耗,故有n i nn a e b φ-'=, n i n n b e a φ'=, 同理,上式可以写成矩阵的形式,表示为00n n i nn n i n n n a a a e b b b e φφ-'⎛⎫⎛⎫⎛⎫⎛⎫=≡ ⎪ ⎪ ⎪ ⎪'⎝⎭⎝⎭⎝⎭⎝⎭Q .矩阵描述了相位变化和功率衰减00n ni i e e φφ-⎛⎫= ⎪⎝⎭Q , 上式中的n φ可表示为n n L n R i R φβπαπ=-,其中n R 为第n 环的半径,0L α=为单位长度上的衰减。

有式()和()可知11n n n n a a b b ++⎛⎫⎛⎫= ⎪ ⎪⎝⎭⎝⎭PQ 。

对于N 个环层叠的微环谐振器有1001122110100n N N N N n a a a A B b b b C D +--+⎛⎫⎛⎫⎛⎫⎛⎫=⋅⋅⋅= ⎪ ⎪ ⎪ ⎪⎝⎭⎝⎭⎝⎭⎝⎭P Q P Q P Q PQ P 当微环谐振器有N 个相同的环组成,如图所示,则有121N -===≡P P P P , 12N ===≡Q Q Q Q故10010100[]n N N n a a a A B b b b C D +-+⎛⎫⎛⎫⎛⎫⎛⎫== ⎪ ⎪ ⎪ ⎪⎝⎭⎝⎭⎝⎭⎝⎭P Q PQ P , 上式中10n a +=为Add ,1n b +为Drop ,0b 为Through ,0a 为输入。

矩阵论课程论文~

矩阵论课程论文~

研究生课程论文西尔维斯特及其矩阵理论课程名称矩阵论姓名郭辉学号1000203040专业检测技术与自动化装置任课教师刘强开课时间2009.09——2010.01教师评阅意见:论文成绩评阅日期课程论文提交时间:10年 3 月 4 日西尔维斯特及其矩阵理论摘要矩阵是伴随着其他理论的研究而产生的,众多数学家为其早期的发展做了大量的工作。

在此基础上,西尔维斯特创用了矩阵一词,引进了与矩阵有关的一些基本概念,给出了矩阵的一些重要结论和著名定理,为矩阵理论的发展做出了重要贡献。

关键词矩阵的早期发展西尔维斯特矩阵名词矩阵理论矩阵思想的萌芽由来已久,早在公元前一世纪,我国最重要的数学经典著作《九章算术》已能够相当成熟地运用矩阵形式解方程组,魏晋时期的数学家刘徽又在《九章算术注》中进一步完善,给出了完整的演算程序[1]。

但那时矩阵概念仅是用来作为线性方程组系数的排列形式解决实际问题,并没有建立独立完善的矩阵理论。

从18世纪早期到19世纪中叶,这种排列形式在线性方程组和行列式计算中应用日益广泛。

在逻辑上,矩阵的概念先于行列式的概念,而在历史上次序正相反[2],因此在矩阵引进的时候它的许多基本性质就已经非常清楚了。

行列式以及代数型的发展为矩阵理论进一步的发展提供了条件。

在矩阵发展的早期,矩形阵列本身并没有引起单独的注意但是,19世纪数学家们在其他数学领域的研究工作导致了矩形阵列更加形式的计算,促进了矩阵理论的诞生。

西尔维斯特在矩阵理论方面的贡献,不仅体现在对矩阵理论内容上的发展,即从不同领域的研究中发展出来的有关矩阵的概念,以及随之引起的相似、对角化和标准型的矩阵分类等等,还有其更深刻的地方:一方面他的工作使得当时比较零散的矩阵知识趋于系统化、理论化,为凯莱创立矩阵理论提供了有利条件;另一方面,西尔维斯特的行列式和矩阵的思想,为代数不变量理论的创立奠定了重要基础。

1.矩阵的早期发展矩阵的早期发展是伴随其他理论的研究而产生的。

矩阵论课程结业论文

矩阵论课程结业论文

浅谈矩阵论的发展在《九章算术》中用矩阵形式解方程组已相当成熟,但那时仅用它作为线性方程组系数的排列形式解决实际问题,并没有建立起独立的矩阵理论。

直到18 世纪末至19 世纪中叶,这种排列形式在线性方程组和行列式计算中应用日益广泛,行列式的发展提供了矩阵发展的条件。

矩阵的早期发展,除了矩阵理论在内容上的发展,即从不同领域的研究中发展出来的有关矩阵的概念,以及随之引起的相似、对角化和标准型的矩阵分类以外,还有矩阵发展中更深刻的一面,即西尔维斯特、凯莱等人在行列式和矩阵理论上的发展及思想,这为代数不变量理论的创立奠定了理论基础。

一、矩阵早期发展的社会与文化背景矩阵是数学中的一个重要的基本概念,是代数学的一个主要研究对象,也是数学研究和应用的一个重要工具。

“矩阵”这个词是由西尔维斯特首先使用的,他是为了将数字的矩形阵列区别于行列式而发明了这个述语。

而实际上,矩阵这个课题在诞生之前就已经发展的很好了。

从行列式的大量工作中明显的表现出来,为了很多目的,不管行列式的值是否与问题有关,方阵本身都可以研究和使用,矩阵的许多基本性质也是在行列式的发展中建立起来的。

在逻辑上,矩阵的概念应先于行列式的概念,然而在历史上次序正好相反。

英国数学家凯莱(A.Cayley,1821-1895) 一般被公认为是矩阵论的创立者,因为他首先把矩阵作为一个独立的数学概念提出来,并首先发表了关于这个题目的一系列文章。

凯莱同研究线性变换下的不变量相结合,首先引进矩阵以简化记号。

1858 年,他发表了关于这一课题的第一篇论文《矩阵论的研究报告》,系统地阐述了关于矩阵的理论。

文中他定义了矩阵的相等、矩阵的运算法则、矩阵的转置以及矩阵的逆等一系列基本概念,指出了矩阵加法的可交换性与可结合性。

另外,凯莱还给出了方阵的特征方程和特征根(特征值)以及有关矩阵的一些基本结果。

凯莱出生于一个古老而有才能的英国家庭,剑桥大学三一学院大学毕业后留校讲授数学,三年后他转从律师职业,工作卓有成效,并利用业余时间研究数学,发表了大量的数学论文。

矩阵秩论文

矩阵秩论文

矩阵的秩的一些结论的证明摘要矩阵是高等代数中主要的一个研究对象,它贯穿着整个高等代数的内容,而矩阵的秩作为矩阵最主要的特征,研究它的结论和性质就变得尤其重要.本文主要从矩阵的秩的结论和矩阵的秩的应用两方面介绍了矩阵的秩,并对矩阵的秩的大量性质进行了研究、证明及应用.其中包括矩阵的秩的求解和矩阵的秩的一些不等式,而且还涉及到了矩阵的秩在求解方程组和向量相关问题上的应用.关键词:矩阵的秩;矩阵的秩的定义;矩阵的秩的结论;矩阵的秩的应用The Conclusion of the Matrix rank’s proofAbstractMatrix is an object in the Advanced Algebra to be studied, which runs through the whole content of the Advanced Algebra, however, the rank of matrix as its main characteristics. The conclusions and the nature’s study become such an important part. The paper is divided into two parts to introduce the matrix, which are the conclusions and the application of the rank. At the same time, the nature of the rank has been studied, proved and used in the paper. Among the applications, including the solution to the rank and some inequality, the paper also includes the application of rank about solving the equations and questions of the vector correlation.Keywords: the rank of the matrix; the definition of the rank; the conclusion of the rank; the applications of the rank.目录引言............................................................................................................................................... - 1 -1.矩阵的秩的两种定义 ............................................................................................................ - 2 -2.引理....................................................................................................................................... - 2 -3.矩阵的秩的一些结论及其证明 ............................................................................................. - 4 -命题1 (4)命题2 (5)命题3 (6)命题4 (6)命题5 (7)命题6 (8)命题7 (9)命题8 (9)命题9 (10)命题10(F ROBENIUS不等式) (10)命题11 (11)命题12 (11)命题13 (12)命题14 (13)3.15命题15 (13)命题16 (14)命题17 (15)3.18命题18 (16)4.矩阵的秩的一些结论的应用............................................................................................... - 17 -总结............................................................................................................................................. - 21 -致谢............................................................................................................................................. - 22 -参考文献 ..................................................................................................................................... - 23 -引言矩阵的秩是高等数学中一个极其重要并广泛应用的概念,是高等代数的一个重要研究对象.因此,矩阵的秩的结论作为高等代数的一个重要工具已经渗透到各章节内容之中,它把高等代数的内容紧紧联系在一起,矩阵的秩作为矩阵的一个重要的本质属性则贯穿矩阵理论的始终.所以对于矩阵的秩的研究不仅能够帮助我们更好的学习矩阵,而且他是我们学习好高等代数各章节的有力保障.矩阵A中的最大阶不为零的子式的阶数就称为矩阵A的矩阵的秩,记为)rank或矩阵的秩)(A.从定义上看, 一个矩阵的秩, 就是一(A个数.事实上,若将矩阵A的每一行看成一个向量,每一列看成一个向量,则行向量组和列向量组中极大无关组中向量的个数是相等的,数量上等与矩阵的秩.若)n)n(mrank≤=,(A)(n(mmrank≤A=,则称A为行满秩的矩阵;若)则称A为列满秩矩阵.n阶方阵的秩等于n时称A为满秩矩阵或可逆矩阵.1. 矩阵的秩的两种定义矩阵的秩是反映矩阵固有特性的一个重要概念,秩是矩阵的一个非常重要的数值特征,是由F.G.Frobenius(1877)提出的.定义1设A 是任意矩阵.若O =A 则说A 的秩为0;若O ≠A 则A 的非零子式的最高阶数就称为A 的秩,记为秩A .定义2设在矩阵A 中有一个不等于O 的r 阶子式D ,则所有1+r 阶子式(如果存在的话)全等于O ,那么D 称为矩阵A 的最高阶非零子式,数r 称为矩阵A 的秩,记作)(A r .并规定零矩阵的秩等于0.2. 引理2.1引理1 A 、B 分别为n m ⨯和s t ⨯矩阵,则)()()(B r A r B A r +=⎥⎦⎤⎢⎣⎡O O 恒成立.[]1证明:设存在可逆矩阵1P ,2P ,1Q ,2Q 使得⎥⎦⎤⎢⎣⎡OO =⎥⎦⎤⎢⎣⎡O O ⎥⎦⎤⎢⎣⎡O O ⎥⎦⎤⎢⎣⎡P O O P s r D D Q Q B A 2121, 其中r D 、s D 分别是由r 个和s 个线性无关的单位向量组成,且[]TrD O 与[]Ts D O是线性无关的向量组,所以s r B r A r D D r s r +=+=⎥⎦⎤⎢⎣⎡OO )()()(,因此得出)()()(B r A r B A r +=⎥⎦⎤⎢⎣⎡O O .引理1中通过分块矩阵构建了秩与两个模块矩阵秩的和相等的矩阵,可以直观方便的通过分块矩阵运算来实现某些性质的证明,有效的简化了证明路径,为以下命题的证明即提供了一种方法,又提供了相应的结论.A 、B 分别为n m ⨯和s t ⨯矩阵,则)()()(B r A r BC A r +≥⎥⎦⎤⎢⎣⎡O 成立. 证明:由引理1得)()()(B r A r B A r +=⎥⎦⎤⎢⎣⎡O O , 因为)()(A r C A r ≥⎥⎦⎤⎢⎣⎡, 所以)()()(B r A r B C A r +≥⎥⎦⎤⎢⎣⎡O .且当)()(A r C A r =⎥⎦⎤⎢⎣⎡时,)()()(B r A r B C A r +=⎥⎦⎤⎢⎣⎡O . 如上证明对引理1做了补充和扩展,对于便于分块的矩阵的秩的确定提供了方法.存在n 阶矩阵A ,T i x x x x ),(21 =为0=Ax 解向量的极大无关组,则n x r A r =+)()(.证明:对方程组,00021212222111211⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡ n nn n n n n x x x a a a a a a a a a化简得()(),000000001001212222111211⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡-⨯-⨯n r n r n x x x k k k k k k得出(),r A r =方程组解为()()()⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡=--⨯-⨯--⨯100010001))((2111211r n r n r n r n r n x x x x x x x ,所以r n x r -=)(,即n x r A r =+)()(.该引理将矩阵秩的性质与方程组解维数联系起来,对于判断方程组解的维数或者通过方程组的解了解相乘矩阵的秩的问题提供了方法.注:引理部分为基础性命题,对以下证明过程起辅助作用,是为了便于以下命题的证明.以上证明过的引理下面的命题均可直接引用.以上命题对矩阵秩的范围,以及矩阵秩与极大线性无关组的关系进行了证明与阐述.3. 矩阵的秩的一些结论及其证明设A 是n 阶方阵,则0≠A 当且仅当n A r =)(.[]2证明:令)()(n r r A r <=,则A 与⎥⎦⎤⎢⎣⎡O O O I r 等价,即存在可逆矩阵P 、Q 使得⎥⎦⎤⎢⎣⎡O O O I =P rAQ , 取其行列式得0=⎥⎦⎤⎢⎣⎡O O O I =P =P r AQ Q A . 所以,当且仅当n A r =)(时,0≠A .该命题是互逆命题,即条件结论可互换,也就是说满秩与行列式非零是等价的,可根据有效条件判断行列式是否等于零或者是否满秩矩阵.矩阵的乘积的秩不超过各因子的秩.即:设A 是n m ⨯矩阵,B 是s n ⨯矩阵,则{})(),(m in )(B r A r B r ≤A .证明:设非零矩阵n m ij a A ⨯=)(,s n ij b B ⨯=)(.AB 可表示为A 的列向量的线性组合,即:⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⋅=ns n n s s n b b b b b b b b b AB 21222211121121),,(ααα, 所以)()(A r AB r ≤.AB 可表示为B 的行向量的线性组合,即:⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⋅⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=n ns n n s s a a a a a a a a a AB βββ 21212222111211,所以)()(B r AB r ≤.可得{})(),(m in )(B r A r B r ≤A .此证明将矩阵分为多个列向量或行向量来处理的,向量组AB 是列向量组A 通过矩阵B 的映射,同时也可说向量组AB 是行向量组B 通过矩阵A 的映射.上述证明说明了映射向量组不能增大基向量组的秩.若可逆矩阵P ,Q 使B PAQ =,则)()(B r A r =.证明:初等行变换与初等列变换不改变矩阵的秩, PAQ 即对矩阵A 进行列变换和行变换,所以)()()(B r A r PAQ r ==.该命题体现了初等变换的性质,以及相似矩阵的特别,对于较为复杂的矩阵的秩的求解起到简化作用,可以通过求解相似或等价矩阵的秩来实现.若2≥n ,则A 的伴随矩阵*A 的秩与A 的秩有如下关系:⎪⎩⎪⎨⎧-≤-===2)(01)(1n )()(*n A r n A r A r n A r 当当当. 证明:当2)(-≤n A r ,O =*A ,所以0)(*=A r ;当1)(-=n A r , 即⎪⎩⎪⎨⎧O ≠O=====-⨯-⨯⨯⨯-⨯⨯)1()1()2(2111n n n n j n j n A A A A A A , 其中,n j ,2,1=.1-≠n j所以1)(*=A r ;当n A r =)(,因为I ⨯=⨯A A A *,所以AA A 1*-=,因为1-A 为满秩矩阵,所以n A r =)(*.伴随矩阵是一特殊矩阵,可用以求解逆矩阵,伴随矩阵的秩与对应矩阵关系如上命题所示,可用以相互求解和验证秩的大小.两矩阵和的秩不超过两矩阵秩的和.即:设 A 、B 均为n m ⨯矩阵,则)()()(B r A r B A r +≤+.[3]证明:由分块矩阵的初等变换①2112⎥⎦⎤⎢⎣⎡+O −−→−⎥⎦⎤⎢⎣⎡O −−→−⎥⎦⎤⎢⎣⎡O O ++B B A AB A A B A bc bc br br ,则)()()(⎥⎦⎤⎢⎣⎡O O =⎥⎦⎤⎢⎣⎡+O ≤⎥⎦⎤⎢⎣⎡O +O O B A r B B A A r B A r , 由引理1得)()()(B r A r B Ar +=⎥⎦⎤⎢⎣⎡OO , 所以)()()(B r A r B A r +≤+.此证明过程用到了分块矩阵,分块矩阵使未知矩阵和方便分块的矩阵的变换变得简单,过程清晰,便于理解.分块矩阵初等变换的规则如下注解所示.上述过程证明了两矩阵和的秩小于两矩阵秩的和,可用于判断和矩阵的范围.注① ⎥⎦⎤⎢⎣⎡O −−→−⎥⎦⎤⎢⎣⎡O O +B A A B A br br 12表示将矩阵⎥⎦⎤⎢⎣⎡O O B A的第一行加到第二行上. ⎥⎦⎤⎢⎣⎡+O −−→−⎥⎦⎤⎢⎣⎡O +B B A AB A A bc bc 21表示将矩阵⎥⎦⎤⎢⎣⎡O B A A 的第二列加到第一列上.设A ,B 均为n m ⨯矩阵,则)()()()()(B r A r B A r B r A r +≤±≤-.证明: 由命题5得)()()(B r A r B A r +≤+,即).()()()()(B r A r B r A r B A r +=-+≤-则由分块矩阵的初等变换⎥⎦⎤⎢⎣⎡−−→−⎥⎦⎤⎢⎣⎡O ±−−→−⎥⎦⎤⎢⎣⎡O O ±+B B B AB B BA B B A bc bc br br 2121, 可得)()()(A r B B B A r B BA r ≥⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡OO ± , 由引理1得),()()(B r B A r B BA r +±=⎥⎦⎤⎢⎣⎡OO ±所以).()()(B r B A r A r +±≤即).()()(B A r B r A r ±≤-因此得出结论).()()()()(B r A r B A r B r A r +≤±≤-上述证明同样运用了分块矩阵初等变换.并进一步求解了和矩阵、差矩阵的秩的范围.设A 为n m ⨯,B 为s n ⨯的矩阵,则.)()()(n B r A r AB r -+≥证明:由分块矩阵的初等变换,①2121⎥⎦⎤⎢⎣⎡I -O −−−→−⎥⎦⎤⎢⎣⎡I O −−−→−⎥⎦⎤⎢⎣⎡I O O ⋅-⋅+n bc B bc n br A br n B A A AB AB得),()()(⎥⎦⎤⎢⎣⎡O O ≥⎥⎦⎤⎢⎣⎡I -O =⎥⎦⎤⎢⎣⎡I O O B A r B A r AB r n n 得),()()(B r A r n AB r +≥+即.)()()(n B r A r AB r -+≥通过分块矩阵的初等变换很方便的求出两矩阵积的秩大于等于两矩阵秩的和减去维数.设A 为n m ⨯,B 为s n ⨯的矩阵,满足O =AB ,则n B r A r ≤+)()(.[4]证明:存在极大无关解向量组T n x x x x ),,,(21 =,使得0=Ax ,由引理3得()n m x r A r ,m in )()(=+,因为O =AB ,所以B 为A 的解向量组T s ),,,(21βββ ,是A 极大无关向量组的线性组合,那么)()(x r B r ≤,得证n B r A r ≤+)()(.该命题用到引理3的结论,是通过构建方程组来确定矩阵的秩的特性,对于该命题还可以通过构矩阵的秩为两矩阵的秩和的分块矩阵来证明,实现较为复杂.设A 是n m ⨯阶矩阵,r A r =)(,从矩阵A 中任取s 行组成矩阵B ,则n s r B r -+≥)(.证明:设t B r =)(,把矩阵B 的t 个无关向量扩充到A 的一个极大无关向量组需要扩展t r -个向量,因为A ,B 不一定为满秩矩阵,所以s n m t r -≤-),min(,即n s r B r -+≥)(.3.10命题10(Frobenius 不等式)设A ,B ,C 分别为l m ⨯,s l ⨯,n s ⨯矩阵,证明)()()()(B r BC r AB r ABC r -+≥.证明: 构造如下矩阵,并进行运算得,,⎥⎦⎤⎢⎣⎡O =⎥⎦⎤⎢⎣⎡I -O I ⎥⎦⎤⎢⎣⎡O BC B AB C B ABC AB 可知)()(⎥⎦⎤⎢⎣⎡O ≥⎥⎦⎤⎢⎣⎡O BC B AB r B ABC AB r . 由引理2得),()()(B r ABC r B ABC AB r +=⎥⎦⎤⎢⎣⎡O ),()()(BC r AB r BC B AB r +=⎥⎦⎤⎢⎣⎡O 所以)()()()(BC r AB r B r ABC r +≥+.即).()()()(B r BC r AB r ABC r -+≥命题11 设k A A A ,,,21 均为n m ⨯矩阵,且1)()()(21====k A r A r A r 则k A A A r k ≤+++)(21 .证明:由命题5得.),()()(),()()(),()(),(21432132121k r r r r r r r r r k k k k =A A +A ≤A A +A +A +A ≤A A +A +A ≤A A +A设A 、B 均为n 阶方阵.则)()()(B r A r B A AB r +≤++.证明:构造如下矩阵并进行运算得:,⎥⎦⎤⎢⎣⎡O +O ++=⎥⎦⎤⎢⎣⎡O I O I +⎥⎦⎤⎢⎣⎡O A AB B A AB B A B A可知)()()()(B r A r A B A r A AB B A AB r +=⎥⎦⎤⎢⎣⎡O ≤⎥⎦⎤⎢⎣⎡O +O ++, 其中),()(B A AB r AAB B A AB r ++≥⎥⎦⎤⎢⎣⎡O +O ++ 所以)()()(B r A r B A AB r +≤++.设A 、C 均为n m ⨯矩阵, B 、D 均为s n ⨯矩阵,则)()()(D B r C A r CD AB r -+-≤-.证明:构造分块矩阵,并进行如下运算,⎥⎦⎤⎢⎣⎡-O --=⎥⎦⎤⎢⎣⎡I O I ⋅⎥⎦⎤⎢⎣⎡-O O -⋅⎥⎦⎤⎢⎣⎡I O I D B CD AB C A B D B C A C s n n m其中⎥⎦⎤⎢⎣⎡I O I n mC、⎥⎦⎤⎢⎣⎡I O I s n B 为可逆矩阵,所以 ),()()()(D B r C A r D B CD AB CA r DB CA r -+-=⎥⎦⎤⎢⎣⎡-O--=⎥⎦⎤⎢⎣⎡-O O -所以).()()()(D B r C A r D B CD AB CA r CD AB r -+-=⎥⎦⎤⎢⎣⎡-O--≤-即)()()(D B r C A r CD AB r -+-≤-.设A 是n 阶方阵,m ,k 为非负整数,则)()()1()(22A ⋅-A ⋅+≥+r m r m A r m .[8]证明:用数学归纳法,当0=m 时显然成立. 由命题10(Frobenius 不等式)).()()()(B r BC r AB r ABC r -+≥得:),()(2)()()()()(23A r A r A r A A r A A r A A A r A r -=-⋅+⋅≥⋅⋅=所以当1=m 时不等式成立. 假设当k m =时不等式成立,即:),()()1()(22A r k A k A r k ⋅-⋅+≥+于是)()1()()2()()()()()1()()()()()(222113A r k A r k A r A r A r k A r k A r A A r A A r A A A r A r k k k ⋅+-⋅+=-+⋅-⋅+≥-⋅+⋅≥⋅⋅=+++,所以当1+=k m 时不等式成立. 故).()()1()(22A r m A r m A r m ⋅-⋅+≥+3.15 命题15 设A 是n 阶方阵,且)()(2A r A r =则对任意自然数k ,有)()(A r A r k =.证明:构造分块矩阵⎥⎦⎤⎢⎣⎡OO 22A A由Frobenius 公式得),()()()()()()(33322222A r A r A A r A A A r A AA r A r A r +=⎥⎦⎤⎢⎣⎡O -O =⎥⎦⎤⎢⎣⎡O -=⎥⎦⎤⎢⎣⎡O ≤+ 由)()(2A r A r =,得),()()()()(2223A r A r A r A r A r =-+≥由定理2得),()()(223A r A A r A r ≤⋅=所以),()(23A r A r =以此类推).()()()(432k A r A r A r A r ==所以得)()(A r A r k =.设A 是非异阵,⎥⎦⎤⎢⎣⎡D C B A是n m ⨯阵,则 ).()()(1B CA D r A r D C B A r --+=⎥⎦⎤⎢⎣⎡[5] 证明:,11⎥⎦⎤⎢⎣⎡-O=⎥⎦⎤⎢⎣⎡⎥⎦⎤⎢⎣⎡I -O I ---B CA D BAD CB ACA r m r而A 是个非异阵,所以)()()()(11B CA D r A r B CA D BA r D CB Ar ---+=⎥⎦⎤⎢⎣⎡-O =⎥⎦⎤⎢⎣⎡. 即).()()(1B CA D r A r D C B A r --+=⎥⎦⎤⎢⎣⎡ 7命题17 设A 是n 阶方阵,且I =2A ,试证.))(())((n A r A r k m =I -+I +其中m 、k 为自然数.[6]证明:因为)(22)()(2I +=I ++=I +⨯I +A A A A A ,所以)()(2I +=I +A r A r .由命题15得)()(I +=I +A r A r m .同理)()(I -=I -A r A r k .0)()(22=I -=I -⨯I +A A A ,由命题8得n A r A r ≤I -+I +)()(.又n A A r A r A r =-I ++I ≥-I ++I )()()(.所以n A r A r =I -+I +)()(,即n A r A r k m =I -+I +))(())((.3.18 命题18 设A 是n 阶方阵,且A A =2,试证n A r A r k m =I -+))(()(,其中m 、k 为自然数.证明:因为A A =2,所以)()(2A r A r =,那么)()(A r A r m =. 因为A A A A A -I =I +-=I -⨯I -2)()(2,所以)())((2I -=I -A r A r .所以)())((I -=I -A r A r k .因为O =-=I -⨯A A A A 2)(,由命题8得n A r A r ≤I -+)()(.又因n r A A r A r A r =I =+-I ≥-I +)()()()(.所以n A r A r =I -+)()(,即n A r A r k m =I -+))(()(.注:以上证明中命题5、6、7、10、12、13、15、16都是采用了分块矩阵,其中包括分块矩阵的和、积以及分块矩阵的初等行列变换,不仅降低了处理矩阵相关问题的难度,还缩减了证明过程,使其过程简明概要,可读性强.分块矩阵对于处理多个矩阵之间的不等式,多个矩阵秩的范围的界定和秩的大小的比较有一定的优越性.命题11、14、15、17、18中都对矩阵的N 次幂进行了秩的运算或比较,用到了归纳、递推、叠代等运算方法,了解到高幂次矩阵的秩的大小或范围,对于处理高幂次矩阵问题,认识高幂次矩阵的性质都十分有用.4. 矩阵的秩的一些结论的应用4.1 例1已知矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---=6116100010A ,求矩阵的秩. 解:存在可逆矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---=P 941321111, 使得⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---=P P -3000200011A , 由命题3可知3)()(1=P P =-A r A r .4.2 例 2 存在矩阵A 、B 分别为n 阶方阵且n B r <)(,试证明0=Ax 的解向量是0=BAx 的解向量的一个子阵.证明:设α为0=A x 的解向量,必然存在0=αBA .由命题2{},)(),(m in )(B r A r B r ≤A得),()(A r BA r ≤则其解向量的秩),()(A r n BA r n -≥-所以0=Ax 的解向量是0=BAx 的解向量的一个子阵.4.3 例3 已知矩阵⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=4041050545450001A ,求其伴随矩阵*A 的秩. 解:对A 进行初等变换,求其秩.,00000500054000010041050005400001004105050545000140410505454500012143124⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡−−−→−⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡−−→−⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡−−→−⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡----r r r c c c c显见3)(=A r .由命题4可知1)(*=A r .4.4 例 4 已知A 是56⨯矩阵, B 为55⨯方矩,4)(=A r ,4)(=B r ,试证明方程组0=ABx 根的个数小于等于3.解:由命题7得,3)()()(=-+≥n B r A r AB r又由命题8,)()(n B r A r ≤+得6)()(≤+x r AB r ,所以3)(=x r ,则命题成立.4.5 例5 已知n 阶方阵A ,证明)(2)2(2A r A A r ≤+. 解:由命题12得),()()(B r A r B A AB r +≤++当A =B 时可写成)()()(A r A r A A AA r +≤++,即)()()(2A r A r A A A r +≤++,因此得).(2)2(2A r A A r ≤+4.6 例6 已知n 阶方阵A 和B ,秩为r B r A r =+)()(,求矩阵⎥⎦⎤⎢⎣⎡I ++n AB B A B A 的秩. 解:由分块矩阵的初等变换,2212⎥⎦⎤⎢⎣⎡I ++O +−−→−⎥⎦⎤⎢⎣⎡I O +−−→−⎥⎦⎤⎢⎣⎡I ++++-n r r n c c n BA AB BA ABBA AB AB B A B A 由命题12的结论)()()(B r A r B A AB r +≤++,再由引理2的结论得)()()(n n r B A r BA AB BA r I ++=⎥⎦⎤⎢⎣⎡I ++O +, 因为n r r B A r n +=I ++)()(,所以n r ABB A BA r n +=⎥⎦⎤⎢⎣⎡I ++)(.4.7 例7 存在矩阵⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=0106011010530202A ,求)(100A r . 解:计算得3)(=A r .⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=104018103116352095420124010601101053020201060110105302022A .得3)(2=A r ,所以)()(2A r A r =.由命题15的结论)()(A r A r k =,所以3)(100=A r .证明过程用到了命题15的结论,对于幂次为100的矩阵秩的求解,可以迅速的实现.4.8 例8 求⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡----=M 1785441610512112012311012的秩.解:由命题16的结论)()()(1B CA D r A r D CB Ar --+=⎥⎦⎤⎢⎣⎡可得)2011102312541011178416512()2312()(1⎥⎦⎤⎢⎣⎡-⋅⎥⎦⎤⎢⎣⎡⋅⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-+⎥⎦⎤⎢⎣⎡=M -r r r )201421128201(2)(⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--+=M r r 再次应用命题16[][][])201282014112()1(2)(1⋅⋅⎥⎦⎤⎢⎣⎡-⎥⎦⎤⎢⎣⎡-++=M -r r r 5212)(=++=M r .证明过程通过两次用到命题16的结论实现了对维数较高矩阵秩的求解,简化了运算,且大大减少了计算量.小结矩阵的秩的内容是非常丰富的,其应用是十分广泛的,证明矩阵秩的有关性质,除了利用分块矩阵以外,在上面还用到了行(列)向量组的极大线性无关组来证,以及矩阵的初等变换来证明,还可以联系到齐次线性方程组的基础解系来证.本文引用到矩阵秩的基本性质及部分定理,对矩阵秩的多条性质进行了证明,并做了相关的应用.其中涉及到了矩阵秩的求解、判断、向量组的相关性、方程组解的情况分析以及秩的不等式、等式等多方面性质.其结论和证明过程可以应用到所涉及的各个领域,包括电子行业,信息处理行业和控制工程域等多个行业.对矩阵秩的多方面的了解对于处理矩阵相关的问题是很有帮助的,例如方程组解的个数问题,模式识别中事物特征的相关性问题等.上述的命题的涉及面广,结论应用性强,所应用的方法较为新颖,希望能对数学及其它领域的发展有所帮助.相信在解决理论研究和解决实际问题上有一定的作用及意义.致谢参考文献[1]杜现昆原永久牛凤文.高等代数[M].北京:高等教育出版社,2006.65-67[2]同济大学数学系编.工程数学.线形代数[M].北京:高等教育出版社,2007.62-65[3]北京大学数学系.高等代数[M].北京:高等教育出版社,1988.87-90[4]张禾瑞.郝炳新.高等代数[M]. 北京:人民教育出版社,1979.76-77[5]张远达.线性代数原理[M]. 上海:上海教育出版社,1982.98[6]北京大学数学系.高等代数[M].北京:高等教育出版社,2003.67[7]赵树媛.线性代数学习与考试指导[M].北京:中国人民大学出版社,1998.56[8]丘维声.高等代数[M].北京:高等教育出版社,2003.87-88[9] 樊恽钱吉林等.代数学辞典[M].武汉:华中师范大学出版社,1994.76[10李书超,蒋君,向世斌等.一类矩阵秩的等式及其推广[J].武汉科技大学学报自然科学版 ,2004 ,27 1 :96-98[11]王松桂,贾忠贞.矩阵不等式[M].合肥:安徽教育出版社,1994.89-90[12]鲍文娣,李维国.关于任意三矩阵秩的一点注记[J].苏州科技学院学报:自然科学版,2005,22(2):39-43。

线性代数论文(矩阵在自己专业中的应用及举例)

线性代数论文(矩阵在自己专业中的应用及举例)

矩阵在自己专业中的应用及举例摘要:I、矩阵是线性代数的基本概念,它在线性代数与数学的许多分支中都有重要的应用,许多实际问题可以用矩阵表达并用相关的理论得到解决。

II、文中介绍了矩阵的概念、基本运算、可逆矩阵、矩阵的秩等内容。

III、矩阵在地理信息系统中也有许多的应用,比如文中重点体现的在计算机图形学中应用。

关键词:矩阵可逆矩阵图形学图形变换正文:第一部分引言在线性代数中,我们主要学习了关于行列式、矩阵、方程、向量等相关性比较强的内容,而这些内容在我们专业的其他一些学科中应用也是比较广泛的,是其它一些学科的很好的辅助学科之一。

因此,能够将我们所学的东西融会贯通是一件非常有意义的事,而且对我们的学习只会有更好的促进作用。

在计算机图形学中矩阵有一些最基本的应有,但是概念已经与线性代数中的有一些不同的意义。

在计算机图形学中,矩阵可以是一个新的额坐标系,也可以是对一些测量点的坐标变换,例如:平移、错切等等。

在后面的文章中,我通过查询一些相关的资料,对其中一些内容作了比较详细的介绍,希望对以后的学习能够有一定的指导作用。

在线性代数中,矩阵也占据着一定的重要地位,与行列式、方程、向量、二次型等内容有着密切的联系,在解决一些问题的思想上是相同的。

尤其他们在作为处理一些实际问题的工具上的时候。

图形变换是计算机图形学领域内的主要内容之一,为方便用户在图形交互式处理过程中度图形进行各种观察,需要对图形实施一系列的变换,计算机图形学主要有以下几种变换:几何变换、坐标变换和观察变换等。

这些变换有着不同的作用,却又紧密联系在一起。

第二部分 研究问题及成果1. 矩阵的概念定义:由n m ⨯个数排列成的m 行n 列的矩阵数表⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡ann an an n a a a n a a a 212222111211 称为一个n m ⨯矩阵,其中an 表示位于数表中第i 行第j 列的数,i=1,2,3,…n ,又称为矩阵的元素。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

矩阵相关文献翻译:Cooperative Spectrum Sensing Using Random Matrix TheoryLeonardo S. Cardoso and Merouane Debbah and Pascal BianchiFROM IEEE字数:5000字基于随机矩阵理论的协作频谱感知摘要本文提出了一种基于随机矩阵理论的协作频谱感知算法,这个算法既适用于AWGN,也适用于衰落信道。

不像先前的研究工作,新算法并不需要噪声统计和方差,并且与随机矩阵的最大和最小特征值有关。

值得注意的是,仿真结果表明,新算法方便随时间变化的拓扑结构,其性能明显优于典型的能量检测算法。

一、前言从美国联邦通信委员会(FCC)频谱政策专责小组[1]的报告中显示,无论是由于稀疏用户访问还是系统的固有缺陷,目前移动通信系统并没有充分利用可用的频谱,这已经成为共识。

可以预见,未来的系统将能够有机会利用这些频谱,通过认知环境的能力的相关知识,以适应相应的无线电参数[2]。

由于微电子和计算机系统的最新进展,这种无线电的时代已经不远,其中最重要的是开发出很好的感知技术。

用最通俗的话来说,频谱检测手段是在一个给定的有噪声的频段下寻找频带中的信号在(也可能包括进行分类的信号)。

这个问题以前得到广泛的研究,如今由于认知无线电研究的部分原因重获关注。

为此,有几个经典的技术,如能量检测(ED)(文献[3] - [5]),匹配滤波器(文献[6])和循环平稳特征检测(文献[7] - [9])。

这些技术有自身的优缺点,而且都是适合于非常特殊的应用场合。

然而,从认知无线电的角度来看,频谱感知有非常严格的要求和限制的问题,例如:•没有信号结构的先验知识(统计、噪音方差值,等等);•在最短的时间内的信号检测;需要具有在严重衰落信道的环境下可靠检测的能力。

Cabric等人的工作[7]、Akyildiz等人的工作[10]、和Haykin[11]提供了从认知网络的角度对这些经典技术进行了汇总。

从这些工作中可以清楚的看到,任何方法都不可能完全应付认知无线电网络的所有需求。

在简单的AWGN(加性高斯白噪声)信道中,经典的方法效果非常好。

然而,在快衰落的情况下,这些技术无法提供满意的解决方案,尤其是隐藏节点问题[12]。

为此,[13]- [16]几部文献已经研究认知无线电的协作频谱感知的情况。

这些工作的目的是通过增加额外的冗余感知方法降低错误概率。

他们还旨在通过减少收集的样本数量,来使用并行测量装置估计次数。

不过,即使人们可以高效的利用空间维度,这些工作也都是是基于相同的基本技术,都需要一个信号的先验信息。

在这项工作中,我们引入一个不需要先验信息的频谱感知方法。

这种方法依赖于多个接收器采用随机矩阵理论(RMT)对接收到的信号进行结构推断。

随机矩阵理论(RMT)是研究大维随机矩阵的经验谱分布函数在一定条件下特殊收敛性质的相关理论,现已被广泛应用于无线通信领域中,如无线信道容量、阵列信号处理、接收机性能分析、通信系统设计等的各个方面。

基于RMT 的频谱感知算法具有能够在噪声统计特性未知的情况下成功感知的优点。

实验表明,我们可以用所收到的少量可靠取样对频谱占用进行估算。

本文的其余部分作如下划分。

在第二节,我们提出了频谱感知的问题。

在第三节中,我们从随机矩阵理论的基础上提出方法。

在第四节中,我们拿出一些实际结果,证实在少量取样下的假设。

然后,在第五节中,我们对所提出的方法显示了性能测试结果。

最后,第六节中,我们得出主要结论,并指出进一步的研究。

二、问题提出频谱感知的基本问题是在有噪声情况下的信号检测。

这是一项艰巨的任务,特别在由于路径损耗或衰减造成接收到的信号功率非常低的情况下,频谱的盲检测是未知的。

这个问题可以构成假设检验,如[3]:其中y(k)是样品在k时刻接收的向量,n(k)是噪声方差(不一定是高斯的),h(k)是衰落分量,s(k)是我们要检测的信号,这样的话,E[s(k)2]=0,上式H0和H1分别只有噪声和信号假设。

我们还假设在N个信道上(K = 1.. N)的h保持不变。

基于能量频谱感知的经典技术检测是将信号能量与已知的阈值相比较,此阈值由V T[3] - [5]中的噪声和信道的统计所得。

下面的式子被公认为判定规则:其中,E[| y(k)|2]是信号能量,V T通常认为是噪声方差。

这种做法有一个缺点,就是噪音/通道的分配和V T都不是已知先验的。

在现实环境中的V T取决于无线的特点,这很难正确估计。

此外,在衰落和路径损耗的情况下,接收到的信号能量可能是噪音序列,使其很难检测,因而使检测到的取样数量N更加有限。

事实上,E[| Y(K)|2]由下式估计得到:这不是一个少量取样情况下的估计公式。

下面,我们在认知网络下提供了一种基于随机矩阵理论的协作频谱感知算法来从原始系统中检测信号而无需噪声方差。

.三、频谱感知的随机矩阵理论考虑图1中描绘的场景,其中用户(用白色表示)与他们的专用(主)基站通信。

二级基站{BS1,BS2,BS3,……,BS k}协作感知信道,以定位白色区域,并利用这些媒介。

在进一步讨论之前,让我们假设以下情况:• K 个基站的二级系统,它们之间的共享信息。

这可以通过有线高速骨干传输实现。

• 各个基站分析的频谱相同的部分。

让我们考虑以下的K×N 矩阵,此矩阵由所有的K 二级基站接收到的取样组成(y i (K )是基站i 在k 时刻的取样)。

随机矩阵理论方法的目标是检验测试不同基站接受到信号之间的独立性。

事实上,在信号的存在的情况下(以H1为例),所有收到的取样之间都是相关的,而当没有信号的情况下(以H0为例的情况下),无论衰落情况如何,取样之间都是不相关的。

因此,在这种情况下,对于固定值K ,当N→∞,样本协方差矩阵N 1YY H 收敛于σ2I 。

然而,在实际情况中,N 可以是与K 有着相同数量级的序列,因此不能直接推断N1YY H 中取样的独立性。

这可以使用随机矩阵理论工具[17]。

在以下渐近随机矩阵理论结果状态的情况下,Y 的各项都是独立的(统一具体的概率分布,与没有信号传输,即H0的情况对应):定理。

考虑K×N 矩阵W ,它的的各项是独立的复杂(或实际)的均值为零且随机变量方差为N 2并且时刻4序列为O (21N )。

当K ,N→∞且K/N→α,根据经验,分布WW H 几乎肯定会收敛到一个非随机的极限分布。

当有趣的是,当没有信号,样本协方差矩阵的特征值(如图2,由MP表示)是有限的,无论是否有噪音存在。

因此,Marchenko-Pastur定理,是在假设矩阵是“所有噪声”的假设下的一个理论预测。

也就是说,从这个特征值分布理论极限的偏差上可以显示出在无噪声的情况下的有关矩阵的信息。

在一信号是(H1)的情况下,矩阵Y可以改写为:其中s(i)和z k(i)=σn k(i)分别是在时刻i和基站k的独立的信号和噪声。

让我们用T表示矩阵:TT H有一个明显特征值λ1= |h i|2+ σ2,,其他所有的等于σ2,。

N1YY H的特征值与大数据量的人口模型样本协方差矩阵的特征值的研究有关[18]。

让我们在这项工作中定义信号信噪比(SNR)ρ :Baik等人最近的工作(文献[18][19])表明:当(这显然满足的假设,当取样数量N 足够高时),N 1YY H 的最大特征值几乎肯定收敛于:这在是优于b =σ2(1+ )2 H0情况下得到的。

因此,不论何时当矩阵N1YY H 特征值的分布不符合背离Marchenko-Pastur 定理(图3)时,探测器知道信号是否存在。

因此,可以使用这个有趣的特征感知频谱。

令λi 为N1YY H 特征值且G= [a ,b],协作感知算法的工作原理如下: A. 噪声分布未知,方差已知。

在这种情况下,用于以下准则:请注意,这种算法的改进(在非渐近情况下考虑出错概率)可以在文献[20]中找到。

这是在H0和H1的情况下,渐近最大的特征值分布的计算结果。

B. 两个噪声分布和方差未知。

请注意,在H0假设的情况下最大和最小特征值的比值不依赖于噪声方差。

因此,为了避免需要噪声的信息,用于以下准则:应该指出,在这种情况下,我们仍然需要较高的取样数量N ,使得条件 式(2)得到满足。

换句话说,取样数的平方与与信噪比的倒数成比例。

此外,请注意测试H1还提供一个很好的估计信噪比ρ。

事实上,N 1YY H 的最大的特征值(b')和最小(a )特征值的比仅与ρ和α有关:据我们所知,这SNR 估计从未在文献中提出。

四、性能分析前面的理论结果显示,我们可以只用取样的协方差矩阵的最大和最小特征值的比值就可以从噪声中区分出信号。

对于有限维数,这样的算法的操作区域仍然是一个问题,这与缩放因子的比例的渐近分布[20]有关。

这节通过分析各种大小不同的矩阵的之间N1YY H 的λmax 和λmin 的比值,说明了在这种情况下的一些特性。

图4和图5显示了在α=1 / 2和α= 1 / 10下,各种矩阵在纯噪声的情况下Y 的λmax/λmin 值。

从这些数字中,我们看到,即使矩阵较小,这两种情况下也可以提供一个很好的近似渐近比例。

如果有需要,例如,N = 100(K= 50,α= 1 / 2;K = 10,α= 1 /10),可以看出,在α=1 / 2和α= 1 / 10的情况下,模拟结果分别等于81%和83%的渐近。

正如预期的那样,对于一个较大的Y 矩阵,图中的比值接近相同的渐近值。

图6和图7显示了α=1 / 2和α= 1 / 10且信号加噪声的情况下,λmax/λmin 的取值。

在这两种情况下,σ2=1/ρ(-5分贝ρ)且 |h i|2=1(运行在表达式(2)的准则下)。

在纯信号的情况下,λmax/λmin=b'/a。

有趣的是,当N = 100(K= 50,α= 1 / 2;K = 10,α= 1 /10),可以看出,在α=1 / 2和α= 1 / 10的情况下,模拟结果分别等于70%和83%的渐近。

正如预期的那样,Y矩阵越大,图中的比值越接近渐近值。

N的一个很好的近似值可低至100。

五、结果基于文献[15][16],我们将随机矩阵探测器的方案和能量协作方案进行了模拟仿真,并对性能进行了对比。

能量探测器的框架在本文的第二部分提到过,用h(k) 瑞利多径衰落方差1/ K进行建模。

造成的差异被归一化,因为考虑到由于路径损耗的缘故,随着基站数量的增加,能量不会没有边界的无限增大。

共对十个二级基站进行了模拟。

基于投票方案,判定准则为以下几点:一是认为整体频谱占用情况判定需要由绝大多数的二级基站选择。

阈值V T为σ2(在噪声方差已知的情况下)。

对于随机矩阵理论为基础的方案,一个固定的总(K = 10)基站被采纳。

相关文档
最新文档