Научная статья на тему 'Моделирование элементов идентификации и контроля параметров системы для интеллектуальной поддержки принятия решений'

Моделирование элементов идентификации и контроля параметров системы для интеллектуальной поддержки принятия решений Текст научной статьи по специальности «Электротехника, электронная техника, информационные технологии»

CC BY
149
36
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
МОЩНОСТЬ / ЭЛЕКТРОЭНЕРГЕТИЧЕСКИЕ СИСТЕМЫ / ФИЛЬТРО-КОМПЕНСИРУЮЩИЕ УСТРОЙСТВА / НАГРУЗКА / ИНТЕЛЛЕКТУАЛЬНЫЕ СИСТЕМЫ / ИДЕНТИФИКАЦИЯ / ОПТИМАЛЬНОЕ УПРАВЛЕНИЕ / POWER / POWER ENERGY SYSTEMS / FILTER-COMPENSATING DEVICES / LOAD / INTELLECTUAL SYSTEMS / IDENTIFICATION / OPTIMAL CONTROL

Аннотация научной статьи по электротехнике, электронной технике, информационным технологиям, автор научной работы — Бордюг Александр Сергеевич, Титов Иван Леонидович, Железняк Александр Александрович

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

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по электротехнике, электронной технике, информационным технологиям , автор научной работы — Бордюг Александр Сергеевич, Титов Иван Леонидович, Железняк Александр Александрович

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

MODELING ELEMENTS FOR IDENTIFICATION AND CONTROL OF THE SYSTEM PARAMETERS FOR INTELLECTUAL DECISION-MAKING SUPPORT

The problem has been reviewed when calculating the nominal power capacity of the filter-compensating device the non-active power consumed by the load from the power source should be taken into account; in the process of operation the compensating device may maintain the required parameter of the electric energy quality regulating the various parameters of quality within the nominal range. There has been considered a way to indentify the circuit current and voltage harmonics in the real time in order to apply in the system of the intellectual decision-making support for optimal governing of the filter-compensating device.

Текст научной работы на тему «Моделирование элементов идентификации и контроля параметров системы для интеллектуальной поддержки принятия решений»

УДК 62-551:51-74

А. С. Бордюг, И. Л. Титов, А. А. Железняк

МОДЕЛИРОВАНИЕ ЭЛЕМЕНТОВ ИДЕНТИФИКАЦИИ И КОНТРОЛЯ ПАРАМЕТРОВ СИСТЕМЫ ДЛЯ ИНТЕЛЛЕКТУАЛЬНОЙ ПОДДЕРЖКИ ПРИНЯТИЯ РЕШЕНИЙ

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

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

Введение

Низкая стоимость и доступность природного газа на буровых платформах позволяют рассмотреть возможность его применения в первичных двигателях электроэнергетических систем (ЭЭС). В первую очередь это связано с наличием нефтяного газа, который добывается вместе с нефтью и может использоваться в качестве топлива, при этом можно сократить долю токсичных выхлопов на 80 % - по сравнению с дизельным топливом, которое сейчас используется в ЭЭС. Из-за отсутствия окиси серы и сажи в продуктах сгорания на элементах двигателя не образуется нагар, масло не загрязняется, в результате чего можно увеличить срок службы приводного двигателя.

Анализируя литературные данные, можно сделать вывод, что применение комбинированных двигателей в ЭЭС буровых платформ является перспективным направлением, в первую очередь для использования попутного газа, что принесет существенную экономию. Тем не менее, существует ряд проблем в использовании данных типов двигателей в многогенераторных ЭЭС. Газодизельные двигатели очень чувствительны к нагрузке из-за их системы автоматического регулирования подачи топлива, в результате чего невозможно обеспечить идентичные динамические характеристики для многогенераторных ЭЭС. Для точной регулировки подачи газа необходимо отсутствие люфтов и малого количества передаточных звеньев. В связи с созданием новой электронной элементной базы средств автоматики, включая микропроцессорные системы, появилась возможность ввести в эксплуатацию комплексную интеллектуальную систему автоматического регулирования генераторным агрегатом, которая также выполняет функцию регулирования частоты вращения приводного двигателя. Таким образом, единая комплексная система повышенной эксплуатационной надежности выполняет все действия по управлению и регулированию контролируемых значений генераторного агрегата и автоматизированной электроэнергетической системы [1-4].

