From f35ea80659b6a49d1bb2ce1d7d002583f3f40947 Mon Sep 17 00:00:00 2001 From: prashantsinalkar Date: Tue, 10 Oct 2017 12:38:01 +0530 Subject: updated the code --- 40/CH8/EX8.24/Exa_8_24.sce | 48 +++++++++++++++++++++++----------------------- 1 file changed, 24 insertions(+), 24 deletions(-) (limited to '40/CH8/EX8.24/Exa_8_24.sce') diff --git a/40/CH8/EX8.24/Exa_8_24.sce b/40/CH8/EX8.24/Exa_8_24.sce index 5ad54e146..4cf5fd699 100755 --- a/40/CH8/EX8.24/Exa_8_24.sce +++ b/40/CH8/EX8.24/Exa_8_24.sce @@ -1,24 +1,24 @@ -//The importance of Periodic extension -//(a)For M=3 -x=[1 2 1]; -XDFT=dft(x,-1) -w=exp(-%i*2*%pi/3); -for i=1:3 - for j=1:3 - WN(i,j)=w^((i-1)*(j-1)); - end -end -WI=WN'; -xn=1/3*WI*XDFT -//The result is periodic with M=3 & 1 period equals x[n] -//(b)For M=4 -y=[1 2 1 0]; -YDFT=dft(y,-1) -w=exp(-%i*%pi/2); -for i=1:4 - for j=1:4 - WN(i,j)=w^((i-1)*(j-1)); - end -end -WI=WN'; -yn=1/4*WI*YDFT +//The importance of Periodic extension +//(a)For M=3 +x=[1 2 1]; +XDFT=fft(x,-1)' +w=exp(-%i*2*%pi/3); +for i=1:3 + for j=1:3 + WN(i,j)=w^((i-1)*(j-1)); + end +end +WI=WN'; +xn=1/3*WI*XDFT +//The result is periodic with M=3 & 1 period equals x[n] +//(b)For M=4 +y=[1 2 1 0]; +YDFT=fft(y,-1)' +w=exp(-%i*%pi/2); +for i=1:4 + for j=1:4 + WN(i,j)=w^((i-1)*(j-1)); + end +end +WI=WN'; +yn=1/4*WI*YDFT \ No newline at end of file -- cgit