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

---
 67/CH3/EX3.31/example331.sce | 42 ++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 42 insertions(+)
 create mode 100755 67/CH3/EX3.31/example331.sce

(limited to '67/CH3/EX3.31/example331.sce')

diff --git a/67/CH3/EX3.31/example331.sce b/67/CH3/EX3.31/example331.sce
new file mode 100755
index 000000000..56b3b546a
--- /dev/null
+++ b/67/CH3/EX3.31/example331.sce
@@ -0,0 +1,42 @@
+//Example 3.31
+clc ;
+clear;
+R=10^3;
+C=10^-3;
+A=1/(R*C);
+Dt=0.005;
+t=0:Dt:10;
+xt=A*exp(-A*t);
+Wmax=2*%pi*1;
+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)];
+figure
+a=gca();
+a.y_location="origin";
+plot(t,xt);
+xlabel('t in sec.');
+ylabel('x(t)');
+title('Continuous Time Signal');
+figure
+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)');
+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('<X(jW)')
+title('Phase Response (CTFT) in Radians');
\ No newline at end of file
-- 
cgit