毕业论文-矩阵的特征值与特征向量的相关研究【范本模板】

本科毕业设计（论文）
( 2015届 )
题目：矩阵的特征值与特征向量的相关研究
学院: 数理与信息工程学院
专业：数学与应用数学
学生姓名：学号:
指导教师：职称：
合作导师: 职称：
完成时间: 201 年月日
成绩:
浙江师范大学本科毕业设计（论文)正文
目录
摘要 (1)
英文摘要 (1)
1 引言 (1)
2 选题背景以及特征值与特征向量的定义与性质 (2)
2. 1 选题背景 (2)
2。

2 特征值与特征向量的定义 (2)
2. 3 特征值与特征向量的性质 (2)
3 矩阵的特征值与特征向量的求解方法 (3)
3。

1 求解数字方阵的特征值与特征向量 (3)
3。

2 已知矩阵A的特征值与特征向量，求与A相关的矩阵的特征值 (7)
4 矩阵的特征值与特征向量的反问题的求解 (7)
4。

1 矩阵的全部特征值与全部特征向量，反求解矩阵A的方法 (7)
4。

2 已知实对称矩阵的全部特征值和部分线性无关的特征向量，反求矩阵A的方法 (9)
5 矩阵的特征值与特征向量的应用 (9)
5. 1 矩阵的特征值与特征向量在线性递推关系上的应用 (9)
5。

2 经济发展和环境污染的增长模型 (14)
6 结论 (16)
参考文献 (16)
矩阵的特征值与特征向量的相关研究
摘要：矩阵的特征值与特征向量占据了高等数学中的一小块, 但是其重要性无可比拟, 它可以应用在数学和生活上，尤其是对现在的科学技术领域, 有着至关重要的作用。

本篇论文主要阐述并归纳了矩阵的特征值与特征向量的概念，性质, 解法以及应用，通过具体的例子, 来体现了矩阵的特征值与特征向量的广泛性和实用性，深刻研究了矩阵的特征值
与特
征向量和它相关的应用。

正文总共分为四个大部分。

第一部分：阐述了它的概念和性质; 第二部分：对于它的求解方法，本篇论文叙述了几种不同的方法，并且有相关例题的作法；第三部分: 关于它的反问题, 本篇论文也有相对应的几种不同的求解方法; 第四部分：关于它在数学领域和生活上的应用.
矩阵; 特征值；特征向量; 反问题；应用关键词:
Correlation matrix eigenvalues and eigenvecto -rs
Mathematical and Information Engineering Mathematics and Applied Mathematics Chen Dong(11170126)
Instructor: Lvjia Feng （Associate Professor）
Abstract: Eigenvalues and eigenvectors occupy the higher mathematics in a small, but its importance is unparalleled，it can be used in mathematics and life，especially in the field of science and technology right now，has a vital role. This paper describes and summarizes the main characteristics and eigenvector matrix concept，nature，solution and applications, through specific examples, to reflect the breadth and practicality matrix eigenvalues and eigenvectors，profound study of matrix eigenvalues and special Eigenvectors and its related applications.
Total body is divided into four parts。

The first part：it describes the concept and nature; Part II：For its solution method, this paper describes several different methods，and relevant examples of practice; Part III：Anti question about it, this papers are also several different corresponding method for solving; part IV: on its application in the field of mathematics and life. Key Words: Matrix; eigenvalues; feature vector; inverse problem；Application
1 引言
在已经有相关深刻探讨的前提下, 本篇论文给出了它的的概念以及它的性质，掌握它的性质是研究其求解方法的前提, 所以要先熟悉它的性质, 再对它的求解方法作详细的步骤和说明. 本篇论文重点介绍了它的求解方法和特它的反问题以及相关应用，展现了它在矩阵运算中的重大作用, 在例题的求解过程中充分运用某些性质，使得问题变得简单, 运算方面上也更简洁, 是简化一些有关矩
阵的比较繁琐问题的一种快捷并且有效的途径. 本篇论文通过一些具体的例题详细说明它的求解方法以及其反问题的求解方法， 并且在数学领域以及生活方面的应用也有其相关的例题来说明矩阵的特征值与特征向量的广泛性以及实用性。

