平均分布,正态分布,一阶滑动和,一阶线性回归 C语言编程.doc

#include stdlib.h #include stdio.h #include time.h #include math.h #includestring.h # define pi 3.1415926 # define sqr 0.707106781//在一阶线性回归出现了参数a double uni[2000]={0};//程序中出现大数组时，很可能导致堆栈溢出，为了避免 double nor[2000]={0};//这个问题，把数组声明为全局变量， double ovlap[1000]; double linreg[1000]; double nor_num[10]; double nor_num_theory[10]={0.0}; double mean( double a[]) { int i; double ever=0.0; for(i=0;i2000;i++) ever+=a[i]/2000.0; return ever; } double std(double a[],double mean) { int i; double stda=0.0; for(i=0;i2000;i++) stda+=(a[i]-mean)*(a[i]-mean)/2000.0; return stda; } double integral(double a,double b) { double i,num=0.0; for(i=a;ib;i+=0.0001) {num+=1/sqrt(2*pi)*exp(-i*i/2)*0.0001;} num=2000*num; return num; } //double B_rela(double a) void main( ) { FILE *fp1=fopen(D:\\data1.txt,w);//用于存放平均分布的相关函数 FILE *fp2=fopen(D:\\data2.txt,w);//用于存放正态分布的相关函数 FILE *fp3=fopen(D:\\data3.txt,w);//用于存放一阶滑动序列的相关函数 FILE *fp4=fopen(D:\\data4.txt,w);//用于存放一阶线性回归的相关函数 FILE *fp=fopen(D:\\data.txt,w); int i,j,k=0,uni_num[10]={0};//检验平均分布 double uni_mean,uni_std; //均匀分布 double nor_mean,nor_std;//正态分布 double ovlap_mean,ovlap_ju,ovlap_std;//一阶滑动序列的平均数，矩，方差 double linreg_mean,linreg_ju,linreg_std;// 一阶线性回归的平均数，矩，方差 double uni_B[21],nor_B[21], ovlap_B[21],linreg_B[21];//相关函数 srand( (unsigned)time( NULL ) ); fprintf(fp,the following are contents of uniform distribution:\n); for( i=0;i2011;i++ ) { uni[i]=rand()/32767.0 ; for(j=0;j=9;j++) { if(i2000(uni[i]=j*0.1)(uni[i](j+1)*0.1)) uni_num[j]++ ; } if(i50) fprintf( fp,%6.4f\t, uni[i]); } fprintf(fp,\n\n); uni_mean=mean(uni); fprintf(fp,the average number of the uniform distribution is:%6.4f\n,uni_mean);//打印平均分布的平均数 uni_std=std(uni,uni_mean); fprintf(fp,the variance of the uniform distribution is :%6.4f\n,uni_std);//打印平均分布的方差 fprintf