* SPICE netlist
R1 Vo 0 {RL}
C1 Vo 0 {Cf}
B1 0 Vo I={ks} * I(V_Its_sns)
B2 0 Vo I={kc} * I(V_Itc_sns)
Vfb N012 0 AC 1
E1 N008 0 V_omega 0 {Irc * Lr}
E2 0 N013 V_omega 0 {Irs * Lr}
L1 N008 N009 {Lr}
B3 N007 N011 V=2* {ks} * V(Vo)
L2 N013 N014 {Lr}
B4 N019 N017 V= 2 * {kc} * V(Vo)
R2 N011 0 {Rs}
R3 0 N017 {Rc}
V_Its_sns N001 N002 0
B5 N005 N009 V={Omega_s * Lr} * I(V_Irc_sns)
B6 N014 N015 V={Omega_s * Lr} * I(V_Irs_sns)
C2 N003 N005 {Cr}
B7 N003 N005 I={Omega_s * Cr } * V(Vrc_sns)
B8 N003 N005 I={Vrc * Cr } * V(V_omega)
C3 N015 N016 {Cr}
B9 N015 N016 I={Omega_s * Cr} * V(Vrs_sns)
B10 N015 N016 I={Vrs * Cr } * V(V_omega)
B11 N002 N007 V={krs} * I(V_Itc_sns)
B12 N023 N019 V= {krc} * I(V_Its_sns)
V_Itc_sns N021 N023 0
E3 Vrs_sns 0 N005 N003 1
E4 Vrc_sns 0 N015 N016 1
L3 N004 N006 {Lm}
B13 N010 N006 V={Omega_s * Lm} * I(V_Imc_sns)
L4 N020 N022 {Lm}
B15 N020 N018 V={Omega_s * Lm} * I(V_Ims_sns)
V_Imc_sns N021 N022 0
V_Ims_sns N001 N004 0
V_Irs_sns N003 N001 0
V_Irc_sns N016 N021 0
B14 0 N010 V={Ims * Lm} * V(V_omega)
B16 N018 0 V={Imc * Lm} * V(V_omega)
E5 V_omega 0 N012 0 {2*pi*F0}
.param Lr = 14u
.param Lm = 60u
.param Cr = 30n
.param n = 4
.param Cf = 660u
.param RL = 2.3
.param Vg = 400
.param F0 = 1/(2*pi*sqrt(Lr*Cr))
.param Z0 = sqrt(Lr/Cr)
.param Req = 8* (n**2) * RL/(pi**2)
.param Qeq = Z0 / Req
.param fs = 1.2* F0
.param Omega_s =2*pi*fs
.param Xeq = Omega_s*Lr - 1/Omega_s/Cr
.param Ln = Lm/Lr
.param fn = fs/F0
.param Its = (2*Vg/pi/Req) * Ln * (fn**2) * (Ln + 1 - 1/(fn**2)) / ((fn**2) * ((Ln + 1 - 1/(fn**2))**2) + (Qeq*Ln*((fn**2) -1))**2)
.param Itc = (2*Vg/pi/Req) * (Ln**2) * fn * Qeq * (1 - (fn**2)) / ((fn**2) * ((Ln + 1 - 1/(fn**2))**2) + (Qeq*Ln*((fn**2) -1))**2)
.param Ims = Itc * Req / Omega_s / Lm
.param Imc = - Its * Req / Omega_s / Lm
.param Irs = Its + Ims
.param Irc = Itc + Imc
.param Vrs = Irc / Omega_s / Cr
.param Vrc = - Irs / Omega_s / Cr
.param ks = (2*n/pi) * Its/sqrt((Its**2) + (Itc**2))
.param kc = (2*n/pi) * Itc/sqrt(Its**2 + Itc**2)
.param Vo = (pi/4/n) * Req * sqrt(Its**2 + Itc**2)
.param Rs = (4*n*Vo/pi) * (Itc**2) / ((Its**2 + Itc**2)**1.5)
.param Rc = -(4*n*Vo/pi) * (Its**2) / ((Its**2 + Itc**2)**1.5)
.param krs = -(4*n*Vo/pi) * (Its * Itc) / ((Its**2 + Itc**2)**1.5)
.param krc = krs
.ac dec 100 10 100k

