способ спектрального анализа периодических многочастотных сигналов, представленных цифровыми отсчетами

Формула изобретения

Способ спектрального анализа многочастотных периодических сигналов, представленных цифровыми отсчетами, заключающийся в том, что для исходного a(t_i) и опорного b₀ (t_i) сигналов, представленных отсчетами мгновенных значений одних и тех же моментов времени t_i=t ₁, t₂, , t_N, находят точки совместного решения а(b ₀) при различных круговых частотах опорного сигнала _0k и фазах опорного сигнала _0k определяют площадь их вольт-амперной характеристики Q_BAXk, отличающийся тем, что исходный a(t_i ) и опорный b₀(t_i) сигналы задают отсчетами мгновенных значений при



где t - шаг дискретизации;

N - количество точек отсчета за время выборки сигнала t_N, для исходного сигнала

a(t_i) определяют производную



одновременно для частот _{01 0m 0n} (Рад/с) с заданным шагом ₀ формируют опорный сигнал



при этом принимают _0m=-180..( _0m+ _0m)..180, определяют площадь вольтамперной характеристики для точек совместного решения



выявляют максимальное значение площади вольт-амперной характеристики при различных фазах max[Q_BAXm( _0m)], выделяют частоты _0k, при которых

Q_BAX( _0m) Q_BAXmax·k,

где Q_BAXmax =max[Q_BAXm( _0m)] - максимальное значение площади вольт-амперных характеристик при различных частотах _0m, k - коэффициент фильтрации, %, для выделенных частот _0k определяют опорный сигнал b₀(t _i) для рассматриваемой частоты _0k



при этом принимают _0k=-180..( _0k+ ₀)..180, после чего определяют площадь вольт-амперной характеристики для точек совместного решения a(b₀)



и определяют амплитуду



где - частота, Гц,

выявляют максимальное значение амплитуды при различных фазах max[Am_k( _0k)], значения А_mk, _k= _0k, _k= _0k, соответствующие максимальному значению амплитуды при различных фазах max[Am_k( _0k)], запоминают, о спектральном составе сигнала судят по значениям Am_k, _k, _k и их либо передают на ЭВМ, либо выводят на дисплей.

Описание изобретения к патенту

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

Известен способ спектрального анализа периодических многочастотных сигналов, представленных цифровыми отсчетами [Патент РФ 2229725, МПК⁷ G01R 23/16, опубл. 12.11.2002], выбранный в качестве прототипа, заключающийся в формировании опорного синусоидального сигнала, который многократно сдвигают по фазе один сигнал относительно другого. Для анализируемого периодического многочастотного сигнала a(t_i) и опорного b₀(t_j)=В_msin( ₀t_j+ ₀) сигнала с амплитудой В_m, представленных отсчетами мгновенных значений для одних и тех же моментов времени t_j=t₁,t₂,..,t_N, где N - число разбиений на периоде Т, находят точки совместного решения а(b₀) при различных круговых частотах опорного сигнала ₀ и фазах опорного сигнала ₀. Строят вольтамперные характеристики и определяют их площадь причем вывод о присутствии гармонической составляющей с круговой частотой _k и фазой _k в анализируемом сигнале a(t_i) делают, исходя из условия F_BAX=0. Затем определяют максимальную площадь вольтамперной характеристики при _k=±90° и находят амплитуду К-й спектральной составляющей по формуле

где

₁ - основная круговая частота.

Далее по значениям _k, _k и судят о спектральном составе анализируемого сигнала a(t _i).

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

Задачей изобретения является расширение арсенала средств аналогичного назначения.

Это достигается тем, что в способе спектрального анализа многочастотных периодических сигналов, представленных цифровыми отсчетами, также как в прототипе, для исходного a(t_i ) и опорного b₀(t_i) сигнала, представленных отсчетами мгновенных значений одних и тех же моментов времени t_i=t₁,t₂,..,t_N, находят точки совместного решения а(b₀) при различных круговых частотах опорного сигнала и фазах опорного сигнала определяют их площадь

Согласно изобретению исходный a(t_i ) и опорный b₀(t_i) сигнал задают отсчетами мгновенных значений, при

