radix 2 disspation in frequency
1 view (last 30 days)
Show older comments
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 Comments
Walter Roberson
on 5 Dec 2012
Are you seeing an error message? If not then what difference do you see between what you are getting and what you expected?
Answers (1)
Azzi Abdelmalek
on 6 Dec 2012
Edited: Azzi Abdelmalek
on 6 Dec 2012
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
0 Comments
See Also
Categories
Find more on Transforms in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!