Алгоритм декомпозиции и сшивания частей радиотехнической системы при подстройке одного из стабилизирующих параметров

Алгоритм декомпозиции и сшивания частей радиотехнической системы при подстройке одного из стабилизирующих параметров

А.В. Скляр

Институт радиотехнических систем и управления Южного федерального университета, г. Таганрог

Аннотация: В данной статье проанализирована возможность декомпозиции радиотехнической системы на части при условии подстройки одного из двух стабилизирующих параметров схемы сшивания. Произведен аналитический расчёт (в общем виде), доказывающий возможность подстройки только одного стабилизирующего параметра в схеме сшивания. Приведён численный пример, подтверждающий вышеизложенную идею.

Ключевые слова: полунатурное моделирование, радиотехническая система, декомпозиция, схема сшивания, стабилизирующий параметр, сшивающий четырехполюсник, идеальная линия задержки, дополнение Шура, правило Крамера, Y-параметры.

Полунатурное моделирование предполагает, что любая радиотехническая система, к примеру - произвольно взятая электрическая цепь, со схемой замещения, приведённой на рис. 1, может быть представлена в виде двух частей - экспериментального макета (аппаратной части) и программной модели (численной части). Такой принцип разделения системы на части называется декомпозицией системы [1].

Щ

а

щ

Рис. 1. - Части А и В исходной системы Этот метод позволяет проверить работоспособность создаваемой системы ещё до окончания её проектирования с целью выявления если не всех, то значительной части её недостатков, которые без ущерба для проекта устраняются разработчиками, непосредственно используя ресурсы этого

метода. В научной литературе такой метод моделирования называется Н1Ь-моделированием [2-9].

Произвольная система, взятая в качестве исходной и показанная на рис. 1, состоит из двух частей А и В, где У1 и У2 - представляют собой выходную и входную эквивалентные проводимости соответственно частей А и В системы. Подразумевается, что одна из этих частей является экспериментальным макетом, а другая - программной моделью. Части А и В системы представлены эквивалентными источниками тока 3$) и З(). Следует отметить, что части А и В представлены в виде схем замещения.

На рис. 2 эта же система представлена в декомпозированном виде и в операторной форме. Разделённая на две части система соединена с помощью схемы сшивания В. Схема сшивания состоит из двух частей, каждая из которых содержит идеализированные источники:

- тока 11(р)ехр(-рт) и 12(р)ехр(-рт);

- напряжения и1(р)ехр(-рт) и и2(р)ехр(-рт) и стабилизирующий элемент г1=1/у1 или г2=1/у2.

Таким образом, схема сшивания представляет собой сшивающий четырёхполюсник, необходимый для согласования аппаратной части с численной частью, и позволяющий обеспечивать информационный обмен между частями декомпозированной системы по методу Якоби. Согласно этому методу, решение уравнений электрического равновесия для каждой части системы (как для экспериментального макета, так и для программной модели) выполняются параллельно друг другу. Затем осуществляется обмен данными между частями декомпозированной системы, и опять происходит решение уравнений для каждой части системы и т.д. Временной интервал между п-ым и п+1-ым обменом данными между частями системы называется шагом интегрирования или просто шагом. Один шаг интегрирования обычно

включает в себя несколько итераций, на каждой из которых решаются уравнения электрического равновесия для декомпозированной системы.

Из рис. 2 видно, что схема сшивания, помимо всего прочего, вносит определённую задержку во взаимодействие частей А и В. Для дальнейшего изложения материала статьи с целью ревизии разбитой на части системы необходимо применить методы анализа радиотехнических цепей [10,11].

Рис. 2. - Система, разбитая на две части с обобщённой схемой сшивания в

операторном виде У параметры сшивающего четырёхполюсника В можно записать в следующем виде:

У =

11 Л

— + — ехр

^ ^_

1- ехр (-2 рт) г 1 1 ^

ехР

' 1 1 ^

— + —

z1 z

ехР

2

+

z1 z

1

1- ехр (-2 рт) 11

— + — 2ехр

2 У

1- ехр

На основе У-параметров для сшивающего четырёхполюсника, имеющего, соответственно, нагрузки У2 и У1, определяются:

