用MATLAB软件解线性方程组.pptxVIP

第1页/共9页3.1矩阵函数 MATLAB提供了许多矩阵函数. 正是因为拥有了为数众多的、完善的矩阵函数，才使得MATLAB具有了强大的功能。 第一页，共10页。第2页/共9页在命令窗口运行帮助命令 ：help elmat ,可以列举出大量的矩阵函数.下面是几个常用的矩阵函数:det计算矩阵的行列式的值inv求矩阵的逆阵rank 求矩阵的秩[V D]=eig(A) 求矩阵A的特征值和特征向量poly 求矩阵的特征多项式rref 用初等变换将矩阵化成行阶梯形null(A,’r’) 给出齐次线性方程组Ax=0 的基础解系fliplr 矩阵左右翻转flipud 矩阵上下翻转trace 求矩阵的迹diag 取得矩阵对角线元素第二页，共10页。第3页/共9页3.2矩阵函数的应用例3–1 设矩阵求A的行列式、秩和逆矩阵。解：A=[3 -4 0; -1 5 2; 4 1 -6]det (A) %求矩阵的行列式的值rank (A) %求矩阵的秩inv (A) %求逆矩阵第三页，共10页。第4页/共9页3.3 解线性方程组3.3.1求线性方程组的唯一解例3–2求线性方程组Ax=B的解，其中：解法1利用矩阵除法：X=A\B 解法2利用求逆矩阵函数 inv：X1=inv(A)*B 比较：解法1比解法2更简便， 解法1 的算法优于解法2 ， 解法1可用于一般矩阵，而解法2只能用于非奇异的方阵 因此，只需运用解法1 .第四页，共10页。第5页/共9页3.3.2求线性方程组的通解例3–3 设求线性方程组Ax=B的通解。第五页，共10页。第6页/共9页解法1： 利用除法 \ 和 null 函数在命令窗口输入以下命令：(注意：这里给出的 A不 是方阵）A=[1 1 -1 -1;2 -5 3 2;7 -7 3 1];B=[5; -4; 7];format ratx1=A\B %求得非齐次方程组Ax=B的一个特解x1Y=null(A,’r’) %求得齐次方程组Ax=0 的基础解系Y输出结果：x1 = 3 2 0 0 Y = 2/7 3/7 5/7 4/7 1 0 0 1 则方程组Ax=B的通解为: x=x1+k1*Y(:,1)+k2*Y(:,2)第六页，共10页。第7页/共9页解法2：利用 rref 函数 在命令窗口输入以下命令：format ratA=[1 1 -1 -1;2 -5 3 2;7 -7 3 1];B=[5; -4; 7];%用初等行变换将增广矩阵 [A B] 化成最简行阶梯形TT=rref([A B])T = 1 0 -2/7 -3/7 3 0 1 -5/7 -4/7 2 0 0 0 0 0 输出结果：于是可得方程组Ax=B的通解为:第七页，共10页。第8页/共9页小 结1.矩阵函数 利用矩阵函数求矩阵的行列式、秩和逆 求方程组的唯一解4. 求方程组的通解第八页，共10页。第9页/共9页感谢您的观看！第九页，共10页。内容总结3.1 矩阵函数。第1页/共9页。第2页/共9页。第3页/共9页。比较：解法1比解法2更简便，。解法1可用于一般矩阵，而解法2只能用于非奇异的方阵。在命令窗口输入以下命令：。第7页/共9页。利用矩阵函数求矩阵的行列式、秩和逆。第8页/共9页。感谢您的观看第十页，共10页。词语：函数 拼音：hánshù 解释： 词语：矩阵 拼音：jǔ zhèn 解释：将m×n个元素排成m行、n列的矩形表，用圆括号(或方括号、双直线)括起来，即a11a12…a1na21a22…a2n……………am1am2…amn称为m行n列矩阵，简称m×n矩阵。当m=n时，称为n阶方阵。对矩阵可按一定法则进行加法、乘法及数与矩阵相乘等运算。 词语：矩阵 拼音：jǔ zhèn 解释：将m×n个元素排成m行、n列的矩形表，用圆括号(或方括号、双直线)括起来，即a11a12…a1na21a22…a2n……………am1am2…amn称为m行n列矩阵，简称m×n矩阵。当m=n时，称为n阶方阵。对矩阵可按一定法则进行加法、乘法及数与矩阵相乘等运算。 词语：解法 拼音：jiě fǎ 解释：解救的方法。元王晔《桃花女》第三折：“你有这许多算法，他可有许多的解法哩。”《平妖传》第十九回：“若有急难时，可解开册子来看便有解法。”鲁迅《南腔北调集·经验》：“我们一向喜欢恭维古圣人，