Киселёв_Статьяx

advertisement
УДК 621.316
Применение быстрого вейвлет преобразования для решения переходных процессов в
электрических цепях
Б.Ю. Киселёв, Д.С. Осипов
Омский государственный технический университет, г. Омск, Россия
Аннотация – В статье рассматривается применение вейвлет-анализа в рекурсивных методах расчета
переходных процессов в цепях переменного тока. Вейвлет-преобразование, имеющее ряд преимуществ перед
преобразованием Фурье для анализа нестационарных процессов в электрических цепях, являются удобным
инструментом для компьютерного моделирования. Предлагаемый алгоритм может быть реализован для расчета
коротких замыканий в электроэнергетических системах.
Ключевые слова: Вейвлет-преобразования, рекурсивный метод, переходные процессы
Быстрое вейвлет преобразование.
В настоящее время вейвлет анализ получил большое распространение во всем мире,
методы вейвлет-преобразования находят широкое применение для решения актуальных
электроэнергетических задач, как за рубежом [5], так и в России [1,3].
Часто при анализе нестационарных сигналов удобно их представление в виде суммы
двух последовательных приближений, аппроксимирующей Am(t) и детализирующей Dm(t)
составляющих [6]:
m
S (t )  Am (t )   D j (t )
(1)
j 1
С последующим их уточнением итерационными методами.
Разложение сигнала на эти самые последовательности позволяет осуществить
Быстрое вейвлет-преобразование (БВП). Оно осуществляет построенный на фильтрации
итерационный алгоритм, число итераций может быть произвольным. Каждая итерация будет
соответствовать все более глубокому уровню разложения [2].
Сигнал, размерностью N можно разложить на 2N / 2 уровней, которые представляют
собой набор аппроксимирующих (Am) и детализирующих (Dm) коэффициентов. Для
нахождения этих коэффициентов в БВП используются специальные коэффициенты фильтра
hl. Эти коэффициенты однозначно определяют функции материнского  (t ) и отцовского
 (t ) вейвлета [6].
Процесс разложения сигнала происходит следующим образом. Исходная функция
(сигнал) v(t) размерностью N интерпретируется как последовательность аппроксимирующих
коэффициентов Am. Далее при помощи коэффициентов фильтра h0 и h1 получаем
последовательность
аппроксимирующих
и
последовательность
детализирующих
коэффициентов на уровень ниже Am-1 и Dm-1. Причем размерность этих сигналов будет равна
N/2. этот процесс соответствует децимации, прореживанию сигнала в два раза. Алгоритм
можно продолжать до тех пор, пока не получим 1 аппроксимирующий коэффициент с00 и
N-1 детализирующих коэффициентов d jk .[3]
Весь процесс разложения сигнала происходит по алгоритму Малла, схема которого
представлена на рисунке 1.
Рис. 1. Схема алгоритма Малла для разложения сигнала
Процесс синтеза или реконструкции сигнала происходит по обратному алгоритму
Малла. Схема алгоритма представлена на рис. 2.
Рис. 2. Схема алгоритма Малла для реконструкции сигнала
Эта процедура использует операции интерполяции и фильтрации фильтрами
реконструкции h2, h3.
Интерполяция – это операция повышения частоты дискретизации в два раза путём
добавления нулевых компонентов вперемежку имеющимся. Далее происходит сложение
пропущенных через фильтры сигналов с jk и d jk и мы получаем на выходе
реконструированный сигнал или набор коэффициентов аппроксимации на уровень выше
с j 1,k . Данный процесс можно продолжать до тех пор пока на выходе мы не получим
исходный сигнал.[6]
Применение вейвлет преобразований для расчета переходных процессов
Ток в катушке индуктивности сразу же после замыкания ключа (выключателя),
должен остаться таким же, каким он был непосредственно до замыкания ключа
(выключателя).
i L (0  )  i L (0  )
(2)
Параметры RL –цепи: R=6.35 Ом, L=0,127 Гн и u (t )  2  220 sin( 2f   ) (источник
переменного тока f=50 Гц.   30 ).
Используя законы Кирхгофа, запишем уравнения напряжения, взяв направление
обхода контура по часовой стрелке Рисунок 3. Мы имеем обыкновенное дифференциальное
уравнение.
u  Ri  L
di
dt
(3)
Рис. 3. Анализируемая схема RL – цепи
Решив уравнение (4) получим:
i L (t )  7.68  sin( 2f  51)  5.97  e
t / 0.02
(4)
В работе [7] предлагается решение обыкновенных дифференциальных уравнений с
помощью вейвлетов – метод вейвлет рекурсии.
Пусть функции u(t ) и iL (t ) представляются векторами, обозначим время как
t  nt , где t -интервал (время) дискретизации.
Дискретное вейвлет преобразование этих векторов будет иметь вид:
(5)
u ( t )    u j , k   j ,k ( t )
j k
(6)
i ( t )    i j ,k   j ,k ( t )
j k
Уравнение обратной конечно-разностной аппроксимации производной в дискретном
времени [7]:
di L i L (t )  i L (t  t )

(7)
dt
t
Уравнение (3) можно переписать, учитывая уравнение (8).
1
t
i L (t )  i L (t  t )  R i L (t )  1 u (t )
L
(8)
L
Подставляя в скалярную версию уравнения (7) дискретные вейвлет серии для u (t ) и
i L (t ) , полученные с помощью уравнений (5),(6) для производного приближения(при малых
значениях ) t получим:
1
k   i  (t )  k   i  (t t )  

