blob: 0bfdbc538705841cdae42ce541693f5039c1c250 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
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)
|