- эквивалентные входная - У11 и выходная - У22 проводимости;

- коэффициенты передачи по току G12, 021 и

- коэффициенты передачи по напряжению К12, К21 в следующем виде:

^ zl

— + 7

7„ =-

V ^

+

--7

V ^

— ехр

72 + — + z

"2

1 - 7

V ^

— + 7

I Г 1 I

7 = -

22

V ^

+

--7

V Z2

ехр

7 + —+

--7

V Z2

ехр

7

^12 =

' 1 1 I

— + —

V ^ Z2 У

ех

Р

— + 7

V Z2

+

1 I 1

--72 — ехР

V Z1

7

^21 =

' 1 1 I

-+-

V ^ ^ у

ехР

— + 7

V ^

+

--7

V ^

—ехр

у ^

К 12 =

' 1 1 I

-+ —

V Z1 Z 2 у

ех

р

72 + — + Z 2

-1 - 7 ^ 2

V ^ У

ехр

к 21 =

' 1 1 I

— + —

V ^ Z 2 у

ех

р

7 +1 +

-1 - 7,

V Z2 У

ехр

Если в выражения К12(р), С12(р), 711(р) подставить значения 1/г1 = 72, а в выражения К21(р), С21(р), 722(р) значения 1/г2 = 71, то они преобразуются к следующему виду:

к12(р) = = ехр(-рт), а12(р) = ^ = -ехр(-рт), 7„(р) = 72,

V

/1

К21 (р) = = ехр(-рт), а,(р) = ^ = -ехр(-рт), 722(р) = 71.

V

/ 2

1

z

z

2

1

z

z

2

2

z

1

z

z

2

1

z

2

z

Анализ полученных выражений позволяет сделать следующие выводы. При выборе параметра г1 = 1/У2 наблюдается совпадение напряжений и токов на входе и выходе четырёхполюсника с учётом задержки на один шаг, при этом эквивалентная входная проводимость схемы сшивания будет равна проводимости нагрузки У2. Аналогичным образом, напряжения и токи на входе и выходе четырёхполюсника совпадают с учётом задержки на один шаг, а выходная проводимость схемы сшивания будет равна У1 , если выбрать параметр 22 = 1/У1. Таким образом, в режиме г1=1/У2 и г2=1/У1 четырёхполюсник В является идеальной линией задержки.

Основной целью статьи будет доказательство того, что сходимости параметров разделённой на части системы к параметрам исходной системы можно добиться и в случае, когда правильно настроен только один из двух стабилизирующих параметров

Если определены матрицы параметров, описывающих части А и В системы (или радиотехнической цепи) и составленные по методу узловых напряжений, то значения проводимостей системы У1 и У2 можно вычислить аналитически. Пусть будут даны матрицы второго порядка, описывающие поведение частей А и В системы в общем виде, как приведено ниже:

У У

1Л\\ Л12

УУ

Л 21 Л 22

и2

4

32

и

УУ

1Б\\ Б12

УУ

Б 21 Б 22

из

и4

Л Л

В данном случае порядок матрицы не имеет значения, главное - это то, что они «пересекаются» одним коэффициентом [10,11,12]. Кроме того, предполагается, что поведение исходной системы или радиотехнической цепи (не разбитой на части) описывается следующей системой уравнений, составленной по методу узловых напряжений:

У 1 Л11 У Л12 0 " и" ' 3 '

У Л21 У + У Л22 1БП У Б12 и2 = 32 + 33 ,

0 У 1 Б 21 У 1Б 22 _ _и4 _ _ 34 _

где и2=и3.

В таком случае система, уже разбитая на части А и В, будет описываться следующей системой уравнений:

7-

у у + Л ^

у

1 А12

— + — ехр

0 0

1 - ехр

0

С 1 1 Л

— + —

V ^ Z2 У

0

ех

Р

С 1 1 >

— + —

V Z1 Z2 у

1 - ехр

1 - ехр(- 2 рт) 0

7В11 +-

— + — ехр(- 2 рт) "2 ^_

1 - ехр(- 2 рт) у

