1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
|
//Perform 16 point DFT of x(n)=(1/3 1/3 1/3) and sketch magnitude and phase spectrum
clc;
clear;
N=16;
xn=zeros(1,16);
xn(1)=1/3;
xn(2)=1/3;
xn(3)=1/3;
for k=0:1:N-1
Xk(k+1)=0;
for n=0:1:N-1
Xk(k+1)=Xk(k+1)+xn(n+1)*exp(-%i*2*%pi*k*n/N);
end
end
wk=0:1:N-1;
disp(Xk,'16 point DFT of x(n)=>');
disp(abs(Xk),'magnitude of 16 point DFT x(n)=>');
disp(atan(imag(Xk),real(Xk)),'phase of 16 point DFT x(n)=>');
subplot(1,2,1)
plot2d3(wk,abs(Xk),2);
xtitle('Magnitude Spectrum','k','|X(k)|');
subplot(1,2,2)
plot2d3(wk,atan(imag(Xk),real(Xk)),2);
xtitle('Phase Spectrum','k','angle(X(k))');
|
