summaryrefslogtreecommitdiff
path: root/716/CH9/EX9.1.III/9_11_ShortQnA_9_1.sce
diff options
context:
space:
mode:
Diffstat (limited to '716/CH9/EX9.1.III/9_11_ShortQnA_9_1.sce')
-rwxr-xr-x716/CH9/EX9.1.III/9_11_ShortQnA_9_1.sce21
1 files changed, 21 insertions, 0 deletions
diff --git a/716/CH9/EX9.1.III/9_11_ShortQnA_9_1.sce b/716/CH9/EX9.1.III/9_11_ShortQnA_9_1.sce
new file mode 100755
index 000000000..d81d75e57
--- /dev/null
+++ b/716/CH9/EX9.1.III/9_11_ShortQnA_9_1.sce
@@ -0,0 +1,21 @@
+//Perform DFT of x(n)=(1 1 -2 -2) and sketch magnitude and phase spectrum
+clc;
+clear;
+N=4;
+xn=[1 1 -2 -2];
+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,'4 point DFT of x(n)=>');
+disp(abs(Xk),'magnitude of 4 point DFT x(n)=>');
+disp(atan(imag(Xk),real(Xk)),'phase of 4 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