diff options
Diffstat (limited to '3812/CH5/EX5.1')
-rw-r--r-- | 3812/CH5/EX5.1/5_1.sce | 45 |
1 files changed, 45 insertions, 0 deletions
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('<X(jW)')
+title('Phase Response(CTFT) in Radians')
|