initial commit / add all books

20 files changed, 329 insertions, 0 deletions

diff --git a/1049/CH4/EX4.10/ch4_10.sce b/1049/CH4/EX4.10/ch4_10.sce
new file mode 100755
index 000000000..fd94c06e2
--- /dev/null
+++ b/1049/CH4/EX4.10/ch4_10.sce
@@ -0,0 +1,24 @@
+clear;

+clc;

+clear

+function [I_TAV]=theta(th)

+    n=360/th;

+    I=1;    //supposition

+    I_av=I/n;

+    I_rms=I/sqrt(n);

+    FF=I_rms/I_av;

+    I_rms=35;

+    I_TAV=I_rms/FF;

+endfunction

+disp("when conduction angle=180");

+th=180;

+I_TAV=theta(th);

+printf("avg on current rating=%.3f A",I_TAV);

+disp("when conduction angle=90");

+th=90;

+I_TAV=theta(th);

+printf("avg on current rating=%.1f A",I_TAV);

+disp("when conduction angle=30");

+th=30;

+I_TAV=theta(th);

+printf("avg on current rating=%.4f A",I_TAV);

diff --git a/1049/CH4/EX4.11/ch4_11.sce b/1049/CH4/EX4.11/ch4_11.sce
new file mode 100755
index 000000000..eb01b8c25
--- /dev/null
+++ b/1049/CH4/EX4.11/ch4_11.sce
@@ -0,0 +1,8 @@
+clear;

+clc;

+f=50;    //Hz

+I_sb=3000;

+t=1/(4*f);

+T=1/(2*f);

+I=sqrt(I_sb^2*t/T);    printf("surge current rating=%.2f A",I);

+r=(I_sb/sqrt(2))^2*T;    printf("\nI^2*t rating=%.0f A^2.s",r);
\ No newline at end of file
diff --git a/1049/CH4/EX4.12/ch4_12.sce b/1049/CH4/EX4.12/ch4_12.sce
new file mode 100755
index 000000000..0944ad47a
--- /dev/null
+++ b/1049/CH4/EX4.12/ch4_12.sce
@@ -0,0 +1,10 @@
+clc

+clear

+V_s=300;//V

+R=60;//ohm

+L=2;//H

+t=40*10^-6;//s

+i_T=(V_s/R)*(1-exp(-R*t/L));

+i=.036;//A

+R1=V_s/(i-i_T);

+printf("maximum value of remedial parameter=%.3f kilo-ohm",R1/1000);
\ No newline at end of file
diff --git a/1049/CH4/EX4.16/ch4_16.sce b/1049/CH4/EX4.16/ch4_16.sce
new file mode 100755
index 000000000..f7468c1d4
--- /dev/null
+++ b/1049/CH4/EX4.16/ch4_16.sce
@@ -0,0 +1,7 @@
+clear;

+clc;

+V_p=230*sqrt(2);

+R=1+((1)^-1+(10)^-1)^-1;

+A=V_p/R;

+s=1;    //s

+t_c=20*A^-2*s;    printf("fault clearance time=%.4f ms",t_c*10^3);
\ No newline at end of file
diff --git a/1049/CH4/EX4.17/ch4_17.sce b/1049/CH4/EX4.17/ch4_17.sce
new file mode 100755
index 000000000..928ca7af0
--- /dev/null
+++ b/1049/CH4/EX4.17/ch4_17.sce
@@ -0,0 +1,22 @@
+clear;

+clc;

+V_s=sqrt(2)*230;    //V

+L=15*10^-6;    //H

+I=V_s/L;    //I=(di/dt)_max

+printf("(di/dt)_max=%.3f A/usec",I/10^6);

+R_s=10;  //ohm

+v=I*R_s;    //v=(dv/dt)_max

+printf("\n(dv/dt)_max=%.2f V/usec",v/10^6);

+

+f=50;    //Hz

+X_L=L*2*%pi*f;

+R=2;

+I_max=V_s/(R+X_L);    printf("\nI_rms=%.3f A",I_max);

+disp("when conduction angle=90");

+FF=%pi/sqrt(2);

+I_TAV=I_max/FF;    printf("I_TAV=%.3f A",I_TAV);

+disp("when conduction angle=30");

