способ определения угла сдвига фаз между двумя сигналами, представленными цифровыми отсчетами

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

Способ определения угла сдвига фаз между двумя сигналами a(t _j) и b(t_j), представленными цифровыми отсчетами с постоянным шагом дискретизации, включающий определение точек перехода сигналов через нулевой уровень от отрицательного значения к положительному в течение одного периода этих сигналов, определение угла сдвига фаз между ними, который пропорционален времени между переходами исследуемых сигналов через нулевой уровень от отрицательных значений к положительным, отличающийся тем, что вычисление углов сдвига фаз производят для сигналов, представленных цифровыми отсчетами, а определение точек перехода сигналов через нулевой уровень от отрицательных значений к положительным осуществляют путем запоминания каждого цифрового отсчета как текущего a(t _j), соответствующего текущему моменту времени t _j, так и предыдущего a(t_j-1), затем производят сравнение значений отсчетов a(t_j-1 ) и a(t_j), на основе которого определяют номера счетных импульсов w и i, соответствующих моментам времени t_w и t_i, при которых происходит первый переход сигнала a(t_j) через нулевой уровень от отрицательного значения к положительному, затем запоминают значение отсчета a(t_p), соответствующего моменту времени t_i+1, далее на основе данных отчетов производят параболическую интерполяцию с помощью квадратного трехчлена

a(t_j )=d·t²+h·t+g, где

а определение момента времени t _n1, соответствующего первому переходу сигнала a(t _j) через нулевой уровень от отрицательного значения к положительному производят следующим образом

затем по аналогичной процедуре определяют момент времени t_n2, соответствующего второму переходу сигнала a(t_j) через нулевой уровень от отрицательного значения к положительному, далее определяют значение фазы _а первого сигнала

одновременно с этим по аналогичной процедуре определяют значение фазы _b второго сигнала, значения угла сдвига фаз между сигналами a(t) и b(t) определяют как разность значений _а и _b.

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

Изобретение относится к измерительной технике.

Известен способ определения угла сдвига фаз между двумя сигналами [Авторское свидетельство СССР №1406511, G01R 25/00, опубл. 30.09.1983], заключающегося в выборе тактовой частоты в течение одного периода входного сигнала при одновременном измерении периода исследуемого сигнала и временного интервала, пропорционального измеряемому сдвигу фаз, и последующем автоматическом преобразовании результатов этих двух измерений в непосредственный отсчет сдвига фаз в градусах. Причем выбор тактовой частоты осуществляется методом подбора частот из исследуемого поддиапазона.

Недостатками способа являются: сложность процедуры измерения и схемной реализации.

Также известен способ определения угла сдвига фаз между двумя сигналами [Метрология и радиоизмерения: Учебник для вузов / В.И. Нефедов, В.И. Хахин, В.К. Битюков и др. / Под ред. профессора В.И. Нефедова. - М.: Высш. шк., 2003. 302-304 с.], выбранный в качестве прототипа, в котором в течение периода Т отслеживают точки пересечения нулевого уровня исследуемыми сигналами. Затем формируют импульс, имеющий длительность t и соответствующий сдвигу фаз между сигналами во времени.

Его заполняют счетными импульсами, количество n которых равно

где - период следования счетных импульсов, m=1, 2, 3....

Угол сдвига фаз между сигналами определяется по соотношению

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

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

Это достигается тем, что в способе определения угла сдвига фаз между двумя сигналами а(t_j ) и b(t_j), представленными цифровыми отсчетами, так же как и в прототипе, включающем определение точек перехода сигналов через нулевой уровень от отрицательных значений к положительным в течение одного периода этих сигналов, и определение угла сдвига фаз между ними, который пропорционален времени между переходами исследуемых сигналов через нулевой уровень от отрицательных значений к положительным.

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

Затем производят сравнение значений отсчетов a(t_j-1) и a(t _j), на основе которого определяют номера счетных импульсов w и i, соответствующих моментам времени t_w и t_i, при которых происходит первый переход сигнала a(t_j) через нулевой уровень от отрицательного значении к положительному. Затем запоминают значение отсчета a(t_p), соответствующего моменту времени t_i+1, далее на основе данных отчетов производят параболическую интерполяцию с помощью квадратного трехчлена

где

Затем производят определение момента времени t_n1, соответствующего первому переходу сигнала а(t_j) через нулевой уровень от отрицательного значения к положительному

Далее по аналогичной процедуре определяют момент времени t_n2, соответствующий второму переходу сигнала а(t_j) через нулевой уровень от отрицательного значения к положительному и затем определяют значение фазы _a первого сигнала

одновременно с этим по аналогичной процедуре определяют значение фазы _b второго сигнала. Значения угла сдвига фаз между сигналами a(t) и b(t) определяют как разность значений _a и _b.

Известно, что угол сдвига фаз между двумя сигналами определяется как разность этих колебаний при прохождении ими определенного (например, нулевого) уровня [Меерсон A.M. Радиоизмерительная техника. - Л.: Энергия, 1978. С.247].

Рассмотрим произвольный сигнал x(t _j), представленный на фиг.1. Для определения фазы данного сигнала производят поиск пар отсчетов, при которых происходит смена полярности сигнала. Затем для повышения точности определения момента перехода сигнала x(t_j) через нулевой уровень от отрицательного значения к положительному производят его интерполяцию на отрезке [t_w; t _р] квадратным трехчленом (3).

Определение момента времени t_n1, соответствующего первому переходу сигнала х(t_j) через нулевой уровень от отрицательного значения к положительному определяют по переходу через нулевой уровень интерполированной функции (3), описывающей сигнал:

