数值实验(城市用水量预测问题)课稿.docxVIP

数值分析实验——城市用水量预测问题插值法Lagrange插值插值模型插值函数：，余项：模型应用记2008-2014年1月份的数据如下表：利用C++语言编程对数据进行Lagrange插值，得到第8年（即2015年）1月份的城市用水量为：1674.15C++代码如下：#includeiostreamusing namespace std;int main(){ double x[7]={1.0,2.0,3.0,4.0,5.0,6.0,7.0}; double y[7]={4032.41,4186.025,4186.025,4374.852,4435.234,4505.427,4505.427};doublen,m; double l[7],l1[7]={1.0,1.0,1.0,1.0,1.0,1.0,1.0},l2[7]={1.0,1.0,1.0,1.0,1.0,1.0,1.0}; n=8.0; m=0.0;inti;for(i=0;i7;i++){int j;for(j=0;j7;j++){if(i==j)continue;l1[i]=l1[i]*(n-x[j]);l2[i]=l2[i]*(x[i]-x[j]); } }for(i=0;i7;i++){l[i]=y[i]*(l1[i]/l2[i]); }for(i=0;i7;i++){ m=m+l[i]; }cout2015年1月用水估计为:mendl;return 1;}Newton插值插值模型插值函数：余项：模型应用同样记2008-2014年1月份的数据如下表：利用C++语言编程对数据进行Newton插值，得到第8年（即2015年）1月份的城市用水量为：1674.15C++代码如下：#includeiostreamusing namespace std;int main(){ double x[7]={1.0,2.0,3.0,4.0,5.0,6.0,7.0},w[8]={1.0,1.0,1.0,1.0,1.0,1.0,1.0};doublen,m; double f[7][7]={4032.41,4186.025,4186.025,4374.852,4435.234,4505.427,4505.427};inti,j; n=8;for(j=1;j7;j++){for(i=0;i7-j;i++){f[j][i]=(f[j-1][i+1]-f[j-1][i])/(x[i+j]-x[i]); } }for(i=1;i8;i++){for(j=0;ji;j++){w[i]=w[i]*(n-x[j]); } } m=f[0][0];for(i=1;i7;i++){ m=m+f[i][0]*w[i]; }cout2015年1月用水估计为:mendl;return 1;}拟合法最小二乘法拟合模型拟合函数：系数求解：其中，，，，为权重系数，这里统一取模型应用同样记2008-2014年1月份的数据如下表：下面分别讨论三种函数形态的拟合线性函数：设，此时，，利用C++语言编程对数据进行一次线性函数的最小二乘法拟合，得到第8年（即2015年）1月份的城市用水量为：4647C++代码如下：#includeiostream#includemath.h#includestdlib.husing namespace std;float **a;float *b;float *x;float p;float w;int n;int c;int k=1;void SOR(float xk[]){//SOR算法求解线性方程组inti,j;float t=0.0;floattt=0.0;float *xl;xl=(float *)malloc(sizeof(float)*(n+1));for(i=1;in+1;i++) { t=0.0;tt=0.0;for(j=1;ji;j++) t=t+a[i][j]*xl[j];for(j=i;jn+1;j++)tt=tt+a[i][j]*xk[j];xl[i]=xk[i]+w*(b[i]-t-tt)/a[i][i]; } t=0.0;for(i=1;in+1;i++) {tt=fabs(xl[i]-xk[i]);tt=tt*tt; t+=tt; } t=sqrt(t);for(i=1;in+1;i++)xk[i]=xl[i];if(k+1=ctp) {k++;SOR(xk); }}int main(){double m; double f[7]={4032.41,4186.025,4186.0