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

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

我才是李波



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

5come5帮你背单词 [ overestimate /'əuvə'restimeit/ vt. & 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 | [楼 主]
huge-stone



性别: 帅哥 状态: 该用户目前不在线
等级: 鹤立鸡群
家族: Westlife Family
发贴: 1265
威望: 0
浮云: 1874
在线等级:
注册时间: 2007-03-12
最后登陆: 2008-06-29

5come5帮你背单词 [ adolescent /ædə'lesənt/ a. 青少年的,青春期的;n. 青少年 ]


这个不是很难哈
顶端 Posted: 2007-12-25 22:39 | [1 楼]
我来我网·5come5 Forum » 学业有成

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