资源描述:
《计算方法源程序》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库。
1、计算方法实验报告实验一舍入误差与数值稳定性班级:计算机cO96姓名:郭睿学号:088115实验一舍入误差与数值稳定性目的与要求:1、通过上机编程,复习巩固以前所学程序设计语言;2、通过上机计算,了解舍入误差所引起的数值不稳定性。3、通过上机计算,了解运算次序对计算结果的影响,从而尽量避免大数吃小数的现象。实验内容:•通过止反两个实例的计算,了解利用计算机进行数值计算中舍入误并所引起的数值不稳定性,深入理解初始小的舍入误差可能造成误差积累从而对计算结果的巨大影响。•通过实际编程,了解运算次序对计算结果的影响,了解实数运算符合的结合律和分配律在计算机里不一定
2、成立。1对n=0,1,2,...,20计算定积分算法1利用递推公式nln6-ln5«0.182322算法2利用递推公式人n=20,19,...,l注意到110511112011nr訂x%*占注和x%算法一:#include#include#include#includeHstdlib.hHvoidmain(void){floaty_0=log(6.0)-log(5.0),y_l;intn=l;system(”cis”);/*清屏*/printf(”y[0]=%・20f”,y_0);while(l){y_
3、l=1.0/n-5*y_0;printf(Ky[%d]=%-20f,y_l);/*输出*/if(n>=20)break;y_0=y_l;n++;if(n%3==0)printf(Hn);}getch();/*保持用户屏幕*/运行结杲:Y[0]=0.182322Y[l]=0.088392Y⑵=0.058039Y[3]=0.043128Y[41=0.034310Y[5]=0.028448Y[6]=0.024428Y[7]=0.020719Y[8]=0.021407Y[9]=0.004076Y[10]=0.079618Y[ll]=-0.307181Y[
4、12]=1.619237Y[13]=-8.019263Y[14]=40.167744Y[15]=-200.772069Y[16]=1003.922729Y[17]=-5019.554688Y[18)=25097.828125Y[19]=-125489.085938Y(20
5、=627445.500000算法二#include#include#include#include"stdlib.hHmain()floaty_0=(l/105.0+l/l26.0)/2,y_l;intn=20;system(ncls
6、");/*清屏*/printf(Hy[201=%-20fy_0);while(l){y_l=l/(5.0*n)-y_0/5.0;printf(Ky[%d]=%-20f-1,y_l);if(n<=l)break;y_O=y_l;n-;if(n%3==0)printf(nn);}getch();运行结杲:Y[20]=0.008730Y[19]=0.008254Y[18]=0.008876Y[17]=0.009336Y[16]=0.009898Y[151=0.010520Y[14]=0.011229Y[13]=0.012040Y[12]=0.0129
7、77Y[ll]=0.014071Y[10]=0.015368Y[9]=0.016926Y⑻=0.018837Y[7]=0.021233Y[6]=0.024325Y[5]=0.028468Y[4]=0.034306Y[3]=0.043139Y(2]=0.058039Y(l]=0.088392Y[0
8、=0.182322说明:从计算结果可以看出,算法1是不稳定的,而算法2是稳定的。计算方法实验报告实验二拉格朗日插值与牛顿插值班级:计算机c096姓名:学号:088115实验二拉格朗日插值与牛顿插值一、目的与要求:熟悉拉榕朗日插值多项式和牛顿插值多项式,注意其不同
9、特点;二、实验内容:通过拉格朗口插值和牛顿插值多项式的两个实例的计算,了解两种求解方法,分析他们的优缺点。三、程序与实例算法1.输入X,.,y;.(i=0,l,2,...,n),令L(xM)=0;2.对=0,1,2,.,„n计算•1/W=nn/=0j右X-Xjxi-XjL,kL“+l,(x)y,程序与实例例1已知函数农Xj0.561600.562800.564010.56521y,-0.827410.826590.825770.82495用三次拉格朗日多项式求x=0.5635的函数近似值。拉格朗日插值:#include#include<
10、conio.h>#include#include"stdli