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 --- 3176/CH4/EX4.4/Ex4_4.sce | 36 ++++++++++++++++++++++++++++++++++++ 1 file changed, 36 insertions(+) create mode 100644 3176/CH4/EX4.4/Ex4_4.sce (limited to '3176/CH4/EX4.4/Ex4_4.sce') diff --git a/3176/CH4/EX4.4/Ex4_4.sce b/3176/CH4/EX4.4/Ex4_4.sce new file mode 100644 index 000000000..352069137 --- /dev/null +++ b/3176/CH4/EX4.4/Ex4_4.sce @@ -0,0 +1,36 @@ +//Ex4_4 +// The Mechanics of Computing the DFT +//Version : Scilab 5.4.1 +// Operating System : Window-xp, Window-7 +//Toolbox: Image Processing Design 8.3.1-1 +//Toolbox: SIVP 0.5.3.1-2 +//Reference book name : Digital Image Processing +//book author: Rafael C. Gonzalez and Richard E. Woods + +clc; +close; +clear; +xdel(winsid())//to close all currently open figure(s). +a=[1 2 4 4]; +//b=fft2(a); +disp('Original Signal:') +disp(a); +M=length(a); +for i=1:M + b(i)=0; + for j=1:M + b(i)=b(i)+(a(j)*exp((-%i*2*%pi*(i-1)*(j-1)/M))); + end +end +disp('DFT of Signal:') +disp(b); + +for i=1:M + d(i)=0; + for j=1:M + d(i)=d(i)+((b(j)*exp((%i*2*%pi*(i-1)*(j-1)/M)))/M); + end +end +disp('IDFT:') +disp(abs(d)); + -- cgit