矩阵在数学计算中的应用

  • 格式:pdf
  • 大小:134.33 KB
  • 文档页数:3

第6卷第3期 2015年3月 黑龙江科学 HEILONGJIANG SCIENCE V01.6 No.3 March 2015 

矩阵在数学计算中的应用 

于荣格,江瑞侠 

(沧州师范学院数学系,河北沧州061001) 

摘要:矩阵是高等代数的一个非常重要的组成部分,通过例题给出了矩阵在数学计算中的一些应用,包 

括算由几个向量所生成的子空间的维数,利用矩阵的秩来判定一个方程组是否有解,通过对二次型的矩 

阵进行合同变换,将二次型化成标准形。 

关键词:子空间:线性方程组:二次型 

中图分类号:O151.21 文献标志码:A 文章编号:1674—8646(2015)03—0014-03 

Applications of matrix in mathematics computation 

YU Rong—ge,JIANG Rui-xia 

(Department of Mathematics,Cangzhou Normal University,Cangzhou 061001,China) 

Abstract:Matrix is a very impo ̄ant part of higher algebra by some application examples ale given in the matrix 

math calculations,including the nlllnber of dimensions considered by several vector subspace generated by the matrix to determine the rank of a system of equations is solvable,by contact quadratic form matrix transformation, 

the quadratic form into standard. 

Keywords:Subspace;lineal equation group;quadratic form 

行列式和矩阵在高等代数中占据着非常重要的地 

位,笔者给出矩阵在高等代数的三类不同计算中的应用。 

1 计算由几个向量所生成的子空间的维数 

等于矩阵 =f 1的秩,所以我们只要求出了矩阵A 1 3 l 

1 2 3 

l 2 l 

1 2 —1 变换第1和嚣4行 

第2行加上第1行 第4行加l卜第l行朵-2 

变换第25u第3行 l 1 

0 2 

O 1 

0 —1 

1 1 

O 1 

O 2 

0 1 

第3杼加L第2行乘 2 第4行加L第2行乘一J 、 

变换第3和第4行 1 1 

0 1 

0 0 

0 O 

1 1 

0 1 

0 O 

0 0 1 1 

.1 l 

O 1 

2 1 

2 .1 

4 2 

2 1 

1 3 

2 .1 

2 l 

4 2 

.1 3 

2 .1 

2 l 

O 0 

1 4 

2 1 

2 l 

1 4 

0 0 2 .1 

2 3 

2 1 

3 1 

矩阵A的秩为3,所以由向量 】, 2,a3,a 4生成的子空 

间维数为3。 

2利用矩阵的秩来判定一个方程组是否有解 

收稿日期:2014—12—29 

作者简介:于荣格(1965一),女,河北东光人,副教授,理学硕士,从事群论研究。 

l4 2 _ O 1 —.....。...。.........。...。.。... L 

= 由[1】中的第四章,我们有下面的定理 

引理1[1]: 

(线性方程组可解的判别法)线性方程组 

j乏。1 X。I+口a:1。2 x:2十+… ̄ ̄ ̄+aCI ln X n =:b  ̄, 

+a…

m2X2…

"-'l-.…

..+C…

lmnX…

n 

有解的充分且必要条件是: 

它的系数矩阵和增广矩阵有相同的秩。 

根据此定理我们可以判断一个方程组是否有解。 

例2: 

解: 判断线性方程组I4x + :一 ,:7是否有解。 

增广矩阵B=r三 ÷二: 

第2行加j二第l行熏.2 第3行加上第1杼乘.4 

鸯珏行加上第2行渠 

秩: 

l l 

l 0 B=秩l0 

厂1 2 f 

2 —3 

=秩f 4 1 2 

—7 

—7 

2 

—7 

0 一l 

3 

3 2 

—1 

一l 

2 —1 2] F1 2 

—7 3 —1 f f 0 ~7 

0 0 0 j=秩l0 0 

所以,由引理1知方程组有解。 

3通过对二次型的矩阵进行合同变换将二次型化成 

标准形 

定义1E ]: 

数域P上n×n矩阵A,B称为合同的,如果有数 

域P上n×n矩阵C,使B=C'AC 

用初等变换把二次型的矩阵化为对角矩阵,为了 

保证所得的矩阵与原矩阵合同,我们要成对地进行初 等变换,即作一次初等行变换后必须作一次相同的列 

变换。 

设二次型厂( ,X2,… )=X'AX。对2n×n矩阵 

[ ]进行初等变换,将A化成对角形,并且保证了合同 

关系。 

如果 

・. 1 l 

P j 

则f(x , 2,… )经线性变换X=PY化成标准形 

dlYl。+ 2Y2 +…+ Y 。。 

