1 files changed, 71 insertions, 0 deletions
diff --git a/1574/CH1/EX1.8/Signals_Ex_1_8.sce b/1574/CH1/EX1.8/Signals_Ex_1_8.sce
new file mode 100755
index 000000000..c3b7fdf63
--- /dev/null
+++ b/1574/CH1/EX1.8/Signals_Ex_1_8.sce
@@ -0,0 +1,71 @@
+clc
+clear
+close
+//Chapter 1 Signals
+//Example 1.8, page no 24
+//given
+
+T=500e-3,w0=2*%pi/T,d=50e-3,A=10
+t=-d/2:0.01:T-d/2
+t1=-d/2:0.01:d/2
+f1=A
+t2=d/2:0.01:T-(d/2)
+f2=0
+a=1
+disp('The fourier series coeff Fn are:')
+for n=-5:5
+    if n==0 then
+       Fr(a)=1,Fi(a)=0
+    else 
+    fa1=f1.*cos(%pi*n*t1/T)
+    fa2=f2.*cos(%pi*n*t2/T)
+     fb1=f1.*sin(%pi*n*t1/T)
+    fb2=f2.*sin(%pi*n*t2/T)
+    end
+    Fr(a)=1/T*(inttrap(t1,fa1)+inttrap(t2,fa2))
+    Fi(a)=1/T*(inttrap(t1,fb1)+inttrap(t2,fb2))
+    mag(a)=abs(Fr(a)+%i*Fi(a))
+    //  disp(mag(a))
+    disp(Fr(a)-%i*Fi(a))
+     x(1,size(t,2))=0
+    x=x+((Fr(a))-%i*Fi(a)).*(cos(%pi*n*t/T)+%i*sin(%pi*n*t/T))
+    a=a+1
+    
+end
+n=-5:5
+subplot(3,1,1),plot(t,f1)
+xlabel("t", "fontsize", 3);
+ylabel("f(t)", "fontsize", 3, "color", "green");
+subplot(3,1,2),plot2d3(n,mag,2,rect=[-11,0,11,1])// expo fourier series coeff
+xlabel("n", "fontsize", 3);
+ylabel("Coeff Magnitude", "fontsize", 3, "color", "red");
+subplot(3,1,3),plot2d(t,x,style=3),plot2d(-t,x,style=3) // one sided spectrum with T=500ms
+xlabel("w", "fontsize", 3);
+ylabel("Fn", "fontsize", 3, "color", "blue");
+legends(['Sampling Function(T=500mS)';'Sampling Function(T=250mS)'],[3,5],opt="lr")
+T1=T/2
+t=-d/2:0.01:T1-d/2
+t1=-d/2:0.01:d/2
+f1=A
+t2=d/2:0.01:T1-(d/2)
+f2=0
+//The Expo fourier series coeff
+for n=-5:5
+    if n==0 then
+       Fr1(a)=1,Fi1(a)=0
+    else 
+    fr1=f1.*cos(%pi*n*t1/T1)
+    fr2=f2.*cos(%pi*n*t2/T1)
+     fi1=f1.*sin(%pi*n*t1/T1)
+    fi2=f2.*sin(%pi*n*t2/T1)
+    end
+    Fr1(a)=1/T1*(inttrap(t1,fr1)+inttrap(t2,fr2))
+    Fi1(a)=1/T1*(inttrap(t1,fi1)+inttrap(t2,fi2))
+    mag(a)=abs(Fr1(a)+%i*Fi1(a))
+    disp(Fr1(a)-%i*Fi1(a))
+     y(1,size(t,2))=0
+    y=y+((Fr1(a))-%i*Fi1(a)).*(cos(%pi*n*t/T1)+%i*sin(%pi*n*t/T1))
+     a=a+1
+end
+plot2d(t,y, style=5),plot2d(-t,y,style=5)// double sided spectrum with T=250ms
+