SXSY_04a-线性代数实验

1、/根标题=数学实验/主标题=线性代数实验/作者=/地址1=电子科技大学数学科学学院/地址2=/地址3=/标题缩减级别=11.1实验基础：线性代数知识11.1.1解线性方程组11.1.2矩阵特征值21.2线性方程组求解21.2.1应用问题：在减肥食谱中的应用21.3矩阵的幂与特征值、特征向量31.4Leslie模型及其应用61.4.1案例1. 简单迁移模型61.5线性变换81.1 实验基础：线性代数知识1.1.1 解线性方程组解线性方程组：Matlab命令：Ab示例：A=rand(3,3);truesol=rand(3,1)b=A*truesolmysol=Ab比较数值结果与真解.1.1.2 矩

2、阵特征值A是n阶方阵,求非零向量和数称为特征向量, 称为特征值.lamda=eig(A) 计算A的特征值,这里lamda是A的全部特征值构成的列向量。P,D=eig(A) 计算出A的全部特征值和对应的特征向量. 其中, D是对角矩阵,保存矩阵A的全部特征值; P是满阵, P的列向量构成对应于D的特征向量组。1.2 线性方程组求解1. 线性方程组的理论应用已经渗透到数学发展的许多分支很多实际问题的处理最后往往归结为比较容易处理的线性方程组的问题2. 线性方程组在工程技术上、空间几何和国民经济的许多领域都有着广泛的应用1.2.1 应用问题：在减肥食谱中的应用 下表是该食谱中的3种食物以及100克每

3、种食物成分含有某些营养素的数量。营养每100克食物所含营养（g）减肥所要求的每日营养量脱脂牛奶大豆面粉乳清蛋白质36511333碳水化合物52347445脂肪071.13如果用这三种食物作为每天的主要食物，那么它们的用量应各取多少才能全面准确地实现这个营养要求？以100克为一个单位，为了保证减肥所要求的每日营养量，设每日需食用的脱脂牛奶个单位，大豆面粉个单位，乳清个单位，则由所给条件得解上方程组得，解为.即为了保证减肥所要求的每日营养量，每日需食用脱脂牛奶27.72克，大豆面粉39.19克，乳清23.32克。求解上述线性方程组的程序：>> A=36 51 13; 52 34 74;

4、 0 7 1.1; b=33;45;3; x=Abx = 0.2772 0.3919 0.23321.3 Leslie模型及其应用科学家LesliePH.于1945年引进一种数学方法,利用某一初始时刻种群的年龄结构现状,动态地预测种群年龄结构及数量随时间的演变过程。1.3.1 案例1. 简单迁移模型每年A镇的人口10%迁往B镇;B镇的人口15%迁往A镇. 假设某年A、B两镇人口各有120人和80人. 问两年后两镇人口数量分布如何?设两镇总人口不变, 人口流动只限于两镇之间. 引入变量:表示 A 镇第 年人口数量;表示 B 镇第 年人口数量.由第 年到第 年两镇人口数量变化规律如下, newpa

5、ge/codebegin%title1:线性代数实验%title2:应用实验%title:简单迁移模型A=0.9,0.15;0.1,0.85; X0=120;80;X2=A2*X0D=eig(A)/codeend, 若，则是特征值1的特征向量.思考：在无限长时间后，两镇人数的变化规律与矩阵A有何关联？1.4 矩阵的幂与特征值、特征向量相关知识点：相似对角化，特征值，特征向量，正交矩阵在许多工程计算中，经常遇到计算矩阵的幂的问题，请问有没有快速的计算方法？概念：定义 设都是阶方阵，若存在可逆矩阵，使，则称是的相似矩阵，并称矩阵与相似对进行运算称为对进行相似变换，称可逆矩阵为相似变换矩阵易证矩阵的

6、相似关系满足：自反性：对于任意阶方阵，有与相似；对称性：若与相似，则与相似；传递性：若与相似，与相似，则与相似定理 若阶矩阵与相似，则与有相同的特征多项式，从而与有相同的特征值相似矩阵的一个重要性质：若阶矩阵与对角阵相似，则即为的个特征值定义 若阶矩阵与对角阵相似，则称可以相似对角化，简称可对角化如果矩阵可以相似对角化，在存在可逆矩阵，满足其中，。由得，可知，为矩阵的特征值，其中为对应于特征值的特征向量。则，如果是正交矩阵，则（正交矩阵）。进而得到.这样就大大简化了计算矩阵的幂的复杂度.需要计算矩阵的特征值、特征向量。以对称矩阵为例说明：利用V,D=eig(A)得：V = 0.8658 0.1

7、968 0.2961 0.3521 -0.4463 0.7064 0.2603 0.4837 0.0822 0.0201 -0.8573 0.5078 -0.2108 -0.6796 0.3311 0.6198D = 1.6310 0 0 0 0 1.8540 0 0 0 0 2.7441 0 0 0 0 6.1709在命令行输入V*V'得：V*V'ans = 1.0000 -0.0000 0.0000 0.0000 -0.0000 1.0000 0.0000 0.0000 0.0000 0.0000 1.0000 00.0000 0.0000 0 1.0000可知eig所计算

8、得到特征向量矩阵V已经是正交矩阵，则计算A的5次幂可用如下命令完成：V*D.5*V'输入A5-V*D.5*V'得：ans = 1.0e-11 * -0.0227 -0.0227 -0.0909 -0.0455 -0.0227 -0.0455 -0.1364 -0.0455 -0.0909 -0.1364 -0.1819 -0.1819 -0.0227 -0.0455 -0.1364 0说明二者相差很小。1.5 特征值的应用newpagenewpageI=imread('cameraman.tif');U,S,V=svd(double(I);s=diag(S);n1=5;Snew=diag(s(1:n1);zeros(size(s,1)-n1,1);figure,ims