ＲＴＲＴ是实验二中那个自己实现ＦＦＴ算法的程序，我的程序是下面这个，但算２５６点的ＤＦＴ很慢．．．．．．．．．．
求强人看看

function xk=myfft1(x,N)
if (N>length(x))
    x=[x,zeros(1,N-length(x))];
else if (N<length(x))
      error('N must not be lower than the length of x');
      N=length(x);
    end
end
l=ceil(log2(N));
  x=[x,zeros(1,(2^l-N))];
  N=2^l;
  W=exp(-i*2*pi/N);
  x1=zeros(1,N/2);
  x2=zeros(1,N/2);
  for p=1:N/2
      x1(p)=x(2*p-1);
      x2(p)=x(2*p);
  end
  if(N>2)
      for k=1:N/2
        xk1=myfft1(x1,N/2);
        xk2=myfft1(x2,N/2);
        xk(k)=xk1(k)+xk2(k)*(W^(k-1));
        xk(k+N/2)=xk1(k)-xk2(k)*(W^(k-1));
      end
  else
      xk(1)=x(1)+x(2);
      xk(2)=x(1)+W*x(2);
  end
         

没写注释，但还只是简单的嵌套而已．．．．．．．．

这个不是很难哈

顺便把我第一个实验的程序分享了，多指教啊

一．
n1=-20:20;
x1=[zeros(1,20),1,zeros(1,20)];
x2=[zeros(1,20),ones(1,21)];
n2=-25:15;
y1=x1;
n3=-12:28;
y2=x2;
subplot(2,2,1);stem(n1,x1);xlabel('n');ylabel('x1');axis([-30,30,-0.5,1.5]);grid;
subplot(2,2,2);stem(n1,x2);xlabel('n');ylabel('x2');axis([-30,30,-0.5,1.5]);grid;
subplot(2,2,3);stem(n2,y1);xlabel('n');ylabel('y1');axis([-30,30,-0.5,1.5]);grid;
subplot(2,2,4);stem(n3,y2);xlabel('n');ylabel('y2');axis([-30,30,-0.5,1.5]);grid;


二．


n=0:100;
x3=sin(pi*n/16);
x4=sin(n/2);
x5=cos(pi*n/12)+cos(3*pi*n/8);
subplot(2,2,1);stem(n,x3);xlabel('n');ylabel('x3');axis([-10,110,-2,2]);grid;
subplot(2,2,2);stem(n,x4);xlabel('n');ylabel('x4');axis([-10,110,-2,2]);grid;
subplot(2,2,3);stem(n,x5);xlabel('n');ylabel('x5');axis([-10,110,-2,2]);grid;

三．


clear all;
n=0:30;
x6=0.2*(0.8.^n);
x7=exp((-1/12+(i*pi/6)).*n);
figure(1);stem(n,x6);xlabel('n');ylabel('x6');grid;
figure(2);subplot(1,2,1);plot(n,abs(x7));xlabel('n');ylabel('abs(x7)');grid;
          subplot(1,2,2);plot(n,angle(x7));xlabel('n');ylabel('angle(x7)');grid;


 

恩，就不知道是怎么的，６４点的ＤＦＴ都可以算，２５６的就卡死了

谁有DSP答案

我的算得很快啊，也没那么复杂

那把你的程序发上来下嘛，谢了

看起来好复杂啊，不懂......我matlab确实糟糕