Решение уравнения (7) находится по соотношению (5).

Затем по аналогичной процедуре определяют момент времени, соответствующий длительности периода сигнала. Далее по соотношению (6) производят определение фазы сигнала х(t_j ).

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

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

На фиг.1 приведена осциллограмма произвольного сигнала x(t_j).

На фиг.2 приведена функциональная блок-схема способа.

Предложенный способ определения угла сдвига фаз между сигналами может быть реализован, например, с помощью функциональной блок-схемы, которая представлена на фиг.2. Она содержит: блок кольцевой памяти 1 (КП); программатор вычисления углов сдвига фаз между сигналами 2 (П).

К входам блока кольцевой памяти 1 (КП) и программатора 2 (П) присоединен источник исследуемых сигналов, представленных в виде цифровых отсчетов. А выходы блока кольцевой памяти 1 (КП) соединены с входами программатора 2 (П).

Выходы программатора 2(П) связаны с входами сегментных индикаторов (не показаны на фиг.2).

Блок кольцевой памяти 1 (КП) может быть реализован на внешней перезаписываемой памяти данных Amtel AT25L256 (32 кБайта). Программатор 2 (П) может быть выполнен на микроконтроллере серии 51 производителя Atmel AT89S53.

Сигналы а(t _j) и b(t_j), представленные в виде цифровых отсчетов, поступают одновременно на вход блока кольцевой памяти 1 (КП) и на вход программатора 2 (П).

В блоке кольцевой памяти 1 (КП) запоминают текущие значения сигналов a(t _j) и b(t_j) для использования их в последующий момент времени.

Затем в программаторе 2 (П) определяют номера счетных импульсов w_a1 и i_a1, при которых происходит первый переход сигнала а(t_j) через нулевой уровень от отрицательного значения к положительному, затем запоминают в блоке кольцевой памяти 1 (КП) значение отсчета соответствующего моменту времени которое также поступает на вход программатора, где определяют коэффициенты интерполированной функции, описывающей сигнал a(t _j) на отрезке

и определяют момент времени соответствующий первому переходу сигнала a(t _j) через нулевой уровень от отрицательного значения к положительному:

Далее определяют номера счетных импульсов w _a2 и i_a2, при которых происходит второй переход сигнала a(t_j) через нулевой уровень от отрицательного значения к положительному, затем запоминают в блоке кольцевой памяти 1 (КП) значение отсчета соответствующего моменту времени которое также поступает на вход программатора.

Затем в программаторе 2 (П) определяют коэффициенты интерполированной функции, описывающей сигнал a(t_j) на отрезке

и определяют момент времени соответствующий второму переходу сигнала а(t _j) через нулевой уровень:

Далее определяют значение фазы _a сигнала а(t_j )

Одновременно с этим в программаторе 2 (П) определяют номера счетных импульсов w_b1 и i _b1, при которых происходит первый переход сигнала a(t _j) через нулевой уровень от отрицательного значения к положительному, затем запоминают в блоке кольцевой памяти 1 (КП) значение отсчета соответствующего моменту времени которое также поступает на вход программатора.

Затем в программаторе 2 (П) определяют по аналогичной процедуре коэффициенты интерполированной функции d_b1, h _b1 и g_b1, описывающей сигнал b(t _j) на отрезке и определяют момент времени соответствующий первому переходу сигнала b(t _j) через нулевой уровень от отрицательного значения к положительному:

Далее определяют номера счетных импульсов w _b2 и i_b2, при которых происходит второй переход сигнала b(t_j) через нулевой уровень от отрицательного значения к положительному, затем запоминают в блоке кольцевой памяти 1 (КП) значение отсчета соответствующего моменту времени которое также поступает на вход программатора.

Затем в программаторе 1 (П) определяют по аналогичной процедуре коэффициенты интерполированной функции d_b2, h _b2 и g_b2, описывающей сигнал b(t _j) на отрезке и определяют момент времени соответствующий второму переходу сигнала b(t _j) через нулевой уровень от отрицательного значения к положительному:

Значение фазы _b сигнала b(t_j ) вычисляют по следующему соотношению

и вычисляют угол сдвига фаз _a,b между сигналами a(t _j) и b(t_j)

С выходов программатора 2 (П) сигнал поступает на соответствующие входы сегментных индикаторов (не показаны на фиг.2).

В качестве примера, возьмем два тестовых сигнала несинусоидальной формы вида:

где =314 с^-1.

При этом период следования счетных импульсов равен 0,625 мс.

В программаторе 2 (П) одновременно вычисляют значения коэффициентов интерполированной функции по соотношениям (8) и (10) и определяют значения _a и _b:

Согласно условию соотношения (13) выбирают значение

Согласно условию соотношения (11) выбирают значение

Затем вычисляют значения угла сдвига фаз _a,b между сигналами a(t _j) и b(t_j)

Относительную погрешность вычисляют по формуле [Бронштейн И.Н., Семендяев К.А. Справочник по математике для инженеров и учащихся Втузов. - М.: Наука, 1980. - 976 с.], в программаторе 2 (П) определяют номера счетных импульсов w _b1 и i_b1, при которых происходит первый переход сигнала а(t_j) через нулевой уровень от отрицательного значения к положительному, затем запоминают в блоке кольцевой памяти 1 (КП) значение отсчета соответствующего моменту времени которое также поступает на вход программатора.

где _T - значение угла сдвига фаз между сигналами, полученное путем решения уравнений

где x= t для сигнала a(t_j) и x ₁= t для сигнала b(t_j), а _T определяется по следующему соотношению

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