Если произвести анализ выражения коэффициента мощности (1), который учитывает показатели мощностей с помощью произведения относительной мощности искажения % и показателя реактивной мощности по первой гармонике cos ф

X = х ^ф, (1)

можно представить обобщенный показатель качества электроэнергии как

U1I1 cos ф

X = хcosф = -

rms rms

где U1 - напряжение сети; и„ - среднеквадратичное значение напряжения; 11 - ток нагрузки; - среднеквадратичное значение тока нагрузки. Используя зависимости

Ц __1

Urms yj1 + СГИ2

u

— = , 1

ф + СГИ ,2 '

где СГИ - суммарные гармонические искажения, получим

= 1__1

+СГИ 2иф+СГИ2

Предполагая, что существует алгоритм управления фильтрокомпенсирующим устройством (ФКУ), при котором достигаются заданные параметры обобщенного показателя качества при наименьшей (или оптимальной) установленной мощности ФКУ, произведем формализацию задачи оптимального управления [3, 5-9].

Материал и результаты исследований

Введем вектор гармоник тока преобразователя в сети:

I = р15 ..., i,..., iN]T, (2)

где ii — i-я гармоника тока; N — количество учитываемых гармоник.

Дополнительно введем вектор гармоник напряжения и вектор неактивной мощности ФКУ:

U = [Ui, ..., ut, ..., Un]t; (3)

X = [Xi, ..., х,, ..., XN]t. (4)

Допустим, что полная установочная мощность ФКУ задается параметром Sopt, а текущая мощность может быть выражена Ss.

Внешними параметрами задачи управления являются векторы I и U, а S0pt — прямое ограничение. Исходя из этого, обобщенный показатель качества X и мощность S2 являются выходными параметрами задачи оптимального управления. Учитывая добавленные векторы H — верхних — и L — нижних границ вектора Х, для i-го элемента векторов справедливо

L < х. < h.. (5)

ill V '

Примем, что даны значения зависимости

X = f (X ); (6)

N.

x

i=1

f (X ) = £ Xi. (7)

Выполним максимизацию X и оптимизацию при принятых ограничениях:

Xщш < Л (X )< Xтах; (8)

Ь < X < Н; (9)

Л (X ) = Зр, (10)

где Xmln, Xmax - минимальное и максимальное значения обобщенного показателя качества.

Создание ограничений для функции (8) вводится для лимитирования выходных параметров задачи управления, а выражение (10) является функциональным лимитом на внутренние управляемые параметры. С физической точки зрения ограничения (10) указывают на то, что при выполнении распределения мощности ФКУ, необходимой для компенсации уровней гармоник и реактивной мощности, текущая суммарная мощность З2 ФКУ равняется установленной мощности Зор1 ФКУ. С математической стороны задачи оптимального управления ФКУ выполнятся системами (2)-(10). Оптимальная точка Х, оптимальное значение целевой функции Х и оптимальное значение целевой функции З^ = /1(Х ) образуют оптимальное решение задачи, которое

может быть локальным и глобальным. Локальное решение представляет собой наименьшее значение целевой функции в ограниченной окрестности точки Х, в то время как глобальное решение дает наименьшее значение целевой функции [10-14].

Чтобы создать оптимальное решение задачи, введем точку Х (оптимальная) - оптимальное значение целевой функции Х , а также оптимальное значение целевой функции ^ = ^(Х). Это решение может быть как локальным, так и глобальным. Локальным решением будет наименьшее значение целевой функции в ограниченной окрестности точки Х, а глобальное решение даст наименьшее значение целевой функции.

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

1) идентифицируются внешние параметры, записываются значения всех постоянных параметров и ограничений (блоки 1, 5, 6);

