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

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

Способ определения действующего значения гармонических составляющих в сигнале, основанный на преобразовании входного сигнала, в соответствии с которым выделяют определенные временные интервалы для анализа, отличающийся тем, что определяют действующее значение входного сигнала, определяют частоту первой гармоники входного сигнала, формируют опорный синусоидальный сигнал с частотой первой гармоники входного сигнала, выбирают первый временной интервал, в котором входной сигнал не изменяет свой знак, определяют середину первого временного интервала, выбирают второй временной интервал, в котором опорный синусоидальный сигнал не изменяет свой знак, определяют середину второго временного интервала, измеряют мгновенные значения входного сигнала в момент времени t₁, а опорного сигнала в момент времени t₂, причем моменты времени t₁ и t₂ выбирают равноотстоящими соответственно от середин первого и второго временных интервалов, определяют модуль отношения мгновенных значений входного и опорного сигналов, многократно определяют модули отношения мгновенных значений входного и опорного сигналов для каждой пары моментов времени, равноотстоящих от середин соответственно первого и второго временных интервалов, определяют усредненное значение модулей отношения мгновенных значений входного и опорного сигналов, а действующее значение гармонических составляющих входного сигнала определяют как корень квадратный из разности квадратов двух величин действующего значения входного сигнала и произведения действующего значения опорного синусоидального сигнала на усредненное значение модулей отношения мгновенных значений входного и опорного сигналов.

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

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

Известен способ спектрального анализа, в соответствии с которым для определения действующего значения высших гармоник в каждом канале сигнал непосредственно фильтруют в узкой полосе частот, детектируют и интегрируют (или возводят в квадрат, детектируют и усредняют) [1.а, б]

Такому способу присущи недостатки громоздкость фильтров не инфранизких частотах, большая погрешность от неравномерности АЧХ фильтров, малое быстродействие.

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

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

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

Недостатки способа на инфранизких частотах аналогичные погрешности АЧХ фильтров и низкое быстродействие.

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

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

Целью изобретения является повышение точности измерений при расширении функциональных возможностей.

Цель в способе определения действующего значения гармонических составляющих в сигнале, основанном на преобразовании входного сигнала, в соответствии с которым выделяют определенные временные интервалы для анализа, достигается тем, что определяют действующее значение входного сигнала, определяют частоту первой гармоники входного сигнала, формируют опорный синусоидальный сигнал с частотой первой гармоники входного сигнала, выбирают первый временной интервал, в котором входной сигнал не изменяет свой знак, определяют середину первого временного интервала, выбирают второй временной интервал, в котором опорный синусоидальный сигнал не изменяет свой знак, определяют середину второго временного интервала, измеряют мгновенные значения входного сигнала в момент времени t₁, а опорного сигнала в момент времени t₂, причем моменты времени t₁ и t₂ выбирают равноотстоящими соответственно от середины первого и второго временных интервалов, определяют модуль отношения мгновенных значений входного и опорного сигналов, многократно определяют модули отношения мгновенных значений входного и опорного сигналов для каждой пары моментов времени, равноотстоящих от середин соответственно первого и второго временных интервалов, определяют усредненное значение модулей отношения мгновенных значений входного и опорного сигналов, а действующее значение гармонических составляющих входного сигнала определяют как корень квадратный из разности квадратов двух величин действующего значения входного сигнала и произведения действующего значения опорного синусоидального сигнала на усредненное значение модулей отношения мгновенных значений входного и опорного сигналов.

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



где U_д действующее значение входного сигнала;

U_од действующее значение опорного синусоидального сигнала;

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

Входной квазисинусоидальный сигнал U_x(t) и опорный синусоидальный сигнал U_y(t) с частотой "" (первой гармоники) представим в виде отдельных функций, рассматриваемых на выбранных временных интервалах bj, не содержащих сигналов, равных нулю:

U_x(t)=U_x(bj), (1)

U_y(t)=U_y(bj), (2)

где t время;