+FF=3.98184;

+I_TAV=I_max/FF;    printf("I_TAV=%.3f A",I_TAV);

+

+printf("\nvoltage rating=%.3f V",2.75*V_s);//rating is taken 2.75 times of peak working voltage unlike 2.5 to 3 times as mentioned int book.
\ No newline at end of file
diff --git a/1049/CH4/EX4.19/ch4_19.sce b/1049/CH4/EX4.19/ch4_19.sce
new file mode 100755
index 000000000..9a97d72ff
--- /dev/null
+++ b/1049/CH4/EX4.19/ch4_19.sce
@@ -0,0 +1,40 @@
+clear;

+clc;

+T_jm=125;

+th_jc=.15;    //degC/W

+th_cs=0.075;    //degC/W

+

+dT=54;    //dT=T_s-T_a

+P_av=120;

+th_sa=dT/P_av;

+T_a=40;    //ambient temp

+P_av=(T_jm-T_a)/(th_sa+th_jc+th_cs);

+if((P_av-120)<1)

+    disp("selection of heat sink is satisfactory");

+end

+dT=58;    //dT=T_s-T_a

+P_av=120;

+th_sa=dT/P_av;

+T_a=40;    //ambient temp

+P_av=(T_jm-T_a)/(th_sa+th_jc+th_cs);

+if((P_av-120)<1)

+    disp("selection of heat sink is satisfactory");

+end

+

+V_m=sqrt(2)*230;

+R=2;

+I_TAV=V_m/(R*%pi);

+P_av=90;

+th_sa=(T_jm-T_a)/P_av-(th_jc+th_cs);

+dT=P_av*th_sa;

+disp("for heat sink");    printf("T_s-T_a=%.2f degC",dT);    printf("\nP_av=%.0f W",P_av);

+P=(V_m/2)^2/R;

+eff=P/(P+P_av);    printf("\nckt efficiency=%.3f pu",eff);

+

+a=60;    //delay angle

+I_TAV=(V_m/(2*%pi*R))*(1+cosd(a));

+printf("\nI_TAV=%.2f A",I_TAV);

+dT=46;

+T_s=dT+T_a;

+T_c=T_s+P_av*th_cs;    printf("\ncase temp=%.2f degC",T_c);

+T_j=T_c+P_av*th_jc;    printf("\njunction temp=%.2f degC",T_j);
\ No newline at end of file
diff --git a/1049/CH4/EX4.2/ch4_2.sce b/1049/CH4/EX4.2/ch4_2.sce
new file mode 100755
index 000000000..37fabe1f1
--- /dev/null
+++ b/1049/CH4/EX4.2/ch4_2.sce
@@ -0,0 +1,12 @@
+clc

+clear

+P=.5;    //P=V_g*I_g

+V=[0.01:0.01:.25];

+plot2d(P./V);

+xlabel('I_g in A');

+ylabel('V_g in V');

+

+t=poly(0,'t');

+P1=P*t/(2*t);

+printf('average gate power dissipation(in watts)');

+disp(P1);
\ No newline at end of file
diff --git a/1049/CH4/EX4.20/ch4_20.sce b/1049/CH4/EX4.20/ch4_20.sce
new file mode 100755
index 000000000..92e8ecc78
--- /dev/null
+++ b/1049/CH4/EX4.20/ch4_20.sce
@@ -0,0 +1,12 @@
+clear;

+clc;

+T_j=125;    //degC

+T_s=70;    //degC

+th_jc=.16;    //degC/W

+th_cs=.08;    //degC/W

+P_av1=(T_j-T_s)/(th_jc+th_cs);    printf("total avg power loss in thristor sink combination=%.2f W",P_av1);

+

+T_s=60;    //degC

+P_av2=(T_j-T_s)/(th_jc+th_cs);    

+

+inc=(sqrt(P_av2)-sqrt(P_av1))*100/sqrt(P_av1);    printf("\npercentage inc in rating=%.2f",inc);
\ No newline at end of file
diff --git a/1049/CH4/EX4.21/ch4_21.sce b/1049/CH4/EX4.21/ch4_21.sce
new file mode 100755
index 000000000..423ff205c
--- /dev/null
+++ b/1049/CH4/EX4.21/ch4_21.sce
@@ -0,0 +1,26 @@
+clear;