2) с учетом прямых ограничений вычисляется вектор Х (блок 2). Здесь же производится вычисление значения целевой функции Sy и ее оценка на оптимальность: если Sy Ф Sopt, то производится корректировка вычисления внутренних управляемых параметров;

3) находится значение обобщенного показателя качества Х, сравнивается со значениями функциональных ограничений Хтт, Хтах (блок 3); если не соблюдаются все условия функционального ограничения на выходные параметры, необходимо произвести корректировку вычислений внутренних управляемых параметров;

4) выполняется вывод координат выходного параметра Х (блок 4).

Рис. 1. Алгоритм решения задачи оптимального управления ФКУ

В данном алгоритме граничные значения (9) задаются для единиц вектора Х, также задается ограничение функций на выходные параметры задачи управления для целевой функции Х. Следует отметить, что в документации (требования ГОСТ и правила Регистра) ограничиваются значения суммарных гармонических искажений (СГИ) напряжения, уровни отдельных гармоник напряжения и cos ф. Имеется в виду, что все граничные значения даны для элементов век-

торов I, U и их функций (СГИ, cos ф). Для работы системы управления ФКУ по алгоритму оптимального управления необходимо наличие информации о гармоническом составе токов и напряжений сети, которая учитывается в векторах I и U. Сигналами на выходе считаются мгновенные значения фазных токов и линейных напряжений сети. Таким образом, при анализе внешних параметров требуется решать задачи идентификации векторов I и U по измеренным токам и напряжениям сети (рис. 2) [4, 14-17].

Рис. 2. Блоки внешних параметров системы оптимального управления ФКУ

Информацию о гармоническом составе сетевых токов и напряжений можно получить от векторов внешних параметров I и U. Чем точнее информация от этих параметров, тем качественнее выходные параметры системы ФКУ, т. к. от частоты, амплитуды и фазы зависит качество подавления той или иной гармонической составляющей при помощи ФКУ. В числе наиболее распространенных способов определения гармонического состава сигнала отметим дискретное преобразование Фурье, фильтр Калмана и искусственные нейронные сети. Основными недостатками данных систем являются характерные особенности применения автономных ЭЭС, например, нестабильная частота напряжения сети. Таким образом, отсутствует возможность использования метода мгновенной мощности для определения реактивной составляющей тока в устройствах компенсации реактивной мощности, следовательно, необходим поиск новых подходов к решению задачи гармонического анализа контролируемого сигнала в реальном времени.

В общем виде входной сигнал y(t) системы может быть аппроксимирован функцией-рядом Фурье:

N

y(t) = £ Zk (t) sin (Ык (t )t + Фк (t));

к=1

N

y(t) = £ [ Ak (t) Sin (кю к (t )t) + Bk (t) cos (кю к (t )t)];

к=1