U_x(bj), U_y(bj) соответствующие сигналы на выбранных временных интервалах bj.

Сигналы U_x(bj) и U_y(bj) на выбранных временных интервалах bj будем аппроксимировать в виде фрагментов синусоид, для которых



где A_x, A_y значения амплитуд аппроксимирующих сигналов;

значения круговой частоты сигналов;

F_x, F_y начальные фазы входного и опорного сигналов.

Разделим мгновенные значения входного и опорного сигналов друг на друга и рассмотрим функцию f(bj):

f(bj)=K_a[sin(t+F_x)]/sin(t+F_y), (5)

где f(bj) функция на интервале времени bj, определяемая отношением двух сигналов U_x(bj) и U_y(bj);

K_a=(A_x/A_y) отношение амплитуд сигналов из (3) и (4).

Для момента времени t₀ на интервале bj, когда значение функции f(bj)= f(t₀)=K_a, должно выполняться условие:

[sin(t_o+F_x)]/[sin(t+F_y)]=1 (6)

Обозначим левую часть уравнения (6) через L, тогда:



Разность фаз F₀ между сигналами U_x(bj) и U_y(bj) равна:

F₀=F_x-F_y, (8)

Если, к примеру, F_x>F_y, то можно принять F_y=0, и после преобразования из выражения (7) получим:

L=cosF_o+(sinF_o)/(tgt) (9)

Если F_x<F, то можно принять F_x=0, и после преобразования из выражения (7) получим:

L=1/[cosF_o+(sinF_octgt)] (10)

Выполнение условия (6) сводится к выполнению условия:

cosF_o+sinF_o/[tg(2/T)t_o]=1, (11)

где t₀ соответствует искомому моменту времени;

T период первой гармоники входного сигнала.

Обозначим (2/T)t_o= значение угла, определяемого положением t₀ на интервале времени периода T. Тогда после перестановок выражение (11) перепишем в следующем виде:

tg=sinF_o/(1-cosF_o). (12)

После преобразований получим:

sinF₀/(1-cosF₀) ctg(F₀/2) (13)

Из (12 и 13) следует:

tg=ctg(F_o/2). (14)

Из выражения (14) получим:

tg=tg[(/2)-(F_o/2)] (15)

После преобразования получим:

tg=tg[(-F_o)/2] (16)

Из равенства (16) получаем выражение для :

b=(-F_o)/2 . (17)

Так как =(2/T) t₀ соответствует моменту времени, когда выполняется условие (6), то из (17) определим положение точки t₀ на интервале bj. Угол соответствует полупериоду, то есть такому интервалу времени, когда значение входного или опорного сигнала имеет один и тот же знак. Положение точки t₀ на интервале bj соответствует середине интервала времени, лежащего внутри одного из полупериодов, из которого исключен интервал времени, соответствующий сдвигу фаз между сигналами. Этот момент времени t₀ находится на одинаковом расстоянии от середин выбранных интервалов времени, где не происходит изменение знака входного или опорного сигнала, см. фиг.1.

Значения синусоидальных сигналов в момент времени t₀ на середине интервала, равного фазовому сдвигу F₀, будут равны соответственно A_xsin(F₀/2) и A_ysin(F₀/2), а отношение мгновенных значений сигналов будет определяться как K_a=A_x/A_y.

Измерение мгновенного значения входного сигнала можно проводить в момент времени t₁, отстоящий от середины выбранного первого временного интервала этого сигнала, к примеру, на интервал времени, соответствующий фазовому сдвигу F₀ или [(/2) F₀] а измерение опорного сигнала в момент времени t₂, отстоящий от середины выбранного второго временного интервала своего сигнала на интервал времени соответственно F₀ или [ (/2) F₀] Для момента времени t₁ мгновенное значение входного сигнала можно представить, как A_xsinF₀= A_xcos[(/2) F₀] а мгновенное значение опорного сигнала в момент времени t₂ можно представить, как A_ysinF₀=A_ycos[(/2) F₀] Модуль отношения мгновенных значений входного и опорного синусоидальных сигналов будет равен K_a=A_x/A_y.