f(xl, 2,x3)= l 一 +2xl 2+2x2 3成标准形。 

二次型I厂( ,X , )的矩阵为 

=㈤ 

1 l 0 

1 0 l 

O 1 —1 

1 O 0 

0 1 0 

0 0 l 第l行柔.1加垂0第2行 

第1列乘.1加到第2N 0 

—1 

1 

一l 

l 

0 0 

l 

—l 

0 

O 

1 l O 

.1 1 

1 .1 

0 O 

1 0 

0 1 

(下转第17页) 

15 1●●,● ●J 

1●,●●●能量,发现掺杂后的ZnO晶胞体积增大。这是因为Pm 

杂质引入使得晶胞内电子间排斥力增加,拉伸了原子 

间距,且Pm原子体积大于Zn原子,因此计算结果发 

现,ZnO晶胞体积增大表明我们的工作是可靠的。 

Pm—ZnO的能带结构如图1所示。纵坐标为能 

量,0点为费米能级。我们测出了纯ZnO禁带宽度为 

0.834eV。这个值与理论值有差距,这是由于GGA近似 

带来的影响,可以用剪刀算符加以修正,然而本文主要 

进行相对值的比较,因此GGA近似并不影响对结果的 

讨论。 

图1 Pm—Zn0能带结构 

Fig.1 Pm—Zn0 energy band structure 

对比电子态密度可知,ZnO价带主要由。一2p态 

构成,导带主要由Zn一3d态构成。Pm掺杂后,在禁带 

区域产生杂质能级。由图1可知,Pm的引入使得原子 

间排斥力增加,使得导带由能量1附近上移至1.5附 

近;价带由0附近下降至一1附近,因此Pm的掺杂使 

得ZnO的原禁带宽度增加。然而Pm的4f电子态在禁 

带区域引入了杂质能级,这使得价带电子可进行分级 

跃迁至导带,且价带已经与杂质能级杂化耦合,因此实 

际禁带宽度仅为0.691eV。所以,Pm的引入有效降低 

了ZnO的禁带宽度,使价带电子更易跃迁至导带。 

2.2吸收光谱分析 对纯ZnO和Pm-ZnO进行吸收光谱分析,如图 

2所示。可以发现,掺杂后的吸收光谱明显发生红移, 

吸收限向红光区移动,这是因为杂质引入使得电子跃 

迁所需的能量减小,使得频率较低的光也可激发电子。 

这一效应有利于光催化材料活性的增强。 

图2吸收光谱图 

Fig.2 Absorption of light spectrum 

3结论 

通过第一性原理,计算了纯ZnO和Pm—ZnO的 

电子结构及吸收光谱,结果表明:掺杂后,Pm原子4f 

电子态产生的杂质能级有效降低了ZnO的禁带宽度, 

使得价带电子可吸收较小能量的光分级跃迁至导带, 

从而使得ZnO的吸收光谱发生红移,提升了其光催化 

活性 

参考文献: 

【1 J XuAW,GaoY,LiuH Q.Preparation,Characterization,andtheir Photocatalytic 

Activities ofRare,一Earth—Doped"rio2 Nanoparticles[J].J.Catal,2002,(2O):207. 12j Wang C,Ao Y H,Wang P F,et a1.Controlled synthesis in large—scale of CdS 

mesospheres and photocatalytic activity[J].Mater Lett,2010,(64):439—441. [3]李聪,孙宵宵,于淼.Pr替位式掺杂GaAs光电性能的第一性原理研究[J]. 

牡丹江教育学院学报,2014,(146):96. 

【4 J P ̄yne M C,Teter M P,Allan D C,et a1.Joannopoulos J D.A based on the 

density functional theory(DFF)and the intermediate neglect of the diferential 

overlap INDO method lJj.Rev.Mod.Phys,1992,(64):1045. 

【5 J Perdew J P,Burke K.Emzerhof M Generalized gradient approximation made simple【Jj.Phys.Rev.Lett,1996,(77):3865. 

1 6 j Vanderbih D.Soft Self-consistent Pseudo Potential in a Generalized Eigenvalue Formalism【J J.Phys.Rev.B,1990,(41):7892. 

(上接第15页) 

第2行加到筘3行 

第2列加到帮3列 l 0 

0 一l 

O 0 

l 一1 

0 l 

O 0 

0 

一l 

O 

一1 

1 

0 0 

O 

O 

一1 

l 

1 令x:[ ]Y,即{兰 l一 

就得/( l, 2, 3)= l 一 2 。 

参考文献: 

[1]张禾瑞,郝钢新.高等代数[M].北京:高等教育出版社,1983. 

[2]北京大学数学系几何与代数教研室代数小组.高等代数[M].北京:高等教 

育出版社,2003. 

17