fft_test4.asv

自己写的matlab程序

fid=fopen('data.txt');
fnum=562;  %采样点数量
flinecnt=fnum*2; %画几条数据线
fstart=fnum;%采样点起始
ffreqz=33;
step_length=2;

%读文件
A=fscanf(fid,'%d %d\n',[8,fnum]);
A=A';
fclose(fid);

%%%%%%%%%第1至7个球的号码统计
B=zeros(33,8);
for index_i=2:8 %第一列是日期
    for index_j=1:fnum
        B(A(fstart*(index_i-1)+index_j),index_i)=B(A(fstart*(index_i-1)+index_j),index_i)+1;
    end
end
%%%%%%%%%

%B1=fft(A(fstart+1:(fstart+fnum)),fnum);
%B1=abs(B1).^2/fnum;
%B2=fft(A(fstart+fnum+1:(fstart+fnum*2)),fnum);
%B2=abs(B2).^2/fnum;
f=ffreqz*(0:fnum-1)/fnum;
f2=ffreqz*(0:32)/33;

%%%%%%%%%% 总体出现情况统计
BTA=zeros(33,1);
BTT=B';
for index_m=1:33
    for index_k=2:7
        BTA(index_m)=BTA(index_m)+ BTT(index_k,index_m);
    end    
end
%%%%%%%%%%
%画第n个球的号码统计曲线
plot(f2,B(33*1+1:33*2),...
    f2,B(33*2+1:33*3),...
    f2,B(33*3+1:33*4),...
    f2,B(33*4+1:33*5),...
    f2,B(33*5+1:33*6),...
    f2,B(33*6+1:33*7),...
    f2,B(33*7+1:33*8)...
    ...
    ...
    );
%画点
BT=B(1:33,2:4);
stem(f2,BT,'*');

%对第n个球的号码统计做FFT分析
BT1=zeros(33,3);
for index_n=1:33*3
    BT1(index_n)=BT(index_n)-16;
end
BT2=fft(BT1,33);
BT2=abs(BT2).^2/33;
plot(f2,BT2(1:33*1),...
     f2,BT2(33*1+1:33*2),...
     f2,BT2(33*2+1:33*3));

 %%%%%%%对各期的7个号码分析
 DQTMP=A';
 DQFX=DQTMP(2:7,1:562);
for index_n=1:33*3
for index_n=1:33*3
    DQFX(index_n)=DQFX(index_n)-16;
end
end
 DQFXFFT=fft(DQFX,6);
 f6=ffreqz*(0:5)/6;
 DQFXFFT=abs(DQFXFFT).^2/6;
 plot(f6,DQFXFFT(1:6*1),...
     f6,DQFXFFT(6*1+1:6*2),...
     f6,DQFXFFT(6*2+1:6*3));
 
axis([0,ffreqz,0,2000]);
grid;