Z(t) = VA2(t) + B2(t); Ф(0 = ctg^,

где Zt(t) - амплитуда и фк(0 - фазовый угол к-й гармоники; N - количество гармоник.

Дискретное во времени представление функции y(t) будет представлено как

N

Ж) = £ [ Ак (tn) Sin (кЮк (tn)) + Bk (tn) cos (кЮк (tn))], (11)

к=1

где Ак(tn) = Zk(tn)cos(фк(tn)); Bk(tn) = 2к(tn)sin(фк(tn)).

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

Xn = [sin<(tn), cos<(tn), sin2<(tn), cos2< (tn), ..., sin N< (tn), cos N< (tn )J (12) и матрицу переменных состояния

Cn =[ A(ín) Bi(tn) ... An (tn) BN (tn)] , (13)

тогда

y = X C . (14)

s n n n V '

Погрешность аппроксимации характеризуется выражением

£n УП УП- (15)

Энергетическая функция ошибки, выражаемая среднеквадратичной погрешностью

1 2

En = 2en - (16)

Функция (16) является функцией m переменных E(cb c2, ..., cm), дифференцируемой по всем аргументам. Запишем метод поиска значения C вектора С, которому соответствует экстремальное (минимальное) значение функции Е:

E (C) > E (C); VE е Zm, (17)

где Zm - область возможных значений вектора C (допустимая область).

В связи с тем, что функция (16) является дифференцируемой, а (14) выражает только линейное преобразование элементов входного вектора (12), для поиска значений матрицы (13) применим градиентный алгоритм. Задачей аппроксимации является такое изменение элементов вектора С, при котором выполняется минимизация функции (16).

Мы подошли к задаче поиска экстремума энергетической функции (16). Воспользуемся методом градиента, при котором вектор (13) изменяют так, что скорость его изменения оказывается пропорциональной градиенту по правилу

dc

d- = ц-grad E (C), (18)

dt ci

где grad E (C) = dE (C) - частная производная функции E по одному из ее аргументов, при

c dci

i = 1, ..., N ; ц - постоянный коэффициент скорости аппроксимации, который при поиске минимума должен удовлетворять условию ц > 0.

Обычно задают вектор C, с помощью которого вычисляют для данной точки градиент энергетической функции (16), затем изменяют каждый компонент вектора С со скоростью, пропорциональной составляющей градиента (18) по этому компоненту. Измерение градиента происходит постоянно в процессе изменения С, таким образом, в каждый момент времени градиент этого момента пропорционален скорости изменения вектора С.

В экстремальной точке C условием экстремума является

grad E (C)

dE dE dE

dc1 dc2 dcN

= 0. (19)

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

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

Сп+1 = Сп + ^^ Е (Сп) (20)

С

или

= дЕ

С, (п+1) = С, (п) + Ц дс '

По правилу (20) вычисления будут продолжаться до тех пор, пока не будет выполнено условие Е(Сп+1) > Е(Сп). Если при расчете имеется дрейф частоты ю основной гармоники (следовательно, происходит пропорциональное изменение частоты всех высших гармоник), произойдет сдвиг экстремума и увеличится ошибка аппроксиматора, которую невозможно компенсировать при постоянной частоте функций элементов входного вектора (12). Чтобы создать наивысшую точность аппроксимации, требуется изменять и частоту в функциях элементов входного вектора (12) таким образом, чтобы она была равна значению угловой частоты ю. Имеется в виду, что дополнительно необходимо следить за частотой аппроксимируемого сигнала.

В случае, когда известно значение выражения (15), при выполнении условия выражения (11) частная производная по частоте ю от среднеквадратичной ошибки (16) даст

N

дю = - ЦюЕп Е кю(*п ) - ЫпВк кю(*п )]> (21)

к=1

где Цга - коэффициент скорости аппроксимации по частоте.

Полученное значение (21) необходимо суммировать с аргументами функций элементов входного вектора (12):

ю^п+1) = ю(*п ) + дю.

Можно сделать вывод, что в функция (16) в области Zm в общем случае может иметь несколько экстремумов. Один является глобальным, а остальные будут локальными (рис. 3). Условие (19) выполнится во всех этих точках, однако неравенство (17) при глобальном экстремуме выполнится во всей области Zm, а при локальном - только в определенной малой области, которая включает точку С . Таким образом, глобальный минимум считается наименьшим из локальных минимумов.

Рис. 3. Наличие глобального и локальных экстремумов многомерной функции

в области ее определения

После завершения итераций алгоритма происходит приведение к локальному экстремуму, который определен к начальному значению С. Но в нашей ситуации экстремумы определяются частотами гармоники и их фазовыми углами (имеется в виду, что существуют ограничения области возможных значений для частоты ю). Это значение и допустимые отклонения содержатся в нормативной документации (требованиях ГОСТа или Регистра), мы точно знаем заданные координаты глобального экстремума С . Гарантировано нахождение функции в области именно глобального экстремума. Помимо этого принимается коэффициент скорости аппроксимации по частоте в несколько раз меньше коэффициента скорости по изменению элементов вектора С (^ю = 0,01^), в результате чего создается сходимость к экстремуму (рис. 4).

Рис. 4. Локализация области определения многомерной функции в окрестности глобального экстремума

Рис. 5 отображает структурную схему аппроксиматора, который выполняет слежение за частотой ю и компонентами Ак, Вк гармоник.

ю

хх(п) = этю/ (п)

