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

Изобретение относится к области систем обработки информации и измерительной технике. Для исходного и опорного сигнала, представленных отсчетами мгновенных значений одних и тех же моментов времени, находят точки совместного решения при различных круговых частотах опорного сигнала. При этом опорный и исходный сигнал задают отсчетами мгновенных значений. Для исходного сигнала определяют производную. Далее для точек совместного решения производной исходного и опорного сигналов определяют площадь вольтамперной характеристики. Выявляют максимальное значение площади вольтамперной характеристики при различных фазах. Также для точек совместного решения исходного и опорного сигналов определяют площадь вольтамперной характеристики. Выявляют максимальное значение амплитуды при различных фазах. О спектральном составе сигнала судят по значениям амплитуд, круговой частоты опорного сигнала и фазы. Эти значения либо передают на ЭВМ, либо выводят на дисплей. Технический результат сводится к расширению арсенала средств аналогичного назначения. 5 ил., 3 табл.

 

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

Известен способ спектрального анализа периодических многочастотных сигналов, представленных цифровыми отсчетами [Патент РФ 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.

Способ спектрального анализа многочастотных периодических сигналов, представленных цифровыми отсчетами, заключающийся в том, что для исходного 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) определяют производную

одновременно для частот ω₀₁…ω_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 и их либо передают на ЭВМ, либо выводят на дисплей.