Matlab在线性代数中的应用.pptVIP

2021/3/14 * 几种变换的行列式与特征值 2021/3/14 * 二维矩阵特征值的几何意义 二维矩阵的特征值表示该变换在原图形的特征向量的方向上的放大量。例如矩阵A1在第一特征向 量 方向的特征值为 ,即横轴 正方向的增益为?1,其结果是把原图中横轴正方向的部分变换到新图的负方向去了; A1在第二特 征向量 的方向的特征值为λ1(2)=1, 即纵轴正方向的增益为1,因而保持了新图和原图在纵轴方向尺度不变。 2021/3/14 * 用eigshow函数看特征值 对于比较复杂的情况,完全凭简单的几何关系去想像是困难的,应当用eigshow函数,联系x和Ax的向量图来思考。 键入eigshow(A4) 。绿色的x表示原坐标系中的单位向量,可以用鼠标左键点住x并拖动它围绕原点转动。图中同时出现以蓝色表示的Ax向量,它表示变换后的新向量。当两个向量处在同一条直线上时（包括同向和反向）,表示两者相位相同,只存在一个（可正可负的）实数乘子λ, Ax?λx 2021/3/14 * Eigshow(A4)产生的图形 2021/3/14 * eigshow([1,2; 2,2])的图形 将eigshow([1,2; 2,2])粘贴到命令窗 2021/3/14 * A是对称实矩阵的情况 特别要注意A是对称实矩阵的情况,所谓对称矩阵是满足AT?A的矩阵。,对2?2矩阵,只要求A(1,2)? A(2,1)。例如令A =[1,2;2,2] 再键入eigshow(A),这时的特点是:Ax?λx出现在Ax椭圆轨迹的主轴上,所以两个特征值分别对应于单位圆映射的椭圆轨迹的长轴和短轴。此时A的特征值为 -0.5616和 3.5616,可以和图形对照起来看。 2021/3/14 * 人口迁徙模型 设在一个大城市中的总人口是固定的。人口的分布则因居民在市区和郊区之间迁徙而变化。每年有6%的市区居民搬到郊区去住,而有2%的郊区居民搬到市区。假如开始时有30%的居民住在市区,70%的居民住在郊区,问10年后市区和郊区的居民人口比例是多少?30年、50年后又如何? 这个问题可以用矩阵乘法来描述。把人口变量用市区和郊区两个分量表示,一年以后,市区人口为xc1? (1?0.06) xc0?0.02xs0,郊区人口xs1? 0.06xc0 ? (1?0.02)xs0, 2021/3/14 * 问题的矩阵描述 用矩阵乘法来描述,可写成: 从初始到k年,此关系保持不变,因此上述算式可扩展为 , 故可用程序ag981n进行计算: A?[0.94,0.02;0.06,0.98], x0?[0.3;0.7] x1?A*x0, x10?A^10*x0, x30?A^30*x0, x50?A^50*x0 得到: 2021/3/14 * 本题特征值和特征向量的意义 无限增加时间k,市区和郊区人口之比将趋向一组常数0.25/0.75。为了弄清为什么这个过程趋向于一个稳态值,我们改变一下坐标系统。在这个坐标系统中可以更清楚地看到乘以矩阵A的效果,先求A的特征值和特征向量,得到 令 它是特征向量的整数化,得到 2021/3/14 * 6 直线和平面的快速绘制程序 平面曲线的快速绘制程序 ezplot(‘ ‘,[a,b]) 引号中函数可以只有一个自变量,代表显函数 ezplot(‘f(x)’, [a,b]) 系统将在 a x b的范围内画出 f = f(x) 引号中的函数若有两个自变量,那就代表隐函数,其典型格式为 ezplot(‘f(x,y)’, [a,b]) 系统将在 a x b的范围内画出 f(x,y)=0。 [a,b]的默认值为[-2π, 2π] 2021/3/14 * 曲线快速绘制举例 ezplot(‘x1+0.2*x2^3+1’) 画出 在x1=[-2π, -2π]的曲线 画多条曲线可按下列方法编程 s1=x1+0.2*x2^3+1 % 方程1 s2=3*x1+2*x2+3 % 方程2 ezplot(s1),hold on % 画方程1,保持 ezplot(s2),grid on % 画方程2,加网格 [x1,x2]=solve(s1,s2) % 解联立方程1,2 2021/3/14 * 解的结果 2021/3/14 * 解的说明 在线性代数中,遇到的都是一次函数,所以不会出现曲线。我们故意用一个三次曲线来说明ezplot的用法,是为了使读者知道,这个命令不限于画直线。 MATLAB用solve命令解题是采用了符号运算工具葙,它的数字精度是32位十进制,而不是一般数值计算时的16位十进制。尽管在本题中有两有效数字就够