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

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

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

t ₁,t₂,...,t_j ,...,t_N;

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

,

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

N - количество точек за время Т,

для последовательности частот ₁, ₂,..., _j,..., _n определяют мгновенную спектральную плотность по выражениям

где S₁( _j) и S₂( _j) - синусная и косинусная составляющая мгновенной спектральной плотности,

затем определяют экстремумы характеристики распределения мгновенной спектральной плотности по частоте S( ), по точкам экстремумов определяют частоты содержащиеся в исходном сигнале a(t_i), для найденных частот определяют амплитуду

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

затем последовательно формируют опорный сигнал

b₀(t_j)=B _msin( ₀t_j+ ₀),

где В_m - амплитуда опорного сигнала,

при частоте опорного сигнала ₀= _k, строят вольтамперную характеристику а(b₀) для исходного сигнала a(t _i), многократно сдвигают по фазе, определяя площадь вольтамперной характеристики F_BAX( ₀=-180°...+180°),

где ₀ - фаза опорного сигнала,

полученные значения F_BAX( ₀=var) запоминают, выбирают значение фазы частотной составляющей _k= ₀, при котором F_BAX 0, значения амплитуды и фазы _k, соответствующие частоте _k, запоминают, данную операцию повторяют для всех найденных значений частоты _k, далее по полученным значениям _k, _k судят о спектральном составе исходного сигнала a(t_i).

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

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

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

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

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

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

Известен способ СА многочастотных периодических сигналов, представленных цифровыми отсчетами [Функциональный контроль и диагностика электротехнических систем и устройств по цифровым отсчетам мгновенных значений тока и напряжения / под редакцией Е.И.Гольдштейна - Томск: Печатная мануфактура, 2003, с.92-94], выбранный в качестве прототипа, заключающийся в том, что для исходного сигнала a(t _i), заданного отсчетами мгновенных значений в моменты времени

t₁, t₂,..., t _j,..., t_N;

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

,

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

N - количество точек за время Т,

для последовательности частот ₁, ₂,..., _j,..., _n определяют мгновенную спектральную плотность (МСП) по выражению

;

;

;

где S₁( _j) и S₂( _j) - синусная и косинусная составляющие мгновенной спектральной плотности.

Частоте _k=f_k·2 , рад (f_k - частота, Гц), содержащейся в исходном сигнале a(t_i), соответствует экстремум функции S( _j). Таким образом определяют последовательность частот ₁, ₂,..., _j,..., _n, содержащихся в сигнале. Для этих частот определяют амплитуду и фазу _k частотной составляющей _k находят по формулам:

.

По полученным значениям _j, _j судят о спектральном составе исходного сигнала a(t_i).

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

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

Это достигается тем, что в способе спектрального анализа многочастотных периодических сигналов, представленных цифровыми отсчетами, так же, как в прототипе, для исходного сигнала a(t_i), заданного отсчетами мгновенных значений в моменты времени

t ₁, t₂,..., t_j ,..., t_N;

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

где t - шаг дискретизации; N - количество точек за время T,

для последовательности частот ₁, ₂,..., _j,..., _n определяют мгновенную спектральную плотность по выражениям

где S₁( _j) и S₂( _j) - синусная и косинусная составляющая мгновенной спектральной плотности. Затем определяют экстремумы характеристики распределения мгновенной спектральной плотности по частоте S( ), по точкам экстремумов определяют частоты _k, содержащиеся в исходном сигнале a(t_i), для найденных частот определяют амплитуду

Согласно изобретению для выявления экстремумов характеристики мгновенной спектральной плотности используют условие

затем последовательно формируют опорный сигнал

b₀(t_i)=B _msin( ₀t_i+ ₀),

где B_m - амплитуда опорного сигнала,

при частоте опорного сигнала ₀= _k строят вольт-амперную характеристику а(b₀) для исходного сигнала a(t _i), многократно сдвигают по фазе, определяя площадь вольт-амперной характеристики F_BAX( ₀=-180°...+180°),

где ₀ - фаза опорного сигнала. Полученные значения F_BAX( ₀=var) запоминают, выбирают значение фазы частотной составляющей _k= ₀, при котором F_BAX 0, значения амплитуды и фазы _k, соответствующие частоте _k, запоминают, данную операцию повторяют для всех найденных значений частоты _k, далее по полученным значениям _k, _k судят о спектральном составе исходного сигнала a(t_i).

Заявленный способ спектрального анализа многочастотных периодических сигналов, представленных цифровыми отсчетами, имеет существенные преимущества, так как повышает точность определения фазы частотной составляющей исходного сигнала за счет того, что фазы всех частотных составляющих сигнала определяют на основе анализа вольт-амперной характеристики a(b₀). Предлагаемый способ позволяет определить спектральный состав сигнала в интересующем диапазоне частот.

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

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

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

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

В табл.2 приведены результаты спектрального анализа по предложенной процедуре для тестового исходного сигнала.

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

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

