数值分析报：Newton插值与romble积分.docVIP

数值分析实验报告 题目：Newton与Romberg联合上机 姓名： 学号： 班级： 信息安全1003班 指导老师： 完成日期： 6月23号 华中科技大学计算机学院 2012年6月目录 一、 实验题目： - 3 - 二、 实验目地： - 3 - 三、 实验要求： - 3 - 四、 主要公式： - 3 - 1. 牛顿插值： - 3 - 2. Romberg数值积分： - 4 - 五、设计流程: - 5 - 七、测试和运行结果： - 10 - 1、测试数据（表7-1）： - 10 - 2、程序运行结果： - 10 - 3、结果分析： - 11 - 八、上机心得体会： - 11 - 1、周： - 11 - 2、姜： - 11 - 实验题目： Newton插值与Romberg积分方法联合上机。 实验目地： 掌握Newton插值和Romberg积分的计算方法； 学会进行误差分析； 能够通过计算机编程解决实际问题； 实验要求： 动态输入节点个数、各节点数据、积分上下限及精度要求； 利用 Newton 插值法计算求积节点的被积函数值 利用 Romberg 算法计算满足精度要求的数值积分近似值（至少计算）； 以二维表形式输出等分数 n、步长 h、Tn、Sn、Cn、Rn以及最终的计算结果 牛顿插值： 其中： Romberg数值积分： 计算顺序如下： 五、设计流程: 程序分为输入部分、求差分表部分、求差值部分、积分部分。 图5-1 程序流程 六、程序实现： // Newton插值与Romberg数值积分程序 // 作者：姜 2012年6月 // GNU General Public License/BSD. // 包含考滤Runge与不考滤Runge两种算法 #include iostream #include cstdio #include cmath #include cstdlib using namespace std; //若定义RUNGE则使用考虑RUNGE的算法否则直接用插值公式 //#define RUNGE 7 //不考虑RUNG将此行注释掉 #define MAX 100 //插分表数组最大值 #define S2N(t2,t1) (t2*4/3 - t1*1/3) #define C2N(s2,s1) (s2*16/15 - s1*1/15) #define R2N(c2,c1) (c2*64/63 - c1*1/63) #define POS(row) (((row+1)*row)1) double val_x[MAX]; //自变量 int num; //记录输入数据个数 double between; //两点间距 //利用一维数组存储二维下三角，采用类实现 class sq_cow{ //插分表列 private: int row; double valy[POS(MAX)];//记录报插值表 friend class sq_row; void set_row(int y){row = y;} public: double operator[](int x){return valy[POS(row)+x];}//重载[] }; class sq_row{ //插分表行 public: sq_cow store; sq_cow operator[](int x){store.set_row(x); return store;} }val_y;//现在可以对val_y用访问二维数组的方式访问 void tab_cha(int num) // 求插分表 { int j; cin val_x[0] val_y[0][0]; for (int i = 1;i num;i++) { cin val_x[i] val_y[i][0]; for (j = 1;j = i;j++) val_y[i][j] = (val_y[i][j - 1] - val_y[i - 1][j - 1]) / (val_x[i] - val_x[i - j]); } } #ifdef RUNGE #define TITLE title 考虑Runge现象 struct refer{