initial commit / add all books

29 files changed, 712 insertions, 0 deletions

diff --git a/671/CH8/EX8.1/8_1.sce b/671/CH8/EX8.1/8_1.sce
new file mode 100755
index 000000000..8ddfe758c
--- /dev/null
+++ b/671/CH8/EX8.1/8_1.sce
@@ -0,0 +1,16 @@
+V1=3300

+f=50

+N1=100

+N2=300

+Z=100+35*%i

+flux_max=V1/(sqrt(2)*%pi*f*N1)

+V2=V1*N2/N1

+I2=V2/Z

+I1=N2/N1*I2

+S=V1*conj(I1)

+Z1=Z*(N1/N2)^2

+

+disp(flux_max)

+disp(I1,I2)

+disp(real(S),imag(S))

+disp(Z1)
\ No newline at end of file
diff --git a/671/CH8/EX8.10/8_10.sce b/671/CH8/EX8.10/8_10.sce
new file mode 100755
index 000000000..41b640a25
--- /dev/null
+++ b/671/CH8/EX8.10/8_10.sce
@@ -0,0 +1,25 @@
+P=50E3

+a=2200/110

+

+/////OC Parameters

+Poc=400

+Ioc=10

+Voc=110

+

+Y0=Ioc/Voc

+Gi=Poc/Voc/Voc

+disp(Gi)

+Bm=sqrt(Y0*Y0-Gi*Gi)

+disp(Bm)

+

+///////////Referred to HV side

+GiHV=Gi/a/a

+disp(GiHV)

+BmHV=Bm/a/a

+disp(BmHV)

+

+//////////Per unit

+GiPU=Gi*P/Voc/Voc

+BmPU=Bm*P/Voc/Voc

+disp(GiPU)

+disp(BmPU)

diff --git a/671/CH8/EX8.11/8_11.sce b/671/CH8/EX8.11/8_11.sce
new file mode 100755
index 000000000..641831991
--- /dev/null
+++ b/671/CH8/EX8.11/8_11.sce
@@ -0,0 +1,40 @@
+P=25000

+V1=2200

+V2=220

+a=V1/V2

+f=50

+r1=2

+r2=0.025

+x1=7

+x2=0.07

+Xm=16000

+

+Im=V2/Xm

+disp("pf=0, 90 degree lag Po=0", Im)

+

+R=r1+a*a*r2

+X=x1+a*a*x2

+Z=sqrt(R*R+X*X)

+Ifl=P/V1

+Vsc=Z*Ifl

+Vsc_percent=Vsc/V1*100

+disp(Vsc,Vsc_percent)

+Im=Vsc/Xm

+Im_percent=Im/Ifl*100

+disp(Im,Im_percent)

+

+Sc_pf=cos(atan(X/R))

+disp(Sc_pf)

+

+////////////Per unit

+ZBpu=1000*2.2*2.2/25

+r1pu=r1/193.6

+r2pu=a*a*r2/193.6

+x1pu=x1/193.6

+x2pu=a*a*x2/193.6

+

+Rpu=r1pu+r2pu

+Xpu=x1pu+x2pu

+

+Xmpu=Xm/193.6

+disp(Rpu,Xpu,Xmpu)

diff --git a/671/CH8/EX8.12/8_12.sce b/671/CH8/EX8.12/8_12.sce
new file mode 100755
index 000000000..3090785ef
--- /dev/null
+++ b/671/CH8/EX8.12/8_12.sce
@@ -0,0 +1,28 @@
+P=10000

+v1=2300

+v2=230

+a=v1/v2

+f=50

+r1=3.96

+r2=0.0396

+x1=15.8

+x2=0.158

+

+R=r1+a*a*r2

+X=x1+a*a*x2

+V2=v2*a

+I=P/V2

+theta=acos(0.80)

+V1=V2+I*(R*cos(theta)+X*sin(theta))

+disp(V1)

+VR=(V1-V2)/V2

+disp(VR)

+

+pf=X/sqrt(R*R+X*X)

+theta2=acos(pf)

+Il=I*(cos(theta)-%i*sin(theta))

+Ic=real(Il*tan(theta2))-imag(Il)

+Rating_Cap=V2*Ic

+disp(Rating_Cap)

+V1=V2

+disp(V1)

diff --git a/671/CH8/EX8.13/8_13.sce b/671/CH8/EX8.13/8_13.sce
new file mode 100755
index 000000000..25e4915c2
--- /dev/null
+++ b/671/CH8/EX8.13/8_13.sce
@@ -0,0 +1,29 @@
+P=100000