* LTspice schematic (.asc)
Version 4
SHEET 1 2864 1232
WIRE 536 -616 136 -616
WIRE 776 -616 616 -616
WIRE -40 -552 -104 -552
WIRE 136 -552 136 -616
WIRE 136 -552 40 -552
WIRE 200 -552 136 -552
WIRE 360 -552 280 -552
WIRE -360 -536 -456 -536
WIRE -232 -536 -280 -536
WIRE -232 -520 -232 -536
WIRE -344 -488 -472 -488
WIRE -200 -488 -296 -488
WIRE -104 -488 -104 -552
WIRE -104 -488 -200 -488
WIRE 360 -488 360 -552
WIRE 776 -456 776 -616
WIRE -472 -424 -472 -488
WIRE -360 -424 -472 -424
WIRE -200 -424 -200 -488
WIRE -200 -424 -296 -424
WIRE 360 -344 360 -408
WIRE 776 -328 776 -376
WIRE -472 -312 -472 -424
WIRE -360 -312 -472 -312
WIRE -200 -312 -200 -424
WIRE -200 -312 -280 -312
WIRE -920 -184 -1016 -184
WIRE -728 -184 -840 -184
WIRE -472 -184 -472 -312
WIRE -472 -184 -648 -184
WIRE -360 -184 -472 -184
WIRE -200 -184 -200 -312
WIRE -200 -184 -280 -184
WIRE 1032 -184 920 -184
WIRE 1272 -184 1032 -184
WIRE 1384 -184 1272 -184
WIRE 1432 -184 1384 -184
WIRE 360 -168 360 -264
WIRE -1016 -136 -1016 -184
WIRE 776 -136 776 -248
WIRE 1032 -136 1032 -184
WIRE -1304 -120 -1432 -120
WIRE -1144 -120 -1256 -120
WIRE -1064 -120 -1144 -120
WIRE 1384 -120 1384 -184
WIRE 920 -104 920 -184
WIRE 1272 -104 1272 -184
WIRE -1432 -72 -1432 -120
WIRE -1256 -72 -1256 -120
WIRE -1064 -72 -1096 -72
WIRE -1304 -56 -1304 -120
WIRE -1016 -24 -1016 -56
WIRE -424 -24 -1016 -24
WIRE 360 -24 360 -88
WIRE 360 -24 -424 -24
WIRE 776 -24 776 -56
WIRE 776 -24 360 -24
WIRE -424 8 -424 -24
WIRE 776 8 776 -24
WIRE -1016 24 -1016 -24
WIRE -1144 40 -1144 -120
WIRE -1064 40 -1144 40
WIRE 920 40 920 -24
WIRE 1032 40 1032 -56
WIRE 1032 40 920 40
WIRE 1272 40 1272 -40
WIRE 1272 40 1032 40
WIRE 1384 40 1384 -40
WIRE 1384 40 1272 40
WIRE 360 72 360 -24
WIRE 920 72 920 40
WIRE -1432 88 -1432 8
WIRE -1304 88 -1304 -8
WIRE -1304 88 -1432 88
WIRE -1256 88 -1256 8
WIRE -1144 88 -1256 88
WIRE -1096 88 -1096 -72
WIRE -1096 88 -1144 88
WIRE -1064 88 -1096 88
WIRE -1432 120 -1432 88
WIRE -1144 120 -1144 88
WIRE -1016 152 -1016 104
WIRE -920 152 -1016 152
WIRE -728 152 -840 152
WIRE -472 152 -648 152
WIRE -376 152 -472 152
WIRE -200 152 -296 152
WIRE 776 184 776 88
WIRE 360 248 360 152
WIRE -472 280 -472 152
WIRE -376 280 -472 280
WIRE -200 280 -200 152
WIRE -200 280 -296 280
WIRE 776 328 776 264
WIRE -472 392 -472 280
WIRE -360 392 -472 392
WIRE -200 392 -200 280
WIRE -200 392 -296 392
WIRE 360 392 360 328
WIRE -360 504 -440 504
WIRE -232 504 -280 504
WIRE -232 520 -232 504
WIRE -472 552 -472 392
WIRE -344 552 -472 552
WIRE -200 552 -200 392
WIRE -200 552 -296 552
WIRE -8 552 -200 552
WIRE 136 552 72 552
WIRE 200 552 136 552
WIRE 360 552 360 472
WIRE 360 552 280 552
WIRE 136 616 136 552
WIRE 536 616 136 616
WIRE 776 616 776 408
WIRE 776 616 616 616
FLAG -424 8 0
FLAG -1432 120 0
FLAG 920 72 0
FLAG 1432 -184 Vo
IOPIN 1392 544 Out
FLAG -232 -520 0
FLAG -456 -536 Vrs_sns
FLAG -232 520 0
FLAG -440 504 Vrc_sns
FLAG -1144 -120 V_omega
IOPIN -1184 608 In
FLAG -1144 120 0
SYMBOL res 1368 -136 R0
SYMATTR InstName R1
SYMATTR Value {RL}
SYMBOL cap 1256 -104 R0
SYMATTR InstName C1
SYMATTR Value {Cf}
SYMBOL bi2 920 -104 R0
SYMATTR InstName B1
SYMATTR Value I={ks} * I(V_Its_sns)
SYMBOL bi2 1032 -136 R0
SYMATTR InstName B2
SYMATTR Value I={kc} * I(V_Itc_sns)
SYMBOL voltage -1432 -88 R0
SYMATTR Value2 AC 1
SYMATTR InstName Vfb
SYMATTR Value ""
SYMBOL e -1016 -152 R0
SYMATTR InstName E1
SYMATTR Value {Irc * Lr}
SYMBOL e -1016 8 R0
SYMATTR InstName E2
SYMATTR Value {Irs * Lr}
SYMBOL ind -936 -168 R270
SYMATTR InstName L1
SYMATTR Value {Lr}
SYMBOL bv 776 -344 R0
SYMATTR InstName B3
SYMATTR Value V=2* {ks} * V(Vo)
SYMBOL ind -936 168 R270
SYMATTR InstName L2
SYMATTR Value {Lr}
SYMBOL bv 776 280 R180
SYMATTR InstName B4
SYMATTR Value V= 2 * {kc} * V(Vo)
SYMBOL res 760 -152 R0
SYMATTR InstName R2
SYMATTR Value {Rs}
SYMBOL res 760 -8 R0
SYMATTR InstName R3
SYMATTR Value {Rc}
SYMBOL voltage 520 -616 R270
SYMATTR InstName V_Its_sns
SYMATTR Value 0
SYMBOL bv -632 -184 R90
SYMATTR InstName B5
SYMATTR Value V={Omega_s * Lr} * I(V_Irc_sns)
SYMBOL bv -744 152 R270
SYMATTR InstName B6
SYMATTR Value V={Omega_s * Lr} * I(V_Irs_sns)
SYMBOL cap -296 -440 R90
SYMATTR InstName C2
SYMATTR Value {Cr}
SYMBOL bi2 -360 -312 R270
SYMATTR InstName B7
SYMATTR Value I={Omega_s * Cr } * V(Vrc_sns)
SYMBOL bi2 -360 -184 R270
SYMATTR InstName B8
SYMATTR Value I={Vrc * Cr } * V(V_omega)
SYMBOL cap -360 408 R270
SYMATTR InstName C3
SYMATTR Value {Cr}
SYMBOL bi2 -296 280 R90
SYMATTR InstName B9
SYMATTR Value I={Omega_s * Cr} * V(Vrs_sns)
SYMBOL bi2 -296 152 R90
SYMATTR InstName B10
SYMATTR Value I={Vrs * Cr } * V(V_omega)
SYMBOL bv 776 -472 R0
SYMATTR InstName B11
SYMATTR Value V={krs} * I(V_Itc_sns)
SYMBOL bv 776 424 R180
SYMATTR InstName B12
SYMATTR Value V= {krc} * I(V_Its_sns)
SYMBOL voltage 520 616 R270
SYMATTR InstName V_Itc_sns
SYMATTR Value 0
SYMBOL e -376 -536 R270
SYMATTR InstName E3
SYMATTR Value 1
SYMBOL e -376 504 R270
SYMATTR InstName E4
SYMATTR Value 1
SYMBOL ind 344 -504 R0
SYMATTR InstName L3
SYMATTR Value {Lm}
SYMBOL bv 360 -248 R180
SYMATTR InstName B13
SYMATTR Value V={Omega_s * Lm} * I(V_Imc_sns)
SYMBOL ind 344 376 R0
SYMATTR InstName L4
SYMATTR Value {Lm}
SYMBOL bv 360 344 R180
SYMATTR InstName B15
SYMATTR Value V={Omega_s * Lm} * I(V_Ims_sns)
SYMBOL voltage 184 552 R270
SYMATTR InstName V_Imc_sns
SYMATTR Value 0
SYMBOL voltage 184 -552 R270
SYMATTR InstName V_Ims_sns
SYMATTR Value 0
SYMBOL voltage -56 -552 R270
SYMATTR InstName V_Irs_sns
SYMATTR Value 0
SYMBOL voltage -24 552 R270
SYMATTR InstName V_Irc_sns
SYMATTR Value 0
SYMBOL bv 360 -72 R180
SYMATTR InstName B14
SYMATTR Value V={Ims * Lm} * V(V_omega)
SYMBOL bv 360 168 R180
SYMATTR InstName B16
SYMATTR Value V={Imc * Lm} * V(V_omega)
SYMBOL e -1256 -88 R0
SYMATTR InstName E5
SYMATTR Value {2*pi*F0}