线性代数几何背景及应用.ppt

线性代数中的几何背景线性代数中的几何背景一、方程及方程组的几何意义一、方程及方程组的几何意义二元一次方程二元一次方程在几何上表示的是一根直线，则两个二元一次方程组在几何上则表示两根直线的位置关系： 相交====〉〉有惟一解 平行====〉〉无解 重合====〉〉无穷多解 例1 求解下列三个线性方程组 （a） （b） （c） 用ezplot(s1),hold on, ezplot(s2),命令 可以解出结果如下图 其中s1和s2分别为方程的字符串表达式 若有三个二元一次方程或更多个数的二元一次方程，代数上称之为 “超定方程”，一般是不相容的和无解的，几何中平面上三根或更多根直线很难交于一点 例2 求解方程组 用图解法解例2  三元一次方程三元一次方程在几何上表示平面，从而两个三元一次方程构成的方程组表示两个平面的交线，三个三元一次方程构成的方程组两两联立求交线，得到两个二元一次方程，对于求得两根交线在xoy面上的投影。

求得两根交线的交点即为方程组的解若三个平面不重合且没有交线或交点，则表示该方程组无解 例3 求解下列线性方程组，并画出三维图形来表示解的情况1） ； （2） ；（3） ； （4）利用MATLAB的M文件编辑器绘图可得： 图3 三元非齐次线性方程组解的几何意义 从图3中可以看出： 方程组（1）的解为三个平面的交点，故该方程组有 唯一解； 方程组（2）的三个平面刚好相交于同一条直线，这个齐次线性方程组有无穷多解，即解空间是一维的 方程组（3）的三个平面没有共同的交点即方程组无解 方程组（4）也无解  推广之后，更多元的线性代数方程组，则表示更高维空间内的方程组，虽然很难想象直观的几何图形，但关于方程的基本概念是一脉相承的，涉及到计算就是从几何概念过渡到代数概念如：阶数、维数等概念二、行列式的几何意义二、行列式的几何意义•二维 已知向量 由向量 和 所构成的 平行四边形的面积为 行列式 的绝对值 •三维 已知三个向量 由这三个向量所构成的平行六面体的体积即为 三阶行列式 的绝对值 如图 三、平面上线性变换的几何意义三、平面上线性变换的几何意义例3 已知向量 ， 矩阵 ， ， ， ， 。

请分析经过线性变换 后，向量 与原向量 的几何关系 •绘制图形如下图所示： 图4 线性变换的几何意义 从图从图4中可以看出：中可以看出：矩阵 对 进行线性变换的结果 为向量 的竖直轴对称向量；矩阵 对 进行线性变换的结果 为向量 的水平轴对称向量；矩阵 对 进行线性变换的结果 为把向量 的横坐标乘以0.5，把 的纵坐标乘以2得到的向量；矩阵 对 进行线性变换的结果 为把向量 按顺时针方向旋转 所得到的向量 例例4：：设x为二维平面上第一象限中的一个单位方块，其四个顶点的数据可写成 把不同的A矩阵作用于此组数据，可以得到多种多样的结果yi=Ai*x用程序ag911进行变换计算，并画出x及yi图形： x[0,1,1,0;0,0,1,1]; subplot(2,3,1), fill([x(1,:),0],[x(2,:),0],'r')A1[1,0;0,1], y1A1*xsubplot(2,3,2), fill([y1(1,:),0],[y1(2,:),0],'g') 绘制几何图形可得： 使用MATLAB时，行列式用Di=det(Ai)求得，特征值和特征向量则用[pi,lamdai]=eig(Ai)计算，算得的结果如下：关于笔算与机算的结合关于笔算与机算的结合 ① 矩阵的赋值和其加、减、乘、除（求逆）命令；② 矩阵化为最简行阶梯型的计算命令；[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, 解高阶线性方程组的方法•解右列方程组•AX=b•可有多种方法,如•(1) X=A\b•(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,3]•b = [2; -3; 2; -2; 1]; •X=inv(A)*b, pause•C=[A,b], [Uc,ip]=rref(C)应用一：线性方程组与矩阵应用一：线性方程组与矩阵1.1 插值多项式例例1 给定t-y平面上的三个点(1,2),(2,3)和(3,6)，求过这三点的二次多项式函数：解：解：本题归结为求a,b,c三个系数，使它们满 足下列各方程 这是典型的三元线性方程组，用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  x矩阵的最后一列即为a，b，c的值，则待求二次多项式为： 例例2 下表给出函数 上4个点的值，试求三次插值多项式 ，并求 的近似值。

ti0123f(ti)30-16解：解：令三次多项式函数过表中已知的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]) 得到x = 1 0 0 0 3 0 1 0 0 -2 0 0 1 0 -2 0 0 0 1 1得到 ，三次多项函数为 , 故 近似等于1.2 平板稳态温度的计算平板稳态温度的计算 在钢板热传导的研究中，常常用节点温度来描述钢板温度的分布例例3 假设图1中钢板已经达到稳态温度分布，上下、左右四个边界的温度值如图所示，而表示钢板内部四个节点的温度。

