《数值分析实验报告.doc》由会员分享,可在线阅读,更多相关《数值分析实验报告.doc(21页珍藏版)》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、 数值分析实验报告计 1310 陈万全学号:132852实验一一实验内容:用两种不同的顺序计算 ,分析其误差的变化。6483.1102n二实验目的:1. 通过上机编程,复习巩固以前所学程序设计语言;2. 通过上机计算,了解舍入误差所引起的数值不稳定性。3. 通过上机计算,了解运算次序对计算结果的影响,从而尽量避免大数吃小数的现象。三算法描述取 float 类型数据从大到和从小到大的累加。然后输入结果四源程序:/* 实验一 */#include using namespace std;float shunxi()float sum=0;for(int i=1;i0;i-)sum=sum+1.0/
2、float(i*i);return sum;int main()cout.precision(10);coutusing namespace std;double inser(double x)double sum1,sum2,sum3,sum;double xi5=-3.0,-1.0,1.0,2.0,3.0;double yi5=1.0,1.5,2.0,2.0,1.0;int i;int a,b,c;for (i=0;ixiib=(i+1)%5;c=(i+2)%5;break;sum1=yia*(x-xib)*(x-xic)/(xia-xib)*(xia-xic);sum2=yib*(x-x
3、ia)*(x-xic)/(xib-xia)*(xib-xic);sum3=yic*(x-xib)*(x-xia)/(xic-xib)*(xic-xia);sum=sum1+sum2+sum3;return sum;int main()coutusing namespace std;#define N 6double fiN;int chazhi();double niudun(double x)double xiN=0.3,0.42,0.5,0.58,0.66,0.72;double yiN=1.04403,1.08462,1.11803,1.5603,1.19817,1.23223;doub
4、le sum=0,sumxN;chazhi();sumx0=1;for (int j=1;j0;j-)for (i=0;iusing namespace std;#include #define A 0#define B 1#define N 8double Yi(double xk)double s;if(xk=0)return 0.0;elses=double(sin(xk)/double(xk);return s;double fx() double sumx=0;double ki;for (int k=0;kB)coutusing namespace std;#include #de
5、fine A 0#define B 1#define N 100double Yi(double xk)double s;if(xk=0)return 0.0;elses=double(sin(xk)/double(xk);return s;double fx() double sumx=0;double ki;for (int k=0;kB)cout“输入的 A 小于 B,错误“endl;return 1;double sumx=fx(),sumy=fy(),sum;double h=abs(A-B);sum=(Yi(A)+4*sumx+2*sumy+Yi(B)*h/double(6*N);cout“sinx/x 在“A“到“B“区间的“N“段分值为:“sumendl;return 0;*/复化辛卜生公式