www.pudn.com > cghost.rar > ZXERCHEN.C

#include "math.h" 
#include "stdio.h" 
void iapcir(x,y,n,a,m,dt) 
int n,m; 
float x[],y[],a[],dt[]; 
{int i,j,k; 
float z,p,c,g,q,d1,d2,s[20],t[20],b[20]; 
for(i=0;i<=m-1;i++)a[i]=0.0; 
if(m>n)m=n; 
if(m>20)m=20; 
z=0.0; 
for(i=0;i<=n-1;i++)z=z+x[i]/(1.0*n); 
b[0]=1.0;d1=1.0*n;p=0.0;c=0.0; 
for(i=0;i<=n-1;i++) 
{p=p+(x[i]-z);c=c+y[i];} 
c=c/d1;p=p/d1; 
a[0]=c*b[0]; 
if(m>1) 
{t[1]=1.0;t[0]=-p; 
d2=0.0;c=0.0;g=0.0; 
for(i=0;i<=n-1;i++) 
{q=x[i]-z-p;d2=d2+q*q; 
c=c+y[i]*q; 
g=g+(x[i]-z)*q*q; 
} 
c=c/d2;p=g/d2;q=d2/d1; 
d1=d2; 
a[1]=c*t[1];a[0]=c*t[0]+a[0]; 
} 
for(j=2;j<=m-1;j++) 
{s[j]=t[j-1]; 
s[j-1]=-p*t[j-1]+t[j-2]; 
if(j>=3) 
for(k=j-2;k>=1;k--) 
s[k]=-p*t[k]+t[k-1]-q*b[k]; 
s[0]=-p*t[0]-q*b[0]; 
d2=0.0;c=0.0;g=0.0; 
for(i=0;i<=n-1;i++) 
{q=s[j]; 
for(k=j-1;k>=0;k--) 
q=q*(x[i]-z)+s[k]; 
d2=d2+q*q;c=c=y[i]*q; 
g=g+(x[i]-z)*q*q; 
} 
c=c/d2;p=g/d2;q=d2/d1; 
d1=d2; 
a[j]=c*s[j];t[j]=s[j]; 
for(k=j-1;k>=0;k--) 
{a[k]=c*s[k]+a[k]; 
b[k]=t[k];t[k]=s[k]; 
} 
} 
dt[0]=0.0;dt[1]=0.0;dt[2]=0.0; 
for(i=0;i<=n-1;i++) 
{q=a[m-1]; 
for(k=m-2;k>=0;k--) 
q=a[k]+q*(x[i]-z); 
p=q-y[i]; 
if(fabs(p)>dt[2])dt[2]=fabs(p); 
dt[0]=dt[0]+p*p; 
dt[1]=dt[1]+fabs(p); 
} 
return; 
} 
main() 
{ 
float y[12]={1,2,3,4,5,6,7,8,9,10,11,12},a[2],dt[3]; 
float x[12]={376.5,378.88,379.89,378.8,384.8,383,402.2,\ 
398.99,405.56,399.88,403.88,399.65}; 
int i; 
for(i=0;i<12;i++) 
{ 
y[i]=250*y[i]/12.0; 
} 
iapcir(x,y,12,a,2,dt); 
printf("a[0]=%f\n",a[0]); 
printf("a[1]=%f\n",a[1]); 
printf("dt[2]=%f\n",dt[2]); 
printf("*******THE END********\n"); 
}