2015计算方法课程设计报告1.docVIP

计算方法课程设计报告 指导教师：何丽丽 学生姓名：李仔麒 学生学号 综合实验一 实验名称 龙贝格外推法 实验目的及要求 求抛物线y=x2与x=y2围成的图形面积，精确到10-4。 算法描述 程序设计源代码 #includestdio.h #includemath.h #define e 1e-6 double fa(double x) {return x*x; } double fb(double z) {return sqrt(z); } int main() {double t[100][100],a,b,s; int i,j,n; while(scanf(%lf%lf,a,b)!=EOF) {t[0][0]=(fa(a)+fb(b))*(b-a)/2; i=1; do {s=0; for(n=1;n=pow(2,i-1);n++) {s+=fa(a+(2*n-1)*(b-a)/pow(2,i)); } t[i][0]=t[i-1][0]+(b-a)*s/pow(2,i); for(j=1;j=i;j++) t[i][j]=((pow(4,j)*t[i][j-1])-t[i-1][i-1])/(pow(4,i)-1); }while(fabs(t[i][i]-t[i-1][i-1])e); printf(%.4lf\n,t[i][i]); } } 结果与分析 综合实验二 实验名称 曲线拟合最小二乘法 实验目的及要求 中国人口拟合问题。据统计，近十年中国人口数据如下（单位：亿） 年 2004 2005 2006 2007 2008 2009 2010 2011 2012 2013 人口 12.9988 13.0756 13.1448 13.2129 13.2802 13.3474 13.4100 13.4735 13.5404 13.6072 根据表中数据，预测2013年末中国人口数。 设人口总数为，根据人口理论的马尔萨斯模型，采用指数函数 对数据进行拟合，其中t是年份。 为了计算方便，将上式两边同取对数，得 ， 令或，变换后的拟合函数为 ， 对人口数据取对数计算，得下表 t 2004 2005 2006 2007 2008 2009 2010 2011 2012 2013 y ln12.9988 ln13.0756 ln13.1448 ln13.2129 ln13.2802 ln13.3474 ln13.4100 ln13.4735 ln13.5404 ln13.6072 由这些数据可确定a，b的值，从而估计的值。 算法描述 程序设计源代码 #includestdio.h #includemath.h #includestdlib.h #define e 2.71828 int main() {int i,j,k,m,n,L; double b[100][100],d[100],x[100],y[100],g[100],c[100][100],q,D,t,sum1,sum,sum2,z; while(scanf(%d%d,m,n)!=EOF) {sum1=sum=sum2=0.0; for(i=0;i=n;i++) scanf(%lf,x[i]); for(j=0;j=n;j++) scanf(%lf,y[j]); for(j=0;j=n;j++) y[j]=log(y[j]); for(i=0;i=m;i++) {for(j=0;j=m;j++) {sum=0; for(k=0;k=n;k++) sum+=pow(x[k],j+i); b[i][j]=sum; } sum1=0; for(k=0;k=n;k++) sum1+=y[k]*pow(x[k],i); d[i]=sum1; } for(i=0;i=m;i++) b[i][m+1]=d[i]; for(k=0;k=m-1;k++) {D=b[k][k]; L=k; for(i=k+1;i=m;i++) {if(fabs(D)fabs(b[i][k])) {D=b[i][k]; L=i; } } if(fabs(D)1e-6) break; else {if(k!=L) {for(j=k;j=m+1;j++)