2 特征值与特征向量的选题背景以及其定义与性质
2. 1 选题背景
随着科技的迅猛发展 ,现在的社会发展的速度日益增加， 高等代数作为一门大学数学的基础学科已经向所有的领域渗透 , 它在所有领域内表现出来的作用已经越来越明显. .
物理、化学、经济等的许多问题在数学上都可以看作是求它的问题. 但是通过特征方程求解它是有一点难度的, 而且在现在的高等数学的教材中用特征方程求它总是要求解带含有参数的行列式， 而且只有先求解出它才能用方程组求解之后的问题。

本篇论文将对它的求解方法、反问题以及相关的应用进行系统性的归纳， 并且有相关的例题给予帮助理解. 2. 2 特征值与特征向量的定义
它在《高等代数》和《线性代数》课程中占据了一席之地， 在大多数的《高等代数》教材中, 把它拉进来就是为了解析线性空间中线性变换A /的， 它的定义如下：
定义 1 设A /是数域P 上的线性空间V 的一个线性变换, 如果对于数域P 中的数λ, 存在一个不是零的向量V ∈ξ, 使得
λξ=A /
那么λ是矩阵A 的一个特征值, 向量x 称作矩阵A 关于特征值λ的特征向量. 在大多数的《线性代数》的教材中， 它的探讨作为矩阵探讨的一个至关重要的组成部分, 它的定义如下所述：
定义2 设 A 是n 阶的方阵, 如果存在数字λ和n 维不是零的向量x ， 使得
x Ax λ=
那么就称λ是A 的特征值， x 是A 的对应特征值λ的特征向量． 2. 3 特征值与特征向量的性质
(1）如果i λ是A 的i r 重的特征值, A 所对应的特征值i λ就会有i S 个线性无关的特征向量．
（2）如果21,x x 都是矩阵A 的属于特征值0λ的特征向量, 那么当21,k k 不全都是零时, 2211x k x k +依然是A 的属于特征值0λ的特征向量．
(3）如果n λλλ,...,,21是矩阵A 的互相不一样的特征值， 而且它所对应的特征向量分别是n x x x ,...,,21， 那么n x x x ,...,,21线性无关．
(4）如果
()n
n ij a A ⨯=的特征值是
n
λλλ,...,,21, 那么
nn n a a a +++=++......221121λλλ， A n =λλλ...21.
（5）实对称矩阵A 的特征值都是实数, 属于不同的特征值的特征向量正交． （6）如果i λ是实对称矩阵A 的i r 重的特征值, 那么所对应特征值i λ刚好有i r 个线性无关的特征向量．
（7）假设λ是矩阵A 的特征值, )(x P 是多项式的函数, 那么)(λP 是矩阵多项式)(A P 的特征值．
3。

矩阵的特征值与特征向量的求解方法
3。

1 求解数字方阵的特征值与特征向量 （1）求解特征多项式()A E f A -=λλ.
（2）特征方程0=-A E λ， 它的全部根n λλλ,...,,21就是A 的全部的特征值． （3）对于任何一个特征值()n i i ≤≤1λ， 求解出齐次的方程组()0=-x A E i λ的一个基础解系ir i i a a a ,...,,21就是A 的属于()n i i ≤≤1λ的线性无关的特征向量。

那么A 的属于i λ的全部的特征向量是ii i i i a k a k a k +++...2211， 其中i k k k ,...,,21是不全都是零的数.
求解特征多项式是解决问题的难度所在， 方法一: 观察特征矩阵的每一行之和, 如果相等而且都是a, 那么将第2列及以后各列都加到第1列， 提取公因子， 再作化简, 而且a 就是其中的一个特征值， ()T
1,...,1,1是A 的属于特征值a
的特征向量．方法二： 将特征矩阵的的两个不是零的常数（不含参数λ）之一化为零, 如果有公因子, 提取出来再作化简.
从上述可以知道， 求解它是相当繁琐的．这里将阐述一个有效的方法， 只是需要对原来的矩阵作行列互逆变换就可以同时求解出它， 所以给出如下定义：
定义: 称矩阵的下列三种变换为行列的互逆变换:
（1）互相更换矩阵的i ，j 两列， 同时互相更换矩阵的i ，j 两行;
（2）矩阵的第i 行乘以不是零的数字k, 同时矩阵的第i 列乘以k
1
；
(3)矩阵的第i 行乘以k 倍加到矩阵的第j 行, 同时第j 列乘以 −k 倍加到矩阵的第i 列.
定理： A 为n 阶的可以对角化的矩阵, 而且()()
T DP −−−−−−→−一系列行列互逆变换n T E A , 其中，
⎥⎥
⎥⎦
⎤⎢⎢⎢⎣⎡=n D λλ 1 ⎪⎪⎪⎭⎫
⎝⎛=n T P ββ 1 ()in i i b b 1=β ()n i ,...,1=,
那么n λλλ,...,,21是A 的全部特征值， T
i i βα=是A 的属于i λ的特征向量. 证明: 因为
D P A P T T T =-1)(
即 D D AP P T ==-1 从而AP=PD 因为
⎥⎥
⎥⎦⎤
⎢⎢⎢⎣⎡=n D λλ 1[]n P αα 1= 所以
()()⎥⎥
⎥⎦
⎤
⎢⎢⎢⎣⎡=n n n A λλαααα 111则 [][]n n n A αλαλαα 111= 所以
n i A i i i ,,1),0( =≠=αλαα
为了运算的简洁， 约定:
（1）i j ka a +表示为矩阵的第i 行乘以k 倍加到第j 行。

（2）i j ka a -表示为矩阵的第i 列乘以-k 倍加到第j 列.
因为用定理求解题目时， 总是会遇到一些类似
⎥⎦⎤⎢⎣⎡=b c a B 0或者⎥⎦
⎤⎢⎣⎡=b c a C 0（b a ≠)形式的矩阵的化对角阵的问题， 所以给出对应的求解方法：
()
⎥
⎦
⎤
⎢⎣⎡−−−−−−−→−⎥⎦⎤⎢⎣⎡=+1001010011221kr -r kr r 2b k a b a c a E B
T
，第二行第一行 或(
)
⎥⎦
⎤⎢⎣⎡-−−−−−−−→−⎥⎦⎤⎢⎣⎡=+-10010100102112kr r kr r 2k b a b c a E C
T
，第二行第一行 其中， )
(b a c k -=
， 所以()T k ,11=α， ()T
1,02=α是B 的分别属于特征值a 和b 的特征向量.
()T 0,11=β， ()T k 1,2-=β是C 的分别属于特征值a 和b 的特征向量。

下面将有3道例题来说明其求解方法, 第一道例题不使用刚才描述的方法，
则后面两道例题运用, 以此来说明这个方法的可操作性以及简便性。

例1： 求解矩阵⎪⎪⎪
⎭⎫ ⎝⎛-----=436102111A 的特征值与特征向量。

解:
()⎪⎪⎪⎪⎪⎪
⎪⎪
⎭
⎫ ⎝⎛---+=⎪⎪⎭
⎫ ⎝⎛10001000143612
11
13λλλλI A []⎪⎪⎪⎪⎪⎪⎪⎪
⎭
⎫ ⎝⎛---+−→
−001010100
63421
111
3,1λλλ []
()[]()⎪⎪⎪⎪
⎪⎪⎪⎪
⎭
⎫
⎝⎛+---+-+--+--−−−→−+--1110101002314111
001211312λλλλλλλλ []()⎪⎪⎪⎪
⎪⎪⎪
⎪
⎭
⎫
⎝⎛--------−−→−+211110100114011001223λλλλλ 所以, 矩阵A 的特征值是1321===λλλ 当1=λ时
()()⎪⎪⎪⎪⎪⎪⎪
⎪⎭
⎫ ⎝⎛---=⎪⎪⎭
⎫ ⎝⎛31111010000
300100111P G
于是, 可以知道属于特征值1=λ的特征向量是()()T
T 3,1,1,1,1,021-=-=ξξ.
例2： 求解 ⎥⎥⎥⎥
⎦
⎤⎢⎢⎢
⎢⎣⎡----=0111101111011110B 的特征值与特征向量. 解
⎥⎥⎥⎥⎦
⎤⎢
⎢⎢
⎢⎣⎡----=10000
111010010110010110100011110)(4E B T
⎥⎥⎥⎥⎦
⎤⎢
⎢⎢
⎢⎣⎡------−−→−++--1100
1020010011100011
2010000110114
3213412r r r r r r r r ⎥
⎥
⎥⎥⎦⎤
⎢
⎢⎢
⎢⎣⎡-----−−→−+-11001020010001101111003000
1
0011244
2r r r r ⎥
⎥
⎥⎥⎦⎤
⎢
⎢⎢
⎢⎣⎡-----−−→−-+-1100102
300100014301111003
000
1
04314
23
21
22
14141r r r r r r ⎥
⎥
⎥⎥⎦⎤
⎢
⎢⎢
⎢⎣⎡-------−−→−+-+21212121100
04143414101001111003
0414
14
14
30
0012
42
31
22
14141r r r r r r ⎥
⎥
⎥⎥⎦
⎤
⎢⎢⎢
⎢⎣⎡-------−→−11113000111101001311001011
130001 所以特征值分别是
;3,14321-====λλλλ
特征向量分别是
()T 1,1,1,31-=α, ()T 1,3,1,12-=α， ()T 1,1,1,13--=α, ()T 1,1,1,14--=α。

下面给出上述定理的推广定理： 定理: A 是任意n 阶的矩阵, 如果
()(
)
T n T
P J
E A
−−−−−−→−一系列行列互逆变换， 其中)(1n r J J J r ≤⎥⎥
⎥⎦
⎤
⎢⎢⎢⎣⎡= 是约当矩阵， ),,1(11r i J i i i =⎥⎥
⎥⎥⎦⎤
⎢⎢⎢
⎢⎣⎡=λλ是约当标准形, ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=i ir i i r T P P P P ββ 11, ),,1(r i =; n r r r r =++ 21所以i λ是A 的特征值， T
ir i i βα=是A 的特征值的特
征向量。

例3： 求解⎥⎥
⎥⎦⎤
⎢⎢⎢⎣⎡--=312130112B 的特征值与特征向量。

解()
⎥⎥
⎥⎦
⎤
⎢⎢⎢⎣⎡--=1003110101310012023E A T
⎥⎥
⎥⎦
⎤
⎢⎢⎢⎣⎡---−−→−+-1004110100311010111
33
1r r r r ⎥⎥
⎥⎦
⎤
⎢⎢⎢⎣⎡---−−−−→−111400111020101012再作一系列变换 所以特征值是
,4,2321====λλλ221===λλ的特征向量()T 1111-=α,43=λ的特征向量
()T 1113-=α.
3. 2 已知矩阵A 的特征值与特征向量， 求与A 相关的矩阵的特征值 此种题目可以运用性质7来求解计算, 用定义就可以求解算得.
4 矩阵的特征值与特征向量的反问题的求解
4. 1 矩阵的全部特征值与全部特征向量， 反过来求解矩阵A 的方法 方法一: 用对角化法求解可逆矩阵P ， 使得B AP P =-1， 那么1-=PBP A 。

方法二: 用对角化法求解正交的矩阵)(1T T T T =-使得B AT T =-1, 所以
T TBT TBT A ==-1。

方法三： 特定元素法
设n 阶矩阵()n n ij a A ⨯=的全部特征值是n λλλ,,,21 , 相应的n 个线性无关的特征向量是n a a a ,,,21 , 所以有
n n n a Aa a Aa a Aa λλλ===,,,222111 ，
从这里可以得到以A 的第1行， 第2行，. . . ,第n 行的元素in i i a a a ,,,21 ),,2,1(n i =当作未知数的n 个非齐次的线性方程组, 求解每个方程
组求出A 中的元素ij a ， 那么就能得到()ij a A =。

例4： 设三阶方阵A 的特征值是,1,0,1321-===λλλ对应的特征向量分别是
()()()T
T T x x x 2,1,2,1,2,2,2,2,1321--=-==， 求解A 。

解： 因为)3,2,1(=i x i 是矩阵A 对应于特征值)3,2,1(=i i λ的特征向量, 所以有
i i x Ax λ=， 令
()⎪⎪⎪
⎭
⎫ ⎝⎛--==212122221,,321x x x P ， 那么
⎪⎪⎪⎭⎫ ⎝⎛---=-212122221
911
P ，
所以有PB AP =, 其中⎪⎪⎪
⎭
⎫ ⎝⎛-=100000001B , 就从上述式子可以得到
⎪⎪⎪⎭
⎫
⎝⎛-==-022********
1PBP A ，
就是问题所要求得的答案.
例5: 设三阶的实对称矩阵A 的特征值是6、3、3, 与特征值6对应的特征向量是()T
p 1,1,11=， 求解A 。

解: 设对应于3的特征向量是()T
x x x X 321,,=. 因为实对称矩阵的不同特征值下
的特征向量正交, 也就是X 的分量满足0321=++x x x ， 又因为特征值3的重数是2， 所以对应于3刚好有2个线性无关的特征向量， 明显0321=++x x x 的基础解系就是对应于3的2个线性无关的特征向量。

从0321=++x x x 得到它的一个基础解系是()T
p 0,1,12-=, ()T
p 1,0,13-=， 令
()⎪⎪⎪
⎭
⎫ ⎝⎛--==101011111,,321p p p P ,
所以可以得到
⎪⎪⎪
⎭
⎫ ⎝⎛==-3000300061B AP P ，
所以， ⎪⎪⎪⎭⎫
⎝⎛==-4111411141
PBP A ，
就是问题所要求得的答案。

4。

2 已经知道实对称矩阵的全部特征值和部分线性无关的特征向量， 反过来求解矩阵A 的方法
从实对称矩阵属于不同的矩阵的特征值的特征向量正交求解出其余的特征向量， 可以运用上述各种的方法求解。

5 矩阵的特征值与特征向量的应用
5. 1 矩阵的特征值与特征向量在线性递推关系上的应用
求解常系数齐次递推关系的方法多种多样， 这里将说明一下如何利用它来求解线性齐次递推关系的一种方法.
设k 阶线性循环数列{}n x 满足递推关系:
),2,1(,2211 ++=+++=---k k n x c x c x c x k n k n n n
其中),,2,1(k i c i =是常数， 并且0≠k c 。

方程组
⎪⎪⎪⎩⎪
⎪⎪⎨⎧===+++=+-+--------1
12
2112211k n k n n n n n k n k n n n x x x x x x x c x c x c x 可以表示为矩阵形式:
⎥
⎥
⎥⎥
⎦⎤
⎢⎢⎢
⎢⎣⎡⎥⎥⎥⎥⎥⎥⎦⎤
⎢
⎢⎢⎢⎢
⎢⎣⎡=⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡----+---k n n n k k k n n n n x x x c c c c x x x x 2112
1121010000100001 （1) ⎥
⎥⎥⎥⎦⎤
⎢⎢⎢⎢⎣⎡=⎥⎥⎥
⎥⎥⎥⎦
⎤
⎢⎢⎢⎢⎢
⎢⎣⎡=⎥⎥⎥⎥⎥⎥⎦
⎤⎢⎢⎢⎢⎢⎢⎣⎡=-----+---+-k n n n k n k k k n n n n k n x x x c c c c A x x x x x 211211211,0100001
00001,α
那么(1）可以写作：
k n k n A -+-=αα1 （2)
由（2）式子递推可以得到1121αααk n k n k n A A ---+-=== . 其中 []T
k k x x x x 1211,,, -=α
所以求解通项n x 就可以归结为求解1+-k n α， 也就是求解k n A -.
如果A 可以对角化, 那么存在可逆矩阵P, 使得A AP P =-1, 所以
1---=P PA A k n k n ， 因为
λλλλ10000
1
01121-------=
--
k k c c c c A E
从第一列开始每一列乘以λ加到后一列上， 就得到如下的矩阵： k k k k k k k k k k k c c c c c c c c c c c ----=⎥⎥
⎥
⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡------------------λλλλλλλλλλλ1111111
2112121010000100001 如果λ是A 的特征值, 明显有()1-=-k A E R λ， 所以线性齐次方程组()0=-X A E λ的基础解系中仅含有一个解向量， 因此当A 有k 个特征值
321,,,λλλ 时, 这k 个特征值对应的特征向量分别是k P P P ,,,21 , 由这k 个特征
向量为列构成的方阵记作P, 那么P 是可逆的， 并且A AP P =-1。

其中
⎥⎥
⎥⎥⎦
⎤⎢
⎢⎢
⎢⎣⎡=n A λλ
λ
000002
1 例6 设数列{}n x 满足递推关系:
)4(22321≥-+=---n x x x x n n n n ,并且3,2,1321=-==x x x , 求解通项n x . 解: {}n x 是三阶循环数列， 将方程组
⎪⎩⎪
⎨⎧==-+=-------22
1
132122n n n n n n n n x
x x x x x x x 用矩阵表示:
⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-----32121010001212n n n n n n x x x x x x 令⎥⎥⎥⎦
⎤
⎢⎢⎢⎣⎡-=010001212A 那么由上式可以递推得到：
⎥⎥⎥⎦
⎤
⎢⎢⎢⎣⎡==⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---------1233432232121x x x A x x x A x x x A x x x n n n n n n n n n n （1） 其中3,2,1321=-==x x x 因为0=-A E λ, 即
022*********=-+-=----λλλλ
λλ, 得到A 的特征值： 2,1,1321=-==λλλ
再从特征方程()()3,2,10==-i X A E i λ解得对应A 的特征值321,,λλλ的特征向量分别是:
⎥⎥⎥⎦⎤
⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=124111,111321P P P []⎥⎥⎥⎦
⎤
⎢⎢⎢⎣⎡-==111211411,,321P P P P 所以
11200010001,20223163361--⎥⎥⎥⎦
⎤
⎢⎢⎢⎣⎡-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---=P P A P 1
3
3
200010001---⎥⎥
⎥⎦
⎤
⎢⎢⎢⎣⎡-=P P A n n
()()()()()()()()()⎥⎥⎥⎦⎤
⎢⎢⎢⎣⎡--+--+-+---+--+-+---+--+-+-=-------------2312312
21233321261332132126133213212613321361n n n n n n n n n n n n n n n 代入（1）式子可以得到：
()()()[]
13
2333)2126()13()213(6
1x x x x n n n n n n --++--++-+-=---
()()⎥⎦
⎤⎢⎣⎡+-+-+-+-=
--+-131
32321611232111961n n n n .
例7 数列()()()()()()()296,185,114,73,42,31,10=======F F F F F F F ,
()477=F , 求解这个数列的通项()n F 。

解: 通过分析这个数列满足条件
()()()),2,1,0(12 =++=+n n F n F n F (1）
根据
()()()
()()
,2,1,01112=⎩⎨
⎧+=+++=+n n F n F n F n F n F 即
),2,1,0(1 ==+n Aa a n n (2）
其中
()()()()()()⎪⎪⎭
⎫
⎝⎛=⎪⎪⎭⎫ ⎝⎛=⎪⎪⎭⎫ ⎝⎛+=⎪⎪⎭⎫ ⎝⎛=⎪⎪⎭⎫ ⎝⎛++=+1301,1,0111,1201F F a n F n F a A n F n F a n n
从（2）式子递推可以得到:
),2,1,0(0 ==n a A a n n （3)
因为
011112
=--=---=
-λλλ
λλA E
得到A 的特征值是
2
5
1,25121-=+=
λλ (4） 对应于21,λλ的特征向量分别是
⎪⎪⎭
⎫ ⎝⎛=⎪⎪⎭⎫ ⎝⎛=1,12211λλX X 取⎪⎪⎭⎫ ⎝⎛=1121λλP ， 所以⎪⎪⎭
⎫ ⎝⎛---=-12211
111λλλλP 那么
⎪⎪⎭
⎫ ⎝⎛-----=⎪⎪⎭⎫ ⎝⎛=++++-n n n n n n n n n n
n
P P A 1221211
1
2121121121121
10
0λλλλλλλλλλλλλλλλ 所以有
⎪⎪⎭
⎫ ⎝⎛-+--+--===⎪⎪⎭⎫ ⎝⎛+++++n n n n n n n n n n a A a n F n F 1221211
1
2121121121033331)()1(λλλλλλλλλλλλλλ 于是
()()
n
n n n n F 1221212
1331λλλλλλλλ-+--=
（5)
把(4）式子代入到(5）式子得到
()1
1
251251++⎪
⎪⎭
⎫
⎝⎛-+⎪
⎪⎭
⎫
⎝⎛+=n n n F
就是题目所要求解的通项.
例8 计算
1
10000110000001110000111000011 =
n D
解: 按照矩阵的第一行展开()321≥==--n D D D n n n （1）
0,121==D D
把（1）变成() ,3,2,112=-=++n D D D n n n 因为
⎩⎨
⎧=-=++++1
112n n n
n n D D D D D 即
() ,3,2,11==+n Aa a n n （2）
其中⎪⎪⎭⎫ ⎝⎛⎪⎪⎭⎫
⎝⎛=⎪⎪⎭⎫ ⎝⎛=⎪⎪⎭⎫
⎝⎛-=⎪⎪⎭⎫ ⎝⎛=++++10,,0111,1211121D D a D D a A D D a n n n n n n 从（2)这个式子递推可以得到
() ,3,2,111==-n a A a n n （3)
因为
012=--=-λλλA E
得到A 的特征值是
2
31,23121i
i -=+=
λλ (4） 对应于21,λλ的特征向量分别是
⎪⎪⎭
⎫
⎝⎛=⎪⎪⎭⎫ ⎝⎛=1,12211λλX X
取⎪⎪⎭
⎫ ⎝⎛=1121λλP , 那么⎪⎪⎭⎫ ⎝⎛---=-12211
111λλλλP 那么
⎪⎪⎭⎫ ⎝⎛-----=⎪⎪⎭⎫ ⎝
⎛=--------11212112111
2212121112111
100n n n n n
n n n n n n P P A
λλλλλλλλλλλλλλλλ 所以就有
⎪⎪⎭⎫ ⎝⎛---===⎪⎪⎭
⎫ ⎝⎛+---1121211
22121111)()1(n n n
n n n a A a n D n D λλλλλλλλλλ 于是
()()
1
2
112
122212
2212111212121111++----=•-=--=n n n n n D λλλλλλλλλλλλλλλλ 5. 2 经济发展和环境污染的增长模型
为了研究某地区的经济发展与环境污染之间的关系, 可以建立如下数学模型：
设00,y x 分别是这个地区目前的环境污染水平和经济发展水平， 11,y x 分别是这个地区若干年后的水平, 而且有下述的关系：
⎩⎨
⎧+=+=0010
01223y x y y x x 令
⎪⎪⎭⎫
⎝⎛=⎪⎪⎭⎫ ⎝⎛=⎪
⎪⎭
⎫ ⎝⎛=2213,,11
1000A y x y x αα 所以上面描述的关系的矩阵形式是01ααA =。

那么经济发展与环境污染的增长模式是
),,2,1(223111
1k i y x y y x x i i i
i i i =⎩⎨
⎧+=+=---- 令⎪⎪⎭
⎫ ⎝⎛=i i i y x α
所以上面描述关系的矩阵形式是k t A i i ,,2,1,1 ==-αα 所以从上述这个形式可以得到：
10
32302120
1αααααααααααi i i A A A A A A A ========-
(*）
下面我们将进行更深一步的讨论： 从矩阵A 的多项式
()()142
21
3
--=----=
-λλλλλA E 得到A 的特征值是1,421==λλ
对于41=λ, 可以求解方程()04=-X A E 得到特征向量⎪⎪⎭
⎫
⎝⎛=111η
对于12=λ, 可以求解方程()0=-X A E 得到特征向量⎪⎪⎭
⎫
⎝⎛-=212η
明显， 21,ηη线性无关 下面分作三种情况分解析:
假设1： ⎪⎪⎭⎫
⎝⎛==1110ηα
从(*）以及它的性质可以知道
⎪⎪⎭
⎫
⎝⎛====1141
110i
i i
i
i A A ηληαα
即⎪⎪⎭⎫ ⎝⎛=⎪
⎪⎭
⎫ ⎝⎛114i i i y x 或者i
i i y x 4== 上面描述的式子表示： 在当前的环境污染水平和经济发展水平的条件下下， i 年后, 当经济发展水平达到相当高的程度时， 环境污染也保持着同步恶化趋势.
假设2： ⎪⎪⎭
⎫
⎝⎛-==2120ηα
因为020<-=y ， 所以不讨论这种情况
假设3： ⎪⎪⎭
⎫
⎝⎛=710α
因为0α不是特征值, 所以不能类似分析， 但是0α可以由21,ηη唯一线性表达出来
21023ηηα-=
由（*）以及特征值与特征向量的性质可以得到:
⎪
⎪⎭
⎫ ⎝⎛+⋅-⋅=⎪⎪⎭⎫ ⎝⎛-⋅-⎪⎪⎭⎫ ⎝⎛⋅=-=-=-==443243211211432323)23(221
12
1210i i i i
i
i
i i i i i A A A A ηληληηηηαα 也就是
⎪⎪⎭⎫
⎝⎛+⋅-⋅=⎪⎪⎭⎫ ⎝⎛443243i i i i y x
443,243+⋅=-⋅=i i i i y x
从上面描述的式子可以预测到这个地区i 年后的水平。

2η因为没有实际的意义所以在假设2中没有作相关的讨论， 但是在假设3中的讨论中起到了至关重要的作用。

6 结论
通过它的概念以及相关的性质学习, 理解了它的各种求解方法, 更是有相关例题求解巩固知识. 而后, 又学习了它的反问题， 也理解了相应的求解方法. 它不仅能应用在数学上, 帮助其更简单的运算, 而且也能应用在生活上, 有效处理生活的各类问题。

学习并且研究数学， 从知识联系到生活， 通过数学的思维或方法来处理某些生活上的问题。

离开数学, 科技无法进步, 生活恐难维持, 所以我们必须热爱数学， 深入探讨数学, 从而促进科技发展, 共创美好的明天。

参考文献
[1] 施劲松，刘剑平。

矩阵特征值、特征向量的确定［J]. 大学数学。

2003. 12,19（6）： 123—126.
［2］ 李迪。

中国数学史简编[M ］. 沈阳： 辽宁人民出版社。

1984,9: 124-125。

［3］ 刘国琪。

。

利用矩阵的初等行变换达到矩阵的特征值与特征向量的同步求解[J]. 江西电力职工大学学报. 1995. 3，8（1)： 30-34。

[4] 陈泽安。

求矩阵特征值与特征向量的新方法[J ］. 长沙通信职业技术学院学报. 2003。

3,2（1）: 66—69.
［5］ 赵文玲,王秀芬。

特征值与特征向量在线性递推关系中的应用［J ］. 山东工程学院学报。

1998. 9，12（3）: 39—44。

[6] 邵丽丽。

矩阵的特征值和特征向量的应用研究［J]。

菏泽学院学报。

2006。

2,28（5): 20—23。

［7］ Roger A. Horn,Charles R 。

Johnson. Matrix Analysis ［M]. POSTS 2005。

9: 1-561. ［8］ 郭华,刘小明. 特征值与特征向量在矩阵运算中的作用[J ］. 渝州大学学报。

2000. 6，17(2） ： 72—75.
［9] CHI BIN ，YE Qing —Kai. COMPUTING THE EIGENVECTORS OF A MATRIX WITH
MULTI
-PLEX EIGENV ALUES BY SVD METHOD[J］。

Published by Shanghai Univsity. 2004. 3，25（3）：10—15 .
[10] 同济大学数学系。

线性代数［M］。

上海:: 上海科学技术出版社,1999：74—76。

［11] 陈公宁. 矩阵理论及其应用［M］。

北京: 科学出版社. 2007,4: 77—80。

［12] 张水坤。

矩阵在计算机科学技术的应用[J]。

科技经济市场，2007,9：1—6. 。

[13］汪祥。

几类结构矩阵的快速算法及其应用［M］。

2004,8: 3-5.
［14] 中外数学简史编写组. 中国数学简史［M]。

山东: 山东教育出版杜，1986，7：16—20。

[15］Bernkopf，Michael。

A history of infinite matrices[M]。

Archive for history of exact sciences。

1968,4：20—25。

.
［16] A。

D. 亚历山大洛夫等著。

数学一它的内容，方法和意义[M]。

北京：科学出版社。

2001,4:
20-21。