华工数学实验七 特征值和特征向量.docVIP

PAGE 2 实验七 特征值与特征向量 地 点： ????计算中心 202房 实验台号： 30 实验日期与时间： 2018年6月6日 评 分： 预习检查纪录： 实验教师： 刘小兰 电子文档存放位置： 电子文档文件名： 信息工程3班-30-邢靖-实验七.docx 批改意见： 实验目的 掌握特征值、特征向量、特征方程、?矩阵的对角化等概念和理论; 掌握将矩阵化为相似对角矩阵的方法; 理解由差分方程xk+1=Axk所描述的动态系统的长期行为或演化; 提高对离散动态系统的理解与分析能力。 问题1 1.当捕食者-被捕食者问题中的捕食参数p是0.125时，试确定该动态系统的演化(给出xk的计算公式).猫头鹰和森林鼠的数量随着时间如何变化?该系统趋向一种被称为不稳定平衡的状态。如果该系统的某个方面(例如出生率或捕食率)有轻微的变动，系统会如何变化? 2.1实验原理 1.特征值与特征向量 2.特征值与特征向量的求法 3.矩阵的对角化 4.离散线性动态系统 5.eig命令 函数: d=eig(A) 功能:求矩阵A的特征值。 说明:返回一列向量d，包含方阵A的所有特征值。 函数: [V,D]=eig(A)或[V,D]=eig(X,nobalance) 功能:求矩阵A的特征值和特征向量。 说明:生成特征值矩阵?D和特征向量构成的矩阵V,使得使得A*V=V*D。矩阵D由A的特征值在主对角线构成的对角矩阵。V是由A的特征向量按列构成的矩阵。[V,D]=eig(A)中，先对A作相似变换再求A的特征值和特征向量;而[V,D]=eig(A,nobalance)中，直接求矩阵A的特征值和特征向量。 2.2算法与编程 % ex1.m求特征值与特征向量 clc A = [0.5 0.4;-0.125 1.1]; [pc,lambda] = eig(A); %求A的特征值和对应的特征向量 [Y,I] = sort(diag(abs(lambda)),descend);%对特征值的绝对值降序排列 temp = diag(lambda); lambda = temp(I) %输出按特征值的绝对值降序排列的特征值 pc = pc(:,I) %与特征值对应的特 %P8_1.m捕食者-被捕食者解的图像表示 % P8_1.m %捕食者-被捕食者解的图像表示 clear, clc a = 0; b = 2000; c = a; d = b; p = 0.1; %确定画图范围 n = 100; %序列迭代次数 xlabel(|\lambda| 1,|u|1) axis([a b c d]),grid on,hold on x = linspace(a,b,30); A = [0.5 0.4;-0.125 1.1]; %特征值绝对值1 [pc,lambda] = eig(A); %求A的特征值和对应的特征向量 [Y,I] = sort(diag(abs(lambda)),descend); %对特征值的绝对值降序排列 temp = diag(lambda); lambda = temp(I) %输出按特征值的绝对值降序排列的特征值 pc = pc(:,I) pc = -pc; z1 = pc(2,1)/pc(1,1)*x; %特征向量v1 z2 = pc(2,2)/pc(1,2)*x; %特征向量v2 h = plot(x,z1),set(h,linewidth,2), text(x(7),z1(7)-100,v1) h = plot(x,z2),set(h,linewidth,2), text(x(20),z2(20)-100,v2) button = 1; while button == 1 [xi yi button] = ginput(1); %用鼠标选初始点 plot(xi,yi,go),hold on X0 = [xi;yi]; X = X0; for i=1:n X = [A*X, X0]; %用这种方式迭代，并画图 h = plot(X(1,1),X(2,1),R.,X(1,1:2),X(2,1:2),r-); hold on text(X0(1,1),X0(2,1),x0)