+clc;

+R=25000;

+I_l1=.021;    //I_l=leakage current

+I_l2=.025;

+I_l3=.018;

+I_l4=.016;

+//V1=(I-I_l1)*R;

+//V2=(I-I_l2)*R;

+//V3=(I-I_l3)*R;

+//V4=(I-I_l4)*R;

+//V=V1+V2+V3+V4

+V=10000;

+I_l=I_l1+I_l2+I_l3+I_l4;

+//after solving

+I=((V/R)+I_l)/4;

+R_c=40;

+V1=(I-I_l1)*R;    printf("voltage across SCR1=%.0f V",V1);

+V2=(I-I_l2)*R;    printf("\nvoltage across SCR2=%.0f V",V2);

+V3=(I-I_l3)*R;    printf("\nvoltage across SCR3=%.0f V",V3);

+V4=(I-I_l4)*R;    printf("\nvoltage across SCR4=%.0f V",V4);

+

+I1=V1/R_c;    printf("\ndischarge current through SCR1=%.3f A",I1);

+I2=V2/R_c;    printf("\ndischarge current through SCR2=%.3f A",I2);

+I3=V3/R_c;    printf("\ndischarge current through SCR3=%.3f A",I3);

+I4=V4/R_c;    printf("\ndischarge current through SCR4=%.3f A",I4);

diff --git a/1049/CH4/EX4.22/ch4_22.sce b/1049/CH4/EX4.22/ch4_22.sce
new file mode 100755
index 000000000..a6185dab9
--- /dev/null
+++ b/1049/CH4/EX4.22/ch4_22.sce
@@ -0,0 +1,14 @@
+clear;

+clc;

+V_r=1000;    //rating of SCR

+I_r=200;    //rating of SCR

+V_s=6000;    //rating of String

+I_s=1000;    //rating of String

+disp("when DRF=.1");

+DRF=.1;

+n_s=V_s/(V_r*(1-DRF));    printf("number of series units=%.0f",ceil(n_s));

+n_p=I_s/(I_r*(1-DRF));    printf("\nnumber of parrallel units=%.0f",ceil(n_p));

+disp("when DRF=.2");

+DRF=.2;

+n_s=V_s/(V_r*(1-DRF));    printf("number of series units=%.0f",ceil(n_s));

+n_p=I_s/(I_r*(1-DRF));    printf("\nnumber of parrallel units=%.0f",ceil(n_p));

diff --git a/1049/CH4/EX4.23/ch4_23.sce b/1049/CH4/EX4.23/ch4_23.sce
new file mode 100755
index 000000000..1a678222a
--- /dev/null
+++ b/1049/CH4/EX4.23/ch4_23.sce
@@ -0,0 +1,14 @@
+clear;

+clc;

+V1=1.6;    //on state voltage drop of SCR1

+V2=1.2;    //on state voltage drop of SCR2

+I1=250;    //current rating of SCR1

+I2=350;    //current rating of SCR2

+R1=V1/I1;

+R2=V2/I2;

+I=600;    //current to be shared

+//for SCR1,     I*(R1+R)/(total resistance)=k*I1                    (1)

+//for SCR2,     I*(R2+R)/(total resistance)=k*I2                    (2)

+//(1)/(2)

+R=(R2*I2-R1*I1)/(I1-I2);

+printf("reqd value of resistance=%.3f ohm",R);
\ No newline at end of file
diff --git a/1049/CH4/EX4.25/ch4_25.sce b/1049/CH4/EX4.25/ch4_25.sce
new file mode 100755
index 000000000..03976fb61
--- /dev/null
+++ b/1049/CH4/EX4.25/ch4_25.sce
@@ -0,0 +1,12 @@
+clear;

+clc;

+f=2000;    //Hz

+C=0.04*10^-6;

+n=.72;

+R=1/(f*C*log(1/(1-n)));    printf("R=%.2f kilo-ohm",R/1000);

+V_p=18;

+V_BB=V_p/n;

+R2=10^4/(n*V_BB);    printf("\nR2=%.2f ohm",R2);

+I=4.2*10^-3;    //leakage current

+R_BB=5000;