у

£12

у

и, и2

и3

и

•1 • 2

•3

Если коэффициенты матриц УА и 7В неизвестны, то необходимо добавить в систему ещё два источника напряжения Е=0 [12], как показано на рис. 3. В этом случае получится выразить токи 1а и 1Ь, которые втекают в В со стороны частей А и В. В этом случае источники напряжения имеют нулевые значения, поэтому они не поменяют общего решения системы.

Е=0

Е=0

-е-

Рис. 3. - Система с добавленными нулевыми источниками напряжения Тогда система примет вид:

у 1 АН у А12 0 0 0 0 0 0

у 1 А21 у А22 0 0 0 0 0 0

0 1 0 -1 0 0 0 0

1 + ± z -

(

0 -1

z

z

2

1 1

-+-

z

Л

1 - г-

С 1 1 Л

— + —

V ^

2 У

-2

0 0 0 0

0 0 0 0

0 0 0 0

V ^

1- г

1+1 -2 -+-Z

00

^ у -1 0 0

-2 - г 1 -2 1- г

0 -1 0 1 0

0 0 1 у у £12

0 0 0 у £21 у £22

и"  • 1

и 2 • 2

1а 0

и 4 0

и 5 0

Ь 0

и 7 • 3

оо _ • 4 _

0

-1

г

0

0

z

После математических преобразований будет получена система, которую необходимо решить:

У 1 А11 У 1 А12 0 0 0 0 0 0 1

У 1 А21 У 1 А22 1 0 0 0 0 0

0 1 0 —1 0 0 0 0

0 0 — 1 1 0 0 0 0

^ 1

0 0 0 0 -1 0 0

0 0 0 0 —1 0 1 0

0 0 0 0 0 1 У Б11 У Б12

0 0 0 0 0 0 У Б21 У 1 Б 22 _

и

и4

и

и7 и

Л

1 и4—2 + Г + ]

22 V 2 2 2 У

из2 + +

V 22 2 У

и-1 -1—

и -1—

4

0 Л 3

Л

0

а

Ъ

Определить значения 21 и системы можно используя следующие выражения:

_ 1/У2 _ УБ 22 22 — УБиУБ 21 У 22 / ^"^С^Б )

¿2 = 1/У = УАП/

В математике такие преобразования называются дополнением Шура. В случае, когда в знаменателе условия получается ноль, то решение резко упрощается. При этом значение стабилизирующих параметров z1 следует выбирать равным 1/Увп, а значение Z2 - равным 1/ УА33 .

При правильном выборе значений стабилизирующих параметров z1 и Z2 система, разбитая на части, решается не итерационным, а прямым методом. Для доказательства этого следует воспользоваться системами (2, 3), выражением (5) и правилом Крамера.

При помощи правила Крамера решение исходной системы (2) можно записать в следующем виде:

(( У — У У + У У )л — У У Т + У У т

и _У В11 Б 22 1Б\21 Б 21 Б 221 А22 Г 1 Б 221 А12^2^ 1 А121 Б12^3

У У У —У У У а-У У У —У У У

Б 22 А11 А 22 А11 Б12 Б 21 А11 Б11 Б 22 Б22 А 21 А

А12

— У У Т + У У Т — У У т

и __ Б 221 А21 1 1Б 221 АП^2 Б12и 3__

2 У У У — У У У + У У У — У У У '

1Б22 А11 А22 1 АН1 Б2Г Б12 """ 1 А1Г БП1 Б22 1 Б221 А21 А12

и — ^2l7A2lJl Y£2l7AllJ2 + (7All7A22 + YAll7£ll 7A2l7Al2 •3

УУУ —V У У -I-V У У —V У У

£ 22 А11 А22 А11 £21 £12 А11 £11 £ 22 £ 22 А21 А12

Теперь, используя правило Крамера, можно найти решение разбитой на части системы (3) с условием, что настраивается только один стабилизирующий параметр. Пусть это будет z1, тогда z1 будет вычисляться в соответствии с выражением (5), а Z2 будет иметь произвольное значение. В этом случае необходимо сравнить результаты решений систем (2) и (3). Решение системы (3) будет в следующем виде:

(у у -у у + у у )• -у у т + у у т р-р

V £11 £22 £12 £21 ^ 1 £221 А22 Г 1 1 £221 А12и 2 ^ 1 А121 £12^3^

и —-

у у у - у у у + у у у - у у у

£ 22 А11 А22 А11 £12 £ 21 А11 £11 £22 £22 А21

А12

-

_ у у т + у у т - у у т р-

и —_ £ 22 А21°1^ £ 22 АП^2 1 АН1 £12^3^__

2 у у у - у у у + у у у - у у у '

£ 22 А11 А 22 А11 £12 £ 21 А11 £11 £ 22 £ 22 А 21 А12

С / \ С у р~2ртЛЛ

(у у у у + у у у у )(1 - Р~2рт)-у у " " " "

V А11 £12 А22 £22 1А2Г А12^ £12^ £ 22 Д1 с / £

'£12

у у - у у + £ 22

1£ 22 £11 £12 £ 21 "

• 3

и —^_V_^_У^ +

^3 у ^

22^ £11 1£121 £ 21 А1 £221 А1Г А22 £2Г £12 """ ^ А1Г £1Г £ 22 ^£22^ А2Г А12/

г2

+

- у у ]р - рт + у у • р ~

