数值分析课程作业.doc

《数值分析课程作业.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、数值分析课程作业（1）解答如下：不稳定的数值计算公式往往会出现“差之毫厘，失之千里”的错误结果。因此，在计算过程中要选用稳定的计算公式。在本题的迭代计算中，由I0往后递推的求解In会使误差快速增大，而先近似计算In然后回代计算In-1…I0会使误差逐渐减小，因此使用此计算方法较为稳定。用C语言实现的代码如下：#include/*计算k的阶乘*/intk_factorial(intk){if(0==k)return1;elsereturnk*k_factorial(k-1);}/*计算1/e，k控制计算精度*/doubleget_e_1(intk){

2、if(0==k)return1;elsereturn1.0*((k%2)?-1:1)/k_factorial(k)+get_e_1(k-1);}/**该函数用于计算公式的近似值。对于给定的n值，首先计算(n+k)值对应的*公式的近似值，然后回代计算n值对应的近似值。*/doubleget_value(intn,intk){doubleres;res=(1+get_e_1(7))/((n+k)+1)/2;while(k>0){res=(1-res)/(n+k);k--;}returnres;}intmain(intargc,char*argv[]){intn;prin

3、tf("请输入一个n值：");scanf("%d",&n);if(n<0)return0;printf("对于给定的n=%d，公式的近似值为：%lf",n,get_value(n,9));return0;}（1）解答如下：考虑拟合函数：，将数据表xx1x2……xmf(x)y1x2……ym中的数据代入，得超定方程（m>n）其系数矩阵为由于多项式q0(x)，q1(x)，q2(x)，……，qn(x)在点集{x1，x2，……，xm}上的正交，所以超定方程组的系数矩阵中不同列的列向量是相互正交的向量组。于是用这一矩阵的转置矩阵去左乘超定方程组左、右两端得正规方程组=>其中

4、，，。因为正规方程组中每一个方程都是一元一次方程可以直接写出原超方程组的最小二乘解，所以拟合函数为这一结果与用次多项式拟合所得结果在理论是完全一样的，只是形式上不同、算法实现上避免了解病态方程组。用C语言实现的代码如下：#include#include#include#includeconstdoubleEPS=1E-10;//运算精度/***读入节点数组x、函数值数组y、权值数组ω及节点数N*/voidreadFile(double*&x,double*&y,double*&omiga,int

5、&N){FILE*fp;fp=fopen(".\ZXECF.txt","r");if(fp==NULL){printf("指定位置的文件不存在，请检查！！");getch();exit(0);}for(N=1;;N++){//计算节点个数Nfscanf(fp,"%*lf%*lf%*lf");if(feof(fp))break;}x=newdouble[N];y=newdouble[N];omiga=newdouble[N];rewind(fp);for(inti=0;i

6、i]);fclose(fp);}/***x数组存放节点的x[i]值，y数组存放节点的y[i]，omiga数组各节点对应权值*a数组为拟和的结果多项式的各系数，m为节点个数，n为多项式最高次数*/voidcurveFitting(double*x,double*y,double*omiga,double*a,intm,intn){double*alpha,*beta,*Q,*b,*t,*s,*d,*q;inti,j,k;doublesum1,sum2;alpha=newdouble[n+1];//α数组beta=newdouble[n+1];//β数组b=newdou

7、ble[n+1];//多项式Q(j-1)[x]的系数t=newdouble[n+1];//多项式Q(j)[x]的系数s=newdouble[n+1];//多项式Q(j+1)[x]的系数Q=newdouble[n+1];//临时值d[j]数组d=newdouble[n+1];//临时值d[j]数组q=newdouble[n+1];//临时值q[j]数组for(i=0;i<=n;i++)//n次多项式的n+1个系数a[i]=0;if(n>m)//拟和多项式的最高次数限定不超过mn=m;//构造Q0(x)Q[0]=b[0]=1;d[0]=sum1=sum2=0;for