* SPICE netlist
R1 Vo 0 {RL}
C1 Vo 0 {Cf}
B1 0 Vo I={ks} * I(V_SNS_IS)
B2 0 Vo I={kc} * I(V_SNS_IC)
Vg Vg 0 AC 1
E1 N005 0 Vg 0 {4/pi}
L1 N005 N006 {L}
B3 N004 N007 V=2* {ks} * V(Vo)
L2 0 N008 {L}
B4 N013 N012 V= 2 * {kc} * V(Vo)
R2 N007 0 {Rs}
R3 0 N012 {Rc}
V_SNS_IS N003 N001 0
C2 N003 N002 {C}
B7 N003 N002 I={Omega_s * C } * V(V_SNS_VC)
C3 N009 N010 {C}
B9 N009 N010 I={Omega_s * C } * V(V_SNS_VS)
B11 N001 N004 V={krs} * I(V_SNS_IC)
B12 N011 N013 V= {krc} * I(V_SNS_IS)
V_SNS_IC N010 N011 0
E3 V_SNS_VS 0 N002 N003 1
E4 V_SNS_VC 0 N009 N010 1
B5 N002 N006 V={Omega_s * L} * I(V_SNS_IC)
B6 N008 N009 V={Omega_s * L} * I(V_SNS_IS)
.param L = 197u
.param C = 51n
.param F0 = 1/(2*pi*sqrt(L*C))
.param Z0 = sqrt(L/C)
.param Cf = 32u
.param Vg = 400
.param RL = 15.5
.param Qs = Z0 / RL
.param fs = 1.2* F0
.param Omega_s =2*pi* fs
.param Req = 8*RL/(pi**2)
.param Xeq = Omega_s * L - 1/Omega_s/C
.param Is = (4*Vg/pi ) * Req/(Req**2 + Xeq**2)
.param Ic = - (4*Vg/pi ) * Xeq/(Req**2 + Xeq**2)
.param Vs = - (4*Vg/pi/Omega_s/C ) * Xeq/(Req**2 + Xeq**2)
.param Vc = - (4*Vg/pi/Omega_s/C ) * Req/(Req**2 + Xeq**2)
.param krs = Xeq*(Req**2)/(Req**2 + Xeq**2)
.param krc = krs
.param ks = 2/pi * Req/sqrt(Req**2 + Xeq**2)
.param kc = -2/pi * Xeq/sqrt(Req**2 + Xeq**2)
.param Rs = Req * (Xeq**2) / (Req**2 + Xeq**2)
.param Rc = Req * (Req**2) / (Req**2 + Xeq**2)
.param Omega_0 = sqrt(1/C)
.param Le = L*(1 + (Omega_0**2)/(Omega_s**2))
.param Rx = Xeq * (Xeq**2)/ (Req**2 + Xeq**2)
.param Ry = Xeq * (Xeq**2 + 2*Req**2) / (Req**2 + Xeq**2)
.param Gs = L*Ic + Vs/Omega_s
.param Gc = L*Is - Vc/Omega_s
.ac dec 20 10 50k