+N1=400

+N2=100

+a=N1/N2

+r1=0.3

+r2=0.015

+x1=1.1

+x2=0.055

+V1=2400

+

+R=r1+a*a*r2

+disp(R)

+X=x1+a*a*x2

+disp(X)

+

+I1=P/V1

+pf=0.8

+theta=acos(pf)

+Vd=I1*(R*cos(theta)+X*sin(theta))

+VR=Vd/V1*100

+V2=(V1-Vd)/a

+disp(VR,V2)

+Vd=I1*(R*cos(theta)-X*sin(theta))

+VR=Vd/V1*100

+V2=(V1-Vd)/a

+disp(VR,V2)

+

+pf=cos(atan(R/X))           /////wrong in the book

+disp("leading",pf)

diff --git a/671/CH8/EX8.14/8_14.sce b/671/CH8/EX8.14/8_14.sce
new file mode 100755
index 000000000..d8dbf19db
--- /dev/null
+++ b/671/CH8/EX8.14/8_14.sce
@@ -0,0 +1,40 @@
+P=10000

+v1=2300

+v2=230

+a=v1/v2

+f=50

+r1=3.96

+r2=0.0396

+x1=15.8

+x2=0.158

+

+R=r1+a*a*r2

+X=x1+a*a*x2

+V2=v2*a

+I=P/V2

+pf=0.8

+theta=acos(pf)

+V1=V2+I*(R*cos(theta)+X*sin(theta))

+Pi=75*V1*V1/V2/V2

+Pc=I*I*R

+Pl=Pi+Pc

+P0=P*pf

+effi=P0/(P0+Pl)*100

+disp(effi)

+

+V1=V2

+Pi=75

+pf2=X/sqrt(R*R+X*X)

+I=P0/V1/pf2

+Pc=I*I*R

+Pl=Pi+Pc

+P0=P*pf

+effi=P0/(P0+Pl)*100

+disp(effi)

+

+I=sqrt(Pi/R)

+Load=V2*I

+P0=Load*pf

+Pl=2*Pi

+effi_max=P0/(P0+Pl)*100

+disp(effi_max)

diff --git a/671/CH8/EX8.15/8_15.sce b/671/CH8/EX8.15/8_15.sce
new file mode 100755
index 000000000..a7bc6ca0d
--- /dev/null
+++ b/671/CH8/EX8.15/8_15.sce
@@ -0,0 +1,62 @@
+P=15E3

+a=3000/250

+

+/////OC Parameters

+Poc=105

+Ioc=0.62

+Voc=250

+

+/////SC Parameters

+Psc=360

+Isc=5.2

+Vsc=157

+

+

+Y0=Ioc/Voc

+Gi=Poc/Voc/Voc

+Bm=sqrt(Y0*Y0-Gi*Gi)

+

+Z=Vsc/Isc

+R=Psc/Isc/Isc

+X=sqrt(Z*Z-R*R)

+

+///////////Referred to HV side

+GiHV=Gi/a/a

+disp(GiHV)

+BmHV=Bm/a/a

+disp(BmHV)

+disp(R)

+disp(X)

+

+//////////Per unit

+ZB=3.12/5.2

+GiPU=Gi*ZB

+BmPU=Y0*ZB

+RPU=R/ZB/1000

+XPU=X/ZB/1000

+disp(GiPU)

+disp(BmPU)

+disp(RPU)

+disp(XPU)

+

+

+pf=0.8

+theta=acos(pf)

+V2=250*a

+I=P/V2

+Vd=I*(R*cos(theta)-X*sin(theta))

+VR=-Vd/V2*100

+disp(VR)

+Pi=105

+Pc=I*I*R

+Pl=Pi+Pc

+P0=P*pf

+effi=P0/(Pl+P0)*100

+

+pf2=cos(atan(R/X))

+disp("leading", pf2)

+

+

+IPU=sqrt(Pi/P/RPU)

+effi_max=IPU

+disp(effi_max)

diff --git a/671/CH8/EX8.16/8_16.sce b/671/CH8/EX8.16/8_16.sce
new file mode 100755
index 000000000..d18e0fd19
--- /dev/null
+++ b/671/CH8/EX8.16/8_16.sce
@@ -0,0 +1,20 @@
+v1=2500

+v2=250

+P=25000

+Pc=130

+Pcfl=320

+pf=0.8

+

+V2=2750

+

