Пассивность: протестируйте, визуализируйте и осуществите пассивность рационального подходящего Выхода
В этом примере показано, как протестировать, визуализируйте и осуществите пассивность выхода от rationalfit функция.
Пассивность данных S-параметра
Анализ временного интервала и симуляция зависят критически от способности преобразовать данные S-параметра частотного диапазона в причинные, устойчивые, и пассивные представления временного интервала. Поскольку rationalfit функционируйте гарантирует, что все полюса находятся в левой полуплоскости, rationalfit выход является и устойчивым и причинным конструкцией. Проблемой является пассивность.
Данные S-параметра N-порта представляют зависимую частотой передаточную функцию H (f). Можно создать объект S-параметров в RF Toolbox™ путем чтения файла Touchstone®, такого как passive.s2p, в sparameters функция.
Можно использовать the ispassive функционируйте, чтобы проверять пассивность данных S-параметра и passivity функционируйте, чтобы построить 2-норму N x N матрицы H (f) на каждой частоте данных.
S = sparameters('passive.s2p');
ispassive(S)ans = logical
1
passivity(S)
Тестирование и визуализация rationalfit Выведите пассивность
rationalfit функция преобразует N-порт sparameter данные, S в матрицу NxN rfmodel.rational объекты. Используя ispassive функция на N x N подходящий выход сообщает об этом даже если входные данные S пассивный элемент, выходная подгонка не является пассивным элементом. Другими словами, норма H (f) больше один на некоторой частоте в области значений [0, Inf].
passivity функционируйте берет N x N подгонка, как введено и строит ее пассивность. Это - график верхней границы нормы (H (f)) на [0, Inf], также известный как H-норму-по-бесконечности.
fit = rationalfit(S); ispassive(fit)
ans = logical
0
passivity(fit)
makepassive функционируйте берет в качестве входа N x N массив подходящих объектов и также исходных данных S-параметра, и производит пассивную подгонку при помощи выпуклых методов оптимизации, чтобы оптимально совпадать с данными входа S S-параметра в то время как удовлетворяющие ограничения пассивности. Остатки C и проходная матрица D выхода pfit изменяются, но полюса выхода pfit идентичны полюсам входной подгонки.
pfit = makepassive(fit,S,'Display','on');
ITER H-INFTY NORM FREQUENCY ERRDB CONSTRAINTS 0 1 + 1.791e-02 17.6816 GHz -40.4702 1 1 + 2.878e-04 275.332 MHz -40.9167 5 2 1 + 9.274e-05 365.612 MHz -40.9092 7 3 1 - 3.458e-07 368.111 MHz -40.906 9
ispassive(pfit)
ans = logical
1
passivity(pfit)
all(vertcat(pfit(:).A) == vertcat(fit(:).A))
ans = logical
1
Запустите makepassive с предписанными полюсами и нулем C и D
Продемонстрировать, что только C и D изменяются makepassive, можно обнулить C и D и повторно выполнить makepassive. Выход, pfit все еще имеет те же полюса как входная подгонка. Различия между pfit и pfit2 возникните из-за различных начальных точек выпуклой оптимизации.
Можно использовать эту функцию makepassive функционируйте, чтобы произвести пассивную подгонку из предписанного набора полюсов без любой идеи запустить C и D.
for k = 1:numel(fit) fit(k).C(:) = 0; fit(k).D(:) = 0; end pfit2 = makepassive(fit,S); passivity(pfit2)
all(vertcat(pfit2(:).A) == vertcat(fit(:).A))
ans = logical
1
Сгенерируйте эквивалентную схему SPICE от пассивной подгонки
generateSPICE функционируйте берет пассивную подгонку и генерирует эквивалентную схему как SPICE subckt файл. Входная подгонка является N x N массив rfmodel.rational объекты, как возвращено rationalfit с S-параметры возражают, как введено. Сгенерированный файл является моделью SPICE, созданной только из пассивного элемента R, L, C элементы и управляемые исходные элементы E, F, G, и H.
generateSPICE(pfit2,'mypassive.ckt') type mypassive.ckt
* Equivalent circuit model for mypassive.ckt .SUBCKT mypassive po1 po2 Vsp1 po1 p1 0 Vsr1 p1 pr1 0 Rp1 pr1 0 50 Ru1 u1 0 50 Fr1 u1 0 Vsr1 -1 Fu1 u1 0 Vsp1 -1 Ry1 y1 0 1 Gy1 p1 0 y1 0 -0.02 Vsp2 po2 p2 0 Vsr2 p2 pr2 0 Rp2 pr2 0 50 Ru2 u2 0 50 Fr2 u2 0 Vsr2 -1 Fu2 u2 0 Vsp2 -1 Ry2 y2 0 1 Gy2 p2 0 y2 0 -0.02 Rx1 x1 0 1 Cx1 x1 0 2.73023889928382e-12 Gx1_1 x1 0 u1 0 -2.05821474710108 Rx2 x2 0 1 Cx2 x2 0 7.77758884882576e-12 Gx2_1 x2 0 u1 0 -2.91590408223011 Rx3 x3 0 1 Cx3 x3 0 2.29141629980399e-11 Gx3_1 x3 0 u1 0 -0.544569559438592 Rx4 x4 0 1 Cx4 x4 0 9.31845201627123e-11 Gx4_1 x4 0 u1 0 -0.654438612763948 Rx5 x5 0 1 Cx5 x5 0 4.89917764982128e-10 Gx5_1 x5 0 u1 0 -0.0811178729493162 Rx6 x6 0 1 Fxc6_7 x6 0 Vx7 18.7776484328739 Cx6 x6 xm6 3.95175907226011e-09 Vx6 xm6 0 0 Gx6_1 x6 0 u1 0 -0.0922164867069764 Rx7 x7 0 1 Fxc7_6 x7 0 Vx6 -0.0836352694088244 Cx7 x7 xm7 3.95175907226011e-09 Vx7 xm7 0 0 Gx7_1 x7 0 u1 0 0.00771255070967324 Rx8 x8 0 1 Cx8 x8 0 1.25490425570277e-08 Gx8_1 x8 0 u1 0 -0.947611765889543 Rx9 x9 0 1 Cx9 x9 0 2.73023889928382e-12 Gx9_2 x9 0 u2 0 -2.08405123616107 Rx10 x10 0 1 Cx10 x10 0 7.77758884882576e-12 Gx10_2 x10 0 u2 0 -2.92719241422735 Rx11 x11 0 1 Cx11 x11 0 2.29141629980399e-11 Gx11_2 x11 0 u2 0 -0.607342254694447 Rx12 x12 0 1 Cx12 x12 0 9.31845201627123e-11 Gx12_2 x12 0 u2 0 -0.692645753148452 Rx13 x13 0 1 Cx13 x13 0 4.89917764982128e-10 Gx13_2 x13 0 u2 0 -0.086070572734409 Rx14 x14 0 1 Fxc14_15 x14 0 Vx15 18.4041527208733 Cx14 x14 xm14 3.95175907226011e-09 Vx14 xm14 0 0 Gx14_2 x14 0 u2 0 -0.0931984769853515 Rx15 x15 0 1 Fxc15_14 x15 0 Vx14 -0.0853325719127744 Cx15 x15 xm15 3.95175907226011e-09 Vx15 xm15 0 0 Gx15_2 x15 0 u2 0 0.00795286573951356 Rx16 x16 0 1 Cx16 x16 0 1.25490425570277e-08 Gx16_2 x16 0 u2 0 -0.948034683426268 Gyc1_1 y1 0 x1 0 -0.141613407560287 Gyc1_2 y1 0 x2 0 -0.0219142462337481 Gyc1_3 y1 0 x3 0 -1 Gyc1_4 y1 0 x4 0 -1 Gyc1_5 y1 0 x5 0 1 Gyc1_6 y1 0 x6 0 -1 Gyc1_7 y1 0 x7 0 -1 Gyc1_8 y1 0 x8 0 0.999868130687755 Gyc1_9 y1 0 x9 0 1 Gyc1_10 y1 0 x10 0 -1 Gyc1_11 y1 0 x11 0 0.809992607413604 Gyc1_12 y1 0 x12 0 0.94190139538071 Gyc1_13 y1 0 x13 0 -0.935564936277693 Gyc1_14 y1 0 x14 0 0.988856148988824 Gyc1_15 y1 0 x15 0 0.957335434638246 Gyc1_16 y1 0 x16 0 -0.999962599647839 Gyd1_1 y1 0 u1 0 0.605803623101776 Gyd1_2 y1 0 u2 0 -0.352481487625827 Gyc2_1 y2 0 x1 0 1 Gyc2_2 y2 0 x2 0 -1 Gyc2_3 y2 0 x3 0 0.899439397281703 Gyc2_4 y2 0 x4 0 0.997165267675473 Gyc2_5 y2 0 x5 0 -0.992305230408251 Gyc2_6 y2 0 x6 0 0.997645853573998 Gyc2_7 y2 0 x7 0 0.965506040538839 Gyc2_8 y2 0 x8 0 -1 Gyc2_9 y2 0 x9 0 -0.26374495608798 Gyc2_10 y2 0 x10 0 0.067827457437524 Gyc2_11 y2 0 x11 0 -1 Gyc2_12 y2 0 x12 0 -1 Gyc2_13 y2 0 x13 0 1 Gyc2_14 y2 0 x14 0 -1 Gyc2_15 y2 0 x15 0 -1 Gyc2_16 y2 0 x16 0 1 Gyd2_1 y2 0 u1 0 -0.336090079737355 Gyd2_2 y2 0 u2 0 0.697953397952857 .ENDS