radix 2 disspation in frequency

1 次查看(过去 30 天)
ayman osama
ayman osama 2012-12-5
for K=0:7
for N=0:3
if mod(K,2)==0
v(N+1)=(x(N+1)+x(N+5))*exp(-1i*N*K*pi/2);
elseif mod(K,2)==1
v(N+1)=((x(N+1)-x(N+5))*exp(-1i*N*pi/4))*exp(-1i*N*K*pi/2);
end
end
Y(K+1)=sum(v);
end
i assume that is radix 2 disspation in frequency what's wron in it
  2 个评论
Walter Roberson
Walter Roberson 2012-12-5
Are you seeing an error message? If not then what difference do you see between what you are getting and what you expected?
ayman osama
ayman osama 2012-12-5
no i didn'y got an error msg i should have fft of x what i see that the first term in fft is right and other terms are wrong i was using dissipation in frequency teq.

请先登录,再进行评论。

回答(1 个)

Azzi Abdelmalek
Azzi Abdelmalek 2012-12-6
编辑:Azzi Abdelmalek 2012-12-6
x=rand(1,8); % Example
xo=x(2:2:end); % odd part
xe=x(1:2:end); % even part
for N=0:3
for k=0:3
Xo(k+1)=xo(k+1)*exp(-j*2*k*N*pi/4)
Xe(k+1)=xe(k+1)*exp(-j*2*k*N*pi/4)
end
Yo(N+1)=sum(Xo)
Ye(N+1)=sum(Xe)
Y(N+1)=Ye(N+1)+exp(-j*2*N*pi/n)*Yo(N+1)
Y(N+5)=Ye(N+1)-exp(-j*2*N*pi/n)*Yo(N+1)
end

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by