数字信号，C++程序，FFT算法圆卷积线卷积
需积分: 15 140 浏览量
2010-06-02
14:11:43
上传
评论
收藏 109KB DOC 举报
数字信号上机作业
1.编程实现 N=2^m,例如 N=256 点 FFT
算法
程序如下：
function X=myt(x) %myt函数 用递归实现
N=length(x);
t=log2(N);
t1=!oor(t);
t2=ceil(t);
if t1~=t2; %若x的长度N不为2的整数次幂，则补0至最接近的2的整数次幂
x=[x zeros(1,2^t2-N)];
N=2^t2;
end
w0=exp(-j*2*pi/N);
X=zeros(1,N);
if N==2
X(1)=x(1)+x(2);
X(2)=x(1)-x(2);
else
n=1:N/2;
xe(n)=x(2*n-1);
xo(n)=x(2*n);
XE=myt(xe); %递归调用
XO=myt(xo);
for n=1:N/2
X(n)=XE(n)+XO(n)*(w0^(n-1));
X(n+N/2)=XE(n)-XO(n)*(w0^(n-1));
end
资源评论
