Построение интегральной модели на примере динамики ветроэнергетической установки

УДК 517.968.2

DOI: 10.14529/ mm p150404

ПОСТРОЕНИЕ ИНТЕГРАЛЬНОЙ МОДЕЛИ НА ПРИМЕРЕ ДИНАМИКИ ВЕТРОЭНЕРГЕТИЧЕСКОЙ УСТАНОВКИ

C.B. Солодуша, Д.О. Герасимов, К.В. Суслов

Данная работа посвящена применению интегро-степенных рядов Вольтерра к описанию нелинейных динамических систем типа «вход - выход». Универсальность данного математического аппарата позволяет создавать программное обеспечение для проведения экспериментов на компьютере. Работа продолжает исследования в области идентификации ядер Вольтерра, начатые в Институте систем энергетики им. Л. А. Ме-лентьева СО РАН. В первой части статьи рассматривается новый алгоритм идентификации полиномов Вольтерра второй степени для систем, допускающих активный эксперимент с помощью тестовых наборов возмущающих воздействий. Во второй части работы приведены результаты численных расчетов для «этадонной» динамической системы. Эталоном послужила модель ветроэнергетической установки с горизонтальной осью вращения. Построены квадратичные полиномы Вольтерра, описывающие нелинейную динамику угловой скорости вращения элементов ветроустановки от угла наклона лопастей и скорости ветра. Практическая идентификация ядер Вольтерра производилась относительно некоторого выбранного стационарного состояния моделируемой системы.

Ключевые слова: нелинейная динамическая система; квадратичный полином Вольтерра; ветроэнергетическая установка с горизонтальной осью вращения.

Введение

В приложениях нередко нелинейные динамические системы (объекты) типа «вход - выход» описываются с помощью аппарата интегростепенных рядов Вольтерра. Назовем лишь некоторые направления исследований, где используется данный математический аппарат: моделирование технических систем [1,2] и электронных устройств [3], нелинейная идентификация каналов связи [4,5] и систем визуализации [6], анализ нестационарных временных рядов |7j.

Применение конечного отрезка ряда Вольтерра [8], дающего представление отклика y

N 1 1 m

y

m=l 0 0 i=l

для моделирования нелинейных динамических объектов самой разной природы, базируется на теореме Фреше [9] и различных ее обобщениях — континуальных аналогах классической теоремы Вейерштрасса об аппроксимации непрерывной функции многочленом.