Х2(П) = соэю/ (п) -

Х2Л-1(п) = этЛю/ (п) -

Х2Л/(п) = соЛю/ (п)

Поиск

Ак Ь к

У Г у Г

Поиск ю

е(п) ▼

У(Ю

Нп)

Рис. 5. Структурная схема аппроксиматора, выполняющего слежение за частотой ю и компонентами Ак, Вк гармоник сигнала у

Рассмотрим входные физические сигналы системы управления ФКУ. В автономных ЭЭС измеряются линейные напряжения и фазные токи сети. Чтобы произвести оценку не только гармонического состава, но и уровня возможных асимметрий (дисбаланса) и других нежелательных характеристик, самым правильным является определение или восстановление напряжений сихронного генератора (СГ) относительно недоступного нулевого провода, т. е. фазных напряжений СГ [16-17]. Таким образом, схема аппроксиматора (рис. 5) примет вид (рис. 6):

Рис. 6. Структурная схема аппроксиматора фазных напряжений сети

Контроль за фазными напряжениями СГ, которые невозможны для физического измерения, можно осуществить, измеряя линейные напряжения сети.

Осуществление гармонического анализа токов выполняется с помощью схемы (рис. 5) с тем отличием, что частота ю определяется один раз для одного из фазных напряжений, и с его помощью анализируются все остальные компоненты токов и напряжений сети.

Результат моделирования процесса аппроксимации синусоидального сигнала приведен на рис. 7.

Ziöl

u I li—

N

-£(ZlJ) 0,002 2 1=

Рис. 7. Результат обучения фильтра аппроксимации синусоидального сигнала

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

Заключение

Определен обобщенный показатель качества электроэнергии в автономной ЭЭС, вводимый как критерий оценки качества функционирования ФКУ.

Произведена формализация задачи оптимального управления ФКУ (обобщенным показателем качества электроэнергии) в автономной ЭЭС.

Решена задача идентификации и контроля внешних параметров системы оптимального управления ФКУ.

СПИСОК ЛИТЕРА ТУРЫ

1. Cooly J. W., Tukey J. W. An Algorithm for Machine Calculation of Complex Fourier Series // Mathematics of Computation. 1965. Vol. 19. P. 297-301.

2. Grady W. M., Heydt G. T. Prediction of Power System Harmonics Due to Gaseous Discharge Lighting // IEEE Transactions on Power Systems. 1985. Vol. 104. P. 554-561.

3. Moreno S. V., Barros G. J. Application of Kalman Filtering for Continuous Real-time Tracking of Power System Harmonics // IEEE Proc. -Gener. Transm. Distrib. 1997. Vol. 144, no. 1. P. 13-20.

4. Черный С. Г., Жиленков А. А. Интеллектуальная поддержка принятия решений при оптимальном управлении для судовых электроэнергетических систем // Вестн. гос. ун-та мор. и реч. флота им. адм. С. О. Макарова. 2014. № 3 (25). С. 68-75.

5. Harris J. On the use of window for harmonic analysis with DFT // Proceedings of the IEEE. 1978. Vol. 66. P. 51-83.

6. Girgis A. A., Ham F. A Quantitative Study of Pitfalls in FFT // IEEE Transactions on Aerospace and Electronic Systems. 1980. Vol. 16, no. 4. P. 434-439.

7. Sharaf M., Dash P. K. A Kalman Filtering Approach for Estimating Power System Harmonics // Proceeding of the Third Intern. Conf. on Power system Harmonics. 1989. P. 34-40.

8. Girgis A. A., Chang W. B., Makram E. B. A Digital Recursive Measurement Sheme for On-Line Tracking of Power System Harmonics // IEEE Transactions on Power Delivery. 1991. Vol. 6, no. 3. P. 1153-1160.

9. Haili M., Girgis A. A. Identification and Tracking of Harmonic Sources in a Power System Using a Kalman Filter // IEEE Transactions on Power Delivery. 1996. Vol. 11, no. 3. P. 1659-1665.

