From b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b Mon Sep 17 00:00:00 2001 From: priyanka Date: Wed, 24 Jun 2015 15:03:17 +0530 Subject: initial commit / add all books --- 716/CH9/EX9.1.III/9_11_ShortQnA_9_1.sce | 21 +++++++++++++++++++++ 1 file changed, 21 insertions(+) create mode 100755 716/CH9/EX9.1.III/9_11_ShortQnA_9_1.sce (limited to '716/CH9/EX9.1.III/9_11_ShortQnA_9_1.sce') 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 -- cgit