Version 4
SHEET 1 1824 1724
WIRE 1152 496 1152 464
WIRE 1152 496 976 496
WIRE 1152 544 1152 496
WIRE -96 608 -96 560
WIRE 32 608 32 560
WIRE 32 608 -96 608
WIRE 208 608 32 608
WIRE 336 608 208 608
WIRE 416 608 336 608
WIRE 496 608 416 608
WIRE 576 608 496 608
WIRE 976 640 976 496
WIRE 208 656 208 608
WIRE 336 656 336 608
WIRE 416 656 416 608
WIRE 720 656 560 656
WIRE 768 656 720 656
WIRE 560 688 560 656
WIRE 720 688 720 656
WIRE 1152 688 1152 624
WIRE 496 704 496 608
WIRE 512 704 496 704
WIRE 336 752 336 720
WIRE 416 752 416 720
WIRE 976 752 976 720
WIRE 336 768 336 752
WIRE 416 768 416 752
WIRE 512 784 512 752
WIRE 560 784 560 768
WIRE 560 784 512 784
WIRE 720 784 720 768
WIRE 720 784 560 784
WIRE 512 800 512 784
WIRE 976 880 976 832
WIRE 1152 880 1152 768
WIRE 1152 880 976 880
WIRE 976 896 976 880
WIRE 1152 960 1152 880
FLAG -96 480 0
FLAG 32 480 0
FLAG 208 720 0
FLAG 336 848 0
FLAG 416 848 0
FLAG 576 608 x
IOPIN 576 608 Out
FLAG 1152 464 Vp
IOPIN 1152 464 BiDir
FLAG 1152 960 Vn
IOPIN 1152 960 BiDir
FLAG 512 800 0
FLAG 768 656 s
IOPIN 768 656 Out
FLAG 336 752 V_x_off
FLAG 416 752 V_x_on
FLAG 0 0 V_x_on
FLAG 976 896 0
SYMBOL bi -96 480 R0
WINDOW 3 -97 153 Left 2
SYMATTR Value I=G_on(I(Vc))*f_on(V(x))
SYMATTR InstName B1
SYMBOL bi 32 480 R0
SYMATTR InstName B2
SYMATTR Value I=G_off(I(Vc))*f_off(V(x))
SYMBOL cap 192 656 R0
SYMATTR InstName C1
SYMATTR Value 1
SYMBOL diode 320 656 R0
SYMATTR InstName D1
SYMBOL diode 432 720 R180
WINDOW 0 24 64 Left 2
WINDOW 3 24 0 Left 2
SYMATTR InstName D2
SYMBOL voltage 336 752 R0
WINDOW 123 0 0 Left 2
WINDOW 39 0 0 Left 2
WINDOW 3 -147 72 Left 2
SYMATTR Value {x_off}
SYMATTR InstName V1
SYMBOL voltage 416 752 R0
WINDOW 3 24 94 Left 2
WINDOW 123 0 0 Left 2
WINDOW 39 0 0 Left 2
SYMATTR Value {x_on}
SYMATTR InstName V2
SYMBOL voltage 1152 672 R0
WINDOW 123 0 0 Left 2
WINDOW 39 0 0 Left 2
SYMATTR InstName Vc
SYMATTR Value 0
SYMBOL bv 1152 528 R0
SYMATTR InstName B4
SYMATTR Value V=cs(I(Vc),V(x))
SYMBOL e 560 672 R0
WINDOW 3 37 48 Left 2
SYMATTR Value {1/x_off}
SYMATTR InstName E1
SYMBOL voltage 976 624 R0
WINDOW 3 24 96 Invisible 2
WINDOW 123 0 0 Left 2
WINDOW 39 0 0 Left 2
SYMATTR Value PULSE(0 1 100n 1f 1f 1u 10u)
SYMATTR InstName V3
SYMBOL res 960 736 R0
SYMATTR InstName R1
SYMATTR Value 1
SYMBOL res 704 672 R0
SYMATTR InstName R2
SYMATTR Value 1
TEXT 24 -32 Left 2 !.param k_off=0.05 i_off=7u alpha_off=3 R_off=100k
TEXT 24 0 Left 2 !.param k_on=-0.05 i_on=-7u alpha_on=3 R_on=1k
TEXT 24 32 Left 2 !.param w_c=100p a_off=2.8n a_on=0.2n
TEXT 24 64 Left 2 !.param x_on=0n x_off=3n
TEXT 16 208 Left 2 !.func f_off(w) {exp(-exp( (w-a_off)/w_c))}
TEXT 16 240 Left 2 !.func f_on(w) {exp(-exp( -(w-a_on)/w_c))}
TEXT 16 280 Left 2 !.func G_off(i) {k_off*(i/i_off-1)**alpha_off*u(i-i_off)}
TEXT 16 312 Left 2 !.func G_on(i) {k_on*(i/i_on-1)**alpha_on*u(i_on-i)}
TEXT 0 -72 Left 2 ;Use fitting parameter of figure 10 and 12
TEXT 24 96 Left 2 !.param lambda=ln(R_off/R_on)
TEXT 16 424 Left 2 !.ic V(x)={x_init}
TEXT 24 128 Left 2 !.param x_init=3n
TEXT 8 952 Left 2 !.model idealD D(Ron=1f Roff=100G Vfwd=1f)
TEXT 416 1016 Left 2 !.tran 0 2u 0 1n
TEXT 16 384 Left 2 !.func cs1(i,x) {i*R_on*exp(lambda*(x_off-x)/(x_off-x_on))}
TEXT 16 352 Left 2 !.func cs(i,w) {i* ( R_on + (R_off-R_on)*(x_off-w)/(x_off-x_on) ) }