1 files changed, 24 insertions, 0 deletions
diff --git a/716/CH9/EX9.2.m/MatLab_Prog9_2.sce b/716/CH9/EX9.2.m/MatLab_Prog9_2.sce
new file mode 100755
index 000000000..93321c76e
--- /dev/null
+++ b/716/CH9/EX9.2.m/MatLab_Prog9_2.sce
@@ -0,0 +1,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))');
+\ No newline at end of file