计算方法上机作业插值与拟合实验报告.docxVIP

计算方法实验题目：班级：学号：姓名：1 实验目的通过编程对拉格朗日插值法、牛顿插值法以及多项式拟合数据的理解观察上述方法的计算稳定性和求解精度并比较各种方法利弊2 实验步骤2.1环境配置：VS2013，C++控制台程序2.2添加头文件#include stdio.h#include stdlib.h#include stdafx.h2.3主要模块程序一共分成三层，最底层是数据结构部分，负责存储数据，第二层是交互部分，即多项式方程部分，负责输入输出获得数据，最上层是核心的算法部分，负责处理已获得的数据。具体功能如下：数据结构部分数据结构部分是整个程序的最底层，负责存储部分。因方程系数作为数据元素插入和删除操作较少，而顺序表空间利用率大且查看方便，故此程序选用顺序表保存系数。数据结构文件中写的是有关顺序表的所有基本操作以供其他文件调用。本次实验使用列主元高斯消元法作为求解方程组的方法，所以也用了二维顺序表存储数组。综上，数据结构部分文件是前两个试验的文件内容和，稍作修改。常系数微分方程部分多项式方程部分是程序的第二层，内容主要是常系数微分方程导数的计算和显示菜单部分。算法部分算法部分分为两个文件，一个是插值部分，一个是拟合部分。插值部分文件负责有关插值的核心算法，处于整个程序最上层部分，负责拉格朗日插值法和牛顿插值法的具体实现过程。调用方程文件的函数，将获得的数据进行处理运算，将结果返回给方程主函数和输出的第二层。每种方法有两个函数，一个为仅仅实现一次插值的算法，另一个是和方程部分联系的函数，负责交互中想实现的整体的算法。拟合部分文件主要负责多项式拟合的算法实现，因为要用到列主元高斯消去法所以也将此部分算法移入其中。主函数负责获取方程系数并显示，算法和方程作为后台程序，顺序表作为存储手段。3 代码3.1主程序部分// Interpolationandfitting.cpp : 定义控制台应用程序的入口点。//#include stdafx.h#includeequation.h#include stdafx.hint _tmain(intargc, _TCHAR* argv[]){GetEquation();while (Exflag){ShowMenu();}return 0;}3.2多项式方程部分方程部分头文件#ifndef _EQUATION_H#define _EQUATION_H#include squencelist.h#include stdio.h#include stdlib.hextern intNumberx;extern intExflag;extern sequenlist *B;extern sequenlist *D;extern sequenlist *L;void GetEquation(void);void ShowMenu(void);void printres(sequenlist *A);void printfunction2(datacoa *A);void printfunctionf(datacoa *A);void Tip(void);#endif方程部分CPP文件#include stdafx.h#include equation.h#include math.h#include alfitting.h#include alinterpolation.h#include squencelist.h#include stdio.h#includeiostream#include iomanip//全局变量intNumberx=0;intExflag = 1;sequenlist *B;sequenlist *D;sequenlist *L;////////////////////////获得给定数据/////////////////////////void GetEquation(void){int j = 0;datatype x = 0;B = InitList();D = InitList();cout 输入给定数据的个数： endl;cin Numberx;cout 从小到大输入x,输入00结束（如y=x^2+2x+1输入1 2 1 00）: endl;cin x;while (x != 00){for (j = 1; j = Numberx; j++){if (!Insert(B, x, j))exit(0);cin x;}}cin.clear();cout 输入f(x),输入00结束（如y=x^2+2x+1输入1 2 1 00）: endl;cin x;while (x != 00){for (j = 1; j = Numberx; j