华工数学实验报告特征值与特征向量.docx

《华工数学实验报告特征值与特征向量.docx》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、《数学实验》报告学院：电子信息学院专业班级：信息工程电联班学号：姓名：实验名称：特征根与特征方程实验日期：2016/05/31特征根与特征方程1.实验目的掌握特征值、特征向量、特征方程、矩阵的对角化等概念和理论；掌握将矩阵化为相似对角矩阵的方法；理解由差分方程xk+1=Axk；提高对离散动态系统的理解与分析能力。2.实验任务1.当捕食者-被捕食者问题中的捕食系数p是0.125时，试确定该动态系统的演化（给出xk的计算公式）。猫头鹰和森林鼠的数量随时间如何变化？该系统趋向一种被称为不稳定平衡的状态。如果该系统的某个方面（例如出生率或捕食

2、率）有轻微的变动，系统如何变化？2.杂交育种的目的是培养优良品种，以提高农作物的产量和质量。如果农作物的三种基因型分别为AA，Aa，aa。其中AA-AAAA-AaAA-aaAa-AaAa-aaaa-aaAA11/201/400Aa01/211/21/20aa0001/41/21AA为优良品种。农场计划采用AA型植物与每种基因型植物相结合的方案培育植物后代，已知双亲基因型与其后代基因型的概率。问经过若干年后三种基因型分布如何？要求：（1）建立代数模型，从理论上说明最终的基因型分布。（2）用MATLAB求解初始分布为0.8，0.2，0时，

3、20年后基因分布，是否已经趋于稳定？概率父体－母体基因型后代基因型1.实验过程1.1实验原理1、特征值与特征向量îùAi›zÈÎrIF'.3ŸWfčn@1FTI•łs:•'lŃ"-Ił"llÍ?4Î91?ÑA%hțÎBÒü-ăñFăfi';?ÑAtfiET’1äÛÍiïGR.łM/Ń.HíÎAHèftiïi•Bfr°èù”xü‹.I=ia.-1-żtči•›g.I2、特征值与特征向量的求法s‹iiI)-.t6.łží[l6Üï:r;IA—MI:i:›gși,ï:r,Itlx—‹ml=a,-«’yüúïa‹,.‹„.‹,.I‹,.+

4、ø»‹n«‹A-‹,ij•-oøë.‹i:«ü.r‹‹źäË‹i«ś.ń3、矩阵的对角化.4、离散线性动态系统.y‹i›i‹Ii&l••I>1.B1‹*•=(k•=).›‹*=fx=y.fs4BoP‹aïIis›aixi>iani

5、.ÒÒ6,Ë‹saaue5、eig命令1.1算法与编程3.2.1clear,clca=-20*100;b=-a;c=a;d=b;p=0.1;n=100;xlabel('

6、lambda

7、>1,

8、u

9、<1')axis([0b0d]),gridon,holdonx=linspace(a,b,30);A=[0.50.4;-0.1251.1];[pc,lambda]=eig(A);[Y,I]=sort(diag(abs(lambda)),'descend');temp=diag(lambda);lambda=temp(I)pc=pc(:,I)

10、pc=-pc;z1=pc(2,1)/pc(1,1)*x;z2=pc(2,2)/pc(1,2)*x;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;whilebutton==1[xiyibutton]=ginput(1);plot(xi,yi,'go'),holdonX0=[xi;yi];X=X0;fori=1:nX=[A*X,X

11、0];h=plot(X(1,1),X(2,1),'R.',X(1,1:2),X(2,1:2),'r-');holdontext(X0(1,1),X0(2,1),'x0')quiver([X(1,2),1]',[X(2,2),1]',[X(1,1)-X(1,2),0]',[X(2,1)-X(2,2),0]',p)set(h,'MarkerSize',6),grid,endend3.2.2clear;A=[10.50;00.51;000];X=[0.8;0.2;0];fori=1:20X=A*X;endX20=XX=[0.8;0.2;0]

12、;C=[111]';n=0;whilenorm(X-C,'fro')>1.0e-16C=X;n=n+1;X=A*X;endformatlong;X,n结果分析.1.2.>>X20=0.9999998092651370.0