From 7f60ea012dd2524dae921a2a35adbf7ef21f2bb6 Mon Sep 17 00:00:00 2001
From: prashantsinalkar
Date: Tue, 10 Oct 2017 12:27:19 +0530
Subject: initial commit / add all books

---
 2279/CH5/EX5.7/Ex5_7.sce | 50 ++++++++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 50 insertions(+)
 create mode 100644 2279/CH5/EX5.7/Ex5_7.sce

(limited to '2279/CH5/EX5.7')

diff --git a/2279/CH5/EX5.7/Ex5_7.sce b/2279/CH5/EX5.7/Ex5_7.sce
new file mode 100644
index 000000000..b1a9f8517
--- /dev/null
+++ b/2279/CH5/EX5.7/Ex5_7.sce
@@ -0,0 +1,50 @@
+//Fourier Series coefficients of half-wave rectifier output
+//Assume the period of the signal T=1
+t=-0.5:0.01:1;
+T = 1;
+for i=1:length(t)
+    if t(i)<T/2 then
+        x(i)=sin(2*%pi*t(i));
+    else
+        x(i)=0;
+    end
+end
+k=-10:10;
+for i=1:length(k)
+    if k(i)==1 then
+        ak(i)=1/(4*%i);
+    elseif k(i)==-1
+        ak(i)=-1/(4*%i);
+    else
+        ak(i)=(cos(k(i)*%pi/2)*exp(-k(i)*%pi/2*-%i))/(%pi-(%pi*k(i)*k(i)));
+    end
+end
+
+
+disp("The fourier series coefficients are...")
+disp(ak)
+disp("magnitude of Fourier series coefficient")
+disp(abs(ak))
+//PLotting frequency spectrum 
+subplot(2,1,1)
+plot(k,abs(ak),'.');
+xtitle("Magnitude Spectrum","k","|ak|");
+for i=1:length(k)
+    if k(i)==0 | k(i)==3 | k(i)==-3 | k(i)==-5 |k(i)==5 then
+        phase(i)=0;
+    elseif k(i)==-1 then
+            phase(i)=%pi/2;
+    elseif k(i)==1  then
+             phase(i)=-%pi/2;
+    elseif k(i)==-2 | k(i)==-4
+          phase(i)=%pi;
+    elseif k(i)==2 | k(i)==4
+            phase(i)=-%pi;
+    else
+            phase(i) =  0;
+    end        
+end
+subplot(2,1,2)
+plot(k,phase,'.');
+xtitle("Phase Spectrum","k","angle(ak)");
+disp(phase)
\ No newline at end of file
-- 
cgit