线性代数 河北工业吹捏学.pptVIP

线性代数MATLAB计算应用;线性方程组求解; 例1. 求解齐次线性方程组的通解。(m1.m);Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd.;2．非齐次线性方程组求解; 例2. 求解非齐次线性方程组的通解。(m2.m);(3) 生成增广矩阵的行最简形R，并将增广矩阵的行最简形中基准元素所在的列号存入向量s中 [R,s]=rref([A,b]);;(4) 对解进行判定并求解 if r1~=r2; disp(’方程组无解’） else if r2==n disp(’方程组有唯一解’） disp (’方程组唯一解为’） x=A\b; else disp (’方程组有无穷解’） disp (’方程组特解为’） x1=A\b disp (’方程组导出组基础解系为’） x0=null(A,’r’) end end ;计算结果显示如下;3．线性方程组求解的几何概念; 例3. 用图形解描述齐次线性方程组的解。 (m3.m);;对上述方程组求解可以发现每个方程都把点限制在某一平面上，但是公共解被限定在一条直线上，显示空面内满足这个方程组的解不唯一，公共点没有被唯一地限位，方程组有无数解。(m4.m);求解发现方程组没有公共解，显示空间三个平面没有公共交点，这种情形方程组无解。(m5.m) ;例4 说明一个桌子分别具有二、三、四条腿的稳定性。; 添加两个桌腿(x1,y1,z1)， (x2,y2,z2)，给平面方程加两个约束;添加三个桌腿(x1,y1,z1)， (x2,y2,z2)， (x3,y3,z3)，给平面方程施加三个约束;四个桌腿，方程组需满足;向量组的极大无关组求法有三种： （1）根据最高阶非零子式求解； （2）初等行变换求行向量组的极大无关组； （3）初等行变换求列向量组的极大无关组;例3：向量组α1=[1，1，0，2，2]T, α2=[3，4，0，8，3]T, α3=[2，3，0，6，1]T, α4=[9，3，2，1，2]T, α5=[6，-2，2，-9，2]T。求出其极大无关组。(m6.m);(3)输入极大无关组 fprintf(’极大无关组为：’) for i=1:r fprintf(‘a%d’ ,s(i)) end for i=1:r A0(:,r)=A(:,s(i)); end A0 ;计算结果显示，极大无关组是 α1, α2，α4，该极大无关组被放入矩阵A0,剩余向量α3, α5被极大无关组线性表示如下;具体计算过程如下：;将其???向量用极大无关组线性表示，求出坐标矩阵K;本节内容涉及的MATLAB 命令;例4：向量组a1=[2，-1，-1，1]T; a2=[-1,2,1,-1] T; a3=[-1,1,2,-1] T;a4=[1,-1,-1,2] T; A=[a1,a2,a3,a4];。 求： ;解：（1）输入向量组，并将其放入矩阵A ;将A中列向量组用P的列向量线性表示（求坐标）;（2）求出A特征值，特征向量，并将特征向量标准正交化；;需要掌握的MATLAB语句;平面运动的概念;研究刚体在平面上投影这个平面图形的运动;平面图形运动分解：;平面图形上的点M(x,y),平移后到达点M1(x1,y1);点M1(x1,y1) 转动φ角后到达点M2(x2,y2);点M(x,y),经过平面运动后到达点M2(x2,y2)的运算关系可以表达为;例5：直角三角形ABC,三点坐标分别为A(0,0),B(5,0),  C(0,8),完成以下图形移动：;（1）向上移动10，向右移动20；（平移）;（3）先顺时针转动π/2，再向下移动15， 向左移动10。 （平面运动）;程序(m9.m);运行结果;常微分方程求解问题; 定义向量;描述(0~10s) 中点的运动轨迹;循环计算可以得到各时间点的位置，速度和加速度;运行结果;The End