MATLAB方程组求解实验报告.doc

《MATLAB方程组求解实验报告.doc》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、《MATLAB应用与仿真教程》实验报告实验序号：1　实验项目名称：方程组求解学　　号姓　　名专业、班实验地点指导教师时间一、实验要求和目的学会用MATLAB语言绘图、学会用MATLAB语言编写程序求解线性、非线性方程组的解。学习线性、非线性方程组求解方法，运用MATLAB语言编写应用程序，完成对线性、非线性方程组的求解。二、主要仪器设备：联想商用电脑、Matlab2010版。三实验任务（1）已知，求A*B，A+B，A-B，A/B，AB（2）求题（1）矩阵A的逆，特征向量、特征值。（注意分别使用函数inv(),[UD]=eig()）（3）已知方程组，其

2、中A为题（1）中矩阵，求该方程组四实验原理（一）线性方程组求解将线性方程的求解分为两类：一类是方程组求唯一解或求特解，另一类是方程组求无穷解即通解。可以通过系数矩阵的秩来判断：若系数矩阵的秩r=n（n为方程组中未知变量的个数），则有唯一解若系数矩阵的秩r

3、1.2利用矩阵的LU、QR和cholesky分解求方程组的解(1).LU分解：LU分解又称Gauss消去分解，可把任意方阵分解为下三角矩阵的基本变换形式（行交换）和上三角矩阵的乘积。即A=LU，L为下三角阵，U为上三角阵。则：A*X=b变成L*U*X=b则X=U(Lb)大大提高运算速度。命令[L，U]=lu(A)(2).Cholesky分解若A为对称正定矩阵，则Cholesky分解可将矩阵A分解成上三角矩阵和其转置的乘积，即：其中R为上三角阵。方程A*X=b变成则命令R=chol(A)(3).QR分解对于任何长方矩阵A，都可以进行QR分解，其中Q为

4、正交矩阵，R为上三角矩阵的初等变换形式，即：A=QR方程A*X=b变形成QRX=b则X=R(Qb)命令[Q,R]=qr(A)2求线性齐次方程组的通解在Matlab中，函数null用来求解零空间，即满足A·X=0的解空间，实际上是求出解空间的一组基（基础解系）。格式：z=null%z的列向量为方程组的正交规范基，满足%z的列向量是方程AX=0的有理基3求非齐次线性方程组的通解非齐次线性方程组需要先判断方程组是否有解，若有解，再去求通解。步骤为：第一步：判断AX=b是否有解，若有解则进行第二步第二步：求AX=b的一个特解第三步：求AX=0的通解第四步：

5、AX=b的通解=AX=0的通解+AX=b的一个特解。五.实验代码（1）已知，求A*B，A+B，A-B，A/B，AB1）a=[211;131;114];b=[12-31;-1307;-3634];>>a*bans=20304369356-1511442）a=[211;131;114];b=[12-31;-1307;-3634];>>a+bans=14-220338-27383)a=[211;131;114];b=[12-31;-1307;-3634];>>a-bans=-10402-27-64-5-304)a=[211;131;114];b=[12-3

6、1;-1307;-3634];>>a/bans=0.17430.04790.01440.09350.10850.00430.11260.02270.10975)a=[211;131;114];b=[12-31;-1307;-3634];>>abans=8.2941-7.9412-4.5882-2.352912.52940.7059-2.23530.35299.4706（2）求题（1）矩阵A的逆，特征向量、特征值。（注意分别使用函数inv(),[UD]=eig()）1)a=[211;131;114];>>b=inv(a)b=0.6471-0.1765-

7、0.1176-0.17650.4118-0.0588-0.1176-0.05880.29412)a=[211;131;114];>>[V,D]=eig(a,'nobalance')V=0.88770.23320.3971-0.42710.73920.5207-0.1721-0.63180.7558D=1.32490002.46080005.2143（1）已知方程组，其中A为题（1）中矩阵，求该方程组方法一a=[211;131;114]a=211131114>>b=[1;-1;3]b=1-13>>[L,U,P]=LU(a)>>[L,U,P]=lu(a)L

8、=1.0000000.50001.000000.50000.20001.0000U=2.000