2 files changed, 34 insertions, 0 deletions
diff --git a/215/CH18/EX18.5/Figure18_5.jpg b/215/CH18/EX18.5/Figure18_5.jpg
new file mode 100755
index 000000000..ef0a5063a
--- /dev/null
+++ b/215/CH18/EX18.5/Figure18_5.jpg
diff --git a/215/CH18/EX18.5/ex18_5.sce b/215/CH18/EX18.5/ex18_5.sce
new file mode 100755
index 000000000..580a1f6aa
--- /dev/null
+++ b/215/CH18/EX18.5/ex18_5.sce
@@ -0,0 +1,34 @@
+clc;
+//Example 18.5
+//Let amplitude be 1 
+A=1;
+Dt=0.01;
+T1=4;
+t=0:Dt:T1/4;
+for i=1:length(t)
+    xt(i)=A
+end
+//Calculate Fourier Transform
+Wmax=2*%pi*1;
+K=4;
+k=-(2*K):(K/1000):(2*K);
+W=k*Wmax/K;
+xt=xt';
+XW=xt*exp(-sqrt(-1)*t'*W)*Dt;
+XW_Mag=real(XW);
+W=[-mtlb_fliplr(W),W(2:1001)];
+XW_Mag=[mtlb_fliplr(XW_Mag),XW_Mag(2:1001)];
+subplot(2,1,1);
+a=gca();
+a.data_bounds=[0,0;1,1.5];
+a.y_location="origin";
+plot(t,xt);
+xlabel('t in sec.');
+title('v(t)vs t');
+subplot(2,1,2);
+a=gca();
+a.y_location="origin";
+plot(W*%pi/2,abs (XW_Mag));
+xlabel('Freq in rad/sec');
+ylabel('|F(jw)|')
+title('|F(jw)| vs t');
+\ No newline at end of file