10. Pileggi D. J., Chandra N. H, Emanuel A. E. Prediction of Harmonic Voltages in Distribution Systems // IEEE Transactions on Power Apparatus and Systems. 2002. Vol. PAS-100, no. 3. P. 1307-1315.

11. Soliman S. A., Christensen G. S., Kelly D. H., El-Naggar К. M. A State Estimation Algorithm for Identification and Measurements of Power System Harmonics // Electric Power System Research Journal. 1990. Vol. 19. P. 195-206.

12. Hartana R. K., Richards G. G. Harmonic Source Monitoring and Identification Using Neural Networks // IEEE Transactions on Power Systems. 1990. Vol. 5, no. 4. P. 1098-1104.

13. Mori H. An Artificial Neural Network Based Method for Power System Voltage Harmonics // IEEE Transactions on Power Delivery. 1992. Vol. 7, no. 1. P. 402-409.

14. Osowski S. Neural Network for Estimation of Harmonic Components in a Power System // IEEE Proceeding С. 1992. Vol. 139, no. 2. P. 129-135.

15. Черный С. Г., Жиленков А. А. Идентификация внешних параметров сигналов для экспертных подсистем в составе устройств судовых электроэнергетических систем // Науч.-техн. ведом. Санкт-Петерб. гос. политехи. ун-та. Информатика. Телекоммуникации. Управление. 2014. № 3 (198). С. 28-36.

16. Zhilenkov A., Chernyi S. Investigation Performance of Marine Equipment with Specialized Information Technology // Procedia Engineering. 2015. Vol. 100. P. 1247-1252.

17. Chernyi S., Zhilenkov A. Modeling of Complex Structures for the Ship's Power Complex Using XILINX System // Transport and Telecommunication. 2015. Vol. 16 (1). P. 73-82.

Статья поступила в редакцию 16.05.2016

ИНФОРМАЦИЯ ОБ АВТОРАХ

Бордюг Александр Сергеевич — Россия, 298309, Керчь; Керченский государственный морской технологический университет; ассистент кафедры электрооборудования судов и автоматизации производства; Laboratory45@yandex.ru.

Титов Иван Леонидович - Россия, 298309, Керчь; Керченский государственный морской технологический университет; ассистент кафедры электрооборудования судов и автоматизации производства; vanya-titov-1993@mail.ru.

Железняк Александр Александрович - Россия, 298309, Керчь; Керченский государственный морской технологический университет; ассистент кафедры электрооборудования судов и автоматизации производства; Zheleznyak13@mail.ru.

A. S. Bordyug, I. L. Titov, A. A. Zheleznyak

MODELING ELEMENTS FOR IDENTIFICATION AND CONTROL OF THE SYSTEM PARAMETERS FOR INTELLECTUAL DECISION-MAKING SUPPORT

Abstract. The problem has been reviewed when calculating the nominal power capacity of the filter-compensating device the non-active power consumed by the load from the power source should be taken into account; in the process of operation the compensating device may maintain the required parameter of the electric energy quality regulating the various parameters of quality within the nominal range. There has been considered a way to indentify the circuit current and voltage harmonics in the real time in order to apply in the system of the intellectual decision-making support for optimal governing of the filter-compensating device.

Key words: power, power energy systems, filter-compensating devices, load, intellectual systems, identification, optimal control.

REFERENCES

1. Cooly J. W., Tukey J. W. An Algorithm for Machine Calculation of Complex Fourier Series. Mathematics of Computation, 1965, vol. 19, pp. 297-301.

2. Grady W. M., Heydt G. T. Prediction of power system harmonics due to Gaseous Discharge Lighting. IEEE Transactions on Power Systems, 1985, vol. 104, pp. 554-561.

3. Moreno S. V., Barros G. J. Application of Kalman Filtering for Continuous Real-time Tracking of Power System Harmonics. IEEE Proc. -Gener. Transm. Distrib., 1997, vol. 144, no. 1, pp. 13-20.

4. Chernyi S. G., Zhilenkov A. A. Intellektual'naia podderzhka priniatiia reshenii pri optimal'nom uprav-lenii dlia sudovykh elektroenergeticheskikh sistem [Intellectual decision-making support under optimal control of marine electric power systems]. Vestnik gosudarstvennogo universiteta morskogo i rechnogo flota im. admirala S. O. Makarova, 2014, no. 3 (25), pp. 68-75.

