线性代数应用例.pptVIP

线性代数应用实例;例 1 平板稳态温度的计算;平板温度计算的模型;A=[4,-1,-1,0; -1,4,0,-1; -1,0,4,-1; 0,-1,-1,4]; b=[30; 50; 60; 80]; U=rref([A,b]);向高阶系统扩展;例 2 交通流的建模;A=[1,-1,0,0;0,1,-1,0;0,0,1,-1;-1,0,0,1]; b=[-100;72;37;-9]; U=rref([A,b]);把上述模型扩展到多个十字路，乃至整个城市，就构成高阶的线性代数方程组。例如下面的 6 节点交通流图，它就要由 6 个方程和 7 个变量来描述。用行最简型方法可以知道，它的解将包括两个自由变量。其物理意义类推。;左图描述了四个城市之间的航空 航线图，其中1、2、3、4 表示四 个城市；带箭头线段表示两个城 市之间的航线。设行号表示起点 城市，列号为到达城市，则 定义邻接矩阵 A 为：;转机航线的数学模型;多次转机到达的城市;例 4 行列式的几何应用;由向量 和 所构成的平行四 边形的面积为行列式 的绝对值。 计算的MATLAB语句为： S=abs(a1*b2-a2*b1) 如果给出的是三角形三个顶点坐标[a1,b1]， [a2,b2]， [a3,b3]，求该三角形面积，则有： MATLAB写成S=abs(det([a2-a1,b2-b1; a3-a1,b3-b1]));多边形可以划分为多个三角形来计算。 先对三角形面积计算构成一个函数程序； 这个子程序名为：cal_area3(A,B,C) A,B,C为三个顶点的二维坐标向量 凸多边形面积只需多次调用这个函数程序； 例如五边形ABCDE，可由 S5= cal_area3(A,B,C)+ cal_area3(A,C,D)+ cal_area3(A,D,E) 求得。也可由多边形面积子程序cal_arean(A)计算。;MATLAB 程序;成药;（1）某医院要购买这 7 种特效药，但药厂的第 3号和第 6 号特效药已经卖完，请问能否用其它特效药配制出这两种脱销的药品。 分析：即 3, 6 向量与其他向量是否线性相关 （2）现在该医院想用这 7 种中草药配制三种新的特效药，下表为新药所需的成分质量 (单位: 克) 。请问如何配制。 分析：这是新药向量与原来药向量是否线性相关的问 题。;;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]; U1=[u1,u2,u3,u4,u5,u6,u7] [V1,r]=rref(U1);运行结果;因此，8，9 两种药可以配出，第 10 种药则不能配出。;假设一个城市的总人口数是固定不变，但人口的分布 情况变化如下：每年都有 5％ 的市区居民搬到郊区； 而有 15％ 的郊区居民搬到市区。若开始有 700000 人 口居住在市区，300000 人口居住在郊区。请分析： （1）10 年后市区和郊区的人口各是多少？ （2）30 年后、50 年后市区和郊区的人口各是多少？ （3）分析（2）中数据相似的原因。;解 这个问题可以用矩阵乘法来描述。令人口变量;A=[0.95,0.15;0.05,0.85]; X0=[700000;300000]; X10=A^10*X0 ; 无限增加时间 n，市区和郊区人口之比将趋向常数 0.75/0.25。为了弄清为什么它趋向于一个稳态值，需要 可以求 Ak, 为此可先将 A 对角化, 然后求其幂。;% 分析 n 年后城市人口分布 clear A=[0.95,0.15; 0.05,0.85]; X0=[700000; 300000]; [P,lambda]=eig(A); syms n % 定义符号变量 n Xn=P*lamda.^n*inv(P)*X0 ;x;其中矩阵：;系数矩阵 A 的行列式为范德蒙 (Vandermonde) 行列式， 且五个坐标点的横坐标各不相同，则该行列式不等于零， 所以方程组有唯一解。; 把五个点的坐标值分别代入二次多项式函数，可以得 到如下线性方程组：;该方程组有三个未知数，但有五个方程，进一步分析 可以得到该方程组无解，即不存在一个二次多项式