快速傅里叶变换算法matlab实现

function y=mefft1（x）
    N=length（x）;
    M=floor（log2（N））;
    if 2^M~=N       %如果N不等于2^M时，需要填补0
        M=M+1;
        N1=2^M;
        x（N+1:N1）=0;
        N=N1;
    end
    b=de2bi（0:N-1M）;%转换为2进制 注：matlab中二进制[0 1 1]=6
    b=b（:end:-1:1）;% 逆序
    adr（:1）=bi2de（b）;%2进制转换为10进制    
    x1=x（adr‘+1）;
    for i=1:N/2