Matlab与线性代数 一.docVIP

Matlab与线性代数 一 Matlab与线性代数(一)2010-06-07 06：30一、线性代数基本方程组 基本方程组： 矩阵表示： 解决问题的视角： 1、解联立方程的视角(行阶梯变换矩阵运算) 着重研究解x，即研究线性方程组的解法。中学里的解方程和MATLAB的矩阵除法就是这样。 要点：矩阵的每一行代表一个方程，m行代表m个线性联立方程。n列代表n个变量。如果m是独立方程数，根据m n、m=n、m n确定方程是欠定、适定还是超定。对这三种情况都会求解了，研究就完成了。必须剔除非独立方程。行阶梯形式、行列式和秩的概念很大程度上为此目的而建立。 2、向量空间中向量的合成的视角(用向量空间解方程组) 把A各列看成n个m维基本向量，线性方程组看成基向量的线性合成： 要点：解x是这些基向量的系数。它可能是常数(适定方程)，也可能成为其中的一个子空间(欠定方程)。要建立其几何概念，并会求解或解空间。 3、线性变换或映射的视角(线性变换及其特征) 把b看成变量y，着重研究把Rn空间的x变换为Rm空间y的效果，就是研究线性变换系数矩阵A的特征对变换的影响。 要点：就是要找到适当的变换，使研究问题的物理意义最为明晰。特征值问题就是一例。 二、线性代数建模与应用概述介绍一些大的系统工程中使用线性代数的情况，使读者知道为什么线性代数在近几十年来变得如此的重要。Leontief教授把美国的经济用500个变量的500个线性方程来描述，在1949年利用当时的计算机解出了42×42的简化模型，使他于1973年获得诺贝尔经济奖，从而大大推动了线性代数的发展。把飞行器的外形分成若干大的部件，每个部件沿着其表面又用三维的细网格划分出许多立方体，这些立方体包括了机身表面以及此表面内外的空气。对每个立方体列写出空气动力学方程，其中包括了与它相邻的立方体的共同边界变量，这些方程通常都已经简化为线性方程。对一个飞行器，小立方体的数目可以多达400,000个，而要解的联立方程可能多达2,000,000个。飞行器的运动要用三个转动和三个平移共六个变量来表示(像在第九章中介绍的那样)。除此以外，为了控制飞行器的三维转动，需要三个控制面，即方向舵(垂直尾翼)、升降舵(水平尾翼)和副翼，它们的偏转角又构成了三个变量。描述飞行器的运动就需要12个变量的组合。而且这已经不单是代数方程，而是微分方程了。卫星上用三种可见光和四种红外光进行摄像，对每一个区域，可以获得七张遥感图象。利用多通道的遥感图可以获取尽可能多的地面信息，因为各种地貌、作物和气象特征可能对不同波段的光敏感。而在实用上应该寻找每一个地方的主因素，成为一张实用的图象。每一个象素上有七个数据，形成一个多元的变量数组，在其中合成并求取主因素的问题，就与线性代数中要讨论的特征值问题有关。在全国设立几十万个观察点，把每一点的经度、纬度和高度三个坐标建立起来。现在对于经度纬度的测量精度要求，已经提高到了若干厘米，对于高度的精度要求更高。在一些边远地区，对于一些特征点的测量，要耗费很大的人力物力。例如对珠穆朗玛峰顶高度的测量，要经过多种方法，取得多种数据，并且用最小二乘法进行误差的处理。 三、行阶梯法解线性方程 1、线性方程的Matlab表示方法 (1)由n个变量组成的m个联立线性代数方程组： 用MATLAB语言表示为： 其中： n是未知数个数，m是独立方程的个数。当m n时，方程组有无数多个解，称为欠定方程；当m n时，方程组无解，称为超定方程；当m=n时，方程组有唯一解，称为适定方程；所以不能简单地看形式上的m和n，还必须剔除其中非独立方程的虚假成分。将要讨论的行阶梯形式、行列式和秩等概念，很大程度上就是为了找到独立方程的数目。 (2)解线性方程组及在matlab中的显示 作图的matlab代码： subplot(2,2,1)%作左上角图，画方程组(a)ezplot(x1-2*x2+1)%画(a)中的x1-2*x2+1=0 hold on ezplot(-x1+3*x2-3)subplot(2,2,2)%作右上角图，画方程组(b)ezplot(x1-2*x2+1)hold on ezplot(-x1+2*x2-3)subplot(2,2,3)%作左下角图，画方程组(c)ezplot(x1-2*x2+1)hold on ezplot(-x1+2*x2-1)subplot(2,2,4)%作右下角图，画方程组(d)ezplot(x1+x2-1)hold on ezplot(x1-x2-3)hold on ezplot(-x1+2*x2+3)使用参考：ezplot 四个方程组在matlab中的作图： 2、初等行变换 将系数矩阵A和B组成增广矩阵： 对增广矩阵的行作以下三种运算不会改变方程组的解，这三种运算组成了矩阵的初等行变换：