计算方法 线性方程组的解法.docVIP

0640303233 张世兵 PAGE 计算机科学与工程学院 实验五 线性方程组的解法 实验目的： 通过本课程的学习，学会编写全主原消去法的计算程序。掌握解线性方程组的最基本算法及其运用，进一步了解该解法的功能、优缺点，领会系数矩阵对解的影响。 实验内容： 高斯消去法： 要求：用随机函数随机产生一个n阶线性方程组，用不换行列的全主元法求解线性方程组，并且对该解进行验证。 方法概要 解线性方程组： 高斯消去法分消元过程和回代过程两个环节。 a. 消元过程 对k=1,2,……,n-1进行下列运算： 1) 选主元，确定r,t,使 若，则系数矩阵为奇异的，停止计算，否则进入下一步。 2) 交换A中的r,t两行及r,t两列，并记下交换的下标j,k. 3)对i=k+1,k+2,…,n;j=k+1,k+2,…,n+1,计算 b. 回代过程： 对k=n,n-1,……,1;按下列公式计算： c. 恢复解的顺序 (2) 计算框图 start start Input A，b，n K=1 选主元素 k=N? 输出 End k=k+1 l=k+1 l=k+1 i=n? i=i+1 从主程序来 d=0? 输出奇异标志 End l=k? = 返回主程序 三 实验结果 运行第一次 运行第二次 (高斯赛德尔迭代失败) 四 心得体会： 此实验在编程时由于循环比较多，所以需要格外注意。 为了顺利的看懂算法，我在做实验前仔细的看了看课本，找了一个例题，根据实验步骤一步一步的解出结果这样有利于理解算法的每一步骤究竟实现什么功能，在后来的试验中证明了此举并不是白做，而是很有成效。其他的同学由于不知每一步的具体功能是干什么的，在调试过程中出现了问题却不知道是哪而出现了问题，不容易改正，建议此方法推广使用。 此实验还有一个相对难实现的功能即将产生的随机方程组以人容易比较理解的方式输出来。 后来通过查阅资料利用两个函数strcat和strvcat，strcat实现将每行的方程连接起来，strvcat实现将每列的方程连接起来成为一个方程组，然后通过set函数输在一个静态文本框中，最后实现功能。 还有一个问题由于方程是随机生成的，所以是对角占优方程的可能性较小，所以就通过一个循环将主对角的元素扩大了5倍，这样对角占优的是可能性就大了。 五 程序代码 function varargout = xianxingfangcheng(varargin) % XIANXINGFANGCHENG M-file for xianxingfangcheng.fig % XIANXINGFANGCHENG, by itself, creates a new XIANXINGFANGCHENG or raises the existing % singleton*. % % H = XIANXINGFANGCHENG returns the handle to a new XIANXINGFANGCHENG or the handle to % the existing singleton*. % % XIANXINGFANGCHENG(CALLBACK,hObject,eventData,handles,...) calls the local % function named CALLBACK in XIANXINGFANGCHENG.M with the given input arguments. % % XIANXINGFANGCHENG(Property,Value,...) creates a new XIANXINGFANGCHENG or raises the % existing singleton*. Starting from the left, property value pairs are % applied to the GUI before xianxingfangcheng_OpeningFunction gets called. An % unrecognized property name or invalid value makes property application % stop. All inputs are passed to xianxingfangcheng_Open