mxulie.zip_m序列_m序列 C++

#include"stdio.h" #include"math.h" void main() { int connections[10]={1,1,1,1,1,1,1,0,0,1}; int registers[10]={0,0,0,0,0,0,0,0,0,1}; int seq[1023]; int m=10; int L=1023; int k,i,j,t; int new_reg_cont[10]; seq[0]=1; for(i=1;i<=L-1;i++) { int temp=0; for(j=0;j<=m-1;j++) { temp=temp+connections[j]*registers[j]; } new_reg_cont[0]=temp%2; for(k=1;k<=m-1;k++) { new_reg_cont[k]=registers[k-1]; } for(t=0;t<=9;t++) registers[t] = new_reg_cont[t]; seq[i]=new_reg_cont[9]; } int save[2045]; for(k=0;k<=1023;k++) { int v=0; for(i=0;i<=k;i++) v+=seq[k-i]*seq[1023-i]; save[k]=v; } for(k=1024;k<=2044;k++) save[k]=save[2046-k]; //} int temp_dx; float PI=3.14; int P,B,K; int temp_re; int temp_im; int FFT_tempnum1r; int FFT_tempnum2r; int FFT_tempnum1i; int FFT_tempnum2i; int LH=512; int Res[1023]; j=LH; for(i=1;i<1023;i++) //倒序最高位加1，逢2向次高位进位 { if(i<j) { temp_dx=seq[i]; seq[i]=seq[j]; seq[j]=temp_dx; } K=512; while((j<K)==0) { j=j-K; K=K/2; } j=j+K; } //对1024点输入序列进行位 for(i=1;i<=10;i++) //第一重循环控制蝶形的级数 { B=pow(2,(i-1)); for(j=0;j<=(B-1);j++) //控制每级的蝶形和旋转因子 { P=pow(2,(i-1))-1; P=pow(2,(10-i))*j; temp_re=cos(2*PI*P/1024); temp_im=sin(2*PI*P/1024); for(k=j;k<=1023;k=k+pow(2,i)) { FFT_tempnum1r=seq[k]; FFT_tempnum2r=seq[k+B]; // FFT_tempnum1i=FFT_Im[k]; // FFT_tempnum2i=FFT_Im[k+B]; seq[k]=FFT_tempnum1r+FFT_tempnum2r*temp_re+FFT_tempnum2i*temp_im; // FFT_Im[k]=FFT_tempnum1i+FFT_tempnum2i*temp_re-FFT_tempnum2r*temp_im; seq[k+B]=FFT_tempnum1r-FFT_tempnum2r*temp_re-FFT_tempnum2i*temp_im; // FFT_Im[k+B]=FFT_tempnum1i-FFT_tempnum2i*temp_re+FFT_tempnum2r*temp_im; } } } for(i=0;i<1024;i++) Res[i]=sqrt(seq[i]*seq[i]); //结果取模 }