13 files changed, 209 insertions, 0 deletions
diff --git a/Working_Examples/215/CH16/EX16.1/ex16_1.sce b/Working_Examples/215/CH16/EX16.1/ex16_1.sce
new file mode 100755
index 0000000..2a112bd
--- /dev/null
+++ b/Working_Examples/215/CH16/EX16.1/ex16_1.sce
@@ -0,0 +1,14 @@
+clc
+//Example 16.1
+disp('Given')
+disp('L=2.5mH Q0=5 C=0.01uF')
+L=2.5*10^-3; Q0=5; C=0.01*10^-6;
+w0=1/sqrt(L*C)
+printf("w0= %3.1f krad/s \n",w0*10^-3);
+f0=w0/(2*%pi)
+alpha=w0/(2*Q0)
+printf("alpha= %3.1f Np/s \n",alpha);
+wd=sqrt(w0^2-alpha^2)
+printf("wd= %3.1f krad/s \n",wd*10^-3);
+R=Q0/(w0*C)
+printf("R= %3.2f ohm \n",R*10^-3);
+\ No newline at end of file
diff --git a/Working_Examples/215/CH16/EX16.10/Figure16_10.jpg b/Working_Examples/215/CH16/EX16.10/Figure16_10.jpg
new file mode 100755
index 0000000..e3cbf61
--- /dev/null
+++ b/Working_Examples/215/CH16/EX16.10/Figure16_10.jpg
diff --git a/Working_Examples/215/CH16/EX16.10/ex16_10.sce b/Working_Examples/215/CH16/EX16.10/ex16_10.sce
new file mode 100755
index 0000000..323d398
--- /dev/null
+++ b/Working_Examples/215/CH16/EX16.10/ex16_10.sce
@@ -0,0 +1,11 @@
+clc
+//Example 16.10
+s=poly(0,'s')
+h=syslin('c',(10*s)/((1+s)*(s^2+20*s+10000)))
+disp(h)
+fmin=0.01
+fmax=10^4
+scf(1);clf;
+//Calculate Bode plot
+bode(h,fmin,fmax)
+
diff --git a/Working_Examples/215/CH16/EX16.11/ex16_11.sce b/Working_Examples/215/CH16/EX16.11/ex16_11.sce
new file mode 100755
index 0000000..e5676bc
--- /dev/null
+++ b/Working_Examples/215/CH16/EX16.11/ex16_11.sce
@@ -0,0 +1,17 @@
+clc
+//Example 16.11
+disp('Given')
+disp('A high pass filter with cutoff frequency of 3k Hz')
+//Cutoff frequency(wc)=1/(R*C)
+//Let us select some standard value of resistor
+disp('Let R=4.7k ohm')
+fc=3*10^3;R=4.7*10^3;
+wc=2*%pi*fc
+C=1/(R*wc)
+printf("\n C=%3.2f nF ",C*10^9);
+s=poly(0,'s')
+h=syslin('c',(R*C*s)/((1+s*R*C)))
+disp(h)
+HW = frmag(h,512);
+w=0: %pi /511: %pi ;
+plot(w,HW)
diff --git a/Working_Examples/215/CH16/EX16.12/ex16_12.sce b/Working_Examples/215/CH16/EX16.12/ex16_12.sce
new file mode 100755
index 0000000..54fcd0b
--- /dev/null
+++ b/Working_Examples/215/CH16/EX16.12/ex16_12.sce
@@ -0,0 +1,26 @@
+clc
+//Example 16.12
+disp('Given')
+disp('Bandwidth = 1M Hz and high frequency cutoff = 1.1M Hz')
+B=10^6;fH=1.1*10^6
+//B=fH-fL
+fL=fH-B
+printf("Low frequency cutoff fL= %d kHz \n",fL*10^-3);
+wL=2*%pi*fL
+printf("wL= %3.2f krad/s \n",wL*10^-3);
+wH=2*%pi*fH
+printf("wH= %3.3f Mrad/s \n",wH*10^-6);
+//Now we need to find values for R,L and C
+//Let X=1/LC
+B=2*%pi*(fH-fL)
+X=(wH-B/2)^2-(B^2/4)
+disp(X)
+disp('Let L=1H')
+L=1;
+C=1/(L*X)
+disp(C,'C=')
+//B=R/L
+R=L*B
+printf("R= %3.3f Mohm \n",R*10^-6);
+
+
diff --git a/Working_Examples/215/CH16/EX16.13/ex16_13.sce b/Working_Examples/215/CH16/EX16.13/ex16_13.sce
new file mode 100755
index 0000000..a46d15d
--- /dev/null
+++ b/Working_Examples/215/CH16/EX16.13/ex16_13.sce
@@ -0,0 +1,26 @@
+clc
+//Example 16.13
+disp('Given')
+disp('Voltage gain = 40dB and cutoff frequency = 10k Hz')
+Av_dB=40
+Av=10^(Av_dB/20)
+f=10*10^3
+B=2*%pi*f
+//From figure 16.41(a)
+disp('1+Rf/R1=100(Gain)')
+//From figure 16.41(b)
+//The transfer function is 
+disp('V+= Vi*(1/sC)/(1+1/sC)')
+//Combining two transfer functions
+disp('V0 = Vi*(1/sC)/(1+1/sC)*(1+Rf/R1)')
+//The maximum value of the combined transfer function is
+disp('Maximum value is V0 = Vi*(1+Rf/R1)')
+disp('Let R1=1k ohm')
+R1=10^3
+Rf=(Av-1)*R1
+printf("Rf= %d kohm \n",Rf*10^-3);
+disp('C=1 uF')
+C=10^-6
+//B=1/(R2*C)
+R2=1/(C*B)
+printf("R2= %3.2f ohm \n",R2);
+\ No newline at end of file
diff --git a/Working_Examples/215/CH16/EX16.2/ex16_2.sce b/Working_Examples/215/CH16/EX16.2/ex16_2.sce
new file mode 100755
index 0000000..40b43f4
--- /dev/null
+++ b/Working_Examples/215/CH16/EX16.2/ex16_2.sce
@@ -0,0 +1,21 @@
+clc
+//Example 16.2
+disp('Given')
+disp('R=40Kohm L=1H C=1/64 uF w=8.2krad/s')
+R=40*10^3; L=1; C=1/64 *10^-6; w=8.2*10^3;
+//The value of Q0 must be at least 5 
+Q0=5;
+w0=1/sqrt(L*C)
+printf("w0= %3.1f krad/s \n",w0*10^-3);
+f0=w0/(2*%pi)
+B=w0/Q0
+printf("Bandwidth= %3.1f krad/s \n",B*10^-3);
+//Number of half bandwidths be N
+N=2*(w-w0)/B
+disp(N)
+//Admittance Y(s)=(1+i*N)/R
+//Finding the magnitude and angle
+magY=sqrt(1+N^2)/R
+angY=atan(N)*(180/%pi)
+disp(angY,'angY=')
+printf("admittance value=%3.2f uS",magY*10^6)
+\ No newline at end of file
diff --git a/Working_Examples/215/CH16/EX16.3/ex16_3.sce b/Working_Examples/215/CH16/EX16.3/ex16_3.sce
new file mode 100755
index 0000000..f81eeff
--- /dev/null
+++ b/Working_Examples/215/CH16/EX16.3/ex16_3.sce
@@ -0,0 +1,24 @@
+clc
+//Example 16.3
+disp('Given')
+disp('R=10 ohm L=2mH C=200 nF w=48 krad/s vs=100*cos(wt) mV')
+R=10; L=2*10^-3; C=200*10^-9; w=48*10^3;
+vsamp=100;
+w0=1/sqrt(L*C)
+printf("w0= %3.1f krad/s \n",w0*10^-3);
+Q0=w0*L/R
+printf("Q0=%d \n",Q0)
+B=w0/Q0
+printf("Bandwidth= %3.1f krad/s \n",B*10^-3);
+//Number of half bandwidths be N
+N=2*(w-w0)/B
+disp(N)
+//Impedance Z(s)=(1+i*N)*R
+//Finding the magnitude and angle
+magZ=sqrt(1+N^2)*R
+angZ=atan(N)*(180/%pi)
+disp(angZ,'angZ=')
+printf("Equivalent impedance value=%3.2f ohm \n",magZ)
+//Approx current magnitude is
+Iamp=vsamp/magZ
+printf("\n Approx current magnitude= %3.2f mA \n",Iamp);
+\ No newline at end of file
diff --git a/Working_Examples/215/CH16/EX16.4/ex16_4.sce b/Working_Examples/215/CH16/EX16.4/ex16_4.sce
new file mode 100755
index 0000000..3f6628e
--- /dev/null
+++ b/Working_Examples/215/CH16/EX16.4/ex16_4.sce
@@ -0,0 +1,17 @@
+clc
+//Example 16.4
+disp('Given')
+disp('R1=2 ohm R2=3 ohm L=1H C=125mF')
+R1=2;R2=3 ; L=1;C=125*10^-3;
+w0=sqrt(1/(L*C)-(R1/L)^2)
+printf("w0=%d rad/s \n",w0)
+//Input admittance is 1/R2+i*w*C+1/(R+I*w*L)
+Y=1/3+%i/4+1/(2+%i*2)
+printf("Y= %3.4f S \n",Y)
+//Now input impedance at resonance 
+Z=1/Y
+printf("Z= %3.4f ohm \n",Z)
+//Resonant frequency f=1/sqrt(L*C)
+f=1/sqrt(L*C)
+printf("f=%3.2f rad/s \n",f);
+
diff --git a/Working_Examples/215/CH16/EX16.5/ex16_5.sce b/Working_Examples/215/CH16/EX16.5/ex16_5.sce
new file mode 100755
index 0000000..b3b172d
--- /dev/null
+++ b/Working_Examples/215/CH16/EX16.5/ex16_5.sce
@@ -0,0 +1,13 @@
+clc
+//Example 16.5
+disp('Given')
+disp('R=5 ohm L=100mH w=100 rad/s')
+Rs=5; Ls=100*10^-3 ;w=100;
+//Let Xs be the capacitive and inductive reactance 
+Xs=w*Ls
+Q=Xs/Rs
+//As Q is greater than 5 we can approximate as
+Rp=Q^2*Rs
+Lp=Ls
+printf("The parallel equivalent is \n");
+printf("Rp= %d ohm \t Lp=%d mH",Rp,Lp*10^3);
+\ No newline at end of file
diff --git a/Working_Examples/215/CH16/EX16.6/ex16_6.sce b/Working_Examples/215/CH16/EX16.6/ex16_6.sce
new file mode 100755
index 0000000..26c1153
--- /dev/null
+++ b/Working_Examples/215/CH16/EX16.6/ex16_6.sce
@@ -0,0 +1,26 @@
+clc
+//Example 16.6
+disp('Given')
+disp('Km=20 Kf=50')
+Km=20; Kf=50;
+s=poly(0,'s')
+//From figure 16.20(a)
+C=0.05; L=0.5;
+//Performing magnitude as well as frequency scaling simultaneously
+Cscaled =C/(Km*Kf)
+Lscaled = L*Km/Kf
+printf("Scaled values are \n")
+printf("Cscaled =%d uF \t Lscaled =%d mH \n",Cscaled*10^6,Lscaled*10^3)
+//Converting the Laplace transform of the circuit
+//From figure 16.20(c)
+disp('Vin=V1+0.5s*(1-0.2*V1)')
+disp('V1=20/s')
+//On substituting V1 in equation of Vin
+
+Zin=(s^2-4*s+40)/(2*s)
+disp(Zin,'Zin=')
+//Now we need to scale Zin
+//We will multiply Zin by Km and replace s by s/Kf
+Zinscaled=horner(Km*Zin,s/Kf)
+disp(Zinscaled,'Zinscaled')
+
diff --git a/Working_Examples/215/CH16/EX16.8/Figure16_8.jpg b/Working_Examples/215/CH16/EX16.8/Figure16_8.jpg
new file mode 100755
index 0000000..5a397bf
--- /dev/null
+++ b/Working_Examples/215/CH16/EX16.8/Figure16_8.jpg
diff --git a/Working_Examples/215/CH16/EX16.8/ex16_8.sce b/Working_Examples/215/CH16/EX16.8/ex16_8.sce
new file mode 100755
index 0000000..0b01847
--- /dev/null
+++ b/Working_Examples/215/CH16/EX16.8/ex16_8.sce
@@ -0,0 +1,14 @@
+clc
+//Example 16.8
+//From figure 16.26
+disp('Writing the expression for voltage gain')
+disp('Vout/Vin=4000*(-1/200)*(5000*10^8/s)/((5000+10^8/s)*(5000+10^6/20s))')
+//On simplification
+s=poly(0,'s')
+h=syslin('c',(-2*s)/((1+s/10)*(1+s/20000)))
+disp(h)
+fmin=0.01
+fmax=10^7
+scf(1);clf;
+bode(h,fmin,fmax)
+