求MATLAB写的FFT源程序！

433309

求MATLAB写的FFT源程序！就是效果相当于调用FFT这个函数。如果能够把仿真模块搭建出来就万分感谢！！

ctwong

直接用matlab里面的fft函数就行了啊

路晨晨

呵呵，谢谢

oyyg

直接是库函数

king7862591

這是 radix-2 FFT/IFFT

1. function  X = radix2(x,inverse)
   
2. x = [1 1 1 1 0 0 0 0 ];
   
3. inverse = 1;
   
4. switch inverse
   
5.     case 1
   
6. xcheck=x;
   
7. p=nextpow2(length(x));
   
8. x=[x zeros(1,(2^p)-length(x))];
   
9. N=length(x);
   
10. M=N/2;                                                                     % Index Helping for Controlling range of indices of the Butterfly at Each Stage (Shrinking)
    
11. for stage=1:log2(N);                                                       % No of times decimation has to occur
    
12.          for index=0:(N/(2^(stage-1))):(N-1);                              % Adjusting the index variations of the butterfly in each stage
    
13.              for n=0:M-1;                                                  % Within a stage, for a given index (single) as reference, develop a local Butterfly Index
    
14.                       a=x(n+index+1)+ x(n+index+M+1);
    
15.                       AA = exp((-j*(2*pi)/N)*(2^(stage-1))*(n));
    
16.                       a = fi(AA);
    
17.                       b=(x(n+index+1)- x(n+index+M+1)).* AA;
    
18.                       x(n+1+index)=a;
    
19.                       x(n+M+1+index)=b;                                    % In place Computation
    
20.              end;
    
21.          end;
    
22. M=M/2;                                                                     % Used for creating Butterfly Pairs (INDEXing the wings) (Shrinkage)
    
23. end;
    
24. X=bitrevorder(x);                                                          % Bit reversing X[k] to obtain X[k] 0<k<N-1
    
25. % Ycheck=fft(xcheck,N)
    
26.                                                                            % Cross Check the answer using inbuilt FFT
    
27.     case 2
    
28. xcheck=x;      
    
29. x=x';
    
30. p=nextpow2(length(x));
    
31. x=[x zeros(1,(2^p)-length(x))];
    
32. N=length(x);
    
33. M=N/2;                                                                     % Index Helping for Controlling range of indices of the Butterfly at Each Stage (Shrinking)
    
34. for stage=1:log2(N);                                                       % No of times decimation has to occur
    
35.          for index=0:(N/(2^(stage-1))):(N-1);                              % Adjusting the index variations of the butterfly in each stage
    
36.              for n=0:M-1;                                                  % Within a stage, for a given index (single) as reference, develop a local Butterfly Index
    
37.                       a=x(n+index+1)+ x(n+index+M+1);
    
38.                       b=(x(n+index+1)- x(n+index+M+1)).*exp((-j*(2*pi)/N)*(2^(stage-1))*(n));
    
39.                       x(n+1+index)=a;
    
40.                       x(n+M+1+index)=b;                                    % In place Computation
    
41.              end;
    
42.          end;
    
43. M=M/2;                                                                     % Used for creating Butterfly Pairs (INDEXing the wings) (Shrinkage)
    
44. end;
    
45. X=bitrevorder(x) ;                                                         % Bit reversing X[k] to obtain X[k] 0<k<N-1
    
46. X = X'/N
    
47. end
    
48. % Ycheck=ifft(xcheck,N)                                                        
    
49.                                                                            % Cross Check the answer using inbuilt IFFT
    
50. end
    
51. 
    

复制代码

victorfan

fft的源程序好难啊

hjy_cool

学习，谢谢！！！！！！！

sunnyffp

回复 5# king7862591


    好人啊，找了很久了

redsox

FFT也可以调用FFTW计算。

redsox

我自己都是直接调用FFTW。