где t - шаг дискретизации;

N - количество точек отсчета за время выборки сигнала t_N.

Для исходного сигнала a(t_i) определяют производную

Одновременно для частот (рад/с) с заданным шагом

₀ формируют опорный сигнал

при этом принимают Определяют площадь вольтамперной характеристики для точек совместного решения

Выявляют максимальное значение площади вольтамперной характеристики при различных фазах выделяют частоты при которых

где - максимальное значение площади вольтамперных характеристикам при различных частотах k - коэффициент фильтрации в процентах. Для выделенных частот определяют опорный сигнал b₀(t_i) для рассматриваемой частоты

при этом принимают После чего определяют площадь вольтамперной характеристики для точек совместного решения а(b₀)

и определяют амплитуду

где - частота в Гц.

Выявляют максимальное значение амплитуды при различных фазах значения соответствующие максимальному значению амплитуды при различных фазах запоминают. О спектральном составе сигнала судят по значениям Am_k, _k, _k и их либо передают на ЭВМ, либо выводят на дисплей.

Заявленный способ спектрального анализа многочастотных периодических сигналов, представленных цифровыми отсчетами, имеет существенные преимущества, поскольку повышает чувствительность выявления частотных составляющих с малыми амплитудами за счет того, что выявляют частоты по спектру производной исходного сигнала, так как известно, что чувствительность - способность выделить и определить параметры спектральной составляющей с малой амплитудой [Функциональный контроль и диагностика электротехнических систем и устройств по цифровым отсчетам мгновенных значений тока и напряжения/ под редакцией Е.И Гольдштейна - Томск: Изд. «Печатная мануфактура», 2003, с.89].

На фиг.1 приведена аппаратная схема устройства, реализующего рассматриваемый способ определения периода.

На фиг.2 приведена осциллограмма тестового исходного сигнала.

На фиг.3 приведена осциллограмма производной тестового исходного сигнала.

На фиг.4 приведена зависимость значения площади вольтамперной характеристики от частоты.

На фиг.5 приведена зависимость амплитуды от частоты.

В табл.1 приведены параметры тестового исходного сигнала.

В табл.2 приведены выявленные частоты и соответствующие значения площади вольтамперной характеристики.

В табл.3 приведены результаты определения параметров тестового исходного сигнала.

Заявленный способ может быть осуществлен с помощью схемы (фиг.1), содержащей датчик анализируемого сигнала 1 (ДАС), который подключен к объекту, сигнал которого исследуется. К датчику анализируемого сигнала 1 (ДАС) последовательно подключены дифференциатор 2 (Д), программатор определения площади вольтамперной характеристики 3 (ПВАХ), первый датчик опорного сигнала 4 (ДОС1). При этом к программатору определения площади вольтамперной характеристики 3 (ПВАХ) последовательно подключены первый программатор 5 (П1), второй датчик опорного сигнала 6 (ДОС2), программатор амплитудного значения по вольтамперной характеристике 7 (АВАХ), второй программатор 8 (П2), который связан с дисплеем или ЭВМ (не показано на фиг.1.). Датчик опорного сигнала 1 (ДАС) соединен с программатором амплитудного значения по вольтамперной характеристике 7 (АВАХ), а первый программатор 5 (П1) соединен со вторым программатором 8 (П2).

В качестве датчика анализируемого сигнала 1 (ДАС) может быть использован датчик тока - промышленный прибор КЭИ-0,1, или датчик напряжения - трансформатор напряжения (220/5В). Программатор площади вольтамперной характеристики 3 (ПВАХ), программатор амплитудного значения 7 (АВАХ), дифференциатор 2 (Д), первый и второй программаторы 5 (П2) и 8 (П2), первый и второй датчики опорного сигнала 4 (ДОС1) и 6 (ДОС2) могут быть выполнены на микроконтроллере серии 51 производителя atmel AT89S53. Для работы пользователя может быть предусмотрена кнопочная клавиатура FT008, имеющая 8 кнопок.

С выхода датчика анализируемого сигнала 1 (ДАС) исходный сигнал a(t,) заданный отсчетами мгновенных значений в моменты времени