+P0=P*pf

+Pl=Pcfl+Pc

+effi=P0/(P0+Pl)*100

+disp(effi)

+

+I2=P/v2

+I1=I2+P/v1

+kVA=I1*v1/1000

+P0=kVA*pf*1000

+effi=P0/(P0+Pl)*100

+disp(effi)

diff --git a/671/CH8/EX8.17/8_17.sce b/671/CH8/EX8.17/8_17.sce
new file mode 100755
index 000000000..0919b7f9d
--- /dev/null
+++ b/671/CH8/EX8.17/8_17.sce
@@ -0,0 +1,14 @@
+a=10

+Load=120000

+V=400

+

+I=Load/V/sqrt(3)

+Vll=a*V/sqrt(3)

+Ill=sqrt(3)*I/a

+all=Vll/V

+disp(Vll,Ill,all)

+

+Vll=a*V*sqrt(3)

+Ill=I/a/sqrt(3)

+all=Vll/V

+disp(Vll,Ill,all)

diff --git a/671/CH8/EX8.18/8_18.sce b/671/CH8/EX8.18/8_18.sce
new file mode 100755
index 000000000..5e8e27aba
--- /dev/null
+++ b/671/CH8/EX8.18/8_18.sce
@@ -0,0 +1,18 @@
+f=50

+N1=500

+

+Pin=60

+Io=0.4

+Vin=220

+r=0.8

+

+Pci=Io*Io*r

+Pi=Pin-Pci

+disp(Pi)

+theta=acos(Pin/Vin/Io)

+Im=Io*sin(theta)

+Xm=Vin/Im

+disp(Xm)

+Iio=Io*cos(theta)

+Ri=Vin/Iio

+disp(Ri)

diff --git a/671/CH8/EX8.19/8_19.sce b/671/CH8/EX8.19/8_19.sce
new file mode 100755
index 000000000..94afd024b
--- /dev/null
+++ b/671/CH8/EX8.19/8_19.sce
@@ -0,0 +1,38 @@
+P=15E3

+v1=2200

+v2=220

+a=v1/v2

+

+/////OC Parameters

+Poc=185

+Ioc=2.72

+Voc=220

+

+/////SC Parameters

+Psc=197

+Isc=6.3

+Vsc=112

+

+Pi=Poc

+disp(Pi)

+

+IHVfl=P/v1

+Pcfl=IHVfl*IHVfl/Isc/Isc*Psc

+disp(Pcfl)

+

+pf=0.85

+Po=P*pf

+Pl=Pi+Pcfl

+effi=Po/(Po+Pl)*100

+disp(effi)

+

+Z=Vsc/Isc

+R=Psc/Isc/Isc

+X=sqrt(Z*Z-R*R)

+pf=0.8

+theta=acos(pf)

+Vd1=IHVfl*(R*cos(theta)+X*sin(theta))

+Vd2=IHVfl*(R*cos(theta)-X*sin(theta))

+VR1=Vd1/v1*100

+VR2=Vd2/v1*100

+disp(VR1, VR2)

diff --git a/671/CH8/EX8.2/8_2.sce b/671/CH8/EX8.2/8_2.sce
new file mode 100755
index 000000000..183e9a44d
--- /dev/null
+++ b/671/CH8/EX8.2/8_2.sce
@@ -0,0 +1,21 @@
+a=2400/240

+

+R=0.2+a*a*2/1000

+X=0.6+a*a*6/1000

+disp(R,X)

+

+Rlv=1/a/a*R

+Xlv=1/a/a*X

+disp(Rlv,Xlv)

+

+I2=150*1000/2400

+Z=R+%i*X

+V=I2*Z

+pcnt=norm(V)/2400*100

+disp(V,pcnt)

+

+Im=2400/1600

+Ii=2400/10000

+Io=Ii-%i*Im

+pf=cos(atan(imag(Io)/real(Io)))

+disp(Io,pf)
\ No newline at end of file
diff --git a/671/CH8/EX8.20/8_20.sce b/671/CH8/EX8.20/8_20.sce
new file mode 100755
index 000000000..147e4c143
--- /dev/null
+++ b/671/CH8/EX8.20/8_20.sce
@@ -0,0 +1,21 @@
+P=50000

+lr=0.9

+

+P0=P*1*0.9

+effi=0.974

+Pl=(1-effi)/effi*P0

+Pi=Pl/2

+Pcfl=Pi/lr/lr

+

+pf=0.8

+P0=P*pf

