1 files changed, 30 insertions, 0 deletions
diff --git a/3523/CH16/EX16.7.6/Ex16_6.sce b/3523/CH16/EX16.7.6/Ex16_6.sce
new file mode 100644
index 000000000..0bfdbc538
--- /dev/null
+++ b/3523/CH16/EX16.7.6/Ex16_6.sce
@@ -0,0 +1,30 @@
+clear all
+clc
+close
+
+C1=0.125*10^-6;//in Farad
+C2=1*10^-9;//in Farad
+R1=360;//in ohms
+R2=544;//in ohms
+theta = sqrt(C1*C2*R1*R2);//in usec
+n = 1/[1+(1+(R1/R2))*(C2/C1)];
+alpha = (R2*C1)/(2*theta*n);
+printf("theta parameter of wave eq %f us \n",theta*10^6)
+printf("n the parameter of circuit eq %f \n",n)
+printf("alpha parameter of circuit eq %f \n",alpha)
+T2 = 10.1*theta;//duration of lightning impulse pulse in us
+T1 = T2/45;//duration of lightning impulse pulse in us
+printf("duration of lightning impulse pulse %f us \n",T2*10^6)
+printf("duration of lightning impulse pulse %f us \n",T1*10^6)
+//answer in the book for T1 is wrong
+
+T = T1/T2;
+printf("generated lighting impulse is %f us \n",T)
+alpha1 = [alpha-sqrt((alpha^2)-1)]/theta;//in us^-1
+alpha2 = [alpha+sqrt((alpha^2)-1)]/theta;//in us^-1
+printf("aplha1 parameter of wave eq is %f us^-1 \n",alpha1*10^-6)
+printf("aplha1 parameter of wave eq is %f us^1 \n",alpha2*10^-6)
+
+//answer in the book is slightly different
+// Now eq of waveform of generated pulse is e(t)=99.75(e^-0.015t - e^-2.77t)
+