t_l,t₂,..,t_i ,.., t_N;

t₂-t₁ =t₃-t₂=t_N-t_N-1= = t;

где t - шаг дискретизации;

N - количество точек отсчета за время выборки сигнала t_N,

поступает на вход программатора определения амплитудного значения по вольтамперной характеристике 2 (АВАХ) и дифференциатора 2 (Д). В дифференциаторе 2 (Д) для исходного сигнала а(t_i ) определяют производную

и передают на программатор площади вольтамперной характеристики 3 (ПВАХ), одновременно с первого датчика анализируемого сигнала 4 (ДОС1) для частот (рад/с) с заданным шагом ₀ формируют опорный сигнал

при этом принимают длина массива N и шаг дискретизации t принимают равными соответствующим параметрам исходного сигнала а(t_i), затем полученный опорный сигнал b ₀(t_i) подают на программатор площади вольтамперной характеристики 3 (ПВАХ), где определяют площадь вольтамперной характеристики для точек совместного решения

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

где - максимальное значение площади вольтамперных характеристикам при различных частотах k - коэффициент фильтрации в процентах. Выделенные частоты передают на второй датчик опорного сигнала 6 (ДОС2) и на второй программатор 8 (П2). Во втором датчике опорного сигнала 6 (ДОС2) для частот определяют опорный сигнал b₀(t_i) для рассматриваемой частоты

при этом принимают длину массива N и шаг дискретизации t принимают равными соответствующим параметрам исходного сигнала a(t_i). Затем определяют площадь вольтамперной характеристики для точек совместного решения а(b₀)

и определяют амплитуду

где - частота в Гц.

В программаторе определения амплитудного значения по вольтамперной характеристике 2 (АВАХ) выявляют максимальное значение амплитуды при различных фазах max[Am_k( _0k)], максимальное значение амплитуды Аm _k и соответствующее значение фазы _k= _0k передают на второй программатор 8 (П2), где их запоминают вместе со значением _k= _0k из первого программатора 5 (П1). После определения амплитуд Am_k и фаз _k для всех значений частоты _0k на выход второго программатора 8 (П2) подают сигнал о готовности результатов спектрального анализа. О спектральном составе сигнала судят по значениям Am_k, _k, _k и их либо передают на ЭВМ, либо выводят на дисплей.

Для проверки работоспособности предложенной процедуры провели спектральный анализ тестового сигнала напряжения (параметры сигнала приведены в табл.1).

a(t _i)=u(t_i)=10·sin(2 ·50·t_i+45)+0,1·sin(2 ·400·t_i+30)+0,05·sin(2 ·800·t_i-20), A.

Для этого рассчитали массив мгновенных значений для N=10000 точек исходного сигнала с шагом дискретизации t=10^-4c(t_N=l c) (осциллограмма тестового исходного сигнала приведена на фиг.2). В дифференциаторе 2 (Д) для исходного сигнала a(t_i) определяют производную (см. фиг.3). Для частот (1,2..900)·2 (рад/с) формируют опорный сигнал и определяют максимальные площади вольтамперной характеристики (см. фиг.4). В программаторе 5 (П1) выделяют частоты например

для _m=400·2 (рад/с) (Q_BAX(400·2 )=49709,83)

Q_BAX(400·2 ) 96138,8561·40%=38455,5424

где

Выделенные частоты и значения Q_BAX представлены в табл.2. Во втором датчике опорного сигнала 6 (ДОС2) для частот определяют опорный сигнал b₀(t_i) для рассматриваемой частоты при этом принимают

Затем определяют площадь вольтамперной характеристики и определяют максимальную амплитуду при различных фазах, например для

Результаты спектрального анализа тестового сигнала представлены в табл.3 и фиг.5. Убедиться в увеличении чувствительности можно по полученным значениям характеристик Q_BAX( _m) и Аm( _m) (см. фиг.4, 5). Таким образом, на фиг.4 видно, что значения Q_BAX( _m) имеют большее значение по сравнению с аналогичными значениями характеристики Аm( _m), причем этот эффект проявляется все больше с ростом задаваемой частоты _m.