* LTspice schematic (.asc)
Version 4
SHEET 1 2176 976
WIRE 432 -488 272 -488
WIRE 432 -440 432 -488
WIRE -192 -392 -288 -392
WIRE -64 -392 -112 -392
WIRE -64 -376 -64 -392
WIRE -176 -344 -304 -344
WIRE -32 -344 -128 -344
WIRE 432 -312 432 -360
WIRE -304 -280 -304 -344
WIRE -192 -280 -304 -280
WIRE -32 -280 -32 -344
WIRE -32 -280 -128 -280
WIRE -784 -168 -880 -168
WIRE -560 -168 -704 -168
WIRE -304 -168 -304 -280
WIRE -304 -168 -480 -168
WIRE -192 -168 -304 -168
WIRE -32 -168 -32 -280
WIRE -32 -168 -112 -168
WIRE 80 -168 -32 -168
WIRE 272 -168 272 -488
WIRE 272 -168 160 -168
WIRE 688 -168 576 -168
WIRE 928 -168 688 -168
WIRE 1040 -168 928 -168
WIRE 1088 -168 1040 -168
WIRE -880 -120 -880 -168
WIRE 432 -120 432 -232
WIRE 688 -120 688 -168
WIRE -928 -104 -1088 -104
WIRE 1040 -104 1040 -168
WIRE 576 -88 576 -168
WIRE 928 -88 928 -168
WIRE -1088 -56 -1088 -104
WIRE -928 -56 -960 -56
WIRE -880 -8 -880 -40
WIRE -256 -8 -880 -8
WIRE 432 -8 432 -40
WIRE 432 -8 -256 -8
WIRE -256 24 -256 -8
WIRE 432 24 432 -8
WIRE 576 56 576 -8
WIRE 688 56 688 -40
WIRE 688 56 576 56
WIRE 928 56 928 -24
WIRE 928 56 688 56
WIRE 1040 56 1040 -24
WIRE 1040 56 928 56
WIRE 576 88 576 56
WIRE -1088 104 -1088 24
WIRE -960 104 -960 -56
WIRE -960 104 -1088 104
WIRE -1088 136 -1088 104
WIRE -880 168 -880 -8
WIRE -784 168 -880 168
WIRE -576 168 -704 168
WIRE -304 168 -496 168
WIRE -208 168 -304 168
WIRE -32 168 -128 168
WIRE 64 168 -32 168
WIRE 272 168 144 168
WIRE 432 200 432 104
WIRE -304 280 -304 168
WIRE -192 280 -304 280
WIRE -32 280 -32 168
WIRE -32 280 -128 280
WIRE 432 344 432 280
WIRE -192 392 -272 392
WIRE -64 392 -112 392
WIRE -64 408 -64 392
WIRE -304 440 -304 280
WIRE -176 440 -304 440
WIRE -32 440 -32 280
WIRE -32 440 -128 440
WIRE 272 488 272 168
WIRE 432 488 432 424
WIRE 432 488 272 488
FLAG -256 24 0
FLAG -1088 136 0
FLAG 576 88 0
FLAG 1088 -168 Vo
IOPIN 1200 -16 Out
FLAG -64 -376 0
FLAG -288 -392 V_SNS_VS
FLAG -64 408 0
FLAG -272 392 V_SNS_VC
FLAG -1088 -104 Vg
IOPIN -976 48 In
SYMBOL res 1024 -120 R0
SYMATTR InstName R1
SYMATTR Value {RL}
SYMBOL cap 912 -88 R0
SYMATTR InstName C1
SYMATTR Value {Cf}
SYMBOL bi2 576 -88 R0
SYMATTR InstName B1
SYMATTR Value I={ks} * I(V_SNS_IS)
SYMBOL bi2 688 -120 R0
SYMATTR InstName B2
SYMATTR Value I={kc} * I(V_SNS_IC)
SYMBOL voltage -1088 -72 R0
SYMATTR Value2 AC 1
SYMATTR InstName Vg
SYMATTR Value ""
SYMBOL e -880 -136 R0
SYMATTR InstName E1
SYMATTR Value {4/pi}
SYMBOL ind -800 -152 R270
SYMATTR InstName L1
SYMATTR Value {L}
SYMBOL bv 432 -328 R0
SYMATTR InstName B3
SYMATTR Value V=2* {ks} * V(Vo)
SYMBOL ind -800 184 R270
SYMATTR InstName L2
SYMATTR Value {L}
SYMBOL bv 432 296 R180
SYMATTR InstName B4
SYMATTR Value V= 2 * {kc} * V(Vo)
SYMBOL res 416 -136 R0
SYMATTR InstName R2
SYMATTR Value {Rs}
SYMBOL res 416 8 R0
SYMATTR InstName R3
SYMATTR Value {Rc}
SYMBOL voltage 64 -168 R270
SYMATTR InstName V_SNS_IS
SYMATTR Value 0
SYMBOL cap -128 -296 R90
SYMATTR InstName C2
SYMATTR Value {C}
SYMBOL bi2 -192 -168 R270
SYMATTR InstName B7
SYMATTR Value I={Omega_s * C } * V(V_SNS_VC)
SYMBOL cap -192 296 R270
SYMATTR InstName C3
SYMATTR Value {C}
SYMBOL bi2 -128 168 R90
SYMATTR InstName B9
SYMATTR Value I={Omega_s * C } * V(V_SNS_VS)
SYMBOL bv 432 -456 R0
SYMATTR InstName B11
SYMATTR Value V={krs} * I(V_SNS_IC)
SYMBOL bv 432 440 R180
SYMATTR InstName B12
SYMATTR Value V= {krc} * I(V_SNS_IS)
SYMBOL voltage 48 168 R270
SYMATTR InstName V_SNS_IC
SYMATTR Value 0
SYMBOL e -208 -392 R270
SYMATTR InstName E3
SYMATTR Value 1
SYMBOL e -208 392 R270
SYMATTR InstName E4
SYMATTR Value 1
SYMBOL bv -464 -168 R90
SYMATTR InstName B5
SYMATTR Value V={Omega_s * L} * I(V_SNS_IC)
SYMBOL bv -592 168 R270
SYMATTR InstName B6
SYMATTR Value V={Omega_s * L} * I(V_SNS_IS)