С выхода датчика анализируемого сигнала 1 (ДАС) исходный сигнал a(t _i) поступает на вход первого программатора 2 (П1) и второго программатора 5 (П2). В первом программаторе 2 (П1) определяют значения мгновенной спектральной плотности S₁ ( _j), S₂( _j) и S( _j) по формулам (1) с задаваемым шагом и диапазоном ₁,..., _n по частоте для всех заданных частот. Затем выявляют экстремумы из условия все амплитуды и фазы частотных составляющих, не удовлетворяющие этому условию, в дальнейших исследованиях не учитываются. Таким образом, с выхода первого программатора 2 (П1) на входы третьего программатора 4 (П3) и датчика опорного сигнала 6 (ДОС) поступают последовательно значения частот _k, содержащиеся в исходном сигнале a(t_i), начиная с минимальной. Сразу после поступления частоты составляющей сигнала _k в датчике опорного сигнала 6 (ДОС) начинают формироваться опорные сигналы b₀ (t_i)=B_msin( ₀t_i+ ₀) с частотой ₀= _k для различных фаз опорного сигнала ₀=-180°...+180°. Одновременно с выхода первого программатора на вход вычислителя 3 (В) поступают последовательно значения мгновенной спектральной плотности S( _k) для определенных частот _k, где определяют амплитуды Am( _k) частотных составляющих _k по формуле (2). С датчика опорного сигнала (ДОС) сигналы b₀(t _i) поступают на вход второго программатора 5 (П2), где строят вольт-амперную характеристику а(b₀ ), определяют значение ее площади F_BAX по формуле

полученные значения F_BAX ( ₀=var) запоминают, выбирают значение _k= ₀ при котором F_BAX 0 и передают на третий программатор 4 (П3). С выхода вычислителя 3 (В) значение передают на вход третьего программатора 4 (П3), где значения амплитуды и фазы _k соответствующие частоте _k, запоминают. Данную операцию повторяют для всех значений, определенных в первом программаторе значений частоты _k. После окончания последней операции на выход программатора 5 (П) подают сигнал о готовности результатов спектрального анализа и либо передают на ЭВМ, либо выводят на дисплей.

Для проверки работоспособности предложенной процедуры провели спектральный анализ тестового сигнала напряжения, содержащего интергармоническую частотную составляющую (параметры сигнала приведены в табл.1) a(t_i)=u(t _i)=10·sin(2 ·50·t_i+130)+4·sin(2 ·100·t_i-26)+1·sin(2 ·170·t_i+150). Для этого рассчитали массив мгновенных значений для N=10000 точек исходного сигнала с шагом дискретизации t=10^-4 с (осциллограмма тестового исходного сигнала приведена на фиг.2) и передали его на первый программатор 2 (П1). В первом программаторе 2 (П1) был задан диапазон частот 1...200 Гц и шаг f=1 Гц. Провели расчет для интергармонической составляющей f=170 Гц

Значения спектральной плотности для различных частот приведены на фиг.3. Далее в первом программаторе 2 (П1) исходя из условия делается вывод о наличии экстремума на рассматриваемой частоте, при этом

S(2 ·170)=5·10³>500,

это говорит о наличии экстремума на частоте _k=2 ·170 Рад, поэтому данную частотную составляющую использовали в дальнейших расчетах. В вычислителе 3 (В) определяли амплитуду по формуле (2)

В датчике опорного сигнала (ДОС) формировали опорный сигнал b₀(t_i) для различных фаз ₀=-180°...+180°, при этом опорный сигнал и значение фазы передавали во второй программатор 5 (П2), где определяли площадь вольт-амперной (a(b ₀)) характеристики F_BAX по формуле (3). Для частоты _k=2 ·170 рад площадь вольт-амперной характеристики будет минимальной F_BAX=1,166·10 ^-3 при ₀=150,726°.

Результаты проверки работоспособности заявленного метода для других заданных частот в тестовом исходном сигнале приведены в табл.2. Сравнив параметры, указанные в табл.1 и табл.2, в частности значения фаз анализируемого сигнала _n, расчетные значения фазы, рассчитанные как в прототипе , и расчетные значения фазы, рассчитанные по предложенному способу _k, можно говорить о точности определения параметров спектрального анализа каждой отдельной частотной составляющей.

О спектральном составе тестового исходного сигнала а(t _i) судят по полученным значениям амплитуды Am _k, частоты _k и фазы _k.

Таблица 1
, В	10	4	1
f n, Гц	50	100	170
n, град	130	-26	150

Таблица 2
, В	10	4	1
S1(f k)	-3,214·104	-1,798·104	-4,33·103
S2(fk)	3,83·104	- 8,767·10 3	2,5·103
S(fk)	5,0·104	2·104	5·10 3
fk , Гц	50	100	170
k, град (как в прототипе)	-50	-26	-30
FВАХ, отн.ед.	1,628·10-13	-3,319·10-13	1,166·10-3
k, град (предлагаемый способ)	130	-26	150,726