基二时间抽取的MATLAB代码

2的整数次幂点数的离散傅里叶变换MATLAB代码

%采用基2时间抽取法编写，N必须为2的整数倍
%function [xn_fft]=fft_N（xn）
clc
clear
N=1024;                                                %采样点数
T=0.8;                                              %数据时常
%fs=（N-1）/T;                                         %采样频率
t=linspace（0TN）;
%f1=50;
%f2=150;
%f3=250;
f1=50;
xn=220*sin（2*pi*f1*t）+50*sin（2*pi*3*f1*t）+20*sin（2*pi*5*f1*t）+15*sin（2*pi*7*f1*t）;
%xn（1:2048）=0;
%xn（2049:6144）=1;
%xn（6145:8192）=0;
WN=exp（-j*pi*2/N）;                              %WN表示底数N=2的变换因子
n=N;
num=0;
while n>1
    n=n/2;
    num=num+1;
end
M=num;                                                     %M表示总共要运算的级数，数据长度对于2的指数
double XK_1[];
double XK_2[];
double xn1[];
int num_new[];

num_new（1）=1;
for num=（M-1）:-1:0
    N_new=length（num_new）;
    num_increase=2^num;
    for num_1=1:N_new
        n

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       3429  2010-06-03 22:49  t_2_fft.m

----------- ---------  ---------- -----  ----

                 3429                    1