5. Harris J. On the use of window for harmonic analysis with DFT. Proceedings ofthe IEEE, 1978, vol. 66, pp. 51-83.

6. Girgis A. A., Ham F. A Quantitative Study of Pitfalls in FFT. IEEE Transactions on Aerospace and Electronic Systems, 1980, vol. AES-16, no. 4, pp. 434-439.

7. Sharaf M., Dash P. K. A Kalman Filtering Approach for estimating power System Harmonics. Proceeding of the Third International Conference on Power system Harmonic, 1989. P. 34-40.

8. Girgis A. A., Chang W. B., Makram E. B. A Digital Recursive Measurement Sheme for On-Line Tracking of Power System Harmonics. IEEE Transactions on Power Delivery, 1991, vol. 6, no. 3, pp. 1153-1160.

9. Haili M., Girgis A. A. Identification and Tracking of Harmonic Sources in a Power System Using a Kalman Filter. IEEE Transactions on Power Delivery, 1996, vol. 11, no. 3, pp. 1659-1665.

10. Pileggi D. J., Chandra N. H, Emanuel A. E. Prediction of Harmonic Voltages in Distribution Systems. IEEE Transactions on Power Apparatus and Systems, 2002, vol. PAS-100, no. 3, pp. 1307-1315.

11. Soliman S. A., Christensen G. S., Kelly D. H., El-Naggar K. M. A State Estimation Algorithm for Identification and Measurements of Power System Harmonics. Electric Power System Research Journal, 1990, vol. 19, pp. 195-206.

12. Hartana R. K., Richards G. G. Harmonic Source Monitoring and Identification Using Neural Networks. IEEE Transactions on Power Systems, 1990, vol. 5, no. 4, pp. 1098-1104.

13. Mori H. An Artificial Neural Network Based Method for Power System Voltage Harmonics. IEEE Transactions on Power Delivery, 1992, vol. 7, no. 1, pp. 402-409.

14. Osowski S. Neural Network for Estimation of Harmonic Components in a Power System. IEEE Proceeding S, 1992, vol. 139, no. 2, pp. 129-135.

15. Chernyi S. G., Zhilenkov A. A. Identifikatsiia vneshnikh parametrov signalov dlia ekspertnykh podsis-tem v sostave ustroistv sudovykh elektroenergeticheskikh sistem [Identification of external signal parameters for expert subsystems in the structure of marine electric power systems]. Nauchno-tekhnicheskie vedomosti Sankt-Peterburgskogo gosudarstvennogo politekhnicheskogo universiteta. Informatika. Telekommunikatsii. Upravlenie, 2014, no. 3 (198), pp. 28-36.

16. Zhilenkov A., Chernyi S. Investigation Performance of Marine Equipment with Specialized Information Technology. Procedia Engineering, 2015, vol. 100, pp. 1247-1252.

17. Chernyi S., Zhilenkov A. Modeling of Complex Structures for the Ship's Power Complex Using XIL-INX System. Transport and Telecommunication, 2015, vol. 16 (1), pp. 73-82.

The article submitted to the editors 16.05.2016

INFORMATION ABOUT THE AUTHORS

Bordyug Alexander Sergeevich - Russia, 298309, Kerch; Kerch State Marine Technological University; Assistant of the Department of Marine Electrical Equipment and Industrial Automation; Laboratory45@yandex.ru.

Titov Ivan Leonidovich - Russia, 298309, Kerch; Kerch State Marine Technological University; Assistant of the Department of Marine Electrical Equipment and Industrial Automation; vanya-titov-1993@mail.ru.

Zheleznyak Alexander Aleksandrovich - Russia, 298309, Kerch; Kerch State Marine Technological University; Assistant of the Department of Marine Electrical Equipment and Industrial Automation; Zheleznyak13@mail.ru.

i Надоели баннеры? Вы всегда можете отключить рекламу.