Мгновенное значение одного сигнала в момент времени будет равно A_xcosF₀, другого сигнала в момент времени будет равно A_ycosF₀. Модуль их отношений будет равен K_a A_x/A_y.

Следовательно, если измерять отношения мгновенных значений сигналов синусоидальной формы в различные моменты времени, равноотстоящие от середин выбранных временных интервалов полуволн своих сигналов, то модули K(t) отношений мгновенных значений входного и опорного сигналов будут равны K(t) K_a=A_x/A_y.

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

На фиг.1 показан пример определений моментов времени t₀, t₁, t₂, равноотстоящих от середин выбранных временных интервалов-полуволн входного и опорного сигналов, для произвольного фазового сдвига между сигналами. Для фазового сдвига F₀ получают несколько пар моментов времени t₁; t₂; и так далее. A_xsinF₀/A_ysinF₀= K_a для измерений в моменты времени t₁ и t₂; A_xcosF₀/A_ycosF₀ для измерений в моменты времени и так далее.

Измерения в момент времени t₀ можно рассматривать как частный случай выбора соответствующей пары, когда t₁=t₂.

Действующее значение U_e гармонических составляющих U_i(t) во входном сигнале U(t) можно определить из соотношений [6]



где U(t) представляет из себя функцию входного сигнала без составляющей первой гармоники входного сигнала.

Запишем входной сигнал напряжения U(t) в виде:

U(t)=sint+g(t) (19)

где амплитуда синусоидального напряжения первой гармоники входного сигнала;

g(t) некоторая функция, значения которой изменяются во времени так, чтобы выполнялось равенство (19).

Предположим, что опорный синусоидальный сигнал U(y)=U"(t) имеет частоту w первой гармоники, тогда:

U(t)=A_osin(t+F_o) (20)

где A₀ амплитуда опорного синусоидального сигнала;

F₀ сдвиги фаз между входным и опорным сигналами.

Для каждого из фазовых сдвигов F₀ будем определять модули K(t) отношения мгновенных значений входного U(t₁) и опорного U"(t₂) сигналов. Тогда входной сигнал U(t) можно представить в виде:

U(t)=[A_osin(t)]K(t) (21)

где K(t) текущие значения модулей отношения мгновенных значений входного U(t) и опорного U"(t) сигналов, вычисленных для различных пар моментов времени t₁ и t₂, равноотстоящих от середин соответственно первого и второго выбранных временных интервалов. При этом можно записать:



Если входной сигнал синусоида, то есть g(t)=0 в (19), то при различных фазовых сдвигах при усреднении значений модулей K(t) отношения мгновенных значений входного и опорного сигналов получают одно и то же значение K_c, которое численно будет равно K_c /A₀, и тогда входной сигнал U(t) представим в виде:

U(t)=(A_oK_c)sin(t)+g(t) (23)

Из (21) и (23) определим выражение для функции g(t):



где K₀ усредненное значение модулей отношения мгновенных значений входного и опорного сигналов.

Подставляя (25) в (19), получим для входного сигнала:

U(t)=A_oK_csin(t)+A_osin(t)[K(t)-K_c] (26)

Из (26) можно определить значение квадрата входного сигнала:



Из [6] получим выражение для квадрата действующего значения входного сигнала в следующем виде:





После преобразования выражение (29) будет иметь вид:



Из (30) среднеквадратическое значение входного сигнала можно представить в следующем виде:



Первое слагаемое правой части выражения (31) представляет собой квадрат произведения действующего значения опорного синусоидального сигнала на усредненное значение K_c модулей отношения мгновенных значений входного и опорного сигналов, что соответствует составляющей первой гармоники входного сигнала, поэтому (31) можно записать в следующем виде:



где U_oд действующее значение опорного синусоидального напряжения;

K(t) -текущие значения модулей отношения мгновенных значений входного и опорного сигналов;

K_с усредненное значение модулей отношения мгновенных значений входного и опорного сигналов.

