diff options
Diffstat (limited to '269/CH16/EX16.1/ex1.sce')
-rw-r--r-- | 269/CH16/EX16.1/ex1.sce | 33 |
1 files changed, 33 insertions, 0 deletions
diff --git a/269/CH16/EX16.1/ex1.sce b/269/CH16/EX16.1/ex1.sce new file mode 100644 index 000000000..561c4399f --- /dev/null +++ b/269/CH16/EX16.1/ex1.sce @@ -0,0 +1,33 @@ +clc;
+//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 |