我来我网
https://5come5.cn
 
您尚未 登录  注册 | 菠菜 | 软件站 | 音乐站 | 邮箱1 | 邮箱2 | 风格选择 | 更多 » 
 

本页主题: 请教DSP实验的编程 显示签名 | 打印 | 加为IE收藏 | 收藏主题 | 上一主题 | 下一主题

我才是李波



性别: 帅哥 状态: 该用户目前不在线
等级: 品行端正
发贴: 244
威望: 0
浮云: 1169
在线等级:
注册时间: 2007-06-20
最后登陆: 2010-04-25

5come5帮你背单词 [ strategy /'strætid3əi/ n. 战略(学),策略,计谋 ]


请教DSP实验的编程

RTRT是实验二中那个自己实现FFT算法的程序,我的程序是下面这个,但算256点的DFT很慢..........
求强人看看

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
         

没写注释,但还只是简单的嵌套而已........
顶端 Posted: 2007-12-25 22:35 | [楼 主]
t-bag



性别: 帅哥 状态: 该用户目前不在线
等级: 人见人爱
家族: 万人坑恋影部落
发贴: 2480
威望: 0
浮云: 1642
在线等级:
注册时间: 2006-12-22
最后登陆: 2011-05-20

5come5帮你背单词 [ festival /'festəvəl/ n. & a. 节日(的) ]


我的算得很快啊,也没那么复杂
顶端 Posted: 2007-12-25 22:52 | [1 楼]
我来我网·5come5 Forum » 学业有成

Total 0.007514(s) query 6, Time now is:12-30 00:12, Gzip enabled
Powered by PHPWind v5.3, Localized by 5come5 Tech Team, 黔ICP备16009856号