[begemot]

В принципе-да, работает.
Правда в моём моделировании показывает примерно в 2 раза меньшее значение. Т.е. не очень точно. Но порядок правильный.

Код:

Fourier components of V(in)
DC component:-7.73628e-014

Harmonic    Frequency     Fourier     Normalized     Phase      Normalized
Number       [Hz]       Component     Component    [degree]    Phase [deg]
    1       2.000e+04    1.000e+00    1.000e+00       -0.00°        0.00°
    2       4.000e+04    1.401e-11    1.401e-11       -0.73°       -0.73°
    3       6.000e+04    2.101e-11    2.101e-11       -0.57°       -0.57°
    4       8.000e+04    2.801e-11    2.801e-11       -0.51°       -0.51°
    5       1.000e+05    3.501e-11    3.501e-11       -0.49°       -0.49°
    6       1.200e+05    4.202e-11    4.202e-11       -0.50°       -0.50°
    7       1.400e+05    4.902e-11    4.902e-11       -0.52°       -0.52°
    8       1.600e+05    5.602e-11    5.602e-11       -0.54°       -0.54°
    9       1.800e+05    6.302e-11    6.302e-11       -0.57°       -0.57°
Total Harmonic Distortion: 0.000000%

N-Period=1
Fourier components of V(ref_out)
DC component:2.23017e-012

Harmonic    Frequency     Fourier     Normalized     Phase      Normalized
Number       [Hz]       Component     Component    [degree]    Phase [deg]
    1       2.000e+04    1.000e+00    1.000e+00      179.98°        0.00°
    2       4.000e+04    5.709e-08    5.709e-08       -0.24°     -180.23°
    3       6.000e+04    8.703e-08    8.704e-08       89.80°      -90.18°
    4       8.000e+04    3.021e-11    3.021e-11        0.78°     -179.21°
    5       1.000e+05    5.689e-11    5.689e-11      127.74°      -52.25°
    6       1.200e+05    4.197e-11    4.197e-11      178.98°       -1.00°
    7       1.400e+05    4.903e-11    4.903e-11      179.03°       -0.96°
    8       1.600e+05    5.603e-11    5.603e-11      179.10°       -0.89°
    9       1.800e+05    6.303e-11    6.303e-11      179.12°       -0.86°
Total Harmonic Distortion: 0.000010%

N-Period=1
Fourier components of V(gen)
DC component:-6.39304e-012

Harmonic    Frequency     Fourier     Normalized     Phase      Normalized
Number       [Hz]       Component     Component    [degree]    Phase [deg]
    1       2.000e+04    1.000e+00    1.000e+00       -0.00°        0.00°
    2       4.000e+04    1.000e-06    1.000e-06      -90.00°      -90.00°
    3       6.000e+04    1.000e-06    1.000e-06     -180.00°     -180.00°
    4       8.000e+04    2.801e-11    2.801e-11       -0.51°       -0.51°
    5       1.000e+05    3.501e-11    3.501e-11       -0.49°       -0.49°
    6       1.200e+05    4.202e-11    4.202e-11       -0.50°       -0.50°
    7       1.400e+05    4.902e-11    4.902e-11       -0.52°       -0.52°
    8       1.600e+05    5.602e-11    5.602e-11       -0.54°       -0.54°
    9       1.800e+05    6.302e-11    6.302e-11       -0.57°       -0.57°
Total Harmonic Distortion: 0.000141%

N-Period=1
Fourier components of V(comp_out)
DC component:1.25397e-012

Harmonic    Frequency     Fourier     Normalized     Phase      Normalized
Number       [Hz]       Component     Component    [degree]    Phase [deg]
    1       2.000e+04    9.999e-03    1.000e+00      179.18°        0.00°
    2       4.000e+04    3.078e-08    3.079e-06       18.71°     -160.47°
    3       6.000e+04    4.476e-08    4.476e-06       76.88°     -102.30°
    4       8.000e+04    2.924e-11    2.924e-09        0.17°     -179.01°
    5       1.000e+05    2.241e-11    2.241e-09       90.65°      -88.53°
    6       1.200e+05    4.420e-13    4.420e-11      153.66°      -25.51°
    7       1.400e+05    5.317e-13    5.318e-11      157.55°      -21.63°
    8       1.600e+05    5.910e-13    5.911e-11      161.88°      -17.29°
    9       1.800e+05    6.561e-13    6.562e-11      164.42°      -14.76°
Total Harmonic Distortion: 0.000543%



comp_out_rms: RMS(v(comp_out))=0.00707034 FROM 0 TO 0.002
gen_rms: RMS(v(gen))=0.707105 FROM 0 TO 0.002
mag: gen_rms/comp_out_rms=100.01

   

Там видно что у генератора V(gen) 2 и 3 гармоники на уровне 1е-6 (-120дБ)

У самого операционника V(ref_out)
2гармоника - 5.709e-08
3гармоника - 8.704e-08

На выходе "моста" V(comp_out)
2гармоника - 3.079e-06
3гармоника - 4.476e-06

Подавление в мосте mag: gen_rms/comp_out_rms=100.01

Т.е, отнормировав с учётом подавления в 100 раз
2гармоника - 3.079e-08
3гармоника - 4.476e-08

Файл моделирования для тех кто как и я сомневался в правомочности метода:

THD_Compensate.zip Размер: 2.95 KB  Загрузок: 32

Правда надо бы оценить что влияет на точность и какие ограничения метода.

Но даже если только порядок совпадает-просто здорово.
Интересно, как бы это проверить более прямыми методами, а-ля режектор.
У меня никогда не получалось намерить на 20Кгц больше -140/145дБ

.