CC BY
0Энергетика. Автоматика
УДК 621.311
Мясоедов Юрий Викторович
Амурский государственный университет г. Благовещенск, Россия e-mail: myv@amursu. ru Савина Наталья Викторовна Амурский государственный университет г. Благовещенск, Россия e-mail: nataly-savina@mail.ru Myasoedov Yuriy Victorovich Amur State University Blagoveshchensk, Russia e-mail: myv@amursu.ru Savina Natalia Victorovna Amur State University Blagoveshchensk, Russia e-mail: nataly-savina@mail.ru
МОДЕЛИРОВАНИЕ ТОКОВ РАСПРЕДЕЛЕННЫХ ИСТОЧНИКОВ ИСКАЖЕНИЯ КАЧЕСТВА ЭЛЕКТРОЭНЕРГИИ ДЛЯ АДАПТИВНОГО УПРАВЛЕНИЯ КАЧЕСТВОМ ЭЛЕКТРОЭНЕРГИИ
MODELING OF CURRENTS OF DISTRIBUTED SOURCES OF POWER QUALITY DISTORTION FOR ADAPTIVE POWER QUALITY CONTROL
Аннотация. В современных энергосистемах и в настоящее время, и при цифровой трансформации и интеллектуализации ухудшается качество электроэнергии. Отличительный признак - увеличение в них доли распределенных источников искажения качества электроэнергии. Выделены инновационные технологии, применение которых сопряжено с появлением новых распределенных источников искажения качества электроэнергии. Предложена обобщенная модель тока распределенных источников искажения качества электроэнергии и ее частные случаи для адаптивного управления качеством электроэнергии при разных видах неопределенности информации. Частные модели разработаны на основе совокупного использования вейвлет-анализа, теории случайного процесса и нечетких множеств.
Abstract. In modern energy systems, both at the present time and with digital transformation and intellectualization, the quality of electricity is deteriorating. A distinctive feature is the increase in the share of distributed sources of electricity quality distortion in them. Innovative technologies are identified, the use of which is associated with the emergence of new distributed sources of electricity quality distortion. The principle of constructing a generalized model is based on the selection of periodic components from a non-stationary random process due to the change in current over time, harmonics of direct and reverse sequence currents, a stationary random process with zero mathematical expectation, differentiable correlation functions and random noise. Particular models are developed on the basis of the combined use of wavelet analysis, random process theory and fuzzy sets.
Ключевые слова: качество электроэнергии, распределенный источник, адаптивное управление, ток, вейвлет-анализ, случайный процесс, нечеткое множество.
Key words: power quality, distributed source, adaptive control, current, wavelet analysis, random process, fuzzy set.
DOI: 10.22250/20730268_2023_103_28
Введение
В современных энергосистемах с нелинейной и несимметричной нагрузкой качество электроэнергии (КЭ), как правило, низкое. Это негативное явление для них, так как снижает функциональную надежность и эффективность работы. Наблюдаются заметный рост потерь электрической энергии, резкое сокращение срока службы электроустановок и кабельных линий, отказы и ложные срабатывания релейной защиты и автоматики, цифровых устройств, используемых как в энергосистемах, так и потребителями. Искажается и теряется информация при диспетчерском и технологическом управлении, становится недостоверным учет электроэнергии.
В то же время в России наблюдается усиление инновационной активности в энергетике, нормативно этот процесс закреплен в [1]. Ее ключевой аспект - цифровая трансформация и интеллектуализация электроэнергетики. Можно выделить инновационные технологии, применение которых повлечет за собой ухудшение КЭ, это возобновляемые источники и накопители энергии; электромобили; распределенная генерация; ряд технологий, придающих электрическим сетям свойства активно-адаптивных сетей. Данный факт подтверждается их принципом действия и конструктивными особенностями [2-6].
Исходя из сказанного, необходимы мониторинг и управление КЭ.
Особенность энергосистем в плане КЭ - множество источников его искажения со случайным изменением параметров режима, распределенных как во времени, так и территориально. Это обусловливает нестабильность показателей КЭ, которые представляют собой многомерные случайные процессы. Несмотря на то, что проблеме КЭ посвящено достаточно большое количество работ, в которых рассматриваются задачи его нормирования, оценка влияния низкого качества электроэнергии на объекты электроэнергетики и потребителей, технические средства и схемы решений, направленные на улучшение КЭ [например, 7-11], вопросы моделирования токов источников искажения КЭ с целью адаптивного управления им требуют развития.
Постановка задачи
Выбор эффективных мероприятий, параметров и законов управления технических средств улучшения КЭ, обеспечение условий электромагнитной совместимости для надежного функционирования электрических сетей энергосистем во многом предопределяются исходной информацией о режимах работы источников искажения КЭ. Различают стационарные и распределенные источники искажения качества электроэнергии (РИ). Стационарные источники не меняют своего территориального положения. Примером служат дуговые сталеплавильные печи, прокатные станы и т.п. Под распределенными понимаются такие источники искажения КЭ, положение которых меняется в пространстве (территориально), а режимы работы - случайным образом во времени. Ярким примером таких источников является электрифицированный железнодорожный транспорт - нечеткое множество электровозов, питание привода которых осуществляется однофазным переменным током, преобразуемым в постоянный посредством выпрямительных устройств, установленных непосредственно на электровозах.
Каждая тяговая подстанция, являясь источником искажения КЭ с постоянно изменяющейся нагрузкой (генератор), в то же время принимает высшие гармонические составляющие тока от соседних тяговых подстанций (рецептор). Все это обусловливает нестабильность показателей КЭ в элек-
трических сетях энергосистем, питающих тяговую нагрузку. При цифровизации и интеллектуализации электрических сетей доля распределенных источников искажения качества электроэнергии увеличится, при этом для моделирования их токов можно в качестве аналога взять электрические сети с тягой переменного тока, так как каждый источник распределенной генерации, являясь генератором высших гармоник тока и интергармоник, токов обратной последовательности, в то же время является рецептором, принимая искажающие токи от других источников распределенной генерации.
Методический подход одинаков. Для определения закономерностей их изменения прежде всего необходимо знать, как меняется ток совокупности распределенных источников искажения качества электроэнергии (РИ).
Режимная информация таких источников зависит от множества факторов: числа поездов на перегоне, графика движения, расположения поездов относительно друг друга, географических особенностей местности, электропотребления по пути движения и т.п. Эта информация обладает неопределенностью. В общем случае для анализа КЭ необходимо учитывать следующие виды неопределенностей: неполнота, недостоверность, нечеткость, неточность, ненадежность. В зависимости от видов неопределенности можно выделить четыре группы информации: детерминированная, вероятностная, нечеткая, интервальная.
Детерминированное представление информации о состоянии РИ для его анализа должно отражать наиболее вероятное (для стохастического описания режимных параметров), наиболее желательное (нечеткое описание) или среднее (интервальное описание) из множества значений, которые может принимать ток РИ. Как правило, детерминированная информация выражается номинальным током тягового трансформатора на подстанции или, в лучшем случае, его долей, идущей на питание контактной сети, которая задается произвольно, исходя из эмпирических наблюдений, или определяется по соотношению электропотребления обмоток средней и низкой стороны трансформатора. Такой подход не раскрывает физической природы режима работы РИ и возможен только при большой степени неопределенности.
Вероятностная информация отражает случайный характер изменения тока распределенных источников, их физическую природу. Исследуя статистику случайных процессов при полном объеме достоверного информационного потока, можно определить закономерности изменения узлового тока тяговой подстанции, найти характеристики случайного процесса и получить модель эквивалентного тока тяговой нагрузки или тока каждого источника распределенной генерации, каждого ВИЭ. При отсутствии достаточного объема статистических выборок, либо при потере информации, либо при ее нечеткости или недостоверности определение вероятностных характеристик случайного процесса затруднительно. Тогда информация является нечеткой. В этом случае модель тока может быть получена сочетанием теории случайных процессов и нечетких множеств.
Если закон распределения случайного процесса неизвестен и найти его не предоставляется возможным из-за неопределенности исходных данных либо их неполноты, связанной с отсутствием соответствующих измерительных средств или системы сбора и обработки информации, то модель тока распределенного источника можно получить, обрабатывая информационный поток с помощью интервального анализа. В данном случае определяются вероятностные интервалы, в которые входят исследуемые токи. При этом учитывается степень зависимостей, - например, между участвующими в суммарной нагрузке подстанции интервальными объектами - электровозами.
Исходя из вышесказанного, можно сделать вывод, что исходная информация, на основе которой принимается решение о КЭ и выборе технических средств, места их установки и режимах работы, является неполной и неточной. Неполнота и недостоверность обусловлены физическими, техническими и эргономическими факторами, что необходимо учитывать при моделировании тока распределенных источников искажения. В зависимости от степени полноты и достоверности информации и должна выбираться эта модель.
Разработка модели тока распределенного источника искажения качества электроэнергии
Обобщенная модель тока РИ должна строиться на основе выделения из нестационарного случайного процесса периодических составляющих, обусловленных изменением тока во времени, гармонических составляющих, вызываемых выпрямлением тока, стационарного случайного процесса с нулевым математическим ожиданием, дисперсией, дифференцируемыми корреляционными функциями и случайной помехи. Она имеет вид:
= 2т=11п=11тп с°5( и>тпЬ + (ртп) + + 2п=2 12пС05(Ш2П1 + <р2п) + 1(1) + 1а() + (1)
где 1тп, штп, <ртп - соответственно амплитуда, частота и фаза m-й периодической составляющей (период во времени) п-й гармоники тока прямой последовательности; ¡2П, ш2п, Ф2п - соответственно амплитуда, частота и фаза п-й гармонической составляющей тока обратной последовательности; т -номер гармоники во времени; п - номер гармоники по частоте; М - число учитываемых гармоник по времени; N - число учитываемых частотных гармоник; I (Ь) - стационарный случайный процесс с нулевым математическим ожиданием; I а¿(Ь) - математическое ожидание тока - среднее значение за рассматриваемый период; % - случайная помеха, представляющая собой шум. Это может быть белый шум, мощность которого не зависит от частоты, а автокорреляционной функцией является 5 - функция Дирака либо Гауссов шум. Согласно центральной теореме, сумма произвольных случайных процессов стремится к Гауссову процессу при возрастании числа слагаемых. При этом сходимость настолько быстра, что при числе слагаемых 5 или 6 результирующий процесс очень близок к Гауссову.
В зависимости от степени неопределенности и группы информационного потока обобщенная модель тока может принимать различные частные формы.
Для получения дискретной модели тока непрерывный процесс его изменения заменяется дискретной последовательностью выборок с некоторым временным интервалом АЬ, по которой с заданной точностью можно восстановить исходный процесс.
При этом целесообразно применять неравномерное квантование, т.е. АЬ - переменная величина. Процесс квантования описывается уравнением:
4( 1) = 1к = Ы(1к)Ла(1-1к), (2)
где ¡к(0 - квантованный параметр; N(1^ - число квантов; Ла(Ь — Ьк) - единичная функция, а оптимальная величина интервала дискретизации определяется с помощью теоремы Котельникова [14, 15].
Тогда токовую нагрузку можно описать дискретно. При этом каждая гармоническая составляющая как прямой, так и обратной последовательности тока представляется дискретной последовательностью:
1к = ¡(к). (3)
На интервалах стационарности корреляционная функция имеет вид:
Як = М1т+к1к = Я, (4)
а спектральная плотность равна
Г(А) = (2л)-121кетК(к), (5)
где А 6 [—л; л].
В этом случае модель тока каждой гармоники прямой и обратной последовательности при ее спектральном представлении имеет вид:
1к = £пе]ЛкГ(А)с1Ш(Л), (6)
а эквивалентное значение тока РИ определяется как среднеквадратичное значение сумм, полученных выше моделей:
т=№=112кп1+1."п=112кп2, (7)
где ¡к - модель тока п-й гармоники прямой последовательности; 1кп2 - то же, но обратной последовательности.
С учетом временных гармонических составляющих ток первой и высших гармоник как прямой, так и обратной последовательностей может быть описан моделью, полученной с помощью формулы Котельникова:
(8)
_ v sinn(t-k)
п(г-к)
корреляционная функция которой имеет вид:
= . (9)
При полном и достоверном информационном потоке частным случаем обобщенной модели (1) является модель, полученная в [16] и представленная стационарным случайным процессом и помехой:
К 0 = 1а + х(1 + КпзтК! + Ът) + (10)
где у - коэффициент вариации случайного процесса; Я(т) - нормированная корреляционная функция; Кп31 - случайный процесс изменения коэффициентов несинусоидальности тока; К21 - случайный процесс изменения коэффициента обратной последовательности тока.
Модели случайных процессов изменения показателей качества электроэнергии, разработанные в [17], полностью подходят для исследования тока распределенных источников.
Возможен другой подход к реализации тока РИ при вероятностном представлении информации - применение вейвлет-анализа в сочетании с теорией случайных процессов. При исследовании нестационарных случайных процессов изменения тока РИ, не имеющих четкого периодического характера (это видно из рис.1), наиболее эффективно использовать не тригонометрические, а некоторые локализованные во времени компактные базисы, коэффициенты разложения по которым сохраняют информацию об изменении параметров аппроксимируемого процесса - вейвлеты.
К(т)
-0,4
х, min
0 2
4
6 8 10 12 14 16 18 20 22 24 26 28 Рис. 1. Корреляционная функция тока.
В отличие от классических преобразований Фурье, вейвлет-анализ имеет базис функций, локализованных и по времени, I, и по частоте /. В результате случайный процесс можно анализировать одновременно в физическом и частотном пространствах. При этом физическое пространство представляется временем, координатой и может отражать как изменение процесса во времени, так и территориально, - например, в фидерной зоне (положение электровозов относительно тяговой подстанции и относительно друг друга). Процесс изменения тока, представленный в системе базисных функций вейвлета, можно записать в виде обобщенного ряда Фурье, подбирая для каждого практического случая наиболее эффективный вейвлет. По своему смыслу вейвлет-преобразование полностью соответствует преобразованию Фурье, но ядром интегрального преобразования вместо экспоненциальной функции е-]т служит вейвлет Ф(1,о), где а - масштаб осцилляции. Использование вейвлетов дает возможность сократить объем информации, необходимой для выявления закономерностей случайно-
го процесса, и снизить погрешность при его аппроксимации за счет учета случайных выбросов и провалов нагрузки.
При плавном изменении процесса на определенном интервале времени вейвлет растягивается по оси t или х (координата) и на нем проводится анализ, и, наоборот, при резком изменении (выброс или провал) - вейвлет сжимается до интервала выброса или провала.
В данном случае модель тока РИ представляется через коэффициенты вейвлет-преобразования и базис-вейвлетов:
1 f
где а = —, fn - частота гармонических составляющих, включая и первую гармонику; х - аналог коор-
Jn
динат времени t, характеризующий изменение тока во времени; Wxa(t) - вейвлет произвольной формы; Wj(x, а) - вейвлет-преобразование - функция двух аргументов - х и а.
При фиксированном значении а распределение Wj (х, а) является сверткой исходного процесса изменения тока i(t).
В качестве примера показан вейвлет-анализ непериодического случайного процесса изменения тока тяговой нагрузки при помощи вейвлета «Сомбреро» (Mexican hat), который является второй производной гауссиана и имеет вид:
X2
W(x) = (x2 - 1)е~. (12)
Как было сказано, переменная х - время.
Для каждой пары х и а структура вычисления вейвлет-преобразования Wj(x,a) следующая: вейвлет W(x) растягивается (сжимается) в а раз по горизонтали и в 1/а по вертикали; полученная функция сдвигается в точку анализа случайного процесса на временной оси, xt - данный вейвлет обозначается Wxa; усредняются значения исследуемого процесса в окрестности точки а (величина, обратная частоте) при помощи вейвлета V'х а; определяется вейвлет- преобразование токовой нагрузки с помощью произвольной функции вейвлета W(x):
W,(x,a)=1f-;™y(i-)I(t)dt. (13)
Распределение значений Wj(x, а) строится в осях х, а: по горизонтали - переменная х (время или координата в пространстве), по вертикали - а. Далее строятся графики вейвлета W(x) при разных значениях х и а, осуществляется их обработка с исследуемым процессом. Методика такого вейвлет-анализа подробно приведена в [18].
Величина любого коэффициента Wj (x, а) показывает, насколько характерный период колебаний а, соответствующий анализируемым частотам, представлен в точке в окрестности исследуемого момента времени х.
Таким образом, вейвлет-преобразование Wj (x, а) содержит информацию о частотных и временных свойствах реализаций случайного процесса одновременно, что позволяет изучать случайные процессы более детально, чем использование Фурье-преобразований. Данный анализ можно реализовать в системе MATLAB.
В общем виде вейвлет-преобразование состоит из двух основных этапов: разложения исследуемого процесса, т.е. нахождение коэффициентов разложения, их анализ, и восстановление процесса по набору коэффициентов вейвлет-разложения.
Схема построения модели тока в пакете прикладных программ (ППП) MATLAB следующая.
1. Выбирается наиболее эффективный тип вейвлета ^o(t) для рассматриваемого процесса, -например, вейвлет Добеши, гауссовский вейвлет, вейвлет "Mexican hat", Морле и др.
При этом вейвлет W(x) Е L2(R) должен удовлетворять условию:
Cv = 2л f+™lul-1i?(u)l2du < (14)
Если W(x) 6 L1(R), то Ч>(0) = ¥(x)dx = 0.
2. Выбираются значение масштабного коэффициента а и параметр смещения (позиции) х.
3. Определяются вейвлет-коэффициенты тока нагрузки, соответствующие выбранным значениям а и х, по формуле (13) по алгоритму [18].
Токовая нагрузка распределенного источника задается выборкой значений {1^}. Объем выборки определяется интервалом дискретизации, найденным с помощью теоремы Котельникова, и периодом наблюдения. В качестве тока берется кусочно-постоянная функция
I( t) = I(kAt) = Ik (15)
при t 6 [kAt, (k + 1) At].
Значения параметра х определяются шагом At: x = 1:11, II = length(I). (16)
Так как аргумент t имеет дискретные значения, определенные шагом A t, то используется вейвлет
V(t) = V0(t/At) (17)
с носителем на промежутке [0, NAt]. Здесь W0(t/At) - выбранный в пункте 1 вейвлет с дискретным представлением аргумента.
Носитель вейвлета W (~~) находится на промежутке [x, x + aNAt], коэффициенты рассчитываются при x = nAt по формуле
CXA = ^Jktn Ik !(Т)М * (лт) dt = ^aA tilT^1 IkCk(a, x), (18)
где для x = nAt коэффициенты Ck (a, x) имеют вид:
Ck(a,x) = v(^)-v(^) (19)
Здесь W(t) - первообразная выбранного вейвлета. Если вейвлет Wg(t) принимается из MATLAB, то его первообразная находится при помощи команды «intwave».
4. Используя формулу обращения, можно получить модель токовой нагрузки:
К0 = Ci-1 ИК2 W,(x, a)W JO-1^, (20)
где Cv = 2л J|^|-1|l//(w)|2 dw вычисляется по вейвлету W(t).
При нечетком задании исходной информации выбор подхода к моделированию токов зависит от вида неопределенности. Например, при ее определении в фидерной зоне количество электровозов можно задавать нечетким множеством. При исследовании нестационарных случайных процессов изменения тока тяговой нагрузки подстанции математическое ожидание каждой реализации тока n-й гармоники, включая первую, можно определять с помощью множественной регрессии с нечеткими коэффициентами. Эти коэффициенты находятся с помощью а-сечений или другим известным из теории нечетких множеств методом [19]:
Icpn(t)=?MZ-ihn. (21)
Если неопределенность обусловлена неполнотой информации, то моделирование тока осуществляется совместным применением теории случайных процессов и нечетких множеств. Сначала осуществляется восстановление информации об изменении тока. Его порядок следующий.
1. Выбираются множества значений трех реализаций случайного процесса: двух ретроспективных (полные потоки информации) и текущей (поток с неполной исходной информацией). Текущее множество нужно восстановить. Ретроспективные множества обозначаются [Р1, /Р2; текущее или восстанавливаемое - /в. Эти множества задаются относительными значениями.
2. Для каждого выбранного множества строится колоколообразная (bell - shapped) функция принадлежности.
Ш = —rrtfB, (22)
I а I
где а, b, с - некоторые числовые параметры, принимающие произвольные действительные значения и упорядоченные отношения: a<b<c, причем b>0, или функция плотности нормального распределения
в предположении, что \ 2по = 1, задаваемая аналитическим выражением
-(1-е)2
Ш = е 2«2 , (23)
где с - математическое ожидание реализации случайного процесса (рассматриваемого множества); а - ее среднее квадратическое отклонение.
Построение графиков функций принадлежности осуществляется в среде fuzzy - TECH MATLAB [20].
3. Рассчитывается относительное расстояние Хемминга по выражениям
eG^i") =1^=1 \ Ц^НО^ )\
S(lP2,Lb) = 1:in=i \ ШП - ш )\ (24)
5(1р2,г)=1ж=1 ¡тп—тт
где 0 < 3([р1,1ь) <1;0< 5(1р2,[ъ) <1;0< 5(1р2,[р1) < 1.
4. В зависимости от значений относительного расстояния Хемминга восстанавливаемое текущее множество заполняется значениями либо первого, либо второго ретроспективного множества. Сравнение обобщенных относительных расстояний Хемминга 8(1р1,1ь), 8(1р2,1ь), 5(1р2,[р1) позволяет определить степень близости одного нечеткого множества к другому. Если 5([р1,1Ь^<5([р2, /й), то 1Ь^1р1.
Нахождение расстояния между 1 значениями тока нечетких множеств [ь и 1р1 для данного случая позволяет выбрать наиболее достоверное восстанавливаемое значение:
а(1в,1Р1) = таХ[Гп(Л)Гп(!!1)] — тт^Ц)^1)]. (25)
Если а(1в,1Р1>)< 0,5, то выбирается значение /(в, в противном случае /Р1.
Полученная совокупность значений тока моделируется рассмотренными выше способами.
Область применения модели тока распределенного источника искажения качества электроэнергии
Учитывая разнообразие распределенных источников искажения качества электроэнергии, стохастический характер их работы, а также высокую степень неопределенности функционирования энергосистем, следует не только осуществлять мониторинг качества электроэнергии, но и адаптивное управление им в электрических сетях.
Под адаптивным управлением качеством электроэнергии понимается такое управление, при котором, используя результаты мониторинга КЭ в адаптивном режиме, автоматически осуществляется управляющее воздействие, направленное на улучшение КЭ, причем доза управляющего воздействия выбирается таким образом, чтобы обеспечить максимальный эффект от ее ввода. Для этого необходимо иметь модели тока РИ, адекватные результатам мониторинга КЭ, т.е. качеству информации о токовой нагрузке РИ.
Предложенные модели тока (обобщенная и частные случаи) позволяют разработать и реализовать алгоритмы адаптивного управления качеством электроэнергии. При этом в алгоритме адаптивного управления КЭ должен присутствовать блок выбора соответствующей модели в зависимости от качества информации по анализируемому РИ, а результаты обработки токов являются основой для выбора узла для ввода управляющего воздействия и дозы управляющего воздействия.
Заключение
1. Показано, что при переходе электроэнергетики на новую технологическую платформу, основанную на цифровизации и интеллектуализации электрических сетей, в энергосистемах увеличится не только доля распределенных источников искажения качества электроэнергии, но и типы таких источников. В то же время для их описания можно использовать единый методический подход. В качестве примера он показан для тяги переменного тока.
2. В зависимости от видов неопределенности выделены четыре группы информации, используемой для моделирования токов распределенных источников качества электроэнергии. Дан их анализ.
3. Получена обобщенная модель тока распределенных источников искажения качества электроэнергии и ее частые случаи при различных группах исходной информации.
4. Показано, как можно моделировать токи гармоник, генерируемых распределенными источниками искажения качества электроэнергии, на основе совокупного применения теории случайных процессов и вейвлет-анализа, теорий случайных процессов и нечетких множеств с использованием ППП МЛТЬЛБ.
1. Энергетическая стратегия РФ на период до 2035 года, утвержденная распоряжением Правительства Российской Федерации от 09.06.2020 г. № 1523- р [Электронный ресурс]: URL: https://www.garant.ru/products/ipo /prime/doc/74148810/ (дата обращения: 15.03.2022 г).
2. Smart Grids - основы и технологии энергосистем будущего / пер. с англ.; под общ. ред. Н.И. Воропая. -М.: Изд. дом МЭИ, 2017. - 461 с.: ил. - ISBN 978-5-383-01228-4.
3. Инновационная электроэнергетика-21 / под ред. В.М. Батенина, В.В. Бушуева, Н.И. Воропая. - М.: ИЦ «Энергия», 2017. - 584 с. ISBN (print): 978-5-98908-457-9.
4. Елистратов, В.В. Возобновляемая энергетика. - Изд. 2-е, доп. - СПб.: Наука, 2013. - 307 с.: ил. - ISBN 978-5-02-038220-6.
5. Велькин, В.И. Методология расчета комплексных систем ВИЭ для использования на автономных объектах - Екатеринбург: Уральский федеральный университет им. Б.Н. Ельцина, 2015. - 226 с. - ISBN 978-5-32102460-7. - EDN WNXMSV.
6. Булатов, Ю.Н. Распределенная генерация и энергетические роутеры в системах электроснабжения железных дорог / Ю.Н. Булатов, А.В. Крюков, Г.О. Арсентьев. - Москва-Берлин: ООО «Директмедиа Паблишинг», 2020. - 172 с. - ISBN 978-5-4499-0083-8. - EDN PHWLQH.
7. Электромагнитная совместимость потребителей / И.В. Жежеленко, А.К. Шидловский, Г.Г. Пивняк [и др.] - М.: Машиностроение, 2012. - 351 с.
8. Высшие гармоники в системах электроснабжения промпредприятий / И. В. Жежеленко. - Изд. 6-е, пере-раб. и доп. - М.: Энергоатомиздат, 2010. - 375 с. - ISBN 978-5-283-033-07-5.
9. Бадер, М.П. Электромагнитная совместимость: учеб. для вузов - М.: УМК МПС, 2002. - 638 с. - ISBN 5-89035-065Х.
10. Управление качеством электрической энергии. Сборник трудов Международной научно-практ. конф. Москва, 23-25 ноября 2016 г. - М.: ООО «Центр полиграфических услуг «РАДУГА», 2017. - 300 с. - ISBN 9785-905486-14-2. - EDN ZSHUCH.
11. Управление качеством электроэнергии / И. И. Карташев [и др.], под ред. Ю.В. Шарова. - М.: Изд. дом МЭИ, 2006. - 319 с. - ISBN 5-903072-13-5. - EDN QMJCZZ.
12. Смирнов, С.С. Высшие гармоники в сетях высокого напряжения - Новосибирск: Академический научно-
изд. производственно-полиграф. и книгораспространительский центр «Наука», 2010. - 326 с. - ISBN 978-5-02023320-1. - EDN QMLHNP.
13. Улучшение качества электроэнергии в системах электроснабжения нетяговых потребителей железных дорог / В.П. Закарюкин, А.В. Крюков, И.А. Любченко, А.В. Черепанов. - Москва-Берлин : Директ-Медиа, 2020.
- 184 с. - ISBN 978-5-4499-1580-1. - EDN KIUCMZ.
14. Баскаков, С.И. Радиотехнические цепи и сигналы. Учеб. для студентов вузов по специальности «Радиотехника» - Изд. 4-е, перераб. и доп. - М.: Высш. шк., 2003 (ГУП ИПК Ульян. Дом печати). - 462 с. - ISBN 5-06003843-2 (в пер.).
15. Нефедов, В.И. Основы радиоэлектроники и связи. Учеб. пособие для студентов высших учебных заведений, обучающихся по специальности «Проектирование и технология радиоэлектронных средств», направления подготовки «Проектирование и технология электронных средств» / В.И. Нефедов, А.С. Сигов; под ред. В.И. Нефедова. - М.: Высш. шк., 2009. - 735 с.: ил. - ISBN 978-5-06-006161-1.
16. Савина, Н.В. Системный анализ потерь электроэнергии в электрических распределительных сетях / отв. ред. Н.И. Воропай; Российская акад. наук, Сибирское отд., Ин-т систем энергетики им. Л.А. Мелентьева, М-во образования и науки Российской Федерации, Амурский гос. ун-т. - Новосибирск: Наука, 2008. - ISBN 978-5-02023222-8.
17. Savina, N.V. Modelling of parameters of the electric power quality characterizing non-sinusoidal and asymmetry of voltage // 6th International conference Electrical power quality and utilization. September 19-21, 2001. - Cracow, Poland. - P.173-178.
18. Смоленцев, Н.К. Введение в теорию вейвлетов. - М.; Ижевск: РХД, 2019. - 292 с. - ISBN 978-5-43440748-9.
19. Борисов, В.В., Федулов, А.С., Зернов, М.М. Основы теории нечетких множеств. Учебное пособие для вузов. Серия «Основы нечеткой математики». - М.: Горячая линия - Телеком, 2014. - ISBN: 978-5-9912037-1-5.
20. Леоненков, А.В. Нечеткое моделирование в среде MATLAB и fuzzyTECH - СПб.: БХВ-Петербург, 2003.
- 719 с.: ил. - ISBN 5941570872.