+Pl=Pi+Pcfl

+effi=P0/(P0+Pl)*100

+disp(effi)

+

+P0=P/2*lr

+Pl=Pi+Pcfl/2/2

+effi=P0/(P0+Pl)*100

+disp(effi)

+

+/////////calculation mistakes in the book

diff --git a/671/CH8/EX8.21/8_21.sce b/671/CH8/EX8.21/8_21.sce
new file mode 100755
index 000000000..e5f2fb0b7
--- /dev/null
+++ b/671/CH8/EX8.21/8_21.sce
@@ -0,0 +1,13 @@
+P=500E3

+effi=0.95

+

+A=[1,1;1,0.6*0.6]

+Pa=inv(A)*[P*(1-effi)/effi;P*0.6*(1-effi)/effi]

+Pi=Pa(1)

+disp(Pi)

+Pc=Pa(2)

+disp(Pc)

+

+Pl=Pi+0.75*0.75*Pc

+effi=P*0.75/(P*0.75+Pl)

+disp(effi)

diff --git a/671/CH8/EX8.22/8_22.sce b/671/CH8/EX8.22/8_22.sce
new file mode 100755
index 000000000..3aabaf878
--- /dev/null
+++ b/671/CH8/EX8.22/8_22.sce
@@ -0,0 +1,14 @@
+v1=2200

+v2=220

+f=50

+emfperturn=12

+

+N1=floor(v1/emfperturn)

+N2=floor(v2/emfperturn)

+disp(N1)

+disp(N2)

+

+fluxmax=emfperturn/4.44/f

+Bmax=1.5

+Ac=fluxmax/Bmax

+disp(Ac)

diff --git a/671/CH8/EX8.23/8_23.sce b/671/CH8/EX8.23/8_23.sce
new file mode 100755
index 000000000..87547d2c0
--- /dev/null
+++ b/671/CH8/EX8.23/8_23.sce
@@ -0,0 +1,24 @@
+v1=3300

+v2=600

+f=50

+Ac=25/10000

+l=1.2

+

+Bmax=1.2

+fluxmax1=Bmax*Ac

+N1=v1/4.44/f/fluxmax1

+N2=v2/4.44/f/fluxmax1

+disp(N1,N2)

+

+I2=20

+pf=0.8

+I1=N2/N1*I2

+I1=I1*(pf-%i*sin(acos(0.8)))

+disp(I1)

+

+Hmax=250

+ATmax=Hmax*l

+immax=ATmax/N1/sqrt(2)

+Ii=0

+I1=-%i*immax+I1

+disp("lagging", real(I1)/norm(I1), norm(I1))

diff --git a/671/CH8/EX8.24/8_24.sce b/671/CH8/EX8.24/8_24.sce
new file mode 100755
index 000000000..8bcc0c9ea
--- /dev/null
+++ b/671/CH8/EX8.24/8_24.sce
@@ -0,0 +1,11 @@
+P=50000

+v1=2400

+v2=240

+f=50

+a=v1/v2

+Pd=375

+pf=0.4

+

+I0=Pd/v1/pf

+I0_prime=I0*a

+disp(I0_prime,pf)

diff --git a/671/CH8/EX8.25/8_25.sce b/671/CH8/EX8.25/8_25.sce
new file mode 100755
index 000000000..e40e6119e
--- /dev/null
+++ b/671/CH8/EX8.25/8_25.sce
@@ -0,0 +1,10 @@
+v1=220

+v2=110

+z1=0.32+%i*0.85

+z2=0.11+%i*0.27

+a=v1/v2

+

+z=z1+a*a*z2

+IHV=v1/norm(z)

+ILV=IHV*a

+disp(IHV,ILV)

diff --git a/671/CH8/EX8.26/8_26.sce b/671/CH8/EX8.26/8_26.sce
new file mode 100755
index 000000000..d3ceb2a7e
--- /dev/null
+++ b/671/CH8/EX8.26/8_26.sce
@@ -0,0 +1,32 @@
+P=1000000

+v1=11000

+v2=230

+f=50

+Vsc=310

+Psc=5210

+

+Isc=P/v1

+Z=Vsc/Isc

+R=Psc/Isc/Isc

+X=sqrt(Z*Z-R*R)

+

+pf=1

+theta=acos(pf)

+Vd=Isc*(R*cos(theta)+X*sin(theta))

+VR=Vd/v1*100

+V1=v1+Vd

+disp(VR,V1)

+

+pf=0.8

+theta=acos(pf)