Разность квадратов двух величин K(t) и K_c в выражении (32) можно представить в виде квадрата корня квадратного из разности квадратов этих же величин, поэтому (32) можно представить в следующем виде:



Первое слагаемое под общим корнем в правой части выражения (33) представляет собой квадрат действующего значения составляющей первой гармоники входного сигнала, а второе слагаемое равно квадрату действующего значения гармонических составляющих входного сигнала.

Следовательно, выражение (33) можно представить в следующем виде:



Из (34) получим окончательное выражение для действующего значения гармонических составляющих U_e входного сигнала:



Таким образом, действующее значение гармонических составляющих входного сигнала равно корню квадратному из разности квадратов двух величин - действующего значения входного сигнала и произведения действующего значения опорного синусоидального сигнала на усредненное значение K_c модулей отношения K(t) мгновенных значений входного и опорного сигналов.

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

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

Устройство содержит формирователь 1, управляемый опорный генератор 2, фазовращатель 3, блок индикации 4, блок 5 для измерения действующего значения входного сигнала. Блоки в устройстве соединены следующим образом. Формирователь 1, управляемый опорный генератор 2 и фазовращатель 3 соединены последовательно. Вход устройства соединен с входом блока 5 для измерения действующего значения входного сигнала, с входом формирователя 1 и с первым входом блока 4 индикации, второй вход которого подключен к выходу фазовращателя 3.

Устройство работает следующим образом. Входной сигнал U_x(t) поступает на блок 5 для измерения действующего значения входного сигнала, на первый вход блока 4 индикации и на вход формирователя 1, на выходе которого формируется напряжение U1, пропорциональное периоду T первой гармоники входного сигнала. Это напряжение U1 поступает на вход управляемого опорного генератора 2, на выходе которого получают напряжение U2 синусоидальной формы, период которых зависит от управляемого напряжения U1 и равен периоду T первой гармоники входного сигнала.

Синусоидальное напряжение U2 амплитудой A_y с выхода опорного генератора 2 поступает на вход фазовращателя 3, на выходе которого получают синусоидальное напряжение U_y(t) той же амплитуды, которая не изменяется при изменениях фазовых сдвигов. Для синусоидального напряжения амплитудой Ay известно его действующее значение, которое равно A_y. Таким образом, на два входа блока индикации двухлучевого осциллографа 4 поступают входные сигналы U_x(t) и опорные синусоидальные сигналы напряжения U_y(t), имеющие между собой фазовой сдвиг, к примеру, как показано на фиг.1. Для каждого фазового сдвига F₀ производят измерения мгновенных значений входного U_x(t₁) и опорного U_y(t₂) сигналов, определяют модули K(t) отношения мгновенных значений при выбранной паре моментов времени t₁ и t₂ соответственно, каждый из которых равноотстоит от середины выбранных временных интервалов, внутри которых входной и опорный сигналы соответственно не изменяют своих знаков. Эти временные интервалы соответствуют полуволнам входного и опорного сигналов. Вычисленные значения модулей K(t) усредняют и получают значения K_c.

Действующее значение гармонических составляющих входного сигнала определяют как корень квадратный из разности квадратов двух величин - действующего значения U5 на выходе блока 5 и произведения действующего значения A_y опорного синусоидального сигнала U_y(t) на усредненное значение K_c модулей Kt отношения мгновенных значений входного и опорного синусоидального сигналов.

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

Используемые источники информации

1. Мирский Г.Я. Аппаратурное определение характеристик случайных процессов. М. Энергия, 1972, а) с.234, б) с. 247.

2. Калинчук Б.А. Пиастро В.П. Анализаторы инфранизких случайных процессов. М. Энергия, 1973, с.116.

3. Авторское свидетельство СССР N 1287034, кл.G01R 23/16, 1987.

4. Авт. свид. СССР N 1113751, кл.G01R, 23/16, 1984 (прототип).

5. Измерения в электронике. Справочник. М. Энергоатом, 1987, с.382 -383.