线性代数几何背景及应用 (2)幻灯片.ppt

线性代数几何背景及应用第1页，共48页，编辑于2022年，星期一一、方程及方程组的几何意义一、方程及方程组的几何意义二元一次方程二元一次方程在几何上表示的是一根直线，则两个二元一次方程组在几何上则表示两根直线的位置关系：相交=有惟一解 平行=无解 重合=无穷多解 例1 求解下列三个线性方程组 （a）（b）（c）第2页，共48页，编辑于2022年，星期一 用ezplot(s1),hold on,ezplot(s2),命令 可以解出结果如下图 其中s1和s2分别为方程的字符串表达式第3页，共48页，编辑于2022年，星期一 若有三个二元一次方程或更多个数的二元一次方程，代数上称之为“超定方程”，一般是不相容的和无解的，几何中平面上三根或更多根直线很难交于一点。例2 求解方程组 用图解法解例2 第4页，共48页，编辑于2022年，星期一 三元一次方程三元一次方程在几何上表示平面，从而两个三元一次方程构成的方程组表示两个平面的交线，三个三元一次方程构成的方程组两两联立求交线，得到两个二元一次方程，对于求得两根交线在xoy面上的投影。求得两根交线的交点即为方程组的解。若三个平面不重合且没有交线或交点，则表示该方程组无解。如下例。第5页，共48页，编辑于2022年，星期一例3 求解下列线性方程组，并画出三维图形来表示解的情况。（1）；（2）；（3）；（4）第6页，共48页，编辑于2022年，星期一利用MATLAB的M文件编辑器绘图可得：图3 三元非齐次线性方程组解的几何意义第7页，共48页，编辑于2022年，星期一 从图3中可以看出：方程组（1）的解为三个平面的交点，故该方程组有 唯一解；方程组（2）的三个平面刚好相交于同一条直线，这个齐次线性方程组有无穷多解，即解空间是一维的。方程组（3）的三个平面没有共同的交点。即方程组无解。方程组（4）也无解。第8页，共48页，编辑于2022年，星期一 推广之后，更多元的线性代数方程组，则表示更高维空间内的方程组，虽然很难想象直观的几何图形，但关于方程的基本概念是一脉相承的，涉及到计算就是从几何概念过渡到代数概念。如：阶数、维数等概念。第9页，共48页，编辑于2022年，星期一二、行列式的几何意义二、行列式的几何意义二维 已知向量 由向量 和 所构成的 平行四边形的面积为 行列式 的绝对值 第10页，共48页，编辑于2022年，星期一三维 已知三个向量 由这三个向量所构成的平行六面体的体积即为 三阶行列式 的绝对值 如图 第11页，共48页，编辑于2022年，星期一三、平面上线性变换的几何意义三、平面上线性变换的几何意义例3 已知向量 ，矩阵 ，。请分析经过线性变换 后，向量 与原向量 的几何关系 。第12页，共48页，编辑于2022年，星期一绘制图形如下图所示：图4 线性变换的几何意义第13页，共48页，编辑于2022年，星期一 从图从图4中可以看出：中可以看出：矩阵 对 进行线性变换的结果 为向量 的竖直轴对称向量；矩阵 对 进行线性变换的结果 为向量 的水平轴对称向量；矩阵 对 进行线性变换的结果 为把向量 的横坐标乘以0.5，把 的纵坐标乘以2得到的向量；矩阵 对 进行线性变换的结果 为把向量 按顺时针方向旋转 所得到的向量。第14页，共48页，编辑于2022年，星期一 例例4：设x为二维平面上第一象限中的一个单位方块，其四个顶点的数据可写成 把不同的A矩阵作用于此组数据，可以得到多种多样的结果yi=Ai*x。用程序ag911进行变换计算，并画出x及yi图形：x0,1,1,0;0,0,1,1;subplot(2,3,1),fill(x(1,:),0,x(2,:),0,r)A11,0;0,1,y1A1*xsubplot(2,3,2),fill(y1(1,:),0,y1(2,:),0,g)第15页，共48页，编辑于2022年，星期一绘制几何图形可得：第16页，共48页，编辑于2022年，星期一 使用MATLAB时，行列式用Di=det(Ai)求得，特征值和特征向量则用pi,lamdai=eig(Ai)计算，算得的结果如下：第17页，共48页，编辑于2022年，星期一关于笔算与机算的结合关于笔算与机算的结合 矩阵的赋值和其加、减、乘、除（求逆）命令；矩阵化为最简行阶梯型的计算命令；U0,ip=rref(A)多元线性方程组MATLAB求解的几种方法；x=inv(A)*b,U=rref(A)行列式的几种计算机求解方法；D=det(A),L,U=lu(A);D=prod(diag(U)n个m维向量组的相关性及其秩的计算方法和命令;r=rank(A),U=rref(A)求线性方程组的基础解系及方程解的MATLAB命令；xb=null(A)矩阵的特征方程、特征根和特征向量的计算命令；f=poly(A);P,D=eig(A)化二次型为标准型的MATLAB命令；yTDy=xTAx;其中y=P-1x,第18页，共48页，编辑于2022年，星期一解高阶线性方程组的方法解右列方程组AX=b可有多种方法,如(1)X=Ab(2)化为行最简型A=3,-4,2,2,-1;0,-6,0,-3,-3;4,-3,4,3,-2;1,2,1,0,-5;-2,6,-2,1,3b=2;-3;2;-2;1;X=inv(A)*b,pauseC=A,b,Uc,ip=rref(C)第19页，共48页，编辑于2022年，星期一应用一：线性方程组与矩阵应用一：线性方程组与矩阵1.1 插值多项式例例1 给定t-y平面上的三个点(1,2),(2,3)和(3,6)，求过这三点的二次多项式函数：解：解：本题归结为求a,b,c三个系数，使它们满 足下列各方程第20页，共48页，编辑于2022年，星期一 这是典型的三元线性方程组，用Matlab时，键入：B=1,1,1,2;1,2,4,3;1,3,9,6;x=rref(B)得到x=1 0 0 3 0 1 0 -2 0 0 1 1 第21页，共48页，编辑于2022年，星期一 x矩阵的最后一列即为a，b，c的值，则待求二次多项式为：例例2 下表给出函数 上4个点的值，试求三次插值多项式 ，并求 的近似值。ti0123f(ti)30-16第22页，共48页，编辑于2022年，星期一解：解：令三次多项式函数过表中已知的4点，可以得到四元线性方程组：应该用计算机求解，键入：A=1,0,0,0;1,1,1,1;1,2,4,8;1,3,9,27,b=3;0;-1;6,s=rref(A,b)第23页，共48页，编辑于2022年，星期一得到x =1 0 0 0 3 0 1 0 0 -2 0 0 1 0 -2 0 0 0 1 1得到 ，三次多项函数为 ,故 近似等于第24页，共48页，编辑于2022年，星期一1.2 平板稳态温度的计算平板稳态温度的计算 在钢板热传导的研究中，常常用节点温度来描述钢板温度的分布。例例3 假设图1中钢板已经达到稳态温度分布，上下、左右四个边界的温度值如图所示，而表示钢板内部四个节点的温度。若忽略垂直于该截面方向的热交换，那么内部某节点的温度值可以近似地等于与它相邻四个节点温度的算术平均值，如。请计算该钢板的温度分布。第25页，共48页，编辑于2022年，星期一图1 平板的温度分布解：根据已知条件可以得到以下线性方程组：化简为标准的矩阵形式如下：第26页，共48页，编辑于2022年，星期一 在MATLAB命令窗口输入：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)第27页，共48页，编辑于2022年，星期一结果为：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中，把钢板内部分成了22个节点，本例把钢板内部分为55个节点，如图2所示。求钢板的稳态温度分布，并绘制温度分布图形。第28页，共48页，编辑于2022年，星期一 钢板的温度分布如图3所示。其中x、y坐标分别表示钢板横、纵方向的节点数，高度表示节点的温度值，该三维图形形象地反映了钢板的温度分布。图3 钢板的温度分布第29页，共48页，编辑于2022年，星期一 1.3 交通流量的分析交通流量的分析例例4 某城市有两组单行道，构成了一个包含四个节点A,B,C,D的十字路口，如图2所示。汽车进出十字路口的流量（每小时的车流数）标于图上。现要求计算每两个节点之间路段上的交通流量 。（假设，针对每个节点，进入和离开的车数相等）图2 单行道4节点交通流图第30页，共48页，编辑于2022年，星期一解：根据已知条件可以得到，四个节点的流 通方程为 节点A：节点B：节点C：节点D：将以上方程组进行整理，得第31页，共48页，编辑于2022年，星期一Matlab程序ea110为 A=1,-1,0,0;0,1,-1,0;0,0,1,-1;-1,0,0,1 b=160;-40;210;-330 U0=rref(A,b)可以得出其最简行阶梯形矩阵 由于U0的最后一行为全零，也就是说，四个方程中实际上只有三个独立方程，所以该方程组为欠定方程，存在无穷多组解。第32页，共48页，编辑于2022年，星期一若以 为自由变量，方程组的解可以表示为：如果有一些车围绕十字路的矩形区反时针绕行，流量 。都会增加，但并不影 响出入十字路的流量。这就是方程组有无穷多解的原因。第33页，共48页，编辑于2022年，星期一人口迁徙问题 例5 假设一个城市的总人口数是固定不变，但人口的分布情况变化如下：每年都有5的市区居民搬到郊区；而有15的郊区居民搬到市区。若开始有700000人口居住在市区，300000人口居住在郊区。请利用分析：（1）10年后市区和郊区的人口各是多少？（2）30年后、50年后市区和郊区的人口各是多少？（3）分析（2）中数据相似的原因。第34页，共48页，编辑于2022年，星期一解：这个问题可以用矩阵乘法来描述。令人口变量 其中 为市区人口所占比例，为郊区人口所占比例。在n+1年的人口分布状态为：用矩阵乘法可写成：第35页，共48页，编辑于2022年，星期一开始市区和郊区的人口数为 可以得到n年后市区和郊区的人口分布：因此10年后的人口可用程序计算如下：A=0.95,0.15;0.05,0.85;X0=700000;300000;X10=A10*X0 程序运行的结果为：市区和郊区人口数约为：744630和255370。第36页，共48页，编辑于2022年，星期一 无限增加时间n，市区和郊区人口之比将趋向一组常数0.25/0.75。为了弄清为什么它趋向于一个稳态值，可以将A对角化。令 ，其中为对角矩阵，则有对角矩阵的幂次可以化为元素的幂次 所以，它就很容易计算。第37页，共48页，编辑于2022年，星期一程序la24%分析n年后城市人口分布A=0.95,0.15;0.05,0.85;X0=700000;300000;P,lamda=eig(A);syms n%定义符号变量nXn=P*lamda.n*inv(P)*X0%.n对矩阵lamda中所有元素进行幂运算计算结果为：随n增大后一项(4/5)n趋近于零。第38页，共48页，编辑于2022年，星期一多项式插值与拟合 例6 下表给出了平面坐标系中五个点的坐标。（1）请过这五个点作一个四次多项式函数，并求x=5时的函数值。用MATLAB绘制多项式函数曲线、通过已知点及插值点。（2）请根据这五个点，拟合一个二次多项式函数，并用MATLAB绘制多项式函数曲线及已知的五个点。x01234y-270210-75第39页，共48页，编辑于2022年，星期一解：（1）根据已知条件，把五个点的坐标值分别代入四次多项式函数，可以得到如下线性方程组：其中矩阵：第40页，共48页，编辑于2022年，星期一系数矩阵A的行列式为范德蒙(Vandermonde)行列式，且五个坐标点的横坐标各不相同，则该行列式不等于零，所以方程组有唯一解。写出程序：x=0;1;2;3;4;%输入已知点坐标y=-27;0;21;0;-75;%构造范德蒙矩阵，也可用内置的vander函数A=x.0,x.1,x.2,x.3,x.4;a=Ay;%得到适定方程组的唯一解a，运行程序，得到a(1)=-27,a(2)=12,a(3)=26,a(4)=-12,a(5)=1第41页，共48页，编辑于2022年，星期一多项式拟合要解一个超定方程把五个点的坐标值分别代入二次多项式函数，可以得到如下线性方程组：其中，第42页，共48页，编辑于2022年，星期一该方程组有三个未知数，但有五个方程，进一步分析可以得到该方程组无解，即不存在一个二次多项式曲线刚好能过已知的五个点。MATLAB软件提供了一个利用最小二乘法解决超定方程组解的方法。求系数的公式也是a=Ay，以找到一条二次曲线来近似地描述已知5点的变化情况。对比插值和拟合的曲线如下图第43页，共48页，编辑于2022年，星期一刚体的平面运动 例7 用平面坐标系中的一个闭合图形来描述刚体，用一个矩阵X来表示它。X的一列表示刚体一个顶点的坐标。为了使图形闭合，X的最后一列和第一列相同；为了实现刚体的平移运算，给矩阵X添加元素值都为1的一行，使矩阵X的形状为3n。若有矩阵：可以证明，矩阵是刚体X沿x轴正方向平移，沿y轴正方向平移后的结果；矩阵是刚体X以坐标原点为中心逆时针转动t弧度的结果。第44页，共48页，编辑于2022年，星期一计算实例设刚体图形的各个顶点坐标如下：构造刚体矩阵X，平移矩阵M及转动矩阵R。x04610853.56.1 6.53.220y014140011664.54.500第45页，共48页，编辑于2022年，星期一X=0,4,6,10,8,5,3.5,6.1,6.5,3.2,2,0;0,14,14,0,0,11,6,6,4.5,4.5,0,0;ones(1,12);%构造刚体矩阵M=1,0,-30;0,1,15;0,0,1;%构造平移矩阵Y1=M*X;%计算平移结果plot(X(1,:),X(2,:);%绘制原来刚体fill(Y1(1,:),Y1(2,:),red);%绘制平移后刚体R=cos(pi/3),-sin(pi/3),0;sin(pi/3),cos(pi/3),0;0,0,1;%构造转动矩阵Y2=R*X;fill(Y2(1,:),Y2(2,:),blue);%绘制转动后刚体Y3=M*R*X;fill(Y3(1,:),Y3(2,:),black);%绘制转动及平移后刚体程序la26的核心语句第46页，共48页，编辑于2022年，星期一程序la26运行的结果如下第47页，共48页，编辑于2022年，星期一请多提宝贵意见！请多提宝贵意见！第48页，共48页，编辑于2022年，星期一