Построить интегральную модель в виде (1) — значит восстановить многомерные переходные характеристики нелинейной динамической системы Km. Сложность практического применения (1) состоит в том, что для однозначного восстановления (i + 1)-мерного ядра Ki(t, s1}..., si), i = 1, m, необходимо задание откликов yi(t,u1, на

¿-параметрическое семейство тестовых возмущений хШ1 (Ь). Это объясняет пристальное внимание многих авторов к проблеме идентификации (см., например, критический обзор в [10]).

Наиболее распространен подход [11], где предложено для восстановления ядер Вольтерра задавать многопараметрическое семейство тестовых сигналов, состоящих из комбинации ^-функций Дирака. Однако такой подход имеет ограниченную область применения [12].

Работы Института систем энергетики им. Л.А. Мелентьева (ИСЭМ) СО РАН в области идентификации ядер Вольтерра были Н&ЧсХТЫ в 90 -X ГОДсХХ (см. обзорную статью [13]). Данная методика [10,14] основана на задании семейства тестовых сигналов в виде специальных ЛИН6ИНЫХ комбинаций функций Хевисайда с отклоняющимся аргументом с амплитудами а^ г = 1,ш. При этом задача идентификации сводится к решению линейных интегральных уравнений Вольтерра I рода, допускающих явные формулы обращения. В [15,16] в терминах амплитуд тестовых сигналов установлены необходимые условия разрешимости соответствующих многомерных интегральных уравнений в определенном классе функций. В то же время, за счет использования полученных условий возможно уменьшение числа тестовых сигналов для разделения отклика системы на составляющие.

Цель данной статьи — на примере «эт^онной> динамической системы рассмотреть новый способ построения интегральных моделей в виде полиномов Вольтерра (1) в наиболее важном для приложений случае N = 2.

1. Новый алгоритм идентификации квадратичного полинома Вольтерра

Опишем основную идею получения линейных многомерных интегральных уравнений Вольтерра I рода, к решению которых сводится задача идентификации ядер Вольтерра. Полагая в (1) N = 2, запишем квадратичную модель

г г г

Уquad(t) = J К1(г,81)х(81№81 + J У К2(г,81,82)х(81)х(82^81¿82, Ь € [0,Т]. (2) 0 0 0

Заметим, что из скалярности х(Ь) вытекает симметричность К2(Ь, 81, 82) в (2) по переменным 81, 82. Если отклик динамической системы уа(Ь,ш1,ш2) (для определенности будем считать, что а > 0) на тестовые сигналы

хаШъШ2

таков, что при 0 < ш1 + с2 < Ь < Т выполняется

а

Ш1 Ш 1 +Ш2

уа(Ь,с1,с2) - а I К1(1,81)й81 + а I К1(1,81)й81'

а2

3уа(Ь,С1, 0)-

0 Ш1

Ш1 Ш1 +Ш2

-3уа

0 Ш1

f2(t,Ui ,U2) = -1

а2

то в обозначениях

Ш1 Ш1 +Ш2

уа(Ь,шх,ш2) — а ! Кх(Ь,вх)ввх + а J Кх(Ь,вх)ввх

0 ш1

справедливо ключевое равенство

¡2^,их,Ш2) = 0) — $¡2^, 0, Шх + Ш2) + ¡2^,Шх + Ш2, —Ш2).

Условие (6) является необходимым для существования решения уравнения

Ш1 Ш1 Ш1 Ш1+Ш2

У У К2(1, вх, 82)ввхвв2 — J К2(Ь, вх, 82)ввхвв2 +

0 0 0 Ш1

Ш1 +Ш2 Ш1+Ш2

+ ! J К2

Ш1 Ш1

в классе непрерывных на = {Ь,вх,в2/0 < вх,в2 < Ь < Т} функций, симметричных по второму и третьему аргументу [17].

Подставляя (3) в (2), приходим к интегральному уравнению

Ш1 Ш1+Ш2 Ш1 Ш1

аJ Кх(Ь, вх)ввх — а ! Кх(Ь, вх)ввх + а2 ^ J К2(Ь, вх, в2)ввхвв2 —

0 Ш1 0 0

Ш1 Ш 1+Ш2 Ш 1+Ш2 Ш 1+Ш2

— J К2

0 Ш1 Ш1 Ш1

Замечание 1. В случае, когда динамическая система стационарна, т.е. ядра Воль-терра Кт в (1) зависят лишь от разности Ь — в^ I = 1,ш, для однозначного восстановления Кх, К2 в модели

г г г

Удиа4(Ь) = ! Кх(вх)х(Ь — вх)ввх + ^ J К<2(вх,в2)х(Ь — вх)х(Ь — в2)ввхвв2 (9) 0 0 0

достаточно уравнений

а

ya(t,ux) - а Kx(sx)dsx

1

~j2

ya(t — ux, —ux) + а Kx(sx)dsx

а2

t — Wl t — Wl

t t t а j Kx(sx)dsx + a2 j j Kj^sx, S2)dsxds2 = ya(t,ux).

t— wi t— wi t— wi

t

t

