线性代数应用举例.pptx

《线性代数应用举例.pptx》由会员分享，可在线阅读，更多相关《线性代数应用举例.pptx（51页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、会计学1线性代数应用举例线性代数应用举例例例 1 平板稳态温度的计算平板稳态温度的计算为了计算平板形导热体的温度分布，将平板划分为许多方格，每一个节点上的稳态温度将等于其周围四个节点温度的平均值。由此可得出阶数与节点数相同的线性方程组，方程的解将取决于平板的边界条件。这个方法可以用来计算飞行器的蒙皮温度等。T1T2T3T4第1页/共51页平板温度计算的模型整理为第2页/共51页A=4,-1,-1,0;-1,4,0,-1;-1,0,4,-1;0,-1,-1,4;A=4,-1,-1,0;-1,4,0,-1;-1,0,4,-1;0,-1,-1,4;b=30;50;60;80;b=30;50;60;8

2、0;U=rref(A,b)U=rref(A,b)MATLAB 程序(ma1)运行结果为：U=1.0000 0 0 0 21.2500 0 1.0000 0 0 26.2500 0 0 1.0000 0 28.7500 0 0 0 1.0000 33.7500第3页/共51页向高阶系统扩展则要解 25 阶的线性方程组。运行书上的程序得温度分布如下 将平板分割得愈细，求出的解就愈精确。如果把上述区域分成 25 个点如右MATLAB 程序程序ma2第4页/共51页例例 2 交通流的建模交通流的建模对于一个有双向车流的十字路口，根据流出流入车数相等的规则，可以列出下列方程组:节点A：x1360 x22

3、60节点B：x2220 x3292节点C：x3320 x4357节点D：x4260 x1251相应的矩阵方程为：第5页/共51页A=1,-1,0,0;0,1,-1,0;0,0,1,-1;-1,0,0,1;A=1,-1,0,0;0,1,-1,0;0,0,1,-1;-1,0,0,1;b=-100;72;37;-9;b=-100;72;37;-9;U=rref(A,b)U=rref(A,b)MATLAB 程序(ma3)运行结果为：U=1 0 0 -1 9 0 1 0 -1 109 0 0 1 -1 37 0 0 0 0 0 由于 U 的最后一行为全零，也就是说，四个方程中实际上只有三个独立。x4 可

4、以任设，因为如果有一些车沿此路口环行，对方程无影响，故方程组的解可如上表示.第6页/共51页n n把上述模型扩展到多个十字路，乃至整个城市，把上述模型扩展到多个十字路，乃至整个城市，就构成高阶的线性代数方程组。例如下面的就构成高阶的线性代数方程组。例如下面的 6 6 节节点交通流图，它就要由点交通流图，它就要由 6 6 个方程和个方程和 7 7 个变量来描个变量来描述。用行最简型方法可以知道，它的解将包括两述。用行最简型方法可以知道，它的解将包括两个自由变量。其物理意义类推。个自由变量。其物理意义类推。向高阶系统扩展第7页/共51页左图描述了四个城市之间的航空航线图，其中1、2、3、4 表示四

5、个城市；带箭头线段表示两个城市之间的航线。设行号表示起点城市，列号为到达城市，则定义邻接矩阵 A 为：例例 3 飞机航线问题飞机航线问题第8页/共51页转机航线的数学模型不难证明：矩阵 A2=A*A 表示一个人连续坐两次航班可以到达的城市，矩阵 A3=A*A*A 表示连续坐三次航班可以到达的城市：其中，第 i 行描述从城市 i 出发，可以到达各个城市的情况，若能到达第 j 个城市，记 A(i,j)=1，否则 A(i,j)=0，规定 A(i,i)=0(其中 i=1,2,3,4)。如第 2 行表示：从城市 2出发可以到达城市 3 和城市 4 而不能到达城市 1 和 2。第9页/共51页多次转机到达

6、的城市分析矩阵 A3 的第二行，可以得出：某人从城市 2 出发，连续坐三次航班可以到达城市 2、3 和城市4，不能到达城市 1，而到达城市 3 和城市 4 的方法各有两种。不难看出，转机两次以下的航线的航路矩阵为At2=A+A2+A3程序为（ma4）A=0,1,1,1;0,0,1,1;0,0,0,0;1,1,0,0;At2=A+A2+A3第10页/共51页例例 4 行列式的几何应用行列式的几何应用二阶行列式的几何意义是两个二维向量构成的平行四边形的面积，三阶行列式的几何意义是三个 3 维向量构成的平行六面体的体积。如下图所示，用 MATLAB 软件来实现面积和体积的运算。第11页/共51页n

