数学实验特征值与特征向量..docVIP

数学实验报告

学院：?

班级：

学号：

姓名:?

完成日期：

实验六矩阵得特征值与特征向量

问题一

一.实验目得

1、掌握特征值、特征向量、特征方程、矩阵得对角化等概念与理论；

２、掌握将矩阵化为相似对角矩阵得方法；

3、理解由差分方程xk+１=Axk所描述得动力系统得长期行为或演化;

4、提高对离散动力系统得理解与分析能力、

二.问题描述

当捕食者－被捕食者问题中得捕食参数p就是0、125时，就是确定该动态系统得演化（给出Xk得计算公式）。猫头鹰与森林树得数量随着时间如何变化？该系统去向一种被称为不稳定平衡得状态。如果该系统得某个方面(例如出生率或者捕食率）有轻微变动，系统会如何变化?

三。问题分析

将线性变换xAxｋ得作用分解为易于理解得成分，其中特征值与特征向量就是分析离散动态系统得关键。

根据已知信息，找到系统对应得差分方程xｋ+１＝Ａxｋ，求出A得特征值与对应得特征向量,再根据不同特征值得个数、绝对值大于1还就是小于1、就是实特征值还就是复数特征值等情形，分析出系统得演化过程。

四．实验过程

问题对应得差分方程为xk+1＝Axk,其中A=0、50、4

-0、１251、1,演化过程求解如下:

第一步:求Ａ得特征值与对应得特征向量.利用如下得代码即可获得:

A=［0、50、4；-0、1251、１］;

［pｃ，ｌamｂda］＝ｅｉg（Ａ)；

［Ｙ,I］＝sort（diag（aｂs（lamｂdａ)），＇descend）；

tｅmp=dｉａg（lambda)；

laｍｂdａ=ｔemp（I）

pc=ｐc(:，I)

运行程序可得A得特征值为lａmbda=

1、0００0

0、6000

A得特征向量ｐc=

—0、6247—0、９７０1

－0、7809-０、242５

显然，这两个特征向量(即pc得第一列与第二列）就是线性无关得,它们构成R２得一组基,为消除小数,选取

V１=4V2=4Ｐ＝44P﹣1ＡP=1、０００

515100、6０

第二步：V１用与V2表示x0与ｘＫ,k=1，2…、因为｛V1，V2｝就是R２得一组基,所以存在系数c1与c2，使得

x0=c１Ｖ1+c2V２、

因为Ｖ1,V2为矩阵A对应于λ=1、0，ｕ=0、6得特征向量,所以ＡV1=λV1，AＶ2=λV2,于就是

Ｘ１=Ax０=A(ｃ１Ｖ1+ｃ2V2)=c1λV1＋c2ｕV2、

Ｘ2＝Ax１=A(c1λV1+ｃ2λV2)=c１λ2Ｖ1+c2u2Ｖ2、

一般地，

Xｋ＝c１λkV1+c2ukV2、

＝c1(1、0）ｋ4＋ｃ2（0、6)k4k=0,1,2,３…、

51

当k趋近于无穷大时，0、6^k趋近于0，假定ｃ1０，则对于所有足够大得k,xk近似地等于ｃ1(１、0)kＶ１，写为

Xｋ≈ｃ1(1、0）k4

5

K越大，近似程度越高，所以对于足够大得k，

Xｋ+１≈c１(1、0)k+14

5

＝Xk

可知猫头鹰与老鼠得数量几乎每月都相当，而且Ｘk约为４

５得倍数，所以每4只猫头鹰对应着5００0只老鼠。

第三步：解得图像表示,见图８—1,其中绿色圆圈代表初始点x0，红色圆点代表迭代序列,箭头代表迭代方向,蓝色直线代表特征向量V1,V2所在得直线.在图８-1中，圆点为鞍点,排斥最快得方向为过圆点与特征向量V１得直线方向。其中V1对应得特征值得绝对值为１、如果x0在这条直线上,则表示c2等于０,且Ｘk始终在原点。吸引最快得方向由特征向量V2决定，其对应得特征值得绝对值大于1、

相应得代码如下：

％P8_1、m

%捕食者—被捕食者解得图像表示

clｅａr，cｌｃ

ａ＝－20*100;ｂ=-ａ；c＝a;d=ｂ；p=0、1；

n=１00；

xｌabel（|\lamｂdａ|=１,｜ｕ｜＜1＇)

axis（［ａｂcd]），ｇridon，hoｌｄｏn

ｘ=ｌinｓｐace（a,b，30)；

A=[0、50、4；－0、1251、1］；

[ｐc,lａｍbｄa］=eig(A）;

[Ｙ，Ｉ]=