若忽略垂直于该截面方向的热交换，那么内部某节点的温度值可以近似地等于与它相邻四个节点温度的算术平均值，如请计算该钢板的温度分布图1 平板的温度分布解：根据已知条件可以得到以下线性方程组：化简为标准的矩阵形式如下： 在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])结果为：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中，把钢板内部分成了2×2个节点，本例把钢板内部分为5×5个节点，如图2所示。

求钢板的稳态温度分布，并绘制温度分布图形 钢板的温度分布如图3所示其中x、y坐标分别表示钢板横、纵方向的节点数，高度表示节点的温度值，该三维图形形象地反映了钢板的温度分布图3 钢板的温度分布 1.3 交通流量的分析交通流量的分析例例4 某城市有两组单行道，构成了一个包含四个节点A,B,C,D的十字路口，如图2所示汽车进出十字路口的流量（每小时的车流数）标于图上现要求计算每两个节点之间路段上的交通流量 假设，针对每个节点，进入和离开的车数相等）图2 单行道4节点交通流图解：根据已知条件可以得到，四个节点的流 通方程为 节点A： 节点B： 节点C： 节点D： 将以上方程组进行整理，得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的最后一行为全零，也就是说，四个方程中实际上只有三个独立方程，所以该方程组为欠定方程，存在无穷多组解。

若以 为自由变量，方程组的解可以表示为： 如果有一些车围绕十字路的矩形区反时针绕行，流量 都会增加，但并不影 响出入十字路的流量这就是方程组有无穷多解的原因 人口迁徙问题 例5 假设一个城市的总人口数是固定不变，但人口的分布情况变化如下：每年都有5％的市区居民搬到郊区；而有15％的郊区居民搬到市区若开始有700000人口居住在市区，300000人口居住在郊区请利用分析：（1）10年后市区和郊区的人口各是多少？（2）30年后、50年后市区和郊区的人口各是多少？（3）分析（2）中数据相似的原因解：这个问题可以用矩阵乘法来描述令人口变量 其中 为市区人口所占比例， 为郊区人口所占比例在n+1年的人口分布状态为：用矩阵乘法可写成： •开始市区和郊区的人口数为 可以得到n年后市区和郊区的人口分布：因此10年后的人口可用程序计算如下：A=[0.95,0.15;0.05,0.85];X0=[700000;300000];X10=A^10*X0 程序运行的结果为：市区和郊区人口数约为：744630和255370  无限增加时间n，市区和郊区人口之比将趋向一组常数0.25/0.75。

为了弄清为什么它趋向于一个稳态值，可以将A对角化令 ，其中Λ为对角矩阵，则有对角矩阵的幂次可以化为元素的幂次 所以，它就很容易计算 程序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趋近于零多项式插值与拟合 例6 下表给出了平面坐标系中五个点的坐标•（1）请过这五个点作一个四次多项式函数，•并求x=5时的函数值用MATLAB绘制多项式函数曲线、通过已知点及插值点•（2）请根据这五个点，拟合一个二次多项式函数，•并用MATLAB绘制多项式函数曲线及已知的五个点x01234y-270210-75解：（1）根据已知条件，把五个点的坐标值分别代入四次多项式函数，可以得到如下线性方程组：•其中矩阵：•系数矩阵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=A\y; % 得到适定方程组的唯一解a，•运行程序，得到•a(1)=-27,a(2)=12,a(3)=26,a(4)=-12,a(5)= 1多项式拟合要解一个超定方程•把五个点的坐标值分别代入二次多项式函数，可以得到如下线性方程组：•其中，•该方程组有三个未知数，但有五个方程，进一步分析可以得到该方程组无解，即不存在一个二次多项式曲线刚好能过已知的五个点MATLAB软件提供了一个利用最小二乘法解决超定方程组解的方法求系数的公式也是a=A\y，以找到一条二次曲线来近似地描述已知5点的变化情况•对比插值和拟合的曲线如下图刚体的平面运动 例7 用平面坐标系中的一个闭合图形来描述刚体，用一个矩阵X来表示它X的一列表示刚体一个顶点的坐标为了使图形闭合，X的最后一列和第一列相同；为了实现刚体的平移运算，给矩阵X添加元素值都为1的一行，使矩阵X的形状为3×n。

•若有矩阵：•可以证明，矩阵是刚体X沿x轴正方向平移，沿y轴正方向平移后的结果；矩阵是刚体X以坐标原点为中心逆时针转动t弧度的结果计算实例•设刚体图形的各个顶点坐标如下：•构造刚体矩阵X，平移矩阵M及转动矩阵R• x04610853.56.1 6.53.220y014140011664.54.500•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的核心语句程序la26运行的结果如下请多提宝贵意见！请多提宝贵意见！。