实验8：矩阵应用.pptVIP

实验八 矩阵应用 数学实验 数学实验 8.1 人口迁移问题 本实验从人口迁移问题出发，讨论了方程组、矩阵乘法、特 征值和特征向量及矩阵对角化问题，并利用MATLAB软件对线性代 数理论进行了验证. 例1 假设一个城市的总人口数是固定不变的，但人口的分布情 况变化如下：每年都有5%的市区居民搬到郊区，而有15%的郊区 居民搬到市区.若开始有700000人口居住在市区，300000人口居 住在郊区，请利用MATLAB软件分析： （1）10年后市区和郊区的人口各是多少？ （2）30年后、50年后市区和郊区的人口各是多少？ （3）分析（2）种数据相似的原因. 数学实验 8.1 人口迁移问题 数学实验 （1）可以得到10年后市区和郊区的人口分布为 A=[0.95,0.15;0.05,0.85]; X0=[700000;300000]; X10=A^10*X0 输出结果： 输入： X10 = 1.0e+005 * 7.4463 2.5537 可以知道，10年后市区和郊区人口数约为744630和255370 数学实验 （2）同理可以得到30年后和50年后的人口分布为 和 A=[0.95,0.15;0.05,0.85]; X0=[700000;300000]; X30=A^30*X0 X50=A^50*X0 输出结果： 输入： X30 = 1.0e+005 * 7.4994 2.5006 X50 = 1.0e+005 * 7.5000 2.5000 从计算结果可以看出，30年后和50年后的市区和郊区人口分布非常接近. 数学实验 （3）若矩阵A可以对角化，则可以利用矩阵对角化的方法来计算矩阵 . A=[0.95,0.15;0.05,0.85]; X0=[700000;300000]; [P,lamda]=eig(A); syms n Xn=P*lamda.^n*inv(P)*X0 输出结果： 输入： Xn = 750000-50000*(4/5)^n 250000+50000*(4/5)^n 从计算结果可以看出，当n越大时，市区和郊区人口就越趋近于750000和250000. 在MATLAB的M文件编辑器中编写程序la24.m，编写完后，在命令窗口输入La24 数学实验 8.2药方配制问题 本实验给出了一个中成药的配制问题，其中涉及到了向量 组的线性相关性、向量组的最大无关组、向量的线性表示及向 量空间等知识.并利用MATLAB软件完成具体计算. 例2 某中药厂用9种中草药（A、B、……、I），根据不同的比 例制成了7种特效药.表1给出了每种特效药每包所需各种成分的 质量（单位：克）,某医院要购买这7种特效药，但药厂的第3号 和第6号特效药已经卖完，请问能否用其他特效药配制出这两种 脱销的药品. 数学实验 1号成药 2号成药 3号成药 4号成药 5号成药 6号成药 7号成药 A 10 2 14 12 20 38 100 B 12 0 12 25 35 60 55 C 5 3 11 0 5 14 0 D 7 9 25 5 15 47 35 E 0 1 2 25 5 33 6 F 25 5 35 5 35 55 50 G 9 4 17 25 2 39 25 H 6 5 16 10 10 35 10 I 8 2 12 0 0 6 20 表1 特效药的成分含量 单位：克 数学实验 u1=[10;12;5;7;0;25;9;6;8]; u2=[2;0;3;9;1;5;4;5;2]; u3=[14;12;11;25;2;35;17;16;12]; u4=[12;25;0;5;25;5;25;10;0]; u5=[20;35;5;15;5;35;2;10;0]; u6=[38;60;14;47;33;55;39;35;6]; u7=[100;55;0;35;6;50;25;10;20]; U=[u1,u2,u3,u4,u5,u6,u7] [U0,r]=rref(U) 输入： 把每一种特效药都看成一个九维列向量，分析向量组 的线性相关性. 若向量组线性无关，则无法配制脱销的特效药；若向量组线性相关，并且能找到 不含 和 的