При этом в очевидных обозначениях справедливо условие [14]: /2(г,ш1) = ¡2(г — и1, необходимое для разрешимости

г г

J У ^

г-ш1 г-ш1

в классе непрерывных на = {81, 82/0 < 81, 82 < Т} симметричных функций.

Замечание 2. В работе [18] описан алгоритм восстановления ядер Кь К2, основанный на сведении исходной задачи к решению системы, состоящей из двух уравнений типа (8) с откликами уа(г,Ш1,Ш2) и у-а(г,Ш1,Ш2) соответственно. Подчеркнем чимость условия (6) при построении интегральной модели в виде (2). Тем самым, интегральных уравнений (4), (8) достаточно для однозначной идентификации ядер

К1 К2

Продемонстрируем эту технику на тестовом примере.

Пример 1. Пусть отклик некоторой математической системы на сигналах (3) имеет следующий вид:

Уег

+ -г4 (484 — 58382 + — 58183 + 82) + у15 ^ — 82) , где 81 = и1, 82 = и1 + и2. Тогда

Уег

Уег

Полагая в (4) и2 = 0 и учитывая (11), (12), приходим к интегральному уравнению Вольтерра I рода

2 4

У К1(г, 8{)с181 = Ьи1 + у г3 и3 + у г5ш5, 0

дифференцируя которое по и1, получаем

2 4

К1

Учитывая (13), из (5) имеем

2 4

¡2(г,81,82 — 81) = у г2( 281 — 82У + у г4 (484 — 583 82 + б^ — 58^ + 84) и, согласуясь с явной формулой обращения (7) [17]:

К (г ) = 1 д212(г,81,82) К2(г,81,82) = — 4 д81д82

находим

1 а2

K2(t, 81,82) = 212 + Y*4(5S1 - 8SlS2 + bs2) , 0 < Si < S2 < t < T.

Итак, квадратичная модель полностью идентифицирована. Легко проверить, что для входов х<а1,ш2 (t) отклики (2) и (10), VqUad(t,M1, ш2) и vtt(t,^1,u2) соответственно, тождественно совпадают.

2. Эталонная динамическая система

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

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

В качестве эталонной динамической системы рассмотрим математическую модель ветроэнергетической установки с горизонтальной осью вращения, представленную с помощью методик [19-21] в следующем виде:

Z

Z(t) = ' z(t) VZ(t) + 0.08b(t) - ¥(t)Ti) ' (14}

Cp

MT t) = Р^Щ^, Ш = MT

2шТ (t) dt J

где шт

— плотность воздуха, S (м2) — ометаемая площадь, R (м) — радиус ветроколеса, b (град) — угол наклона лопастей по нормали от направления ветра, V (м/с) — скорость ветра; безразмерные величины

: Cp — коэффициент использования энергии ветра, Z быстроходность, z — текущее значение быстроходности. Численное решение системы (14) - (16) выполнено с помощью метода Рунге -Кутта четвертого порядка. Отклики Ашт(t) = шт (t) — шТо на тестовые воздействия Ab(t) = b(t)-b0 вида, (3) и V(t) = V0 (AV(t) = V(t)-V0n b(t) = b0) для t e [0, 20] (с) на равномерной сетке с шагом h = 1 (с) использовались для построения квадратичных полиномов Вольтерра

t t t АшТ(t) = j K1 (t,s1)Ab(s1)ds1 + У J Kl(t,s1,s2)Ab(s1)Ab(s2)ds1ds2, (17) 0 00

г г

Аш2т (г)= К2(г,в1)АУ (в^вг + К1(г,в1,в2)АУ (в')АУ (в2 ¿82

о о

соответственно. Такой выбор шага связан с экспериментальными данными, полученными при решении (14) - (16).

Для обеспечения лучшей точности амплитуда пробных сигналов, используемых при нахождении ядер Вольтерра в (17), (18), была согласована с величиной действующих возмущений. Заметим, что модель, построенную только по одной группе пробных сигналов, нельзя рассматривать как «всережимную», т.е. одинаково пригодную для расчетов во всем диапазоне допустимых изменений входных воздействий. Практическая идентификация ядер проводилась относительно некоторого выбранного стационарного состояния моделируемой системы. С целью повышения точности моделирования были введены опорные начальные режимы, для которых выполнено построение интегральных моделей. При этом все отклонения входных и выходных параметров, принимаемые в качестве воздействующих сигналов и реакций на них, вычислялись относительно избранного опорного режима. Работа продолжает исследования, начатые в [22]. Результатам исследования моделей вида (17), (18) посвящен следующий раздел.

3. Вычислительный эксперимент

Так как априори информация о том, является ли моделируемая динамическая система стационарной или нет, вообще говоря, неизвестна, проведем апостериорный анализ. Для этого рассмотрим двумерный отклик системы на пробные сигналы хШ1 (£) = е(Ь) — е(Ь — шг). Если окажется, что К'г = —К' , то это означает, что К1(1,ш1) = К'(Ь — шг) и, следовательно, система стационарна.

На рис. 1 приведены графики е = \К 11г\ — \К1 ' \ пр и £ = Ш'. Графи к е' получен в случае, когда ядра Вольтерра из (17) построены для У0 = 8 (м/с), а е2 — для У0 = 10 (м/с) при АЬ = 10 (град). На рис. 2 приведены графики е = \К2гг\ — \К21ш1 \ при £ = График ег получен в случае, когда ядра Вольтерра из (18) построены для Ь0 = 10 (град), е2 — дл я Ь0 = 20 (град), ез — дл я Ь0 = 0 (град) при АУ = 5 (м/с).

/ / г-"1

• / / / 1 1 1

1 1 1 е--< >—о— >--в--< >~в— ь-в—1

1 I 1 1 1 I / ! / / /

1 К

II II II а в —в— е

Г| г > -Е > > > > > [>

г

1 >1 к у 1--в--< >--в—« >--в—< 1—э--<>

--в--' Е1 - е2

* Ц/ | Г/ 1 V 1

2 4 6 8 10 12 14 16 18 20

ее

г

е

Е

2

4

6

8

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

Заключение

В данной статье рассмотрен новый способ построения квадратичного полинома Вольтерра для описания нестационарной динамической системы типа «вход-выход:». Основная идея состоит в модификации метода разделения отклика интегральной модели на составляющие, которая основана на использовании необходимых условий разрешимости соответствующих многомерных интегральных уравнений Вольтерра I рода. Приведены результаты численных расчетов для отгонной» модели ветроэнергетической установки с горизонтальной осью вращения, иллюстрирующие нестационарность исследуемого динамического объекта. Дальнейшее развитие работы связано с использованием кубичного полинома Вольтерра.

Работа проводилась при частичной финансовой поддержке гранта РФФИ, проект № 15-01-01J/25a.

Литература

1. Веников, В.А. Кибернетические модели электрических систем / В.А. Веников, О.А. Суханов. - М.: Энергоиздат, 1982. - 328 с.

2. Пупков, К.А. Функциональные ряды в теории нелинейных систем / К.А. Пупков, В.И. Капалин, А.С. Ющенко. - М.: Наука, 1976. - 448 с.

3. Stegmayer, G. Comparison of Volterra Models Extracted from a Neural Network for Nonlinear Systems Modeling / G. Stegmayer // Lecture Notes in Computer Science. - 2005. - V. 3697. -P. 457-463.

4. Tong Zhou, G. Nonlinear Channel Identification and Performance Analysis with PSK Inputs / G. Tong Zhou, G.B. Giannakis // First IEEE Signal Processing Workshop on Signal Processing Advances in Wireless Communications, New York. - 1997. - P. 337-340.

5. Chen, C.H. Fifth-Order Volterra Kernel Estimation for a Nonlinear Communication Channel with PSK and QAM Inputs / C.H. Chen, E.J. Powers // Ninth IEEE Signal Processing Workshop on Statistical Signal and Array Processing, New York. - 1998. - P. 435-438.

6. Lin, J.N. 2-D Adaptive Nonlinear Equalizers / J.N. Lin, R. Unbehauen // Proceedings of EUSIPCO-92, Sixth European Signal Processing Conference, Brussels. - V. 1. - 1992. -P. 135-138.

7. Minu, K.K. Volterra Kernel Identification by Wavelet Networks and Its Applications to Nonlinear Nonstationary Time Series / K.K. Minu, C. Jessy John // Journal of Information and Data Management. - 2012. - V. 1, № 1. - P. 4-9.

8. Вольтерра, В. Теория функционалов, интегральных и интегро-дифференциальных уравнений / В. Вольтерра. - М.: Наука, 1982. - 304 с.

9. Frechet, М. Sur les Funktionnoles Continues / M. Frechet // Ann. de PEcole Normale Sup. -1910. - V. 27. - P. 193-216._

. Bulletin of the South Ural State University. Ser. Mathematical Modelling, Programming

& Computer Software (Bulletin SUSU MMCS), 2015, vol. 8, no. 4, pp. 40-49

10. Апарцин, А.С. Неклассические уравнения Вольтерра I рода в интегральных моделях динамических систем: теория, численные методы, приложение: дис. ... д-ра физ.-мат. наук / А.С. Апарцин. - Иркутск: Иркутский гос. ун-т, 2000.

11. Данилов, Л.В. Теория нелинейных динамических цепей / Л.В. Данилов, Л.Н. Матханов,

B.С. Филиппов. - М.: Энергоиздат, 1990. - 256 с.

12. Льюнг, Л. Идентификация систем. Теория для пользователя / Л. Льюнг. - М.: Наука, 1991. - 432 с.

13. Apartsyn, A.S. Modeling of Nonlinear Dynamic Systems with Volterra Polynomials: Elements of Theory and Applications / A.S. Apartsyn, S.V. Solodusha, V.A. Spiryaev // International Journal of Energy Optimization and Engineering. - 2013. - V. 2, № 4. - P. 16-43.

14. Апарцин, А.С. Неклассические уравнения Вольтерра I рода: теория и численные методы / А.С. Апарцин. - Н.: Наука, 1999. - 193 с.

15. Апарцин, А.С. Об оптимизации амплитуд тестовых сигналов при идентификации ядер Вольтерра / А.С. Апарцин, С.В. Солодуша // Автоматика и телемеханика. - 2004. -№ 3. - С. 116-124.

16. Солодуша, С.В. Численное моделирование динамики теплообмена модифицированным квадратичным полиномом Вольтерра / С.В. Солодуша // Вычислительные технологии. - 2013. - Т. 18, № 2. - С. 83-94.

17. Апарцин, А.С. К идентификации нелинейных нестационарных динамических систем / А.С. Апарцин // Краевые - Иркутск: ИГУ, 1997. - С. 91-99.

18. Сидоров, Д.Н. Методы анализа интегральных динамических моделей: теория и приложения / Д.Н. Сидоров. - Иркутск: ИГУ, 2013. - 293 с.

19. Пронин, Н.В. Модель ветрогенератора ВЭУ-3 в пакете MATLAB / Н.В. Пронин, А.С. Мартьянов // Вестник ЮУрГУ. Серия: Энергетика. - 2012. - № 37 (296). -

C. 143-145.

20. Perdana, A. Dynamic Response of Grid-Connected Wind Turbine with Doubly Fed Induction Generator During Disturbances / A. Perdana, O. Carlson, J. Persson // Proc. of IEEE Nordic Workshop on Power and Industrial Electronics. - Trondheim, 2004.

21. Sedaghat, A. Aerodynamic Design of a 300 kW Horizontal Axis Wind Turbine for Province of Semnan / A. Sedaghat, M. Mirhosseini // Energy Conversión and Management. - 2012. -V. 63. - P. 87-94.

22. Суслов, К.В. Повышение качества электрической энергии при управлении элементами активно-адаптивных сетей / К.В. Суслов, Д.О. Герасимов, С.В. Солодуша // Управление качеством электрической энергии: сб. тр. междунар. науч.-практ. конф. - М.: Нац. исслед. ун-т МЭИ, 2014. - С. 191-198.

Светлана Витальевна Солодуша,

кандидат физико-математических наук, заведующий лабораторией, лаборатория «Неустойчивые задачи вычислительной математики:», Институт систем энергетики им. Л.А. Мелентьева СО РАН (г. Иркутск, Российская Федерация), solodusha@isem.irk.ru.

Дмитрий Олегович Герасимов, доцент, кафедра «Электроснабжение и электро:

(г. Иркутск, Российская Федерация), gerasimovdo@mail.ru. Константин Витальевич Суслов,

кандидат технических наук, .доцент. кафедра «Электроснабжение и электротехника», Иркутский национальный исследовательский технический университет (г. Иркутск, Российская Федерация), souslov@istu.edu.

Поступила в редакцию 25 января 2015 г.

MSC 93A30, 45D05 DOI: 10.14529/mmp150404

Construction of an Integral Model

by the Example of Wind Turbine Dynamics

S. V. Solodusha, Melentiev Energy Systems Institute Siberian Branch of the Russian

Academy of Sciences, Irkutsk, Russian Federation, solodusha@isem.irk.ru,

D. 0. Gerasimov, Irkutsk National Research Technical University, Irkutsk,

Russian Federation, gerasimovdo@mail.ru,

K. V. Suslov, Irkutsk National Research Technical University, Irkutsk, Russian

Federation, souslov@istu.edu

This study addresses the application of Volterra integral-power series to describe the nonlinear dynamic "input-output" systems. The universality of this mathematical tool makes it possible to create a software for computer experiments. This study is a continuation of the research on the identification of the Volterra kernels, which was started at the Energy Systems Institute, Siberian Branch of the Russian Academy of Sciences. The first part of the paper presents a new algorithm for the identification of the second-degree Volterra polynomials for the systems which can be used for an experiment based on the test sets of disturbances. In the second part the numerical calculation results for a "reference" dynamic system are given. The reference system is represented by a model of a horizontal-axis wind turbine. Quadratic Volterra polynomials are constructed. They describe the nonlinear dynamics of the angular velocity of the wind turbine components, depending on a blade lean angle and wind speed. The Volterra kernels were practically identified with respect to some chosen stationary state of the simulated system.

Keywords: nonlinear dynamic system; quadratic Volterra polynomial; horizontal - axis wind turbine.

References

1. Venikov V.A., Sukhanov O.A. Kiberneticheskie rnodeli elektricheskikh sistem [Cybernetic Models of Electric Power Systems]. Moscow, Energoizdat, 1982. 328 p.

2. Pupkov K.A., Kapalin V.I., Yushenko A.S. Funktsional'nye ryady v teorii nelineynykh sistem [Functional Series in the Theory of Non-Linear Systems]. Moscow, Nauka, 1976. 448 p.

3. Stegmayer G. Comparison of Volterra Models Extracted from a Neural Network for Nonlinear Systems Modelling. Lecture Notes in Computer Science, 2005, vol. 3697, pp. 457-463. DOI: 10.1007/11550907,72

4. Tong Zhou G., Giannakis G.B. Nonlinear Channel Identification and Performance Analysis with PSK Inputs. First IEEE Signal Processing Workshop on Signal Processing Advances in Wireless Communications, New York, 1997, pp. 337-340.

5. Chen C.H., Powers E.J. Fifth-Order Volterra Kernel Estimation for a Nonlinear Communication Channel with PSK and QAM Inputs. Ninth IEEE Signal Processing Workshop on Statistical Signal and Array Processing, New York, 1998, pp. 435-438.

6. Lin J.N., Unbehauen R. 2-D Adaptive Nonlinear Equalizers. Proceedings of EUSIPCO-92, Sixth European Signal Processing Conference, vol. 1, Brussels, 1992, pp. 135-138.

7. Minu K.K., Jessy John C. Volterra Kernel Identification by Wavelet Networks and Its Applications to Nonlinear Nonstationary Time Series. Journal of Information and Data Management, 2012, vol. 1, no. 1, pp. 4-9.

8. Volterra V. Teoriya funktsionalov, integral'nykh i integro-differentsial'nykh uravneni [A Theory of Functionals, Integral and Integro-Differential Equations]. Moscow, Nauka, 1982.

9. Frechet M. Sur les Funktionnoles Continues. Ann. de l'Ecole Nórmale Sup, 1910, vol. 27, pp.193-216.

10. Apartsyn A.S. [Nonclassical Volterra Equations of the First Kind in Integral Models of Dynamic Systems: Theory, Numerical Methods, Application. Dissertation of the Doctor of Physical and Mathematical Sciences]. Irkutsk, Irkutsk State University, 2000.

11. Danilov L.V., Matkhanov L.N., Filippov V.S. Teoriya nelineynykh dinamicheskikh tsepey [A Theory of Non-Linear Dynamic Circuits]. Moscow, Energoizdat, 1990.

12. Ljung L. Identifikatsiya sistem. Teoriya pol'zovatelya [System Identification. Theory for the User]. Moscow, Nauka, 1991.

13. Apartsyn A.S., Solodusha S.V., Spiryaev V.A. Modelling of Nonlinear Dynamic Systems with Volterra Polynomials: Elements of Theory and Applications. International Journal of Energy Optimization and Engineering, 2013, vol. 2, no 4, pp. 16-43. DOI: 10.4018/ijeoe.2013100102

14. Apartsyn A.S. Nonclassical Linear Volterra Equations of the First Kind. Utrecht, Boston, VSP, 2003.

15. Apartsyn A.S., Solodusha S.V. Test Signal Amplitude Optimization for Identification of the Volterra Kernels. Automation and Remote Control, 2004, vol. 65, no 3, pp. 464-471. DOI: 10.1023/B:AURC.0000019379.43119.dO

16. Solodusha S.V. [Numerical Modelling of Heat Exchange Dynamics by Modified Quadratic Volterra Polynomial]. Vychislitel'nye tekhnologii [Computational Technologies], 2013, vol. 18, no 2, pp. 83-94.

17. Apartsyn A.S. [By the Identification of Nonlinear Nonstationary Dynamic Systems]. Kraevye zadachi [Boundary Value Problems], Irkutsk, Irkutsk State University, 1997, pp. 91-99.

18. Sidorov D.N. Medody analiza integralnykh dinamicheskikh modeley: teoriya i prilozheniya [Methods of Analysis of Integrated Dynamic Models: Theory and Applications]. Irkutsk, Irkutsk State University, 2013.

19. Pronin N.V., Martyanov A.S. Model of Wind Turbine VEU-3 in the Package MATLAB. Bulletin of the South Ural State University. Series: Power Engeneering, 2012, no. 37

20. Perdana A., Carlson O., Persson J. Dynamic Response of Grid-Connected WTind Turbine with Doubly Fed Induction Generator During Disturbances. Proc. of IEEE Nordic Workshop on Power and Industrial Electronics, Trondheim, 2004.

21. Sedaghat A., Mirhosseini M. Aerodynamic Design of a 300 kW Horizontal Axis WTind Turbine for Province of Semnan. Energy Conversion and Management, 2012, vol. 63, pp. 87-94.

22. Suslov K.V., Gerasimov D.O., Solodusha S.V. [Increasing Power Quality During Control of Smart Grid Elements]. Upravlenie kachestvom elektricheskoy energii [Proceedings of the International Conference Power Quality Management], Moscow, National Research University "MPEI", 2014, pp. 191-198.

Received January 25, 2015