£22 А21 1 £22 А11и 2е"

у у у - у у у + у у у - у у у '

1 £ 22 АП^ А22 -'АП^ £12^ £21 ^ 1 АП^ £1И £ 22 ^£22^ А21^ А12

СС у л

у у | £22 £11 £22 """

УаП^ £11 АП^ А22 А21 А12/

• 3

и —_^_^_У_+

4 у ^

22 £11 1 £121 £ 21 ^ А1 £ 221 АП^ А22 £12^ £ 21 ^ ^ АП^ £11^ £ 22 ^£22^ А21Г А12/

г2

(-у у (у у + у у р~2рт - у у р~2рт) + 7All7£l27£21(р р 1))•

V £12 £ 21Ч7 АП7 £11 ^ 7 АП7 А22е 1 А211 А12^ / """

+-—-+

у

У1 £ 22 £11 £12 £ 21 ^ /V"1 £ 221 АП^ А22 £12^ £ 21 ^ ^ АП^ £1Г £ 22 1 £ 221 А2Г А12/

г2

у у •р-рт- у у • р~рт

+ _ А21 £211 А11 £2^ 2

ууу - у у у + у у у - у у у

£22 А11 А22 А11 £12 £21 А11 £11 £ 22 £ 22 А21 А12

Полученные выражения (7) показывают, что если выполнить три итерации без изменения значения вектора свободных членов, т.е.: -т) = -2т); -т) • -2т);

-т) =3з(г -2т) или ерт=1, е-2рт=1, то выражения (7) совпадут с выражениями (6). А это означает, что если правильно подстроен хотя бы один из стабилизирующих параметров и выполнено три итерации, то решение системы (3) всегда совпадает с решением системы (2).

Подводя итог, можно сделать вывод, что при правильном выборе значений двух стабилизирующих параметров схема сшивания обеспечивает устойчивость и сходимость итерационного процесса за две итерации, а если правильно выбрано значение только одного параметра, то за три итерации.

В качестве простого примера, поясняющего вышеизложенные расчёты, будет рассмотрена простейшая радиотехническая цепь, состоящая из источников тока и проводимостей, как показано на рисунке ниже:

\А (1) ¿д (2) }

Мр) Ул

Ф Р

