close

Вход

Забыли?

вход по аккаунту

?

Патент BY16511

код для вставкиСкачать
ОПИСАНИЕ
ИЗОБРЕТЕНИЯ
К ПАТЕНТУ
РЕСПУБЛИКА БЕЛАРУСЬ
(46) 2012.10.30
(12)
(51) МПК
НАЦИОНАЛЬНЫЙ ЦЕНТР
ИНТЕЛЛЕКТУАЛЬНОЙ
СОБСТВЕННОСТИ
(54)
BY (11) 16511
(13) C1
(19)
G 01R 23/16 (2006.01)
СПОСОБ СПЕКТРАЛЬНОГО АНАЛИЗА
ЭЛЕКТРИЧЕСКОГО СИГНАЛА
(21) Номер заявки: a 20100064
(22) 2010.01.20
(43) 2011.08.30
(71) Заявитель: Желудкевич Михаил
Ефимович (BY)
(72) Автор: Желудкевич Михаил Ефимович (BY)
(73) Патентообладатель: Желудкевич Михаил Ефимович (BY)
(56) SU 1663571 A1, 1991.
RU 2363005 C1, 2009.
RU 2001103694 A, 2003.
SU 1120252 A, 1984.
CN 101261292 A, 2008.
DE 4134472 A1, 1993.
CN 101368870 A, 2009.
US 6408696 B1, 2002.
(57)
Способ спектрального анализа электрического сигнала, при котором посредством
фильтрации исследуемого сигнала выделяют косинусные и синусные гармонические составляющие ряда Фурье по отсчетам i = 1, N + 1 , где N - количество отсчетов, соответствующим середине интервалов разбиения диапазона изменения исследуемого сигнала
180°
Y(n), где n = 1, N , при одинаковом шаге между отсчетами, равном
, причем значения
N
Y(n), соответствующие середине n-го интервала, принимают равными полусумме значений на его концах, запоминают упомянутые значения, центрируют их и по центрирован-
BY 16511 C1 2012.10.30
°
ным значениям отсчетов сигнала Y(n ) , где n = 1, N , определяют два ряда Фурье: первый
ряд косинусных составляющих сигнала в базисе ортогональных тригонометрических дискретных многочленов Чебышева первого рода и второй ряд синусных составляющих сигнала в базисе ортогональных тригонометрических дискретных функций Чебышева
второго рода, определяют соответствующие амплитуды ak, bk и bN этих составляющих из
выражений:
Фиг. 1
Фиг. 2
BY 16511 C1 2012.10.30
2 N °
π⋅k
(n − 0,5) ,
a k = ∑ Y(n ) cos
N
N n =1
где k = 1, N − 1,
°
1 N
Y(n ) = Y(n ) − ∑ Y(n ),
N n =1
N °
2
π⋅k
(n − 0,5),
b k = ∑ Y(n )sin
N n =1
N
1 N °
n +1
b N = ∑ Y(n )(− 1) ,
N n =1
запоминают ak, bk и bN, по которым определяют N гармонических компонентов
 π⋅k
(n − 0,5) + ψ k , где k = 1, N , спектра исследуемого сигнала, параметры Ak,Ψk и
A k sin 
 N

