diff options
Diffstat (limited to '67/CH3/EX3.8')
-rwxr-xr-x | 67/CH3/EX3.8/example38.sce | 43 |
1 files changed, 43 insertions, 0 deletions
diff --git a/67/CH3/EX3.8/example38.sce b/67/CH3/EX3.8/example38.sce new file mode 100755 index 000000000..78121960d --- /dev/null +++ b/67/CH3/EX3.8/example38.sce @@ -0,0 +1,43 @@ +clc ;
+close ;
+// Analog S i g n a l
+A =1; // Ampl i tude
+Dt = 0.005;
+t = 0: Dt :10;
+xt = exp(-A*t);
+Wmax = 2* %pi *1; // Analog Fr equency = 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)]; // Omega f rom Wmax to Wmax
+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
+figure
+a = gca ();
+a.y_location ="origin";
+plot (t,xt);
+xlabel ( ' t in sec. ' );
+ylabel ( ' x ( t ) ' )
+title ( ' Continuous Time Signal ' )
+figure
+// Pl o t t i n g Magni tude Re spons e o f CTS
+subplot (2 ,1 ,1);
+a = gca ();
+a.y_location ="origin";
+plot (W, XW_Mag );
+xlabel ( ' Fr equency i n Radians / Seconds> W' );
+ylabel ( ' abs (X(jW) ) ' )
+title ( 'Magni tude Re spons e (CTFT) ' )
+// Pl o t t i n g Phase Reponse o f CTS
+subplot (2 ,1 ,2);
+a = gca ();
+a.y_location = "origin";
+a.x_location = "origin";
+plot (W, XW_Phase *%pi /180) ;
+xlabel ( ' Fr equency in Radians / Seconds> W' );
+ylabel ( '<X(jW) ' )
+title ( ' Phase Re spons e (CTFT) i n Radians ' )
\ No newline at end of file |