+R1=(V_BB/I)-R2-R_BB;    printf("\nR1=%.0f ohm",R1);
\ No newline at end of file
diff --git a/1049/CH4/EX4.26/ch4_26.sce b/1049/CH4/EX4.26/ch4_26.sce
new file mode 100755
index 000000000..497433665
--- /dev/null
+++ b/1049/CH4/EX4.26/ch4_26.sce
@@ -0,0 +1,15 @@
+clear;

+clc;

+V_p=18;

+n=.72;

+V_BB=V_p/n;

+I_p=.6*10^-3;

+I_v=2.5*10^-3;

+V_v=1;

+R_max=V_BB*(1-n)/I_p;    printf("R_max=%.2f kilo-ohm",R_max/1000);

+R_min=(V_BB-V_v)/I_v;    printf("\nR_min=%.2f kilo-ohm",R_min/1000);

+

+C=.04*10^-6;

+f_min=1/(R_max*C*log(1/(1-n)));    printf("\nf_min=%.3f kHz",f_min/1000);

+f_max=1/(R_min*C*log(1/(1-n)));    printf("\nf_max=%.2f kHz",f_max/1000);

+

diff --git a/1049/CH4/EX4.27/ch4_27.sce b/1049/CH4/EX4.27/ch4_27.sce
new file mode 100755
index 000000000..1d1b4e777
--- /dev/null
+++ b/1049/CH4/EX4.27/ch4_27.sce
@@ -0,0 +1,28 @@
+clear;

+clc;

+R1=1000;

+C=.5*10^-6;

+f=50;

+w=2*%pi*f;

+V_s=230;

+X_c=1/(w*C);

+v_c=30;

+R=0;

+Z=sqrt((R+R1)^2+X_c^2);

+phi=atand(X_c/(R+R1));

+I1=V_s/(Z*complex(cosd(-phi),sind(-phi)));

+V_c=I1*X_c*complex(cosd(-90),sind(-90));

+a=abs(V_c);    //magnitude of V_c

+b=-atand(imag(V_c)/real(V_c));    //argument of V_c

+//v_c=sqrt(2)*a*sind(a1-b)

+a1=asind(v_c/(sqrt(2)*a))+b;    printf("min angle=%.1f deg",a1);

+

+R=25000;

+Z=sqrt((R+R1)^2+X_c^2);

+phi=atand(X_c/(R+R1));

+I1=V_s/(Z*complex(cosd(-phi),sind(-phi)));

+V_c=I1*X_c*complex(cosd(-90),sind(-90));

+a=abs(V_c);    //magnitude of V_c

+b=-atand(imag(V_c)/real(V_c));    //argument of V_c

+//v_c=sqrt(2)*a*sind(a2-b)

+a2=asind(v_c/(sqrt(2)*a))+b;    printf("\nmax angle=%.2f deg",a2);

diff --git a/1049/CH4/EX4.3/ch4_3.sce b/1049/CH4/EX4.3/ch4_3.sce
new file mode 100755
index 000000000..dc52af731
--- /dev/null
+++ b/1049/CH4/EX4.3/ch4_3.sce
@@ -0,0 +1,9 @@
+clear;

+clc;

+P=.5;    //P=V_g*I_g

+s=130;    //s=V_g/I_g

+I_g=sqrt(P/s);

+V_g=s*I_g;

+E=15;

+R_s=(E-V_g)/I_g;    printf("gate source resistance=%.2f ohm",R_s);

+//Answers have small variations from that in the book due to difference in the rounding off of digits.
\ No newline at end of file
diff --git a/1049/CH4/EX4.4/ch4_4.sce b/1049/CH4/EX4.4/ch4_4.sce
new file mode 100755
index 000000000..5540acace
--- /dev/null
+++ b/1049/CH4/EX4.4/ch4_4.sce
@@ -0,0 +1,16 @@
+clear;

+clc;

+R_s=120;    //slope of load line is -120V/A. This gives gate source resistance

+printf("gate source resistance=%.0f ohm",R_s);

+

+P=.4;    //P=V_g*I_g

+E_s=15;

+//E_s=I_g*R_s+V_g,    after solving this

+//120*I_g^2-15*I_g+0.4=0    so

