Mathematica(Gauss消元法)解线性方程组-高斯消元法-列主元法详解
高斯列主元消元法解线性方程组
高斯列主元消元法解线性方程组 一、题目：用Gauss 列主元消去法解线性方程组Ax b =，其中， A=17.031 -0.615 -2.991 1.007 -1.006 0.000-1.000 34.211 -1.000 -2.100 0.300 -1.7000.000 0.500 13.000 -0.500 1.000 -1.5004.501 3.110 -3.907 -61.705 12.170 8.9990.101 -8.012 -0.017 -0.910 4.918 0.1001.000 2.000 3.000 4.500 5.000 21.803?? ? ? ? ? ? ? ? ??? 0.230 -52.322 54.000 240.236 29.304 -117.818b ?? ? ? ?= ? ? ? ? ??? T X=(0.907099 -1.961798 3.293738 -4.500708 3.029344 -5.255068) 二、原理及步骤分析 设 n n ij R a A ?∈=][)1(，n n R b b b b ∈=],,,[)1()2(2)1(1 。若约化主元素 ),,2,1(0)(n k a k kk =≠，则通过高斯消元法将方程b AX =约化为三角形方程组求解。 如果在消元过程中发现某个约化主元0) (=k kk a , 则第K 次消元就无法进行。此外，即 使所有约化主元全不为零，虽然可以完成方程组的求解，但也无法保证结果的可靠性，因为计算过程中存在舍入误差。 为减少计算过程中的舍入误差对解的影响，在每次消元前，应先选择绝对值尽可能大的元作为约元的主元，如果在子块的第一列中选取主元，则相应方法称为列主元消元法。相应过程为： （1）选主元：在子块的第一列中选择一个元) (k k i k a 使) (max k ik n i k k k i a a k ≤≤= 并将第k 行元与第k i 行元互换。 （2）消元计算：对k=1,2,……n-1依次计算 ()()()?? ?? ?????++=-=++=-=++==++n k k i b m b b n k k j i a m a a n k k i a a m k k ik k i k i k kj ik k ij k ij k kk k ik k ik ,,2,1,,2,1,,,2,1) ()()1() ()()1()() ()( （3）回代求解
用高斯消元法求解线性代数方程组
用高斯消元法求解线性代数方程组 1234111 5 -413-2823113-2104151 3-21719x x x x ??????????????????=?????? ?????? ?????? 1111X *??????=?????? （X*是方程组的精确解） 1 高斯消去法 1.1 基本思想及计算过程 高斯（Gauss ）消去法是解线性方程组最常用的方法之一，它的基本思想是通过逐步消元，把方程组化为系数矩阵为三角形矩阵的同解方程组，然后用回代法解此三角形方程组得原方程组的解。 为便于叙述，先以一个三阶线性方程组为例来说明高斯消去法的基本思想。 ??? ??=++II =++I =++III) (323034)(5 253)(6432321 321321x x x x x x x x x 把方程（I ）乘（2 3 - ）后加到方程（II ）上去，把方程（I ）乘（2 4- ）后加到方程（III ）上 去，即可消去方程（II ）、（III ）中的x 1，得同解方程组 ?? ? ??=+-II -=-I =++III) (20 223)(445.0)(6 4323232321x x x x x x x 将方程（II ）乘（ 5 .03 ）后加于方程（III ），得同解方程组： ?? ? ??-=-II -=-I =++III) (42)(445.0)(6432332321x x x x x x 由回代公式（3.5）得x 3 = 2，x 2 = 8，x 1 = -13。 下面考察一般形式的线性方程组的解法，为叙述问题方便，将b i 写成a i , n +1，i = 1, 2,…,n 。
高斯消元法(完整)
高斯消元法解线性方程组 在工程技术和工程管理中有许多问题经常可以归结为线性方程组类型的数学模型，这些模型中方程和未知量个数常常有多个，而且方程个数与未知量个数也不一定相同。那么这样的线性方程组是否有解呢？如果有解，解是否唯一？若解不唯一，解的结构如何呢？这就是下面要讨论的问题。 一、线性方程组 设含有n 个未知量、有m 个方程式组成的方程组 a x a x a x b a x a x a x b a x a x a x b n n n n m m mn n m 11112211211222221122+++=+++=+++=???????ΛΛΛΛΛΛΛΛΛ （3.1） 其中系数a ij ，常数b j 都是已知数，x i 是未知量（也称为未知数）。当右端常数项b 1, b 2, …, b m 不全为0时，称方程组（3.1）为非齐次线性方程组；当b 1=b 2= … =b m = 0时，即 a x a x a x a x a x a x a x a x a x n n n n m m mn n 111122121122221122000 +++=+++=+++=???????ΛΛΛΛΛΛΛΛΛ （3.2） 称为齐次线性方程组。 由n 个数k 1, k 2, …, k n 组成的一个有序数组（k 1, k 2, …, k n ），如果将它们依次代入方程组（3.1）中的x 1, x 2, …, x n 后，（3.1）中的每个方程都变成恒等式，则称这个有序数组（k 1, k 2, …, k n ）为方程组（3.1）的一个解。显然由x 1=0, x 2=0, …, x n =0组成的有序数组（0, 0, …, 0）是齐次线性方程组（3.2）的一个解，称之为齐次线性方程组（3.2）的零解，而当齐次线性方程组的未知量取值不全为零时，称之为非零解。 （利用矩阵来讨论线性方程组的解的情况或求线性方程组的解是很方便的。因此，我们先给出线性方程组的矩阵表示形式。） 非齐次线性方程组（3.1）的矩阵表示形式为： AX = B 其中 A = ????????????mn m m n n a a a a a a a a a ΛΛΛΛΛΛΛ212222111211，X = ????????????n x x x M 21， B = ????? ???????n b b b M 21 称A 为方程组（3.1）的系数矩阵，X 为未知矩阵，B 为常数矩阵。将系数矩阵A 和常数矩阵B 放在一起构成的矩阵
列主元消去法
实验一 列主元消去法 【实验内容】1. 掌握列主元消去法的基本思路和迭代步骤 2. 并能够利用列主元的高斯消去法解任意阶数的线性方程组； 【实验方法与步骤】列主元消去法编写程序 1．列主元消去法基本思路 设有线性方程组Ax b =，设A 是可逆矩阵。列主元消去法的基本思想就是通过列主元的选取将初等行变换作用于方程组的增广矩阵[]|B A b =，将其中的A 变换成一个上三角矩阵，然后求解这个三角形方程组。 2．列主元高斯消去法算法描述 将方程组用增广矩阵[]()(1)|ij n n B A b a ?+==表示。 步骤1：消元过程，对1,2,,1k n =- （1） 选主元，找{},1,,k i k k n ∈+ 使得 ,max k i k ik k i n a a ≤≤= （2） 如果,0k i k a =，则矩阵A 奇异，程序结束；否则执行（3）； （3） 如果k i k ≠，则交换第k 行与第k i 行对应元素位置，k kj i j a a ?， ,,1j k n =+ ； （4） 消元，对,,i k n = ，计算/,ik ik kk l a a =对1,,1j k n =++ ，计算 .ij ij ik kj a a l a =- 步骤 2：回代过程： （1） 若0,nn a =则矩阵奇异，程序结束；否则执行（2）； （2） ,1/;n n n nn x a a +=对1,,2,1i n =- ，计算 ,11/n i i n ij j ii j i x a a x a +=+??=- ??? ∑ 习题3第一题程序如下
#include 《Gauss列主元消去法》实验报告 实验名称：Gauss列主元消去法程序设计？？？成绩：_________ 专业班级：数学与应用数学1202班？姓名：王晓阳？??学号： 实？验？日？期：？2014?年11月10日 实验报告日期：？2014年？11月10日 一.实验目的 1. 学习Gauss消去法的基本思路和迭代步骤. 2. 学会运用matlab编写高斯消去法和列主元消去法程序，求解线性方程组. 3. 当绝对值较小时，采用高斯列主元消去法? 4. 培养编程与上机调试能力. 二、实验内容 用消去法解线性方程组的基本思想是用逐次消去未知数的方法把原线性方程组Ax二b 化为与其等价的三角形线性方程组，而求解三角形线性方程组可用回代的方法求解 1. 求解一般线性方程组的高斯消去法? (1) 消元过程： 设a kk k-0 ,第i个方程减去第k个方程的m ik Tk k倍，("k 1^1, n)，得到 A k1x=b k1. 经过n-1次消元，可把方程组A1^b1化为上三角方程组A n x=b n. ⑵回代过程: 以解如下线性方程组为例测试结果 2. 列主元消去法 由高斯消去法可知，在消元过程中可能出现a kk k =0的情况，这是消去法将无法进行, 即使主元素a kk k-0但很小时，用其作除数，会导致其他元素数量级的严重增长和舍入误差的扩散，最后也使得计算解不可靠.这时就需要选取主元素，假定线性方程组的系数矩阵A是菲奇异的. (1)消元过程： 对于k =1,2,川,n -1，进行如下步骤: 1) 按列选主元，记 2) 交换增广阵A的p,k两行的元素 A(k,j)=A(p,j) ( j=k,…,n +1) 3) 交换常数项b的p,k两行的元素。 b(k)=b(p) 4) 计算消元 (2) 回代过程 (3) 以解如下线性方程组为例测试结果 三、实验环境 MATLAB R2014a 四、实验步骤 计算方法实验报告1 【课题名称】 用列主元高斯消去法和列主元三角分解法解线性方程 【目的和意义】 高斯消去法是一个古老的求解线性方程组的方法，但由它改进得到的选主元的高斯消去法则是目前计算机上常用的解低阶稠密矩阵方程组的有效方法。 用高斯消去法解线性方程组的基本思想时用矩阵行的初等变换将系数矩阵A 约化为具有简单形式的矩阵（上三角矩阵、单位矩阵等），而三角形方程组则可以直接回带求解 用高斯消去法解线性方程组b Ax =（其中A ∈Rn ×n ）的计算量为：乘除法运算步骤为 32(1)(1)(21)(1)(1)262233n n n n n n n n n n n MD n ----+= +++=+-，加减运算步骤为 (1)(21)(1)(1)(1)(25) 6226n n n n n n n n n n AS -----+= ++= 。相比之下，传统的克莱姆 法则则较为繁琐，如求解20阶线性方程组，克莱姆法则大约要19 510?次乘法，而用高斯消去法只需要3060次乘除法。 在高斯消去法运算的过程中，如果出现abs(A(i,i))等于零或过小的情况，则会导致矩阵元素数量级严重增长和舍入误差的扩散，使得最后的计算结果不可靠，所以目前计算机上常用的解低阶稠密矩阵方程的快速有效的方法时列主元高斯消去法，从而使计算结果更加精确。 2、列主元三角分解法 高斯消去法的消去过程，实质上是将A 分解为两个三角矩阵的乘积A=LU ，并求解Ly=b 的过程。回带过程就是求解上三角方程组Ux=y 。所以在实际的运算中，矩阵L 和U 可以直接计算出，而不需要任何中间步骤，从而在计算过程中将高斯消去法的步骤进行了进一步的简略，大大提高了运算速度，这就是三角分解法 采用选主元的方式与列主元高斯消去法一样，也是为了避免除数过小，从而保证了计算的精确度 【计算公式】 1、 列主元高斯消去法 设有线性方程组Ax=b ，其中设A 为非奇异矩阵。方程组的增广矩阵为 第1步（k=1）：首先在A 的第一列中选取绝对值最大的元素 1 l a ，作为第一步的主元素: 111211212222112[,]n n n l n nn n a a a a b a a a b a a a b ?????? ?? =?????? ?? ????a b 消元法解线性方程组 学校：青海师范大学 院系：数学系 专业：数学与应用数学 班级：10B 指导教师：邓红梅 学号：20101611218 姓名：梅增旺 摘要：线性方程组在数学的各个分支，在自然科学，工程技术，生产实际中经常遇到，而且未知元的个数及方程的个数可达成百上千，因此它的理论是很重要的，其应用也很广泛。本篇将就解线性方程组在此做一浅谈，以消元法为主要方法。消元法是解一般线性方程组行之有效的方法，早在中学大家都已经有接触，消元法的基本思想是通消元变形把方程组化成容易求解的同解方程组进行求解。 关键字：线性方程组消元法求解 Abstract: linear equations in various branches of mathematics, natural science,engineering technology, often encountered in actual production, and the unknown element number and the number of equations can be hundreds, so itis important in the theory, its application is very extensive. This article on thesolution of linear equations based on a discussion, mainly by means ofelimination method. Elimination method is the general linear equations ofeffective early in high school, everyone has a contact, the basic idea ofelimination method is through the elimination of the equations of deformationinto easy to solve with the solution of equations. Keywords:elimination method for solving linear equations #include "Stdio.h" #include "Conio.h" /*L是矩阵的行减1,从程序上看是最外层循环的次数 N 对应矩阵的行数，M对应矩阵的列数 可以通过改变L、N、M来控制矩的阶数 */ #define L 3 #define N 4 #define M 5 void gauss(double a[N][M],double x[N]) {int i,j,l,n,m,k=0; double temp[N]; /*第一个do-while是将增广矩阵消成上三角形式*/ do{n=0; for(l=k;l Matlab之Gauss消元法解线性方程组 1.Gauss消元法 function x=DelGauss(a,b) %Gauss消去法 [n,m]=size(a); nb=length(b); det=1;%存储行列式值 x=zeros(n,1); for k=1:n-1 for i=k+1:n if a(k,k)==0 return end m=a(i,k)/a(k,k); for j=k+1:n a(i,j)=a(i,j)-m*a(k,j); end b(i)=b(i)-m*b(k); end det=det*a(k,k);%计算行列式 end det=det*a(n,n); for k=n:-1:1%回代求解 for j=k+1:n b(k)=b(k)-a(k,j)*x(j); end x(k)=b(k)/a(k,k); end Example: >>A=[1.0170-0.00920.0095;-0.00920.99030.0136;0.00950.0136 0.9898]; >>b=[101]'; >>x=DelGauss(A,b) x= 0.9739 -0.0047 1.0010 2.列主元Gauss消去法： function x=detGauss(a,b) %Gauss列主元消去法 [n,m]=size(a); nb=length(b); det=1;%存储行列式值 x=zeros(n,1); for k=1:n-1 amax=0;%选主元 for i=k:n if abs(a(i,k))>amax amax=abs(a(i,k));r=i; end end if amax<1e-10 return; 《数值分析》实验报告 一、实验目的与要求 1．掌握高斯消去法的基本思路和迭代步骤； 2．培养编程与上机调试能力。 二、实验内容 1．编写用高斯消元法解线性方程组的MATLAB程序，并求解下面的线性方程组，然后用逆矩阵解方程组的方法验证. 5x?2x?x?80.101x?2.304x?3.555x?1.183??312312??（1）（2） 21x?8x?32x?2.137x?3.712x?4.623?1.347x???312312??1x?3x?6x??2.835x?1.072x?5.643x?3.035??132 312 2．编写用列主元高斯消元法解线性方程组的MATLAB程序，并求解下面的线性方程组，然后用逆矩阵解方程组的方法验证. 5x?2x?x?80.101x?2.304x?3.555x?1.183??312312??（1）（2） 2x?8x?3x?212.137?4.6231.347?x?3.712x?x??321321??1x?3x?6x??2.835x?1.072x?5.643x?3.035??132 312三．MATLAB计算源程序 AX?b MATLAB1. 程序用高斯消元法解线性方程组的b；输入的量：系数矩阵和常系数向量A RA,RB, n方程组中未知量的个数的秩输出的量：系数矩阵和增广矩阵BA.及其解的信息和有关方程组解X gaus(A,b) function [RA,RB,n,X]=B=[A b]; n=length(b); RA=rank(A); RB=rank(B);zhica=RB-RA; if zhica>0, disp('RA~=RB.') ，所以此方程组无解请注意：因为return end if RA==RB if RA==n disp('RA=RB=n.') ，所以此方程组有唯一解请注意：因为X=zeros(n,1); C=zeros(1,n+1); for p= 1:n-1 for k=p+1:n m= B(k,p)/ B(p,p); B(k,p:n+1)= B(k,p:n+1)-m* B(p,p:n+1); 列主元高斯消去法求逆矩阵程序代码： #include M[i][j] = M[i][j] - M[t][j]*m; E[i][j] = E[i][j] - E[t][j]*m; } } } void HuiDai(int t) { int i,j; double max; max=M[t][t]; for(i=t;i 数值分析计算实习题 第5章解线性方程组的直接方法 【选题 列主元高斯消去法解线性方程组。 书上的计算实习题1、2、3都要求用列主元高斯消去法解线性方程组，所以考虑写一个普适的程序来实现。 对于线性方程组Ax二b,程序允许用户从文件读入矩阵数据或直接在屏幕输入数据。 文件输入格式要求: (1)第一行为一个整数n (2<=n<=100)，表示矩阵阶数。 (2)第2~n+l行为矩阵A各行列的值。 (3)第n+2~n+n+2行为矩阵b各行的值。 屏幕输入：按提示输入各个数据。 输出：A. b、det(A).列主元高斯消去计算过程、解向量X。 【算法说明】 设有线性方程组Ax=b,其中设A为非奇异矩阵。方程组的增广矩阵为 ?12 ?21 [Nb] = 第1步（k=l ）：首先在A的第一列中选取绝对值最大的元素?I,作为第一步的主元素: ?|| H0 然后交换（A, b）的第1行与第I行元素，再进行消元计算。 设列主元素消去法已经完成第1步到第k?l步的按列选主元，交换两行，消元计算得到与原方程组等价的方程组 A（k）x=b（k） 4? …4;) …唸) ? 忒 ? ? 輕 ■ [A.b]T[A ⑹,b")] = ??■ 咲■ ■ ■ ■ ■ * *■ 〃伏) ?? - % ■ 第k步计算如下: 对于 k=l, 2, ?…，0-1 （1）按列选主元：即确定t使 （2）如果tHk,则交换［A, b］第t行与第k行元素。（3）消元计算 5 4* J 叫=一鱼(=^ + 1,…,H) % 吗 <-?y + 〃如伽 (fJ = R + l,…/) b- <-勺+加汝仇， (i = /c + l,…,《) 消元乘数mik 满足: n (%-D 内) X1 < ------ -- ---- 9(j = ? 一 1,?一2■…J)tk M 1,(,=斤 +1, ???,?) fet e (4)回代求解 例：用Gauss 列主元消去法、QR 方法求解如下方程组： 12342212141 312.4201123230x x x x ?????? ? ? ?- ? ? ?= ? ? ?-- ? ? ??????? 1. 1）Gauss 列主元法源程序： function x=Gauss(A,b) [m,n]=size(A); if m~=n error('矩阵不是方阵') return end B=[A,b]; n=length(A); for j=1:n-1 q=[zeros(j-1,1);B(j:n,j)]; [c,r]=max(abs(q)); %c 为列主元，r 为所在行 if r~=j temp=B(j,:); %交换两行 B(j,:)=B(r,:); B(r,:)=temp; end for i=j+1:n B(i,:)=B(i,:)-B(j,:)*(B(i,j)/c); end end x(n)=B(n,n+1)/B(n,n); for i=n-1:-1:1 for j=i:n-1 B(i,n+1)=B(i,n+1)-B(i,j+1)*x(j+1); end x(i)=B(i,n+1)/B(i,i); end 2)在命令窗口输入A ，b,得到x 的近似解： >> A=[2,2,1,2;4,1,3,-1;-4,-2,0,1;2,3,2,3]; >> b=[1;2;1;0]; >> x=Gauss(A,b) x = 1.5417 -2.7500 0.0833 1.6667 2. 1）QR方法源程序： function [Q,R,X]=qrfj(A,b) [m,n]=size(A); if m 解线性方程组的消元法及其应用 （朱立平 曲小刚） ● 教学目标与要求 通过本节的学习，使学生熟练掌握一种求解方程组的比较简便且实用的方法—高斯消元法，并能够熟练应用消元法将矩阵化为阶梯形矩阵和求矩阵的逆矩阵. ● 教学重点与难点 教学重点：解线性方程组的高斯消元法，利用消元法求逆矩阵. 教学难点：高斯消元法，利用消元法求逆矩阵. ● 教学方法与建议 先向学生说明由于运算量的庞大，克莱姆法则在实际应用中是很麻烦的，然后通过解具体的方程组，让学生自己归纳出在解方程组的时候需要做的三种变换，从而引出解高阶方程组比较简便的一种方法—高斯消元法，其三种变换的实质就是对增广矩阵的初等行变换，最后介绍利用消元法可以将矩阵化为阶梯形矩阵以及求矩阵的逆。 ● 教学过程设计 1.问题的提出 由前面第二章的知识，我们知道当方程组的解唯一的时候，可以利用克莱姆法则求出方程组的解，但随着方程组阶数的增高，需要计算的行列式的阶数和个数也增多，从而运算量也越来越大，因此在实际求解中该方法是很麻烦的. 引例 解线性方程组 ??? ??=+-=+=++132724524321 21321x x x x x x x x )3()2()1( 解 （1）???→??)2()1(?????=+-=++=+13245247 232132121x x x x x x x x )3()2()1(????→?+-?+-?) 3()2()1()2()4()1(?????-=+-=+=+133524567232 3221x x x x x x )3()2()1( ????→?+-?)3()65 ()2(??????? =--=+=+76 724567233221x x x x x )3()2()1( 用回代的方法求出解即可. 问题：观察解此方程组的过程，我们总共作了三种变换：（1）交换方程次序，（2）以不等于零的数乘某个方程，（3）一个方程加上另一个方程的k 倍.那么对于高阶方程组来说，是否也可以考虑用此方法. 2.矩阵的初等变换 定义1 阶梯形矩阵是指每一非零行第一个非零元素前的零元素个数随行序数的增加而增加的矩阵. 定义2 下面的三种变换统称为矩阵的初等行变换： i. 互换矩阵的两行（例如第i 行与第j 行，记作j i r r ?）， ii. 用数0≠k 乘矩阵的某行的所有元素（例如第i 行乘k ，记作i kr ）， iii. 把矩阵某行的所有元素的k 倍加到另一行的对应元素上去（例如第j 行的k 倍加到第i 行上，记作j i kr r +）. 同理可以定义矩阵的初等列变换. 定义 3 如果矩阵A 经过有限次初等变换变为矩阵B ，则称矩阵A 与B 等价，记作 A ~ B . 注：任意一个矩阵总可以经过初等变换化为阶梯形矩阵. 3. 高斯消元法 对于一般的n 阶线性方程组 ?????? ?=++=+++=+++n n nn n n n n n n b x a x a x a b x a x a x a b x a x a x a ΛΛΛΛΛΛΛΛΛΛΛΛΛΛΛ22112 22221211 1212111 )()2()1(n （3.1） 若系数行列式0det ≠A ，即方程组有唯一解，则其消元过程如下： 第一步，设方程（1）中1x 的系数01≠l a 将方程)(l 与（1）对调，使对调后的第一个方程1x 的系数不为零.作)1(11 1 a a i i - ),3,2(n i Λ=，得到同解方程组 ?? ? ????=++=++=+++)1()1(2)1(2) 1(2 )1(22)1(22)0(1)0(12)0(121)0(11n n nn n n n n n b x a x a b x a x a b x a x a x a ΛΛΛΛΛΛΛΛΛΛΛΛΛ （3.2） 第二步，设0) 1(22≠a ，保留第二个方程，消去它以下方程中的含2x 的项，得 河北工业大学计算机软件技术基础（VC）课程设计报告 学院管理班级管理104班姓名杨立宝 __ 学号 101707____ 成绩 __ ____ 一、题目： 求线性代数方程组的解（高斯消去法）（C13） 二、设计思路 1、总体设计 1）分析程序的功能 第一：编写输入程序，通过键盘先输入对应的已知量及函数的大小n和系数a[i]和得数b[i]。 第二：编写中间程序，通过函数的调用先定义线性代数方程，然后通过程序求出方程的梯形矩阵系数，并最终得出结果。 第三编写输出程序，输出最终结果。 2）系统总体结构：设计程序的组成模块，简述各模块功能。 模块一：各函数的具体内容 A：三个输入函数，分别输入n，一维数组，二维数组。即输入已知量。 B：中间运算函数，计算是使得方程系数所成的矩阵成梯形矩阵，未知数的结果。即计算中间变量及结果。 C：最后输出函数，输出最后计算结果。 模块二：各函数原型的声明 a写头文件。 b变量声明:存放输入数据的数组的声明，存放中间变量的数组的声明，存放运算结果的数组的声明。分别存放对应数据。 c输入有关操作的文字 d函数调用，在运算中自动调用对应的函数解决对应问题。 模块三:主函数 2、各功能模块的设计：说明各功能模块的实现方法 模块一：各个函数的声明，直接声明。 模块二：各函数都通过for循环来实现各个数组之间的基本运算。 3、设计中的主要困难及解决方案 在这部分论述设计中遇到的主要困难及解决方案。 1）困难1 函数调用是怎么用？ 解决方案： 仔细阅读课本，以及同学之间的讨论，和老师的帮助。 4、你所设计的程序最终完成的功能 1）说明你编制的程序能完成的功能 输入线性代数的系数后，运行程序即可得到梯形矩阵和结果。 2）准备的测试数据及运行结果 三、程序清单 如果是使用一个文件完成的程序，只需列出程序代码。 如果是使用多文件完成的程序，首先说明程序中的代码存放在哪些文件中，说明文件名（例如：本程序包含first.cpp、second.cpp、third.cpp和all.h四个文件）；然后依次给出每个文件名及该文件清单，例如： #include 贵州师范大学数学与计算机科学学院学生实验报告 课程名称： 数值分析 班级： 数本（一）班 实验日期： 年 月 日 学 号： 090704020098（81） 姓名： 吴胜 指导教师： 杨一都 实验成绩： 一、实验名称 实验五:线性方程组的数值解法 二、实验目的及要求 1. 让学生掌握用列主元gauss 消去法、超松弛迭代法求解线性方程组. 2. 培养Matlab 编程与上机调试能力. 三、实验环境 每人一台计算机，要求安装Windows XP 操作系统，Microsoft office2003、MATLAB6.5(或7.0). 四、实验内容 1. 编制逐次超松弛迭代(SOR 迭代)函数(子程序),并用于求解方程组 ????? ??=-++=+-+=++-=+++-1 4141 4144321 432143214321x x x x x x x x x x x x x x x x 取初始向量T x )1,1,1,1()0(=,迭代控制条件为 5 )1()(10 2 1||||--?≤ -k k x x 请绘制出迭代次数与松弛因子关系的函数曲线,给出最佳松弛因子.SOR 迭代 的收敛速度是否一定比Gauss-Seidel 迭代快? 2. 编制列主元 Gauss 消去法函数(子程序),并用于解 ??? ??=++-=-+-=+-6 15318153312321 321321x x x x x x x x x 要求输出方程组的解和消元后的增广矩阵. 注：题2必须写实验报告 五、算法描述及实验步骤 Gauss 消去法： 功能 解方程组b Ax = . 输入 n ,n n ij a A ?=)(,T n b b b b ),,,(21 =. 输出 方程组的解T n x x x x ),,,(21 =或失败信息. 高斯消元法简介 一，教学目标 知识与技能：了解高斯消元法 过程与方法：直接演示说明，学习做简单练习 情感，态度和价值观：进一步体会解方程组的根本思想消元，通过高斯消元的学习增强学习数学的能力 二，重点与难点：高斯消元法 三，课型：新授课 四，教学过程： 1.在前面的几节课，已经用加减消元和代入消元法求解二元或者三元一次方程组，其基本的思想就是从已知的方程导出未知数较少的方程组，直到最后得到一个一元一次方程，这种做法可适用于一般的n 元线性方程组（线性方程组），但是由于未知数的增加，我们希望我们的消元是有规律的，以避免混乱，下面介绍高斯消元法 2.例1：解方程组 1234123412341234251027612632517315292763 x x x x x x x x x x x x x x x x ---=?? -++-=?? ---=??--++=-? 解：把第一个方程的2倍，-3倍，5倍分别加到第2，3，4个方程上，可以消去2，3，4个 方程的未知数1x 12342342342342510 522226 2 1 7213 x x x x x x x x x x x x x ---=?? +-=?? +-=??--+=-? 为了使以后少出现分数运算，交换第二，三个方程的位置 12342342342342510 2 1 522226 7213 x x x x x x x x x x x x x ---=?? +-=?? +-=??--+=-? 把第2个方程的-5倍，7倍分别加到第3，4个方程，可以消去第3，4个方程未知数2x 123423434342510 2 1 31221 6126 x x x x x x x x x x x ---=?? +-=?? --=??-=-? 整理一下方程，第3个方程的左右两边乘以13 - ，第4个方程左右两边乘以1 6 123423434342510 2 1 47 21 x x x x x x x x x x x ---=?? +-=?? +=-??-=-? 用高斯消元法和列主元消去法求解线性代数方程组 （X*是方程组的精确解） 1 高斯消去法 1.1 基本思想及计算过程 高斯（Gauss ）消去法是解线性方程组最常用的方法之一，它的基本思想是通过逐步消元，把方程组化为系数矩阵为三角形矩阵的同解方程组，然后用回代法解此三角形方程组得原方程组的解。 为便于叙述，先以一个三阶线性方程组为例来说明高斯消去法的基本思想。 ??? ??=++II =++I =++III) (323034)(5 253)(6 432321 321321x x x x x x x x x 把方程（I ）乘（2 3 - ）后加到方程（II ）上去，把方程（I ）乘（2 4- ）后加到方程（III ）上 去，即可消去方程（II ）、（III ）中的x 1，得同解方程组 ?? ? ??=+-II -=-I =++III) (20 223)(445.0)(6 4323232321x x x x x x x 将方程（II ）乘（ 5 .03 ）后加于方程（III ），得同解方程组： ?? ? ??-=-II -=-I =++III) (42)(445.0)(6432332321x x x x x x 由回代公式（3.5）得x 3 = 2，x 2 = 8，x 1 = -13。 下面考察一般形式的线性方程组的解法，为叙述问题方便，将b i 写成a i , n +1，i = 1, 2,…,n 。 ??? ?? ??=++++=++++=+++++++1,3322111 ,223232221211,11313212111n n n nn n n n n n n n n n a x a x a x a x a a x a x a x a x a a x a x a x a x a （1-1） 如果a 11 ≠ 0，将第一个方程中x 1的系数化为1，得 ) 1(1,1)1(12)1(121+=+++n n n a x a x a x 其中)0(11 ) 0()1(1a a a ij j = ， j = 1, …, n + 1（记ij ij a a =) 0(，i = 1, 2, …, n ; j = 1, 2, …, n + 1） 从其它n –1个方程中消x 1，使它变成如下形式 ?? ? ????=++=++=++++++)1(1,)1(2)1(2) 1(1 ,2)1(22)1(22) 1(1,1)1(12)1(121n n n nn n n n n n n n a x a x a a x a x a a x a x a x （1-2） 其中n i a m a a ij i ij ij ,,2)1(1) 1( =?-=，1,,3,211 )1(1 1+== n j a a m i i 由方程（1-1）到（1-2）的过程中，元素11a 起着重要的作用，特别地，把11a 称为主元素。 如果（1-2）中0) 1(22≠a ，则以) 1(22a 为主元素，又可以把方程组（1-2）化为： ?? ? ??????=++=++=+++=+++++++)2(1 ,)2(3)2(3) 3(1,3)2(33)2(33) 2(1 ,2)2(23)2(232) 1(1,1)1(12)1(121 n n n nn n n n n n n n n n n a x a x a a x a x a a x a x a x a x a x a x （1-3） 针对（1-3） 继续消元，重复同样的手段，第k 步所要加工的方程组是： ?? ?? ?? ?? ? ????=++=++=+++=+++=++++-+---+---+-----++) 1(1,)1()1() 1(1,)1()1() 1(1,1)1()1(11) 2(1 ,2)2(23)2(232) 1(1,1)1(13)1(132)1(121 k n n n k nn k k nk k n k n k nn k k kk k n k n k kn k k k k n n n n n n a x a x a a x a x a a x a x a x a x a x a x a x a x a x a xGauss列主元消去法程序设计
列主元高斯消去法和列主元三角分解法解线性方程
消元法解线性方程组
高斯消元法讲解
MATLAB之GAUSS消元法解线性方程组
完整版高斯消元法MATLAB实现
列主元高斯消去法求逆矩阵
数值分析计算实习题列主元高斯消去法解线性方程组
Gauss列主元消去法、QR(MATLAB)
(完整版)解线性方程组的消元法及其应用
C++课程设计高斯消元法求线性代数方程组的解
Gauss列主元消去法
高斯消元法简介
作业一 高斯消元法和列主元消元法
- 数值分析列主元消去法的实验报告
- Mathematica(Gauss消元法)解线性方程组-高斯消元法-列主元法详解
- 高斯列主元消元法解线性方程组
- 用高斯列主元消元法解线性方程组
- 数值分析计算实习题列主元高斯消去法解线性方程组
- 高斯消元法MATLAB实现
- Gauss列主元消去法程序设计
- 高斯列主元消元法解线性方程组
- 列主元高斯消去法和列主元三角分解法解线性方程
- 高斯列主元消元法解线性方程组
- Matlab程序Gauss列主元消去法
- 第3章3-01高斯消元法,列主元法
- 全主元高斯消去法求矩阵的
- 高斯法和列主元高斯消去法解线性方程组(MATLAB版)
- 列主元高斯消去法求逆矩阵
- 高斯列主元消元法解方程组
- (完整版)高斯消元法MATLAB实现
- 列主元高斯消元法
- Gauss列主元消去法程序设计
- 用高斯列主元消元法解线性方程组