matlab方程组求解实验报告

《matlab方程组求解实验报告》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、《MATLAB应用与仿真教程》实验报告实验序号：1实验项目名称：方程组求解学号姓名专业、班实验地点指导教师时间一、实验要求和目的学会用MATLAB语言绘图、学会用MATLAB语言编写程序求解线性、非线性方程组的解。学>』线性、非线性方程组求解方法，运用MATLAB语言编写应用程序，完成对线性、非线性方程组的求解。二、主要仪器设备：联想商用电脑、Matlab2010版。三实验任务<211、<12-31、(1)已知力=131B=-1307<11勺「3634?求A*B,A+B，A-B,A/B,AB(2)求题(1)矩阵A的逆，特征

2、向量、特征值。(注意分别使用函数iiw()，fUD]=eig())f1(3)己知方程组/^=氏/7=-1，其中A为题(1)中矩阵，求该方程组ft.3,四实验原理(一)线性方程组求解将线性方程的求解分为两类：一类是方程组求唯一解或求特解，另一类是方程组求无穷解即通解。可以通过系数矩阵的秩來判断：若系数矩阵的秩r=n(n为方程组中未知变量的个数)，则冇唯一解若系数矩阵的秩r

3、要用于解低阶稠密矩阵——直接法；另一类是解大型稀疏矩阵——迭代法。1.1利用矩阵除法求线性方程组的特解(或一个解)方程：AX=b解法：X=Ab1.2利用矩阵的LU、QR和cholesky分解求方程组的解(l).LU分解:LU分解又称Gauss消去分解，可把任意方阵分解为K三角矩阵的基本变换形式(行交换)和上三角矩阵的乘积。即A=LU，L为下三角阵，U为上三角阵。贝ij:A*X=b变成L*U*X=b则X=U(Lb)大大提高运算速度。命令[L，U]=lu(A)(2).Cholesky分解若A为对称正定矩阵，则Cholesky

4、分解可将矩阵A分解成上三角矩阵和其转置的乘积，即：A=R'*R其屮R为上三角阵。方程A*X=b变成R'*R*X=b则X=R(R、b)命令R=chol(A)(3).QR分解对于任何松方矩阵A,都可以进行QR分解，其屮卩为正交矩阵，R为上三角矩阵的初等变挽形式，B卩：A=QR方程A*X=b变形成QRX=b贝ijX=R(Qb)命令[Q，R]=qr(A)2求线性齐次方程组的通解在Matlab中，函数null用来求解零空间，即满足A•X=0的解空间，实际上是求出解空间的_组基(基础解系)。格式：z=null%z的列向量为方程组的正

5、交规范基，满足Z'xZ=Iz=null(A,'r')%z的列向量是方程AX=0的有理基3求非齐次线性方程组的通解非齐次线性方程组需要先判断方程组是否有解，若有解，再去求通解。步骤为：第一步：判断AX=b是否有解，若有解则进行第二步第二步：求AX=b的一个特解第三步：求AX=0的通解第四步：AX=b的通解=AX=0的通解+AX=b的一个特解。五.实验代码<211、<12-31>(1)已知4=131，B=-1307，求A*B,A+B,A-B,A/B,AB<11幻「3634;1)a=[211;131;114];b=[12-31;-

6、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-30a=[211;131;114];b=[12-31;»a/bans=-1307;-3634];0.17430.04790.01440.09350.10850.00430.11265)0.02270.1

7、097a:[211;131:114];b=[12-31;»abans=-1307;-3634];8.2941-7.9412-4.5882-2.352912.52940.7059-2.23530.35299.4706(2)求题(1)矩阵A的逆，特征向量、特扯值。(注意分别使用函数inv()，[UD]=eig())1)a=[211;1»b=inv(a)b=31:114];0.6471-0.1765-0.1176-0.17650.4118-0.0588-0.11762)-0.05880.2941a=[211;131;114];»[

8、V,D]=eiV=g(a，’nobalance’)0.88770.23320.3971-0.42710.73920.5207-0.1721-0.63180.7558D=1.32490002.46080005.214311(2)已知方程组/1%=/7，/?=-1，其中A为题(1