+a=120;    b=-15;    c=0.4;

+D=sqrt((b^2)-4*a*c);

+I_g=(-b+D)/(2*a);    V_g=P/I_g;

+printf("\ntrigger current=%.2f mA",I_g*10^3);    printf("\nthen trigger voltage=%.3f V",V_g);

+

+I_g=(-b-D)/(2*a);    V_g=P/I_g;

+printf("\ntrigger current=%.2f mA",I_g*10^3);    printf("\nthen trigger voltage=%.2f V",V_g);

diff --git a/1049/CH4/EX4.5/ch4_5.sce b/1049/CH4/EX4.5/ch4_5.sce
new file mode 100755
index 000000000..dad080ca7
--- /dev/null
+++ b/1049/CH4/EX4.5/ch4_5.sce
@@ -0,0 +1,16 @@
+clear;

+clc;

+//V_g=1+10*I_g

+P_gm=5;    //P_gm=V_g*I_g

+//after solving, eqn becomes 10*I_g^2+I_g-5=0

+a=10;    b=1;    c=-5;

+I_g=(-b+sqrt(b^2-4*a*c))/(2*a);

+E_s=15;

+//using E_s=R_s*I_g+V_g

+R_s=(E_s-1)/I_g-10;    printf("reistance=%.3f ohm",R_s);

+

+P_gav=.3;    //W

+T=20*10^-6;

+f=P_gav/(P_gm*T);    printf("\ntriggering freq=%.0f kHz",f/1000);

+

+dl=f*T;    printf("\nduty cycle=%.2f",dl);
\ No newline at end of file
diff --git a/1049/CH4/EX4.6/ch4_6.sce b/1049/CH4/EX4.6/ch4_6.sce
new file mode 100755
index 000000000..1b37acf0a
--- /dev/null
+++ b/1049/CH4/EX4.6/ch4_6.sce
@@ -0,0 +1,13 @@
+clear;

+clc;

+I=.1;

+E=200;

+disp("in case load consists of (a)L=.2H");

+L=.2;

+t=I*L/E;    printf("min gate pulse width=%.0f us",t*10^6);

+disp("(b)R=20ohm in series with L=.2H");

+R=20;

+t=(-L/R)*log(1-(R*I/E));    printf("min gate pulse width=%.3f us",t*10^6);

+disp("(c)R=20ohm in series with L=2H");

+L=2;

+t=(-L/R)*log(1-(R*I/E));    printf("min gate pulse width=%.2f us",t*10^6);

diff --git a/1049/CH4/EX4.8/ch4_8.sce b/1049/CH4/EX4.8/ch4_8.sce
new file mode 100755
index 000000000..807690b24
--- /dev/null
+++ b/1049/CH4/EX4.8/ch4_8.sce
@@ -0,0 +1,9 @@
+clc

+clear

+E_s=16;

+R_s=128;

+P=.5;

+y=poly([P -E_s R_s],'i','coeff');

+a=roots(y);

+printf('trigger current=%.1f mA',a(1)*1000);

+printf('\ntrigger voltage=%.0f V',P/a(1));

diff --git a/1049/CH4/EX4.9/ch4_9.sce b/1049/CH4/EX4.9/ch4_9.sce
new file mode 100755
index 000000000..aa1afc0cb
--- /dev/null
+++ b/1049/CH4/EX4.9/ch4_9.sce
@@ -0,0 +1,22 @@
+clear;

+clc;

+function [I_TAV]=theta(th)

+    I_m=1;    //supposition

+    I_av=(I_m/(2*%pi))*(1+cosd(th));

+    I_rms=sqrt((I_m/(2*%pi))*((180-th)*%pi/360+.25*sind(2*th)));

+    FF=I_rms/I_av;

+    I_rms=35;

+    I_TAV=I_rms/FF;

+endfunction

+disp("when conduction angle=180");

+th=0;

+I_TAV=theta(th);

+printf("avg on current rating=%.3f A",I_TAV);

+disp("when conduction angle=90");

+th=90;

+I_TAV=theta(th);

+printf("avg on current rating=%.3f A",I_TAV);

+disp("when conduction angle=30");

+th=150;

+I_TAV=theta(th);

+printf("avg on current rating=%.3f A",I_TAV);