[>

Рис. 4. Исходная радиотехническая цепь Данная цепь будет подвержена декомпозиции и разделена на 2 части причем эти части будут соединены схемой сшивания (см. рисунок ниже):

(4)___________(52

Рис. 5. Радиотехническая цепь после её декомпозиции Номиналы проводимостей будут иметь следующие значения: УА1=5 См, уд2=8 См, УВ1=4 См, Ув2=6 См. Ток источника 11 будет постоянным и равным 0.1 А, а ток источника 12 будет изменяться по закону вт(4*л*п*^) А, где п - номер текущего шага интегрирования, а N - общее количество шагов

интегрирования. В данном эксперименте значение N будет выбрано равным 50.

Матрицы проводимостей уА и ув в таком случае будут иметь вид:

Га =

5 - 5 - 5 13

и Г£ —

10 - 6 - 6 6

Тогда система, состоящая из этих двух частей, будет иметь следующий вид:

5 -5 0 " и"  Л"

-5 13 +10 - 6 и2 — 0

0 -6 6 и _ _ Л _

Также пусть имеется система, составленная на основе вышеприведённой системы, согласно выражению (4):

5 - 5 - 5 13 0 1

0 1

0 0

0

0

0 0 0

0 -1 0 -1 1

0 0 0 0 0 0 0 0 0 0 0 0

0 0 0 0

0

0 0 0

1

00

0 0 — -1

0 0 -1 0 1 0

0 0 0 1 10 - 6

0 0 0 0 - 6 6

"и/ и2

и4

и5

и7

и

- - и;2 + г2

-—и;2 + 2

•1 0

0

11

-+ —

V г2 г1 у

и- -1- 2

г

11

— + —

V 2 2 г1 у 0

0

и -1-

4 2Ь

Целью эксперимента является расчёт и сравнение решений этих двух систем. В исходной и декомпозированной системах будут рассчитываться значения напряжений между нулевым и первым узлами. Пусть такое напряжение для исходной системы будет обозначено как Ш, а для разбитой на части системы как ЦШ. Расчёт будет производиться в программе на МАТЬАВ.

Пусть значение стабилизирующего параметра г1 равно 100 Ом, а стабилизирующий параметр г2 настроен согласно правилу (5) и равен 0.125 Ом. В таком случае будут получены следующие результаты:

а

Ь

Сопоставление результатов моделирования

Рис. 6. Напряжение между узлами 1 и 0 в исходной (Ц1) и декомпозированной (ЦЦ1) системах при z1=100 Ом На рис. 6-10 сплошной линией обозначена зависимость напряжения Ц1 между узлами 1 и 0 исходной системы от номера шага интегрирования п. «Звёздочками» обозначена зависимость напряжения ЦЦ1 между узлами 1 и 0 декомпозированной системы от номера шага интегрирования п.

Теперь следует изменить значение стабилизирующего параметра z1. Пусть значение параметра z1 равно 1 МОм, а значение параметра z2 также остаётся равным 0.125 Ом, как и в предыдущем случае. Тогда будут получены следующие результаты:

Сопоставление результатов моделирования

Рис. 7. - Напряжение между узлами 1 и 0 в исходной (Ц1) и декомпозированной (ЦЦ1) системах при z1= 1 МОм

Пусть значение стабилизирующего параметра 22 равно 100 Ом, а стабилизирующий параметр 21 настроен согласно правилу (5) и равен 0.25 Ом. В таком случае будут получены следующие результаты:

Сопоставление результатов моделирования

Рис. 8. - Напряжение между узлами 1 и 0 в исходной (Ш) и декомпозированной (ЦШ) системах при 22=100 Ом Если значение параметра 22 равно 1 МОм, а значение параметра 21 так же остаётся равным 0.25 Ом, то будут получены следующие результаты:

Сопоставление результатов моделирования

Рис. 9. - Напряжение между узлами 1 и 0 в исходной (Ц1) и декомпозированной (ЦЦ1) системах при 22= 1 МОм Если же настраивать стабилизирующие параметры 21 и 22 без использования дополнения Шура, а задаться произвольными значениями

(например, z1 = ЮОкОм и z2=10 Ом), то будут получены результаты, приведенные на рис. 10:

Сопоставление результатов моделирования

О 5 10 15 20 25 30 35 40 45 50

Рис. 10. - Напряжение между узлами 1 и 0 в исходной (Ш) и декомпозированной (ЦШ) системах при значениях стабилизирующих параметров z1 =100 кОм и z2 =10 Ом Как видно из рис. 6 - 9, напряжение между узлами 0 и 1 декомпозированной системы (ЦЦ1) совпадает с напряжением между узлами 0 и 1 исходной системы (Ц1), если корректно настроен хотя бы один из стабилизирующих параметров ^ или z2). Если же оба стабилизирующих параметра настроены неверно (рис. 10), то значения напряжений между узлами 0 и 1 декомпозированной (ЦЦ1) и исходной (Ц1) систем не совпадают. На основе полученных данных можно сделать следующий вывод: если правильно настроен только один из стабилизирующих параметров, то параметры разделённой на части системы все равно сойдутся к параметрам исходной системы, как и было заявлено ранее. В таком случае, нет необходимости подстраивать оба стабилизирующих параметра. Это существенно упростит процесс РН1Ь-моделирования, т.к. подстройка стабилизирующего параметра в аппаратной части является далеко не простой задачей. Можно задаться каким-то определённым значением стабилизирующего параметра в аппаратной части, а подстроить только

стабилизирующий параметр в программной модели (численной части) и в результате получить сходимость параметров разбитой на части системы к параметрам исходной системы.

Литература

1. Копысов С.П. Методы декомпозиции и параллельные распределённые технологии для адаптивных версий метода конечных элементов: дис. ... док. физ.-мат. наук: 05.13.18. Ижевск, 2006. 404 с.

2. Ren, W., 2007. Accuracy Evaluation of Power Hardware-in-the-Loop (PHIL) Simulation, PhD thesis, Florida State University, Tallahassee, pp: 1-68.

3. Avras, A., A.J. Roscoe and G.M. Burt, 2014. Scalable Real-Time Controller Hardware-In-the-Loop Testing for Multiple Interconnected Converters. The UPEC 2014 conference, IEEE Date Views 02.06.2017 URL: dx.doi.org/10.1109/UPEC.2014.6934620.

4. Скляр А.В. Полунатурное моделирование многофункциональных модульных систем. Условия выбора модели. Инженерный вестник Дона, 2017, №2 URL: ivdon.ru/ru/magazine/archive/N2y2017/4236

5. Шайкин А.С., Шайкина Е.В. Применение комплекса полунатурного моделирования в процессе проектирования информационно-измерительных и управляющих систем. Инженерный вестник Дона, 2014, №1. URL: ivdon.ru/ru/magazine/archive/n1y2014/2248

6. Скляр А.В., Мережин Н.И. Полунатурное моделирование комплексных систем // Компьютерные и информационные технологии в науке, инженерии и управлении "КомТех-2017": материалы Всероссийской научно-технической конференции с международным участием. Таганрог: ЮФУ, 2017.С. 21-24.

7. Попов В.П., Максимов М.Н., Мережин Н.И. Об устойчивости и сходимости моделирования по частям // Вестник Южного научного центра РАН. Т.1. № 3. 2005. С. 11-21.

8. Мережин Н.И. Стенд для аналого-цифрового моделирования с использованием адаптивных схем сшивания // Материалы международной научной конференции «Методы и алгоритмы принятия эффективных решений». Ч.2. - Таганрог: ТТИ ЮФУ, 2009. С. 47-51.

9. Максимов М.Н., Мережин Н.И., Скляр А.В., Мережин Д.Н. Использование оператора Пуанкаре-Стеклова для обеспечения устойчивости PHIL моделирования // Сотрудничество стран БРИКС для устойчивого развития: материалы Международной научно-практической конференции молодых ученых стран БРИКС. Ростов-на-Дону: ЮФУ, 2015. Т.2. С. 81-82.

10. Максимов М.Н., Мережин Н.И., Федосов В.П., Лабынцев А.В., Максимов А.А. Эквивалентная схема сшивающего четырехполюсника // Радиотехника и электроника. 2016. Т.61. №2. С. 162-169.

11. Maksimov M., Merezhin N., Lyashev V., Sinyutin S. Poincare-Steklov filter in hardware-in-the-loop modeling. 2017 International Siberian Conference on Control and Communications (SIBCON 2017). Date Views 15.12.2018 URL: doi.org/10.1109/SIBC0N.2017.7998531.

12. Скляр А.В., Максимов М.Н. Расщепление системной матрицы с помощью обобщённой схемы сшивания при полунатурном моделировании радиотехнических цепей // Материалы Всероссийской научно-технической конференции «Современные проблемы радиоэлектроники». Красноярск: СФУ, 2018. С. 41-45.

References

1. Kopysov S.P. Metody dekompozitsii i parallel'nye raspredelennye tekhnologii dlya adaptivnykh versiy metoda konechnykh elementov [Decomposition methods and parallel-distributed technologies for adaptive versions of the finite element method]: dis. Ph.D.: 05.13.18. Izhevsk, 2006. 404 p.

2. Ren, W., 2007. Accuracy Evaluation of Power Hardware-in-the-Loop (PHIL) Simulation, PhD thesis, Florida State University, Tallahassee. pp: 1-68.

3. Avras, A., A.J. Roscoe and G.M. Burt, 2014. Scalable Real-Time Controller Hardware-In-the-Loop Testing for Multiple Interconnected Converters. The UPEC 2014 conference, IEEE Date Views 02.06.2017. URL: dx.doi.org/10.1109/UPEC.2014.6934620.

4. Sklyar A.V. Inzenernyj vestnik Dona (Rus), 2017, №2. URL: ivdon.ru/ru/magazine/archive/ N2y2017/4236

5. Shaykin A.S., Shaykina E.V. Inzenernyj vestnik Dona (Rus), 2014, №1. URL: ivdon.ru/ru/magazine/archive/n1y2014/2248

6. Sklyar A.V., Merezhin N.I. Polunaturnoye modelirovaniye kompleksnykh sistem [PHIL simulation of complex systems]. Komp'yuternyye i informatsionnyye tekhnologii v nauke, inzhenerii i upravlenii "KomTekh-2017": materialy Vserossiyskoy nauchno-tekhnicheskoy konferentsii s mezhdunarodnym uchastiyem. Taganrog: YUFU, 2017. pp. 21-24.

7. Popov V.P., Maksimov M.N., Merezhin N.I. Vestnik Yuzhnogo nauchnogo tsentra RAN. 2005. Vol.1. № 3. pp. 11-21.

8. Merezhin N.I. Stend dlya analogo-tsifrovogo modelirovaniya s ispol'zovaniem adaptivnykh skhem sshivaniya [Stand for PHIL simulation using adaptive coupling schemes]. Materialy mezhdunarodnoy nauchnoy konferentsii «Metody i algoritmy prinyatiya effektivnykh resheniy». Taganrog:TTI YuFU, 2009. Vol.2. pp. 47-51.

9. Maksimov M.N., Merezhin N.I., Sklyar A.V., Merezhin D.N. Ispol'zovanie operatora Puankare-Steklova dlya obespecheniya ustoychivosti PHIL modelirovaniya [The using of the Poincare-Steklov operator to ensure the stability of the PHIL simulation]. Sotrudnichestvo stran BRIKS dlya ustoychivogo razvitiya: materialy Mezhdunarodnoy nauchno-prakticheskoy konferentsii

molodykh uchenykh stran BRIKS. Rostov-na-Donu: YuFU, 2015. Vol.2. pp. 8182.

10. Maksimov M.N., Merezhin N.I., Fedosov V.P., Labyntsev A.V., Maksimov A.A. Radiotekhnika i elektronika. 2016. Vol.61. №2. pp. 162-169.

11. Maksimov M., Merezhin N., Lyashev V., Sinyutin S. Poincare-Steklov filter in hardware-in-the-loop modeling. 2017 International Siberian Conference on Control and Communications (SIBCON 2017). Date Views 15.12.2018 URL: doi.org/10.1109/SIBC0N.2017.7998531.

12. Sklyar A.V., Maksimov M.N. Rasshcheplenie sistemnoy matritsy s pomoshch'yu obobshchennoy skhemy sshivaniya pri polunaturnom modelirovanii radiotekhnicheskikh tsepey [Splitting of the system matrix using a generalized coupling scheme for PHIL simulation of radio engineering circuits]. Materialy Vserossiyskoy nauchno-tekhnicheskoy konferentsii «Sovremennye problemy radioelektroniki». Krasnoyarsk: SFU, 2018. pp. 41-45.