main.cpp

使用Visual C++生成Frank,cubic,chu序列

#include "Frank.cpp"
#include <stdio.h>
#include <iostream>

void main()
{
  int i,j,p,m,mSeqLength,sequence[1000];
  int length;
  int num1,num2,num3;

  //complex<double> complex1;
  
  FILE *fp;
  char filename[10];
 

 cout<<"Please input the square root of sequence length:";
 cin>>length;
 cout<<"Please input ACF num:";
 cin>>num1;
 cout<<"Please input CCF num1:";
 cin>>num2;
 cout<<"Please input CCF num2:";
 cin>>num3;
 cout<<"please input the filename that we can write data on:\n";
 scanf("%s",filename);

 if ((fp=fopen(filename,"w"))==NULL)
 {
  printf("cannot open file\n");
  exit(0);
 } 

 Frank(length,sequence);

 long lNumber;
 int *a=new int[length];
 m=coprime(length,a);
 lNumber = (long)m;
 long SeqPhase;  
 SeqPhase = length;


 float fPacf[1000];
 float fPccf[1000];
	//for(i=0;i<length;i++)
	//{
	//	fPacf[i] = i;
	//   fprintf(fp,"%5f",fPacf[i]);
	//}
 length = length*length;
 ACF(num1,length,lNumber,sequence, SeqPhase,fPacf);
 CCF(num2,num3,length,lNumber,sequence, SeqPhase,fPccf);
 fprintf(fp,"Length=%d,the Frank sequence as following:\n",length);

 for (i=0;i<length*m;i++)
 {
	   fprintf(fp,"%3d",sequence[i]);
	 if ((i+1) %length == 0)
	 {
	 fprintf(fp,"\n");
	 }
 }


 fprintf(fp,"\n");
 fprintf(fp,"Periodic ACF of No.%d code is:\n",num1);
 for (i=0;i<length;i++)
	   fprintf(fp,"%10f",fPacf[i]);
 fprintf(fp,"\n");
 fprintf(fp,"Periodic CCF of No.%d and No.%d code is:\n",num2,num3);
 for (i=0;i<length;i++)
	   fprintf(fp,"%10f",fPccf[i]);
 fclose(fp);
}