- 1、本文档共12页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
数学实验报告
学 院:
班 级:
学 号:
姓 名:
完成日期:
实验六 矩阵的特征值与特征向量
问题一
一.实验目的
1.掌握特征值、特征向量、特征方程、矩阵的对角化等概念和理论;
2.掌握将矩阵化为相似对角矩阵的方法;
3.理解由差分方程xk+1 = Axk所描述的动力系统的长期行为或演化;
4.提高对离散动力系统的理解与分析能力.
二.问题描述
当捕食者-被捕食者问题中的捕食参数p是0.125时,是确定该动态系统的演化(给出Xk的计算公式)。猫头鹰和森林树的数量随着时间如何变化?该系统去向一种被称为不稳定平衡的状态。如果该系统的某个方面(例如出生率或者捕食率)有轻微变动,系统会如何变化?
三.问题分析
将线性变换x Axk的作用分解为易于理解的成分,其中特征值与特征向量是分析离散动态系统的关键。
根据已知信息,找到系统对应的差分方程xk+1 = Axk,求出A的特征值和对应的特征向量,再根据不同特征值的个数、绝对值大于1还是小于1、是实特征值还是复数特征值等情形,分析出系统的演化过程。
四.实验过程
问题对应的差分方程为xk+1 = Axk,其中A= 0.5 0.4
-0.125 1.1 , 演化过程求解如下:
第一步:求A的特征值和对应的特征向量。利用如下的代码即可获得:
A=[0.5 0.4;-0.125 1.1];
[pc,lambda]=eig(A);
[Y,I]=sort(diag(abs(lambda)),descend);
temp=diag(lambda);
lambda=temp(I)
pc=pc(:,I)
运行程序可得A的特征值为lambda =
1.0000
0.6000
A 的特征向量pc =
-0.6247 -0.9701
-0.7809 -0.2425
显然,这两个特征向量(即pc的第一列和第二列)是线性无关的,它们构成R2的一组基,为消除小数,选取
V1= 4 V2= 4 P= 4 4 P﹣1AP = 1.00 0
5 1 5 1 0 0.60
第二步:V1用和V2表示x0和xK ,k=1,2….因为{ V1,V2}是R2的一组基,所以存在系数c1和c2,使得
x0= c1 V1+ c2 V2.
因为V1,V2为矩阵A对应于λ=1.0,u=0.6的特征向量,所以AV1=λV1,A V2=λV2,于是
X1=Ax0=A(c1 V1+ c2 V2)= c1λV1+ c2uV2.
X2=Ax1=A(c1λV1+ c2λV2)= c1λ2V1+ c2u2V2.
一般地,
Xk= c1λkV1+ c2ukV2.
= c1 (1.0)k 4 + c2 (0.6)k 4 k=0,1,2,3 ….
5 1
当k趋近于无穷大时,0.6^k 趋近于0,假定c10,则对于所有足够大的k,xk近似地等于c1 (1.0)k V1,写为
Xk≈c1(1.0)k 4
5
K越大,近似程度越高,所以对于足够大的k,
Xk+1≈c1(1.0)k+1 4
5
= Xk
可知猫头鹰和老鼠的数量几乎每月都相当,而且 Xk约为 4
5 的倍数,所以每4只猫头鹰对应着5000只老鼠。
第三步:解的图像表示,见图8-1,其中绿色圆圈代表初始点x0,红色圆点代表迭代序列,箭头代表迭代方向,蓝色直线代表特征向量V1,V2所在的直线。在图8-1中,圆点为鞍点,排斥最快的方向为过圆点和特征向量V1的直线方向。其中V1对应的特征值得绝对值为1.如果x0在这条直线上,则表示c2等于0,且Xk始终在原点。吸引最快的方向由特征向量V2决定,其对应的特征值的绝对值大于1.
相应的代码如下:
% P8_1.m
% 捕食者-被捕食者解的图像表示
clear,clc
a=-20*100;b=-a;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];
[pc,lambda]=eig(A);
[Y,I]=sort(diag(ab
您可能关注的文档
- 数学建模——食品质量安全抽检数据分析.doc
- 社团成立答辩.ppt
- 数学建模收入分配定位.doc
- 社戏(景物描写).ppt
- 数学建模贪心算法.ppt
- 数学建模统计建模论文.pdf
- 射频识别技术(RFID)基础理论.ppt
- 射频微波器件测试解决方案.pdf
- 数学建模一等奖——二孩问题.pdf
- 摄影测量学 第4章 解析空中三角测量.ppt
- 护理学相关知识复习测试卷共三套.doc
- 护理学相关知识复习试题含答案(3套).doc
- 2025届高考语文复习:补写句子 课件.pptx
- 气压带和风带对气候的影响(第1课时)(教学设计).docx
- 气压带和风带对气候影响教学设计2024-2025学年高中地理人教版(2019)选择性必修1.docx
- 《故都的秋》课件 2024-2025学年统编版高中语文必修上册.pptx
- 《屈原列传》课件 2024-2025学年统编版高中语文选择性必修中册.pptx
- 《巫溪家乡文化》课件-2024-2025学年高一语文同步备课课件(统编版必修上册).pptx
- 《苏武传》课件 2023-2024学年统编版高中语文选择性必修中册.pptx
- 郑州中控ZKTime8.3 WEB考勤软件培训文档.pptx
文档评论(0)