解线性方程组的直接法的matlab实现1

《解线性方程组的直接法的matlab实现1》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。

1、赤峰学院本科毕业论文（设计）学校编码：******分类号密级学号：*********UDC本科毕业论文(设计)解线性方程组的直接法的Matlab实现学生姓名：***********************所属院部：***********************专业：***********************指导教师：***********************2014年5月14日-17-赤峰学院本科毕业论文（设计）摘要：给出用MATLAB解线性方程组的各种方法,用MATLAB直接操作,不用编程,便可立即求出线性方程组的解.方法直

2、观、简便、速度快,具有较强的实用性,另外提供了Jacobi迭代法程序.关键字：线性方程组数值解程序设计MATLABJacobi迭代法数据结构1引言在自然科学与社会科学的研究中，常常需要求解线性代数方程组，如实验数据的曲线、曲面的拟合和用差分法或有限元法解偏微分方程等都要用到线性代数方程组的求解。由于从不同的问题导出的线性代数方程组的系数矩阵不同，比如：矩阵阶数的大小、矩阵中的非零元稠密情况等，粗略地，系数矩阵可以分为低阶稠密矩阵和大型稀疏矩阵。关于线性代数程组的求解，主要分为直接法和迭代法两大类，在理论上，用直接法可以通过有限步的计算

3、得到精确解，而迭代法是通过逐次迭代逼近来求得近似解，实际上，由于舍入误差的影响，由直接法得到的解也不精确，因此，在某些需要高精度解的问题中，常常把由直接法得到的解再运用迭代法迭代若干步，以提高解的精度，一般地说，对于低阶稠密的线性代数方程组以及大型带形方程组的求解，采用直接法比较有效，而对于大型稀疏（非带形）方程组则用迭代法求解比较有利。当然，采用直接法，还是迭代法，还是直接法与迭代法交替运用，要根据具体情况确定。本章主要讨论一些最基本的直接法，并在此基础上讨论它的变形情况，大量的科技与工程实际问题，常常归结为解线性代数方程组，有关线

4、性方程组解的存在性和唯一性在“线性代数”理论中已经作过详细介绍，本章的主要任务是讨论系数行列式不为零的n阶非齐次线性方程组Ax=b的两类主要求方法：直接法（精确法）和迭代法。对它的解法我们最熟悉的就是主元消去法,但它只是适用于A是低价稠密的矩阵,对于由工程技术中产生的大型稀疏矩阵方程组(即A的阶数n很大,但零元素较多,例如求某些偏微分方程数值解所产生的线性方程组,n≥104),还需利用迭代法求解。如在计算机内存和运算两方面,都可以根据A中有大量零元素的特点采用迭代法。本文将介绍两种常见的迭代:Jacobi迭代法和Gauss-Seide

5、l迭代,并用迭代法在数学软件Matlab上实现线性方程组的解。1迭代法的基本思想迭代法是按照某种规则构造一个向量序列{x(k)},使其极限向量x*是Ax=b的精确解。因此,对迭代法来说一般有下面几个问题:(1)如何构造迭代序列?(2)构造的迭代法序列是否收敛?在什么情况下收敛?(3)如果收敛,收敛的速度如何?我们应该给予量的刻划,用以比较各种迭代法收敛的快慢。2相关知识线性方程组的概念及分类线性方程组的一般形式为a11x1+a12x2+…+a1nxn=b1a21x1+a22x2+…-17-赤峰学院本科毕业论文（设计）+a2nxn=b2

6、am1x1+am2x2+…+amnxn=b{n(1)若记X=x1x2(…xn)T,b=b1b2(…bn)TA=a11a12…a1na21a22…a2n…am1am2…amn则线性方程组(1)记为AX=b.(2)若b的元素不全为零,则称方程组(1)为非齐次线性方程组;若b的元素全为零,即b1=b2=…=bn=0,则AX=0.(3)并称方程组(3)为齐次线性方程组,也称作方程组(2)的导出方程组,称(Ab)=a11a12…a1n…b1a21a22…a2n…b2…am1am2…amn…bn为线性方程组(1)的增广矩阵,记作A.

7、若在方程组(1)中,当mn,3、算法用高斯消元法解线性方程组bAX的MATLAB程序输入的量：系数矩阵A和常系数向量b；输出的量：系数矩阵A和增广矩阵B的秩RA,RB,方程组中未知量的个数n和有关方程组解X及其解的信息.function[RA,RB,n,X]=gaus(A,b)B=[Ab];n=length(b);RA=rank(A);RB=rank(B);zhica=RB-RA;ifzhica>0,disp('请注意：因为RA~=RB，所以此方程组无解.')returnendifRA==RBifRA==ndisp('请注意：因为RA

8、=RB=n，所以此方程组有唯一解.')X=zeros(n,1);C=zeros(1,n+1);forp=1:n-1（1）LU分解法lu分解法解线性方程组functionx=luxiaoyuan(A,b)[m,n]=siz