性方程组数值解法.docxVIP

计算方法实验题目：班级：学号：姓名：1 实验目的通过编程加深对列主元高斯消去法、LU三角分解法和雅克比迭代法等求解多项式方程方法的理解观察上述三种方法的计算稳定性和求解精度并比较各种方法利弊2 实验步骤2.1环境配置：VS2013，C++控制台程序2.2添加头文件#include stdio.h#include stdlib.h#include stdafx.h#includeiostream2.3主要模块程序一共分成三层，最底层是数据结构部分，负责存储数据，第二层是交互部分，即多项式方程部分，负责输入输出获得数据，最上层是核心的算法部分，负责处理已获得的数据。具体功能如下：数据结构部分数据结构部分是整个程序的最底层，负责存储部分。因数组作为数据元素插入和删除操作较少，而顺序表空间利用率大且查看方便，故此程序选用二维顺序表保存系数。数据结构文件中写的是有关其的所有基本操作以供其他文件调用。多项式方程部分多项式方程部分是程序的第二层，内容是有关方程组的所有函数、构建方程、输出方程等等，同时在此文件中获得方程系数并储存，同时此文件还负责显示菜单部分。算法部分此文件负责核心算法，处于整个程序最上层部分，负责列主元高斯消去法、LU三角分解法和雅克比迭代法的具体实现过程。通过调用方程文件的函数，将获得的数据进行处理运算，可以得到结果返回给方程主函数和输出的第二层。总结：主函数负责获取方程系数并显示，算法和方程作为后台程序，顺序表作为存储手段。3 代码3.1主程序部分// Solutionoflinearquations.cpp : 定义控制台应用程序的入口点。//#include stdio.h#include stdlib.h#include stdafx.h#include squencelist.h#include equation.h#include algorithm.h#includeiostreamint _tmain(int argc, _TCHAR* argv[]){while (Exflag){GetEquation();ShowMenu();}return 0;}3.2多项式方程部分方程部分头文件#ifndef _EQUATION_H#define _EQUATION_H#include stdio.h#include stdlib.h#include squencelist.hextern int Xnumbers;extern int Fnumber;extern int Exflag;extern datacoa *A;void GetEquation(void);void ShowMenu(void);void printfunction(datacoa *A);void printresx(datacoa *A);void Tip(void);#endif方程部分CPP文件#include stdafx.h#include equation.h#include math.h#include algorithm.h#include stdio.h#include stdlib.h#includeiostream#include iomanipusing namespace std;//全局变量int Xnumbers = 0;int Fnumber = 0;int Exflag = 1;datacoa *A;////////////////////////多项式函数系数/////////////////////////void GetEquation(void){int i, j,flag=1;float x;A = InitStruct();while (flag){cout 方程未知量和解总个数： endl;cin Xnumbers;cout 方程个数： endl;cin Fnumber;cout 输入方程系数,输入00结束（如输入 2 1 5 或者2 1 5 3 4 1 00 endl;cout 3 4 1 00 ）: endl;cin x;while (x != 00){for (i = 1; i = Fnumber; i++){for (j = 1; j = Xnumbers; j++){if (!Insert(A, x, i, j))exit(0);cin x;}}}j = 1;printfunction(A);if (Xnumbers == Fnumber+1)flag = 0;else{cout 方程可能无解 endl;A-m = 0;A-n = 0;}}