怎樣用matlab寫Radix-2 FFT?

wccs

怎樣用matlab寫Radix-2 FFT?

aoe_liuc

顶一下!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

Griswold

挺简单的吧

sigmadelta2

谢谢谢谢谢谢

华夏甲骨文

心机版

路晨晨

呵呵，用下

king7862591

希望對板主有用。

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. 
    

复制代码

hjy_cool

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

nafic

感谢分享学习资料

sunnyffp

感谢分享！