+Vd=Isc*(R*cos(theta)+X*sin(theta))

+VR=Vd/v1*100

+V1=v1+Vd

+disp(VR,V1)

+

+pf=0.8

+theta=acos(pf)

+Vd=Isc*(R*cos(theta)-X*sin(theta))    //calculation mistake in the book at this point

+VR=Vd/v1*100

+V1=v1+Vd

+disp(VR,V1)

diff --git a/671/CH8/EX8.27/8_27.sce b/671/CH8/EX8.27/8_27.sce
new file mode 100755
index 000000000..dc6161c75
--- /dev/null
+++ b/671/CH8/EX8.27/8_27.sce
@@ -0,0 +1,25 @@
+P=10000

+v1=2200

+v2=220

+r1=4

+r2=0.04

+x1=5

+x2=0.05

+a=v1/v2

+R=r1+a*a*r2

+X=x1+a*a*x2

+I1=P/v1

+pf=0.8

+V2=v2*a

+theta=acos(pf)

+Vd=I1*(R*cos(theta)+X*sin(theta))

+VR=Vd/v1*100

+V1=v1+Vd

+disp(VR,V1)

+

+pf=cos(atan(X/R))

+disp("leading",pf)

+

+Pr=P*pf

+Q=-P*sin(acos(pf))

+disp(Q,Pr)

diff --git a/671/CH8/EX8.28/8_28.sce b/671/CH8/EX8.28/8_28.sce
new file mode 100755
index 000000000..01c8f91e2
--- /dev/null
+++ b/671/CH8/EX8.28/8_28.sce
@@ -0,0 +1,25 @@
+P=20000

+v1=200

+v2=400

+f=50

+V1=600

+V2=200

+

+a=V1/V2

+disp(a)

+

+I1=P/v2

+VA=V1*I1

+disp(VA)

+

+disp(I1)

+I2=30000/V2

+Is=I2-I1

+disp(I2)

+disp(Is)

+

+VAtrans=v2*I1

+VAcond=30000-VAtrans

+disp(VAtrans)

+disp(VAcond)

+

diff --git a/671/CH8/EX8.29/8_29.sce b/671/CH8/EX8.29/8_29.sce
new file mode 100755
index 000000000..086952265
--- /dev/null
+++ b/671/CH8/EX8.29/8_29.sce
@@ -0,0 +1,17 @@
+P=100000

+v1=11000

+v2=400

+

+Vl=400

+Vp=Vl/sqrt(3)

+

+a=ceil(v1/Vp)

+disp(a)

+

+Il=P/sqrt(3)/Vl

+Ip=Il

+disp(Il,Ip)

+

+Il=P/sqrt(3)/v1

+Ip=Il/sqrt(3)

+disp(Il,Ip)

diff --git a/671/CH8/EX8.3/8_3.sce b/671/CH8/EX8.3/8_3.sce
new file mode 100755
index 000000000..93802f8a9
--- /dev/null
+++ b/671/CH8/EX8.3/8_3.sce
@@ -0,0 +1,19 @@
+P=200E3
+V1=11000
+V2=415
+f=50
+N2=80
+
+N1=V1/V2*N2
+a=N1/N2
+disp(a)
+
+I2=P/V2
+disp(I2)
+I1=I2/a
+disp(I1)
+
+Z2=V2*V2/P
+disp(Z2)
+Z21=Z2*a*a
+disp(Z21)
diff --git a/671/CH8/EX8.4/8_4.sce b/671/CH8/EX8.4/8_4.sce
new file mode 100755
index 000000000..979d53302
--- /dev/null
+++ b/671/CH8/EX8.4/8_4.sce
@@ -0,0 +1,25 @@
+V1=200

+f=50

+N1=150

+Ac=10*5/10000

+

+

+Flux_max=V1/(sqrt(2)*%pi*f*N1)

+Bmax=Flux_max/Ac

+Hmax=250

+l=2*(25+10)+2*(20+10)

+l=l/100

+AT_max=Hmax*l

+Im_max=AT_max/150

+Im_rms=Im_max/sqrt(2)

+disp(Im_rms)

+

+Fe_loss=23000

+

+Cv = 2*(25+2*10)*10*5 + 2*20*10*5

+Cv=Cv/1000000

+Cl=Fe_loss*Cv

+disp(Cl)

+

+Ii=150/V1

+disp(Ii)

