应用MATLAB解决高等代数问题.pptxVIP

第四讲：应用MATLAB解决高等代数问题矩阵的初等变换与方程的MATLAB求解交换矩阵中的两个行向量的位置；用一个非零数乘以矩阵的某一行向量把矩阵的某一个行向量乘以实数并加到矩阵的另一行上

例：求解线性方程组线性代数方法用增广矩阵初等变换即消元法过程

经过初等行变换将矩阵01变为矩阵02这时矩阵对应的方程组03此方程组的解为04

A=[3-153;1-121;1-2-12]%输入矩阵的数据01A([13],:)=A([31],:)%交换第一行和第三行数据02A(2,:)=A(2,:)-A(1,:)%将第一行乘以-1加到第二行03A(3,:)=A(3,:)-3*A(1,:)%将第一行乘以-3加到第三行04A(3,:)=A(3,:)-5*A(2,:)%将第二行乘以-5加到第三行05方法之一：初等变换法

A=[3-153;1-121;1-2-12]%输入矩阵的数据formatrat%分数数据格式rref(A)%化简矩阵方法之二：Cramer法则A=[3-15;1-12;1-2-1]%输入矩阵的数据B=[312]‘;%输入线性方程组的常数项S=[000]‘;%给解向量赋初值fori=1:3%for循环

C=A;%将矩阵A赋给临时矩阵CC(:,i)=B;%将常数项赋给矩阵C的第i列即求AiS(i)=det(C)/det(A);%求xiendformatrat%数据格式说明为分数形式S%显示S

方法之三：利用矩阵的左除“\”A=[3-15;1-12;1-2-1]；b=[312]‘；x=A\bx=10/7-1/7-2/7二、线性方程组的解结构1。齐次方程组的解结构AX=0求解方程过程如下

根据最简行阶梯形矩阵写出简化方程组确定自由求知量整理方程组为向量形式量提取方程组右端各自由求知量的系数形成的向量组即为基础解系将系数矩阵化为最简行阶梯形矩阵

例：解线性方程组：01A=[1–11–1;1–1–11;1–1–22]%输入矩阵rref(A)%将矩阵化为最简阶梯形矩阵03应用MATLAB计算过程如下：02

-11-1-1-11-1-22A=-10001-1000ans=运行结果为：

取x2，x4为自由求知量，得方程组的解的向量形式为由运行结果知化简的等价方程组为：

所以基础解系为：所以齐次方程组的通解为

2.非齐次方程组的解的结构求解非齐次线性方程组的通解的步骤如下：写出非次方程组的增广矩阵；将增广矩阵化为最简行阶梯形矩阵；观察增广矩阵与系数矩阵的秩是否相等，若相等则方程组有解，若不相等则方程组无解；写出对应的简化的线性方程组；确定自由求知量整理方程组为向量形式。

例：求解下列非齐线性方程组

A=[1231;1452;2983;3772];b=[3;2;7;12]；01计算结果如下03formatratc=[Ab]；rref(c);02在MATLAB中输入的命令如下

ans=100-1/231/601002/30011/2-7/600000所以简化方程组为：

取x4为自由求知量所以原线性方程组的通解为：

010203040506向量组线性相关与线性无关的定义：01如果存在m个不全为零的一组数k1,k2,…,km使02成立，则称向量组03是线性相关的。04如果仅当k1=k2