7、n 实例实例实例实例n n（I I I I）已知三角形已知三角形ABCABC三个顶点的坐标分别为：三个顶点的坐标分别为：(1,2),(3,3),(4,1)(1,2),(3,3),(4,1)，计算该三角形的面积；，计算该三角形的面积；n n（IIIIIIII）已知凸九边形九个顶点的坐标分别为：已知凸九边形九个顶点的坐标分别为：(0,8.5),(3,7),(6,0),(3,-4),(1,-5),(-5,-3),(0,8.5),(3,7),(6,0),(3,-4),(1,-5),(-5,-3),(-7,0),(-5,6),(-3,8),(-7,0),(-5,6),(-3,8),计算该九边形的面积。计

8、算该九边形的面积。n n（IIIIIIIIIIII）在平面坐标系中画出以上三角形和九边在平面坐标系中画出以上三角形和九边形。形。平行四边形面积计算平行四边形面积计算第12页/共51页解解解解：(I)(I)如图所示，三角形如图所示，三角形 ABC ABC 的面积就等于向量的面积就等于向量ABAB和向量和向量ACAC所构成平行四边形面积的一半。所构成平行四边形面积的一半。其中：其中：第13页/共51页n n由向量由向量 和和 所构成的平行四所构成的平行四 边形的面积为行列式边形的面积为行列式 的绝对值。的绝对值。计算的计算的MATLABMATLAB语句为：语句为：S=abs(a1*b2-a2*b1

9、)S=abs(a1*b2-a2*b1)实例实例实例实例给出的是三角形三个顶点坐标给出的是三角形三个顶点坐标aa1 1,b,b1 1，aa2 2,b,b2 2，aa3 3,b,b3 3，求该三角形面积，则有：，求该三角形面积，则有：MATLABMATLAB写成写成S=abs(det(a2-a1,b2-b1;a3-a1,b3-b1)S=abs(det(a2-a1,b2-b1;a3-a1,b3-b1)第14页/共51页n n多边形可以划分为多个三角形来计算。多边形可以划分为多个三角形来计算。n n先对三角形面积计算构成一个函数程序；先对三角形面积计算构成一个函数程序；n n这个子程序名为：这个子程序

10、名为：cal_area3(A,B,C)cal_area3(A,B,C)n nA,B,CA,B,C为三个顶点的二维坐标向量为三个顶点的二维坐标向量n n凸多边形面积只需多次调用这个函数程序；凸多边形面积只需多次调用这个函数程序；n n例如五边形例如五边形ABCDEABCDE，可由，可由S5=cal_area3(A,B,C)+cal_area3(A,C,D)+S5=cal_area3(A,B,C)+cal_area3(A,C,D)+cal_area3(A,D,E)cal_area3(A,D,E)求得。（求得。（MATLABMATLAB程序程序ma4ma4）也可由多边形面积子程序也可由多边形面积子程

11、序cal_arean(A)cal_arean(A)计算。计算。扩展至多边形面积计算第15页/共51页解：解：(II)(II)如图所示，凸九边形面积是由如图所示，凸九边形面积是由9-2=79-2=7个三个三角形面积组成。角形面积组成。第16页/共51页 在在MATLABMATLAB命令窗口运行程序命令窗口运行程序ma5.mma5.m，即可以算出，即可以算出三角形和九边形面积，同时可以得到图形：三角形和九边形面积，同时可以得到图形：第17页/共51页MATLAB 程序程序function s=cal_area3(a,b,c)%a,b,c 应为同形的 2 维行向量或列向量，%格式检验语句略去 ab=

12、b-a;%计算向量ABac=c-a;%计算向量ACif size(ab)=1,2%判读向量AB是否为行向量 A=ab;ac;%构造矩阵Aelse A=ab,ac;ends=abs(det(A)/2;%根据公式计算三角形面积第18页/共51页成药1号成药2号成药3号成药4号成药5号成药6号成药7号成药A10214122038100B1201225356055C531105140D79255154735E012255336F255355355550G94172523925H651610103510I821200620例例 5 药方配置问题药方配置问题第19页/共51页n n（1 1）某医院要购买这

13、）某医院要购买这 7 7 种特效药，但药厂的第种特效药，但药厂的第 3 3号和第号和第 6 6 号特效药已经卖完，请问能否用其它特号特效药已经卖完，请问能否用其它特效药配制出这两种脱销的药品。效药配制出这两种脱销的药品。分析：即分析：即 3,6 3,6 向量与其他向量是否线性相关向量与其他向量是否线性相关n n（2 2）现在该医院想用这）现在该医院想用这 7 7 种中草药配制三种新的种中草药配制三种新的特效药，下表为新药所需的成分质量特效药，下表为新药所需的成分质量 (单位单位:克克)。请问如何配制。请问如何配制。分析：这是新药向量与原来药向量是否线性相关的问分析：这是新药向量与原来药向量是否

14、线性相关的问 题。题。问题及分析思路第20页/共51页1号新药2号新药3号新药A4016288B6214167C14278D4410251E53607F5015580G7111838H416821I145230新药的成分要求新药的成分要求第21页/共51页n nu1=10;12;5;7;0;25;9;6;8;u1=10;12;5;7;0;25;9;6;8;n nu2=2;0;3;9;1;5;4;5;2;u2=2;0;3;9;1;5;4;5;2;n nu3=14;12;11;25;2;35;17;16;12;u3=14;12;11;25;2;35;17;16;12;n nu4=12;25;0;5

15、;25;5;25;10;0;u4=12;25;0;5;25;5;25;10;0;n nu5=20;35;5;15;5;35;2;10;0;u5=20;35;5;15;5;35;2;10;0;n nu6=38;60;14;47;33;55;39;35;6;u6=38;60;14;47;33;55;39;35;6;n nu7=100;55;0;35;6;50;25;10;20;u7=100;55;0;35;6;50;25;10;20;n nU U1 1=u1,u2,u3,u4,u5,u6,u7=u1,u2,u3,u4,u5,u6,u7n n V1V1,r=rref(U,r=rref(U1 1)问题

16、(1)的 MATLAB 程序（ma6)第22页/共51页运行结果运行结果V1=1 0 1 0 0 0 0 0 1 2 0 0 3 0 0 0 0 1 0 1 0 0 0 0 0 1 1 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0r=1 2 4 5 7可见这七种特效药是“相关的”,3、6 两种药可用其它 5种药线性配制出来,但第1、2、4、5、7 种药“无关”。第23页/共51页因此，8，9 两种药可以配出，第 10 种药则不能配出。V2=1 0 1 0 0 0 0 1 3 0 0 1 2 0 0

17、 3 0 3 4 0 0 0 0 1 0 1 0 2 2 0 0 0 0 0 1 1 0 0 0 0 0 0 0 0 0 0 1 0 1 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0s=1 2 4 5 7 10为求第二个问题,把 3 种新药与 7 种原药组成矩阵 U2，求 rref，得：第24页/共51页n ns1=40 62 14 44 53 50 71 41 14s1=40 62 14 44 53 50 71 41 14n ns2=162 141 27 102 60 155

18、 118 68 52s2=162 141 27 102 60 155 118 68 52n ns3=88 67 8 51 7 80 38 21 30s3=88 67 8 51 7 80 38 21 30n nU2=U1,s1,s2,s3U2=U1,s1,s2,s3n nV2 r=rref(U2)V2 r=rref(U2)问题(2)的 MATLAB 程序第25页/共51页假设一个城市的总人口数是固定不变，但人口的分假设一个城市的总人口数是固定不变，但人口的分布布情况变化如下：每年都有情况变化如下：每年都有 5 5 的市区居民搬到郊区；的市区居民搬到郊区；而有而有 1515 的郊区居民搬到市区。若

19、开始有的郊区居民搬到市区。若开始有 700000 700000 人人口居住在市区，口居住在市区，300000 300000 人口居住在郊区。请分析：人口居住在郊区。请分析：（1 1）10 10 年后市区和郊区的人口各是多少？年后市区和郊区的人口各是多少？（2 2）30 30 年后、年后、50 50 年后市区和郊区的人口各是多少年后市区和郊区的人口各是多少？（3 3）分析（）分析（2 2）中数据相似的原因。）中数据相似的原因。例例 6 人口迁徙问题人口迁徙问题第26页/共51页解 这个问题可以用矩阵乘法来描述。令人口变量其中 xn 为市区人口所占比例，yn 为郊区人口所占比例。在 n+1年的人口

20、分布状态为：用矩阵乘法可写成：第27页/共51页A=0.95,0.15;0.05,0.85;X0=700000;300000;X10=A10*X0 开始市区和郊区的人口数为 可以得到 n 年后市区和郊区的人口分布：因此(1)10 年后的人口可用程序计算如下：运行结果为：故市区和郊区人口数约为：744630和255370。第28页/共51页 无限增加时间 n，市区和郊区人口之比将趋向常数0.75/0.25。为了弄清为什么它趋向于一个稳态值，需要求Ak,为此可先将 A 对角化,再求其 k 次幂。对角矩阵的幂次可以化为元素的幂次 余下很容易计算。令其中为对角矩阵，则有第29页/共51页%分析分析 n

21、 n 年后城市人口分布年后城市人口分布(ma7)(ma7)clearclearA=0.95,0.15;A=0.95,0.15;0.05,0.850.05,0.85;X0=700000;X0=700000;300000300000;P,lambda=eig(A);P,lambda=eig(A);syms syms n n%定义符号变量定义符号变量 n nX Xn n=P*lamda.=P*lamda.n n*inv(P)*X0*inv(P)*X0 MATLAB 程序显然,随 n 增大(4/5)n 趋近于零,而 Xn 趋于运行结果为：750000250000第30页/共51页x01234y-270

22、210-75例例 7 多项式插值与拟合多项式插值与拟合求:(1)过这五个点作一个四次多项式函数(2)请根据这五个点，拟合一个二次多项式函数下表给出了平面坐标系中五个点的坐标。并求 x=5 时的函数值 p4(5)。用 MATLAB 绘制多项式函数 p4(x)的曲线、已知点及插值点(5,p4(5)。并用 MATLAB 绘制 p2(x)的曲线及已知的五个点。第31页/共51页其中矩阵：解：(1)根据已知条件，把五个点的坐标值分别代入四次多项式函数，可以得到如下线性方程组：第32页/共51页系数矩阵 A 的行列式为范德蒙(Vandermonde)行列式，且五个坐标点的横坐标各不相同，则该行列式不等于零

23、，所以方程组有唯一解。MATLAB 程序:(ma8)x=0;1;2;3;4;%输入已知点坐标y=-27;0;21;0;-75;A=x.0,x.1,x.2,x.3,x.4;%构造 vandermonde 矩阵a=Ay;%得到适定方程组的唯一解 a运行程序，得到a(1)=-27,a(2)=12,a(3)=26,a(4)=-12,a(5)=1.第33页/共51页 把五个点的坐标值分别代入二次多项式函数，可以得到如下线性方程组：其中，(2)多项式拟合要解一个超定方程第34页/共51页该方程组有三个未知数，但有五个方程，进一步分析可以得到该方程组无解，即不存在一个二次多项式曲线刚好能过已知的五个点。MA

24、TLAB 软件提供了一个利用最小二乘法解决超定方程组解的方法。求系数的公式也是 a=Ay，以找到一条二次曲线来近似地描述已知 5 个点的变化情况。对比插值和拟合的曲线如下图第35页/共51页用平面坐标系中的一个闭合图形来描述刚体，用一个矩阵 X 来表示它。X 的一列表示刚体一个顶点的坐标。为了使图形闭合，X 的最后一列和第一列相同；为了实现刚体的平移运算，给矩阵 X 添加元素值都为 1 的一行，使 X 为 3n 矩阵。若有矩阵：则可以证明，矩阵 Y1 是刚体 X 沿 x 轴正方向平移 c1，沿 y 轴正方向平移 c2 后的结果；矩阵 Y2 是刚体 X 以坐标原点为中心逆时针转动 t 弧度的结果

25、。例例 8 刚体的平面运动刚体的平面运动第36页/共51页x04610853.56.1 6.53.220y014140011664.54.500实 例用下列数据表示字母 A:对 A 进行以下平面运动,并绘制移动前后的图形。(1)向上移动 15,向左移动 30;(2)逆时针转动/3;(3)先逆时针转动3/4,然后向上平移 30,向右平移 20。第37页/共51页解 构造刚体矩阵 X，平移矩阵及转动矩阵。M1=100010-30201M2=10001020301R1=00100R2=00100MATLAB 程序（ma9)第38页/共51页 X=0,4,6,10,8,5,3.5,6.1,6.5,3.

26、2,2,0;0,14,14,0,0,11,6,6,4.5,4.5,0,0;ones(1,12);%构造刚体矩阵 X M1=1,0,-30;0,1,15;0,0,1;%构造平移矩阵 M1 Y1=M1*X;%计算平移结果 fill(Y1(1,:),Y1(2,:),red);%绘制平移后刚体 Y1 plot(X(1,:),X(2,:);%绘制原来刚体 X hold on axis equal R1=cos(pi/3),-sin(pi/3),0;sin(pi/3),cos(pi/3),0;0,0,1;%构造转动矩阵 R1 Y2=R1*X;%计算旋转结果 fill(Y2(1,:),Y2(2,:),blu

27、e);%绘制转动后刚体 Y2 fill(Y3(1,:),Y3(2,:),black);%绘制转动及平移后刚体 grid on hold off M2=1,0,20;0,1,30;0,0,1;%构造转动矩阵 M2 R2=cos(3*pi/4),-sin(3*pi/4),0;%构造转动矩阵 R2 sin(3*pi/4),cos(3*pi/4),0;0,0,1;Y3=M2*R2*X;%旋转平移后的结果 Y3第39页/共51页运行的结果如下第40页/共51页线性变换的线性变换的几何含义几何含义实例实例1 1已知向量已知向量 。请分析经过线性变换。请分析经过线性变换 后，向量后，向量 与向量与向量 的几

28、何关系。其中的几何关系。其中 分别为：分别为：第41页/共51页 在在MATLABMATLAB命令窗口运行程序命令窗口运行程序ma10.mma10.m，可以得到图形：，可以得到图形：第42页/共51页线性变换的线性变换的几何含义几何含义实例实例2 2 已知矩阵已知矩阵MATLABMATLAB分析特征向量的几何含义。分析特征向量的几何含义。，求它们的特征值和特征向量，并用，求它们的特征值和特征向量，并用第43页/共51页解：用解：用解：用解：用MATLABMATLABMATLABMATLAB求矩阵特征值和特征向量的方法为：求矩阵特征值和特征向量的方法为：求矩阵特征值和特征向量的方法为：求矩阵特征

29、值和特征向量的方法为：用用MATLABMATLAB演示矩阵演示矩阵A A的特征向量几何含义的命令为：的特征向量几何含义的命令为：（1 1）r=eig(A),r=eig(A),列向量列向量r r为矩阵为矩阵A A的特征值的特征值（2 2）V,D=eig(A)V,D=eig(A)，对角矩阵，对角矩阵D D的对角线元素为矩阵的对角线元素为矩阵A A的特的特征值，矩阵征值，矩阵V V的列向量为矩阵的列向量为矩阵A A的特征向量。的特征向量。eigshow(A)eigshow(A)第44页/共51页 在在MATLABMATLAB命令窗口运行程序命令窗口运行程序ma11.mma11.m，可以分别得到图形：

30、，可以分别得到图形：第45页/共51页二次型的二次型的几何含义几何含义1 1、利用正交变换化二次型为标准形的几何含义。、利用正交变换化二次型为标准形的几何含义。例例1 1 用正交变换，把下列二次型化为标准形，并讨论变用正交变换，把下列二次型化为标准形，并讨论变换前后所对应的二次曲线换前后所对应的二次曲线 及及 。解：用解：用MATLABMATLAB命令命令eigeig可以算出二次型矩阵的特征值分别可以算出二次型矩阵的特征值分别为：为：4.38204.3820，6.6180 6.6180 和和3.70163.7016，-2.7016-2.7016。第46页/共51页 在在MATLABMATLAB命令窗口运行程序命令窗口运行程序ma12.mma12.m，可以得到图形：，可以得到图形：第47页/共51页二次型的二次型的几何含义几何含义2 2、二次型正定、负定的几何含义。、二次型正定、负定的几何含义。例例2 2 分析下列二次型的正定性，并画出对应的二次曲面分析下列二次型的正定性，并画出对应的二次曲面 。第48页/共51页 在在MATLABMATLAB命令窗口运行程序命令窗口运行程序ma13.mma13.m，可以得到图形：，可以得到图形：第49页/共51页谢谢！谢谢！第50页/共51页