j k Ljk jk
 j k Ljk jk
t 

R

1

1
 

(9)
k   i  (t ) 
k   u  (t )
 j k Ljk jk  L  j k jk jk 
L
Где
 jk (t ) 
 jk (t  t ) 
1
a0
a0
t

 b0 k 
a

 t  t
 a
 
(10)


 b0 k 
(11)
БВП имеет быструю рекурсивную форму, что наиболее легко реализовать с
использованием теории группы фильтров [8].
Используя теорию КМА описываемой Малла в [9] Вейвлет разложение функции (в
данном случаи функции тока на катушке индуктивности и напряжения источника) будут
иметь вид:

 
A
D
(12)
i L ( t )   i L  j ,k ( t )    i L  j ,k ( t )
j
,
k
j  J k  j ,k
k 

 
A
D
 u j ,k  j ,k ( t )    v j ,k  j ,k ( t )
(13)
j  J k 
k 
В уравнениях (6,9-6.10) индексы D и A представляют собой детализирующие и
аппроксимирующие коэффициенты.
Ниже будут представлены численные результаты эксперимента. Простое
дифференциальное уравнение (3) будет решено при помощи предложенного в [7] метода
формулы правой обратной разности (конечно – разностные аппроксимации производных).
Этапы процедуры анализа переходного процесса при помощи метода вейвлет
рекурсии:
Шаг 1. Производим БДП используя коэффициентов фильтра, выбранного типа
вейвлета. Тем самым находим коэффициенты вейвлет преобразования. (уравнения 5-6)
Шаг 2. Решаем систему уравнений, составленную в вейвлет области.
1
1
R A
1 A
A
A
(14)
j i Ljk  j  i Ljk 1  i Ljk  u jk
L
L
a 0 t
a 0 t
1 D
1
R D
1 D
D
(15)
i


i

i

u jk
Ljk
Ljk 1
Ljk
j
j
L
L
a 0 t
a 0 t
u (t ) 
В этом численном
4
t  10 (с)  100 мс.
эксперименте
j=1,
a0j  21  2
,
временной
шаг
A
D
Выразив из уравнений i Ljk и i Ljk получим:
A
i Ljk
 u Ajk

L




  5  10 3  i A
Ljk 1


D
i Ljk
 u Djk

L




  5  10 3  i D
Ljk 1


3
 L
(16)
5  10  R
 L
3
5  10  R
(17)
Шаг 3. Выполняем обратное вейвлет преобразование (реконструкцию сигнала),
A
D
последовательности коэффициентов i
и i Ljk преобразуются в единый дискретный
Ljk
сигнал.
Рисунок 4. иллюстрирует результаты рекурсивного решения системы уравнений
(результаты после 600 итераций).
15
i(t),A
10
5
0
-5
-10
0
0.01
0.02
0.03
0.04
0.05 0.06
time(sec)
0.07
0.08
0.09
0.1
Рис. 4. Результаты рекурсии после 600 итераций
Рассмотренный метод, называемый методом вейвлет рекурсии, использует
относительно новый математический аппарат, вейвлет-преобразование, которое имеет
определенные преимущества по сравнению с методами, такими, как анализ Фурье при
анализе переходных и нестационарных сигналов. Преимущества и недостатков метода
вейвлет рекурсии: БВП производит разрежение представленных сигналов, таким образом, к
данному методу могут быть применены редкие матричные приёмы, позволяющие увеличить
скорость решения. Во многих случаях, метод вейвлет рекурсия сходится относительно
быстро Программное обеспечение вейвлет-анализа легко доступны и адаптированы к
различным приложениям, в том числе к энергетике.
Библиографический список
1.
Аббакумов, А. А. Разработка методики и алгоритмов идентификации
отклонений от нормативов параметров качества электроэнергии в системах
электроснабжения : дис. … канд. техн. наук : 05.13.18 / А. А. Аббакумов / Мордовский гос.
ун-т им. Н.П. Огарева – Саранск, 2005. – 180 с.
2.
Дьяконов, В. П. Вейвлеты. От теории к практике / В. П. Дьяконов.– М. : СолонР, 2002. – 448 с.
3.
Карпенко, С. В. Математическое моделирование нестационарных
электрических процессов в электротехнических системах на основе численных методов
вейвлет-анализа : дис. … канд. техн. наук : 05.13.18 / С. В. Карпенко. – Новокузнецк, 2006. –
164 с.
4.
Малла, С. Вейвлеты в обработке сигналов : пер. с англ. / С. Малла. – М. : Мир,
2005. – 671 с.
5.
Мисриханов, А. М. Применение методов вейвлет-преобразования в
электроэнергетике / А. М. Мисриханов // Автоматика и телемеханика. – 2006. – № 5. – С. 5–
23.
6.
Яковлев, А. Н. Введение в вейвлет-преобразования : учеб. пособие / А. Н.
Яковлев. – Новосибирск : НГТУ, 2003. – 104 с.
7.
Calli, A. W. "Analysis of Electrical Transients In Power System Via A Novel
Wavelet Recursion Method" / A. W. Calli. – Ph.D. Dissertation, Purdue University, 1997. – 143 p.
8.
Strang, G. Wavelets and Filter Banks /
G. Strang, T. Nguyen / WelleslyCambridge Press. – Wellesly. MA, 1996. – 527 p.
9.
Mallat, S. A wavelet tour of signal processing / S. Mallat / Academic Press. – New
York, 1998. – 805 p.
Download