diff --git a/671/CH8/EX8.5/8_5.sce b/671/CH8/EX8.5/8_5.sce
new file mode 100755
index 000000000..7df0ae10f
--- /dev/null
+++ b/671/CH8/EX8.5/8_5.sce
@@ -0,0 +1,18 @@
+function [x,y]=polar_to_cart(r,theta)

+    theta=theta/180*%pi

+    x=r*cos(theta)

+    y=r*sin(theta)

+endfunction

+

+[Zr,Zc]=polar_to_cart(5,30)

+Z2=complex(Zr,Zc)

+

+V=200

+N2=75

+N1=100

+

+I2=V/Z2

+I1=N2/N1*I2

+

+disp(norm(I1))

+disp(" lagging" , real(I1)/norm(I1))

diff --git a/671/CH8/EX8.6/8_6.sce b/671/CH8/EX8.6/8_6.sce
new file mode 100755
index 000000000..15c1d6123
--- /dev/null
+++ b/671/CH8/EX8.6/8_6.sce
@@ -0,0 +1,14 @@
+R=60

+V=6

+Rs=2400

+

+a=sqrt(Rs/R)

+disp(a)

+

+Pl_max=1/2*(R*R)/(Rs+Rs)

+I1=V/2/Rs

+

+Il=I1*a

+disp(Il)

+Vl=V/2/a

+disp(Vl)

diff --git a/671/CH8/EX8.7/8_7.sce b/671/CH8/EX8.7/8_7.sce
new file mode 100755
index 000000000..2b44d1c0d
--- /dev/null
+++ b/671/CH8/EX8.7/8_7.sce
@@ -0,0 +1,19 @@
+f=50

+a=1100/220

+Rhv=0.125

+Xhv=0.625

+Rlv=0.005

+Xlv=0.025

+

+Zhv=Rhv+%i*Xhv

+Zlv=Rlv+%i*Xlv

+

+Z1=Zhv+a*a*Zlv

+disp(Z1)

+Z2=Z1/a/a

+disp(Z2)

+

+Zpu1=Z1*0.05/1.1/1.1

+disp(Zpu1)

+Zpu2=Z2*0.05/0.22/0.22

+disp(Zpu2)

diff --git a/671/CH8/EX8.8/8_8.sce b/671/CH8/EX8.8/8_8.sce
new file mode 100755
index 000000000..7e05fe177
--- /dev/null
+++ b/671/CH8/EX8.8/8_8.sce
@@ -0,0 +1,26 @@
+P=600E3

+a=2400/600

+r1=0.05

+r2=0.004

+x1=0.025

+x2=0.016

+R1=1667

+Xm=417

+

+/////As seen from the LV side

+

+Zlv=r2+%i*x2+1/a/a*(r1+%i*x1)

+disp(Zlv)

+

+RiLV=R1/a/a

+disp(RiLV)

+XmLV=Xm/a/a

+disp(XmLV)

+

+

+Zpu=Zlv*0.6/0.6/0.6

+disp(Zpu)

+Ri=RiLV*0.6/0.6/0.6

+disp(Ri)

+Xm=XmLV*0.6/0.6/0.6

+disp(Xm)

diff --git a/671/CH8/EX8.9/8_9.sce b/671/CH8/EX8.9/8_9.sce
new file mode 100755
index 000000000..de3559235
--- /dev/null
+++ b/671/CH8/EX8.9/8_9.sce
@@ -0,0 +1,48 @@
+P=50E3

+a=2200/220

+

+/////OC Parameters

+Poc=405

+Ioc=5

+Voc=220

+

+/////SC Parameters

+Psc=805

+Isc=20.2

+Vsc=95

+

+

+Y0=Ioc/Voc

+Gi=Poc/Voc/Voc

+Bm=sqrt(Y0*Y0-Gi*Gi)

+

+Z=Vsc/Isc

+R=Psc/Isc/Isc

+X=sqrt(Z*Z-R*R)

+

+///////////Referred to HV side

+GiHV=Gi/a/a

+disp(GiHV)

+BmHV=Bm/a/a

+disp(BmHV)

+disp(R)

+disp(X)

+

+////////////Referred to LV side

+disp(Gi)

+disp(Bm)

+RLV=R/a/a

+disp(RLV)

+XLV=X/a/a

+disp(XLV)

+

+

+//////////Per unit

+GiPU=GiHV/0.0103

+BmPU=BmHV/0.0103

+RPU=R/96.8

+XPU=X/96.8

+disp(GiPU)

+disp(BmPU)

+disp(RPU)

+disp(XPU)