ИНФОРМАЦИОННО-УПРАВЛЯЮЩИЕ СИСТЕМЫ ^
УДК 681.5
с1о1:10.15217/1ззп1684-8853.2017.5.62
НЕЧЕТКОЕ УПРАВЛЕНИЕ СОЛНЕЧНОЙ БАТАРЕЕЙ
М. В. Буракова, канд. техн. наук, доцент
В. Ф. Шишлаков3, доктор техн. наук, профессор
аСанкт-Петербургский государственный университет аэрокосмического приборостроения,
Санкт-Петербург, РФ
Постановка проблемы: солнечная батарея является одним из основных инструментов для создания возобновляемых источников энергии. Нелинейные характеристики солнечной батареи как объекта управления требуют совершенствования алгоритмов управления для эффективной выработки энергии в условиях меняющейся внешней среды. Цель исследования: разработка имитационной модели солнечной батареи для синтеза и исследования нечеткого алгоритма управления точкой максимальной мощности. Результаты: рассмотрены варианты моделирования солнечной батареи с помощью библиотеки физического моделирования в 51тзсаре и с помощью графического программирования в Э1тиНпк. Модели были использованы при синтезе системы управления точкой максимальной мощности солнечного модуля. Исследованы два варианта управления: алгоритм возмущения и наблюдения и нечеткое логическое управление. Показано, что нелинейный закон управления, обеспечиваемый нечетким регулятором, позволяет сократить время переходного процесса, исключив колебания выходной мощности. Практическая значимость: применение нечеткого управления может оказаться полезным при проектировании систем управления широким кругом объектов солнечной энергетики.
Ключевые слова — солнечная батарея, моделирование, метод возмущения и наблюдения, нечеткое управление.
Введение
Использование возобновляемых источников электрической энергии является перспективным направлением современной энергетики [1, 2]. В рамках этого направления солнечная энергия занимает ключевые позиции в силу своей практической неисчерпаемости и доступности во всех точках поверхности Земли. Солнечная энергетика стремительно завоевывает позиции на мировом энергетическом рынке. Фотоэлектрические системы способны непосредственно преобразовывать энергию солнечной радиации в электрическую энергию без выброса вредных веществ в атмосферу. Они отличаются высокой надежностью, поскольку не имеют подвижных частей. Солнечные батареи пригодны как для выработки энергии в больших масштабах, так и для решения проблем энергообеспечения отдельных зданий, энергоснабжения автономных подвижных объектов, при создании солнечных холодильных систем и во многих других приложениях [3-5]. Вариации солнечной активности могут быть достаточно точно прогнозируемы, но главный недостаток солнечной энергетики заключается в зависимости от погодных и климатических условий. Для обеспечения стабильной выработки электроэнергии необходимо управлять работой солнечной батареи, поскольку мощность, генерируемая солнечной батареей, зависит от ее напряжения, значение которого на выходе постоянно изменяется в соответствии с освещенностью, временем суток, температурой панели. Контроллер солнечной батареи должен работать в режиме обеспечения максимальной выходной мощности
при вариациях температуры и солнечной радиации.
В настоящее время известны различные методы отслеживания точки максимальной мощности. Наиболее распространенным является метод возмущения и наблюдения [6, 7]. Перспективным для решения этой задачи является использование нечеткой логики [8, 9] и искусственных нейронных сетей [10]. В настоящей работе выполнено сравнение метода возмущения и наблюдения с алгоритмом нечеткого управления точкой максимальной мощности. При анализе методов управления используются имитационные компьютерные модели, разработанные средствами пакета Ма^аЬ.
Математическое моделирование солнечной батареи
Фотоэлектрический элемент (ФЭ) может быть представлен в виде функционального блока, имеющего внешние, внутренние и выходные параметры (рис. 1). К внешним параметрам относят освещенность и температуру ФЭ (в и Т). К вну-
G Ux.x* 1к.з U
I
T
P
Рис. 1. Параметры ФЭ Fig. 1. Photocell parameters
тренним параметрам относят напряжение холостого хода (Ux x) и ток короткого замыкания (1к з). Выходные параметры — выходное напряжение (U), ток нагрузки (I) и выходная мощность (P).
Основной характеристикой ФЭ является вольт-амперная характеристика — зависимость между током нагрузки и напряжением при постоянных значениях температуры и интенсивности поступающего солнечного излучения.
Имитационное моделирование солнечных батарей необходимо для синтеза систем управления, позволяющих эффективно вырабатывать энергию в условиях меняющейся внешней среды. Пакет MatLab обладает широкими возможностями для моделирования систем солнечной энергетики. Решить поставленную задачу можно двумя способами: с использованием Simulink MatLab или с помощью Simscape MatLab [11, 12].
Simulink MatLab является графическим языком программирования, который позволяет представить динамическую систему в виде набора блоков, описываемых определенными математическими зависимостями. Библиотека Simscape
,rO
+V
CD
rÙ Solar Cell ■
Solar Celli
■
rÙ Solar Cell ■
Solar Cell3
■
Solar Cell4
rÙ Solar Cell5 ■
-V 3
Рис. 2. Объединение ФЭ в Simscape Fig. 2. Combining photocells in Simscape
позволяет использовать физический язык моделирования, в котором имитационная модель объекта формируется с помощью библиотеки физических компонентов, которые взаимодействуют друг с другом путем обмена энергией через двунаправленные порты связи. Одновременное использование Simulink и Simscape позволяет выполнить верификацию имитационной модели солнечной батареи.
Для получения большого выходного напряжения ФЭ подключаются последовательно, а для достижения большого выходного тока ячейки соединяются параллельно. Модель ФЭ в библиотеке Simscape имеет три контакта: уровень солнечной радиации и токовые вход и выход. Последовательное соединение шести ФЭ показано на рис. 2.
Аналогично можно сформировать массив из 36 ФЭ, входящих в состав стандартного солнечного модуля (рис. 3).
Схема исследования солнечной батареи в Simscape, с помощью которой можно получить вольт-амперную и мощностную характеристики солнечного модуля, представлена на рис. 4.
Полученные в результате вычислительного эксперимента зависимости показаны на рис. 5, а и б.
Полученные характеристики являются основными для солнечного модуля, однако могут представлять интерес также температурные зависимости, которые нельзя получить, используя мо-
О
+V
+V
6 cells1
+V
6 cells3
-V
А
+V
-V
6 cells5
+V
Ir
-V
6 cel ls2
+V
Ir
-V
6 cel ls4
+V
Ir
-V
6 cells6
Рис. З. Солнечный модуль Fig. 3. The solar module
- Electrical Reference
Рис. 4. Блок-схема эксперимента в Simscape Fig. 4. Block diagram of the experiment in Simscape
а)
5 4,5 4
3,5 I, А 3
2,5
2
1,5
1
0,5
10
U, В
15
20
б)
P, Вт
90 80 70 60 50 40 30 20 10 0
1
2 —\ Л
3 \
4 " Л1 \\
>
и11
10
U, В
15
20
Рис. 5. Вольт-амперная (а) и мощностная (б) характеристика солнечного модуля при различных значениях освещенности: 1 — 1200 Вт/м2; 2 — 1000 Вт/м2; 3 — 800 Вт/м2; 4 — 600 Вт/м2
Fig. 5. Volt-ampere (а) and power (б) characteristic of the solar module for various values of illumination: 1 — 1200 W/m2; 2 — 1000 W/m2; 3 — 800 W/m2; 4 — 600 W/m2
Rn I
G
Рис. 6. Электрическая схема замещения ФЭ: G — солнечное излучение; Rn и Вш — последовательное и параллельное (шунтирующее) сопротивление (паразитные параметры ФЭ)
Fig. 6. Electrical circuit of photocell replacement: G — Solar radiation; Вп, Вш — series and shunt resistance
дель ФЭ в Ятвсаре. Для исследования влияния температуры на параметры солнечного модуля можно использовать язык графического моделирования, опирающийся на аналитическое описание системы.
Математическая модель ФЭ [13, 14] основана на анализе структуры, представленной на рис. 6.
В соответствии с первым законом Кирхгофа можно записать
1 ~ 1ф 1д
(1)
где 1ф — фототок; 1д — ток диода; 1ш — ток шунтирующего сопротивления. По второму закону Кирхгофа
=
U + IRa
(2)
Величина тока, протекающего через р-п-пере-ход (через диод), зависит от приложенного к нему напряжения согласно формуле
С Л
'д = J0
exp
q(U + IRz) AkT
-1
(3)
где 10 — обратный ток насыщения; q — заряд электрона; и — выходное напряжение; А — диодный фактор (константа идеальности); I — выходной ток; к — постоянная Больцмана; Т — абсолютная температура ФЭ, К.
Величина фототока описывается формулой
О
Ль Рн ) - ^г" ^к.з (гн );
G.
/ф(Т) = /ф(Тн)(1+Tk(TB -т)),
(4)
где Тн — номинальная температура, 25 °C; G — освещенность солнечного элемента, Вт/м2; Gн —
5
5
Таблица 1. Параметры фотоэлектрического модуля Table 1. Parameters of the photoelectric module
Параметр Значение
Последовательное сопротивление Дп 0,001 Ом
Заряд электрона д 1,6 • 1019 Кл
Постоянная Больцмана к 1,38 • 10-23 Дж/К
Фактор идеальности диода А 1,6
Температурный коэффициент Тк 0,68 мА
Параллельное сопротивление Вш 1 Ом
Номинальная освещенность 1000 Вт/м2
Номинальная температура Тн 25 °С
Количество последовательно соединенных ячеек N 36
а) 5 4
I, А
10
15
20
U, В
б) 80 70 60 50
P Вт 40 30 20
10 0
10 15 U,B
20
2 3 1
4 \
25
1 2
32
4 \
\\
у
/
25
Рис. 7. Вольт-амперная характеристика (а) и изменение мощности (б) солнечной батареи при различных значениях температуры: 1 — 25 °С; 2 — 35 °С; 3 — 45 °С; 4 — 55 °С Fig. 7. Volt-ampere characteristic (а) and change in the power (б) of a solar battery for different temperature values: 1 — 25 °С; 2 — 35 °С; 3 — 45 °С; 4 — 55 °С
номинальная освещенность, Вт/м2; Т — температура окружающей среды; Тк — температурный коэффициент фототока.
Температурный коэффициент вычисляется по формуле
™ _ ^к.з (Т2 ) ~ ^к.з №)
1ь —-.
й Т2 - Т-у
Подставляя (2)-(4) в (1), получаем выражение для описания тока солнечного элемента
I = ^- (1+ЩТ - Тн ))-
- In
exp
q(U + IRa)
AkT
-1
U + IRn Rm
(5)
Можно считать, что Дш стремится к бесконечности, тогда в (5) можно пренебречь последним слагаемым:
I = *к.з — (1+Щт - Тв ))-
- 1Г
exp
q(U + IRU)
AkT
Л ^ -1
(6)
где I0 — обратный ток насыщения, описываемый формулой
i m Л3 (iV Л
Iq = Io (Тж)
т
v н у
exp
qEL Ak
; j
(7)
где — полоса пропускания диода.
Обратный ток насыщения при номинальной температуре рассчитывается по формуле
1о (Тж) =-
-^к.з (Г=)
exp
f g^x.x (Г= ) ^ AkT„
(8)
-1
При моделировании были приняты параметры, соответствующие солнечному модулю ОТ650К (табл. 1).
Формулы (1)-(8) позволяют собрать в 81ти1тк модель солнечного модуля и выполнить вычислительные эксперименты при любых заданных условиях.
Влияние температуры на вольт-амперную и мощностную характеристики солнечной батареи исследовано на рис. 7, а и б.
Результаты моделирования показывают, что характеристики солнечного модуля существенно нелинейны и обнаруживают сильную зависимость от внешних условий — солнечной радиации и температуры окружающей среды.
0
5
0
5
Управление точкой максимальной мощности
Принцип управления точкой максимальной мощности солнечной батареи иллюстрирует рис. 8.
На выходе устройства управления располагается широтно-импульсный модулятор, который изменяет скважность импульсов, поступающих на преобразователь напряжения.
Задача обеспечения максимальной мощности может решаться с помощью метода возмущения и наблюдения, который предполагает уменьшение или увеличение скважности управляющих импульсов Б с последующей оценкой изменения выходной мощности [8, 9]. Если в момент времени к мощность Р(к) и напряжение и(к) больше, чем в предыдущий момент времени, то Б увеличивается на ДБ, иначе уменьшается на ДБ. Величина ДБ подбирается экспериментально.
Закон управления можно описать с помощью правил:
если ДР < 0 и и(к) > и(к - 1), то Б(к + 1) = = Б(к) - ДБ;
если ДР < 0 и и(к) < и(к - 1), то Б(к + 1) = = Б(к) + ДБ;
■ Рис. 8. Система управления солнечной батареей
■ Fig. 8. Solar battery control system
если ДР > 0 и U(k) < U(k - 1), то D(k + 1) = = D(k) - ДБ;
если ДР > 0 и U(k) > U(k - 1), то D(k + 1) = = D(k) + ДБ.
Реализация алгоритма возмущения и наблюдения в MatLab представлена на рис. 9.
Нетрудно заметить, что метод возмущения и наблюдения является, по существу, разновидностью метода подъема по градиенту. Скважность изменяется дискретно, поэтому выбор ДБ оказывает заметное влияние на переходные процессы. Уменьшение значения ДБ приводит к росту длительности переходного процесса, увеличение ДБ вызывает колебания тока и напряжения (рис. 10).
Добиться переменного шага при изменении скважности импульсов позволяет использование нелинейного закона управления, реализуемого с помощью нечеткого логического регулятора (НЛР).
Основой нечеткого логического управления являются продукционные правила, связывающие нечеткие посылки и заключения. Для многих объектов эти правила неочевидны, и для их синтеза требуется применять процедуры обучения [15-17]. Однако в рассматриваемой задаче правила управления могут быть сформулированы эвристически.
В задаче нечеткого управления точкой максимальной мощности солнечной батареи можно выбрать такие же входные и выходные переменные, что и в методе возмущения и наблюдения: входами НЛР являются приращение напряжения AU(k) и приращение мощности AP(k), а выходом — приращение скважности импульсов АБ(к).
Для описания входов и выходов НЛР будем использовать лингвистические переменные с термами: {NB, NS, Z, PS, PB}, где NB — «negative
d> P
_2_ U
Л
lTs
Subtract
Constant2
Memory1
П
Product Sign
Memory
Subtract1
0.001
Constant5
3
Switch1
Product1
Constant3
Subtract2
Л
№
Saturation2
D
0 J Switch2
Memory2
Constant4
Рис. 9. Блок-схема алгоритма возмущения и наблюдения в Simulink MatLab Fig. 9. Block diagram of the perturb and observe algorithm in Simulink MatLab
U, В
2G 18 1б 14 12 1G 8 б 4 2 G
2
1.......7...............
/1
G,5
1,б
2,б
t, c
Рис. 10. Изменение выходного напряжения при AD = 0,01 (линия 1) и AD = 0,1 (линия 2) Fig. 10. Change in the output voltage at AD = = 0,01 (line 1) and at AD =0,1 (line 2)
g G,8 с о
ьн M £ G,6
е
S I
и р
н
-G,1 -G,G6 -G,G2 G,G2 G,G6 Вход 1, вход 2, выход
G,1
Рис. 11. Лингвистическое описание входов и выхода НЛР
Fig. 11. Linguistic description of inputs and output of fuzzy logic controller
big», NS — «negative small», Z — «zero», PB — «positive big», PS — «positive small» (рис. 11).
Структура НЛР представлена на рис. 12.
При составлении управляющих правил используется принцип изменения скважности импульсов с переменным шагом: при малых изменениях мощности изменение скважности замедляется, а при больших — растет.
Например:
1. Если напряжение сильно увеличилось и мощность тоже сильно увеличилась, то следует продолжать сильно увеличивать напряжение, т. е. сильно уменьшить скважность.
2. Если напряжение сильно увеличилось, а мощность слегка увеличилась, то следует слегка увеличивать напряжение, т. е. слегка уменьшить скважность.
Управляющие правила могут быть сведены в табл. 2.
Схема сравнения алгоритмов слежения за точкой максимальной мощности представлена на рис. 13. Будем считать, что время квантования существенно выше периода широтно-импульс-ной модуляции, тогда преобразователь может быть представлен эквивалентным сопротивлением. Инерционность преобразователя может быть учтена введением апериодического звена на его выходе.
Таблица 2. Управление скважностью импульсов Table 2. Pulse width control
AP/AU NB NS Z PS PB
NB PB PS NS NS NB
NS PS PS PS NS NS
Z PS PS Z NS NS
PS NS NS PS PS PB
PB NB NS NS PS PB
G
1
2
1
G
Œ>
P
l_T
Add1
JV
Zero-Order Saturation1 Hold1
Transport Delay
Œ>
U
Add3
JV
h
Fuzzy Logic Controller
Zero-Order Hold
Saturation
+ +
w
Add4 Saturation2
1 z
Out1
Unit Delay2
Transport Delay1
Рис. 12. Блок-схема НЛР в Simulink MatLab
Fig. 12. Block diagram of the fuzzy logic controller in Simulink MatLab
Subsystem2
+v
_^ S PS Ir
—► -V
Simulink-PS
Converter1
Signal 1
Signal Builder
f(x)=0
Solver Configuration
FLC
•©r
PS S —►
Current Sensor
S PS
Simulink-PS Converter Variable Resistor
PS-Simulink Converter
[I]
Goto2
PS S —►
Voltage Sensor
PS-Simulink Converter1
Scope
К [U]
Goto1
Scope1
Product
[P] X
— Electrical Reference
Goto
P dU U
P&O
U Out1 P
—
Manual Switch
[U]
From2
О
Scope2
In1 1-1
Out1 1 w
In2 Scope4
Subsystem
Scope3
Рис. 13. Схема сравнения алгоритмов управления точкой максимальной мощности Fig. 13. Comparison of control algorithms for the maximum power point tracking
а)
4,1 4,05
I, А 4 3,95 3,9 3,85 3,8
Y 2
V
i
б)
P, Вт
0 0,05 0,1 0,15 0,2 0,25 0,3 0,35 0,4 t, c
0 0,05 0,1 0,15 0,2 0,25 0,3 0,35 0,4 t, c
Рис. 14. Изменение тока нагрузки (а) и мощности (б) при управлении с помощью НЛР (линия 1) и с помощью
алгоритма возмущения и наблюдения (линия 2) Fig. 14. Change of the load (а) and power (б) current during control by means of fuzzy logic controller (line 1) and with the help of the perturb and observe algorithm (line 2)
Переходные процессы при скачкообразном изменении освещенности солнечной батареи при разных вариантах реализации алгоритма управления показаны на рис. 14, а и б. Из анализа сравнения графиков видно, что использование НЛР обеспечивает значительное преимущество по сравнению с алгоритмом возмущения и наблюдения, позволяя сократить время переходного процесса, избежав при этом колебаний выходной мощности.
Заключение
Приведенные результаты вычислительных экспериментов соответствуют техническим характеристикам солнечного модуля МР650К, па-
Литература
1. Messenger R. A., Ventre J. Photovoltaic Systems Engineering. — CRC PRESS, 2004. — 435 p.
2. Markvart T. Solar Electricity. — Wiley, 2000. — 280 p.
3. Winter C. J., Sizmann R. L., Vant-Hull L. L. Solar Power Plants: Fundamentals, Technology, Systems, Economics. — Springer-Verlag, 1991. — 425 p.
4. Garg H. P. Solar Energy: Fundamentals and Applications. — Tata McGraw-Hill Education, 2000. — 434 p.
5. Kalogirou S. A. Solar Energy Engineering: Processes and Systems. — Academic Press, 2013. — 840 p.
6. Ropp M. E., Hohm D. P. Comparative Study of Maximum Power Point Tracking Algorithms // Progress in Photovoltaics: Research and Applications. 2003. Vol. 11(1). P. 47-62.
7. Esram T., Chapman P. L. Comparison of Photovoltaic Array Maximum Power Point Tracking Techniques // IEEE Trans. Energy Convers. 2007. Vol. 22. N 2. P. 439-449.
8. Cheikh M. S., Larbes C., Tchoketch G. F. Maximum Power Point Tracking using a Fuzzy Logic Control Scheme // Revue des Energies Renouvelables. 2007. Vol. 10. N 3. P. 387-395.
9. Aredes M. A., França B. W., Aredes M. Fuzzy Adaptive P&O Control for MPPT of a Photovoltaic
раметры которого использовались при моделировании. Модель, полученная с помощью физического моделирования в Simscape, показала параметры, практически идентичные параметрам модели, собранной с помощью графического программирования Simulink. Обе модели были использованы при проектировании системы управления точкой максимальной мощности солнечного модуля. Были исследованы два варианта управления: алгоритм возмущения и наблюдения и нечеткое логическое управление. Показано, что нелинейный закон управления, обеспечиваемый НЛР, позволяет значительно улучшить качество переходных процессов. Применение нечеткого управления может оказаться полезным при проектировании систем управления широким кругом объектов солнечной энергетики.
Module // Journal of Power and Energy Engineering.
2014. N 2. Р. 120-129.
10. Hadjab M., Berrah S., Abid H. Neural Network for Modeling Solar Panel // International Journal of Energy. 2012. Iss. 1. Vol. 6. Р. 9-16.
11. Черных И. В. Моделирование электротехнических устройств в Matlab, SimPowerSystems и Simulink. — M.: ДМК Пресс, 2008. — 288 с.
12. Дьяконов В. П. MATLAB R2007/2008/2009 для радиоинженеров. — M.: ДМК Пресс, 2010. — 976 с.
13. Васильев А. М., Ландсман А. П. Полупроводниковые фотопреобразователи. — М.: Сов. радио, 1976. — 246 с.
14. Раушенбах Г. Справочник по проектированию солнечных батарей. — М.: Энергоатомиздат, 1983. — 360 с.
15. Бураков М. В., Кирпичников А. П. Нечеткий регулятор ПИД-типа для нелинейного объекта // Вестник Казанского технологического университета.
2015. Т. 18. № 4. C. 242-244.
16. Бураков М. В., Яковец О. Б. Нечеткое управление силовым гироскопическим прибором // Изв. вузов. Приборостроение. 2015. Т. 58. № 10. С. 157-166.
17. Бураков М. В., Коновалов А. С., Яковец О. Б. Эволюционный синтез нечетких регуляторов // Информационно-управляющие системы. 2015. № 6. С. 28-33. doi:10.15217/issn1684-8853.2015.6.28
UDC 681.5
doi:10.15217/issn1684-8853.2017.5.62 Fuzzy Control of Solar Power Plant
Burakov M. V.a, PhD, Tech., Associate Professor, bmv@sknt.ru Shishlakov V. F.a, Dr. Sc., Tech., Professor, svfmail@yandex.ru
aSaint-Petersburg State University of Aerospace Instrumentation, 67, B. Morskaya St., 190000, Saint-Petersburg, Russian Federation
Introduction: A solar battery is one of the main tools for the creation of renewable energy sources. Nonlinear characteristics of a solar battery as a control object require the improvement of control algorithms for the efficient generation of energy in a changing environment. Purpose: Our goal is to develop a simulation model of a solar battery for the synthesis and study of a fuzzy algorithm for the control of the maximum power point. Results: Some options have been considered for modeling a solar battery using a physical modeling library in Simscape, and using graphical programming in Simulink. These models were used in the synthesis of a system for the control of the maximum power point of a solar module. Two control options were studied: perturb and observe algorithm, and fuzzy logic control. It is shown that the nonlinear control law provided by the fuzzy controller allows you to reduce the time of the transient process by eliminating the oscillations of the output power. Practical relevance: The application of fuzzy control can be useful in the design of control systems for a wide range of solar energy facilities.
Keywords — Photovoltaic Module, Simulation, Perturb and Observe Method, Fuzzy Control.
References
1. Messenger R. A., Ventre J. Photovoltaic Systems Engineering. CRC PRESS, 2004. 435 p.
2. Markvart T. Solar Electricity. Wiley, 2000. 280 p.
3. Winter C. J., Sizmann R. L., Vant-Hull L. L. Solar Power Plants: Fundamentals, Technology, Systems, Economics. Springer-Verlag, 1991. 425 p.
4. Garg H. P. Solar Energy: Fundamentals and Applications. McGraw-Hill Education, 2000. 434 p.
5. Kalogirou S. A. Solar Energy Engineering: Processes and Systems. Academic Press, 2013. 840 p.
6. Ropp M. E., Hohm D. P. Comparative Study of Maximum Power Point Tracking Algorithms. Progress in Photovoltaics: Research and Applications, 2003, vol. 11(1), pp. 47-62.
7. Esram T., Chapman P. L. Comparison of Photovoltaic Array Maximum Power Point Tracking Techniques. IEEE Trans. Energy Convers., 2007, vol. 22, no. 2, pp. 439-449.
8. Cheikh M. S., Larbes C., Tchoketch G. F. Maximum Power Point Tracking using a Fuzzy Logic Control Scheme. Revue des Energies Renouvelables, 2007, vol. 10(3), pp. 387-395.
9. Aredes M. A., França B. W., Aredes M. Fuzzy Adaptive P&O Control for MPPT of a Photovoltaic Module. Journal of Power and Energy Engineering, 2014, no. 2, pp. 120-129.
10. Hadjab M., Berrah S., Abid H. Neural Network for Modeling Solar Panel. International Journal of Energy, 2012, iss. 1, vol. 6, pp. 9-16.
11. Chernykh I. V. Modelirovanie elektrotekhnicheskikh ustroistv v Matlab, SimPowerSystems i Simulink [Simulation of Electrical Devices in Matlab, SimPowerSystems and Simu-link]. Moscow, DMK Press Publ., 2008. 288 p. (In Russian).
12. Dyakonov V. P. MATLAB R2007/2008/2009 dlia radio-inzhenerov [MATLAB R2007/2008/2009 for Radio Engineers]. Moscow, DMK Press Publ., 2010. 976 p. (In Russian).
13. Vasilyev A. M., Landsman A. P. Poluprovodnikovye foto-preobrazovateli [Semiconductor Photovoltaics]. Moscow, Sovetskoe radio Publ., 1976. 246 p. (In Russian).
14. Raushenbakh G. Spravochnik po proektirovaniiu sol-nechnykh batarei [Guide to the Design of Solar Panels]. Moscow, Energoatomizdat Publ., 1983. 360 p. (In Russian).
15. Burakov M. V., Kirpichnikov A. P. Fuzzy PID Controller for Nonlinear Plant. Vestnik Kazanskogo tekhnologicheskogo universiteta, 2015, vol. 18, no. 4, pp. 242-244 (In Russian).
16. Burakov M. V., Yakovets 0. B. Fuzzy Control of Power Gyroscopic Unit. Izvestiia vuzov. Priborostroenie, 2015, vol. 58, no. 10, pp. 157-166 (In Russian).
17. Burakov M. V., Konovalov A. S., Yakovets O. B. Evolutionary Synthesis of Fuzzy Controllers. Informatsionno-uprav-liaiushchie sistemy [Information and Control Systems], 2015, no. 6, pp. 28-33 (In Russian). doi:10.15217/issn1684-8853.2015.6.28