AN, ΨN которых определяют из выражений:
1 2
Ak =
a k + b 2k ,
2
где k = 1, N − 1,
a
ψ k = arctg k ,
bk
1
AN = bN ,
2
ΨN = 0.
Изобретение относится к измерительной технике и может использоваться для выделения гармонических составляющих электрического сигнала в информационно-измерительных устройствах.
Известен способ спектрального анализа дискретного сигнала (SU 1702323 A1, G 01R
23/16, 1991), в котором для определения гармоник спектра исследуемого сигнала используют последовательность низкодобротных цифровых фильтров. Частотная характеристика
первого фильтра перекрывает весь анализируемый частотный диапазон сигнала, эквивалентная частотная характеристика всей последовательности является самой узкой и определяет разрешающую способность анализатора. Эквивалентные частотные характеристики
соседних последовательностей отличаются одна от другой на величину разрешающей
способности. На первом этапе фильтрации на выходе последнего фильтра выделяют одну
гармонику, соответствующую разрешающей способности, и определяют ее параметры. На
втором этапе на выходе предпоследнего фильтра выделяют уже две гармоники, одну известную, параметры которой определены на первом этапе, а параметры второй определяют.
Аналогично реализуется алгоритм фильтрации и на последующих этапах. Недостатком
указанного способа спектрального анализа является необходимость знать частотный диапазон исследуемого сигнала и значение разрешающей способности.
Известен способ графоаналитического определения спектра сигнала (Бессонов Л.А.
"Теоретические основы электротехники". М.: Высшая школа, 1978. - с. 162-165), основанный на использовании разложения исследуемого сигнала в ряд Фурье с использованием
"правила трапеций", при переходе от интеграла Фурье к суммированию, что обеспечивает
более точное определение параметров косинусно-синусных составляющих по сравнению
с традиционно используемым "правилом прямоугольника".
2
BY 16511 C1 2012.10.30
Наиболее близким к предлагаемому способу определения спектра сигнала является
способ определения спектра сигнала, реализованный в устройстве по А.с. SU 1663571 G
01R 13/16, 1991. По этому способу записанный в устройстве массив отсчетов дискретного
сигнала фильтруют последовательностью нерекурсивных фильтров, реализующих скользящее усреднение. Подбором параметров цифровых фильтров их порядка уменьшают
массив отсчетов, подлежащих ДПФ, до величины, обеспечивающей устойчивое выделения заданного числа гармоник и определение их амплитуд.
Недостатки аналога-прототипа заключаются в следующем.
Применение цифровых фильтров на основе скользящего усреднения эффективно для
стационарных сигналов. Для определения числа усредняемых отсчетов, количества циклов фильтрации (порядка фильтра), числа каскадов фильтров необходимы априорные
данные о частотных свойствах исследуемого сигнала.
N
Применение ДПФ позволяет выделить только
гармоник исследуемого сигнала при
2
N −1
- при нечетном N.
четком числе отчетов N или
2
Задачей предлагаемого решения с использованием чебышевских разложений является
увеличение количества выделяемых гармоник спектра по сравнению с прототипом, где
используется ДПФ.
Поставленная задача достигается тем, что в способе спектрального анализа электрического сигнала, при котором посредством фильтрации исследуемого сигнала выделяют косинусные и синусные гармонические составляющие ряда Фурье по отсчетам i = 1, N + 1 ,
где N - количество отсчетов, соответствующим середине интервалов разбиения диапазона
изменения исследуемого сигнала Y(n), где n = 1, N , при одинаковом шаге между отсчета180 o
, причем значения Y(n), соответствующие середине n-го интервала, приN
нимают равными полусумме значений на его концах, запоминают упомянутые значения,
ми, равном
°
центрируют их и по центрированным значениям отсчетов сигнала Y(n ) , где n = 1, N ,
определяют два ряда Фурье: первый ряд косинусных составляющих сигнала в базисе ортогональных тригонометрических дискретных многочленов Чебышева первого рода и
второй ряд синусных составляющих сигнала в базисе ортогональных тригонометрических
дискретных функций Чебышева второго рода, определяют соответствующие амплитуды
ak, bk и bN этих составляющих из выражений:
2 N °
π⋅k
(n − 0,5) ,
a k = ∑ Y(n ) cos
N
N n =1
где k = 1, N – 1, k = 1, N – 1
°
1 N
Y(n ) = Y(n ) − ∑ Y(n ),
N i =1
N °
π⋅k
2
(n − 0,5),
b k = ∑ Y(n )sin
N
N n =1
1 N °
π⋅k
(n − 0,5),
b N = ∑ Y(n )sin
N
N n =1
запоминают ak, bk и bN, по которым определяют N гармонических компонентов
 π⋅k
(n − 0,5) + ψ k , где k = 1, N , спектра исследуемого сигнала, параметры Ak,Ψk и
A k sin 
 N

AN, ΨN которых определяют из выражений:
3
BY 16511 C1 2012.10.30
Ak =
1 2
a k + b 2k ,
2
где k = 1, N − 1,
ψ k = arctg
ak
,
bk
1
bN ,
2
ΨN = 0.
Заявленный способ спектрального анализа электрического дискретного сигнала с использованием чебышевских разложений имеет существенные преимущества, так как позволяет увеличить число выделяемых гармоник спектра по сравнению с ДПФ, повышает
точность оценок параметров гармоник с использованием усреднения по формуле:
1
Y(n ) = [Y(i − 1) + Y(i)] ,
2
что соответствует использованию "правила трапеций" более точного, чем "правила прямоугольников", и при этом снижается уровень априорных сведений о частотных свойствах используемого сигнала. Использование упрощенной формулы сглаживания
повышает быстродействие в определении спектральных составляющих и обеспечивает
универсальность предлагаемого способа.
Для подтверждения технического результата предлагаемого решения был проведен
вычислительный эксперимент, в ходе которого реальный производственный электрический сигнал с 13 исходными цифровыми отсчетами был подвергнут ДПФ (фиг. 1) и чебышевским разложениям (фиг. 2).
На фиг. 1 представлены амплитуды 6 гармоник спектра исследуемого сигнала, которые были получены с использованием ДПФ, на фиг. 2 представлены амплитуды 12 гармоник исследуемого сигнала с использованием предлагаемого решения, однако значение
амплитуд четных гармоник в два раза меньше, чем значения, полученные с использованием ДПФ. Таким образом, при увеличении числа выделяемых гармоник в два раза их амплитуда соответственно уменьшается.
AN =
Национальный центр интеллектуальной собственности.
220034, г. Минск, ул. Козлова, 20.
4
Документ
Категория
Без категории
Просмотров
0
Размер файла
160 Кб
Теги
by16511, патент
1/--страниц
Пожаловаться на содержимое документа