From 7f60ea012dd2524dae921a2a35adbf7ef21f2bb6 Mon Sep 17 00:00:00 2001 From: prashantsinalkar Date: Tue, 10 Oct 2017 12:27:19 +0530 Subject: initial commit / add all books --- 3812/CH5/EX5.1/5_1.sce | 45 +++++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 45 insertions(+) create mode 100644 3812/CH5/EX5.1/5_1.sce (limited to '3812/CH5/EX5.1/5_1.sce') diff --git a/3812/CH5/EX5.1/5_1.sce b/3812/CH5/EX5.1/5_1.sce new file mode 100644 index 000000000..b8e7ad45e --- /dev/null +++ b/3812/CH5/EX5.1/5_1.sce @@ -0,0 +1,45 @@ +//Example 5.1: +//Continuous Time Fourier Transform of a +clear; +clc; +close; +// Analog Signal +A =1;//Amplitude +Dt = 0.005; +t = 0:Dt:10; +xt=exp(-A*t); +//Continuous-time Fourier Transform +Wmax=2*%pi*1; //Analog Frequency = 1Hz +K=4; +k=0:(K/1000):K; +W=k*Wmax/K; +XW=xt*exp(-sqrt(-1)*t'*W)*Dt; +XW_Mag= abs(XW); +W=[-mtlb_fliplr(W), W(2:1001)]; +XW_Mag=[mtlb_fliplr(XW_Mag),XW_Mag(2:1001)]; +[XW_Phase,db]=phasemag(XW); +XW_Phase=[-mtlb_fliplr(XW_Phase),XW_Phase(2:1001)]; +//Plotting Continuous Time Signal +a=gca(); +a.y_location="origin"; +plot(t,xt); +xlabel('t in sec.'); +ylabel('x(t)') +title('Continuous Time Signal') +//Plotting Magnitude Response of CTS +subplot(2,1,1); +a = gca(); +a.y_location = "origin"; +plot(W,XW_Mag); +xlabel('Frequency in Radians/Seconds---> W'); +ylabel('abs(X(jW))') +title('Magnitude Response (CTFT)') +//Plotting Phase Reponse of CTS +subplot(2,1,2); +a = gca(); +a.y_location = "origin"; +a.x_location = "origin"; +plot(W,XW_Phase*%pi/180); +xlabel('Frequency in Radians/Seconds---> W'); +ylabel('