(21), (22) Заявка: 2007111016/28, 26.03.2007
(24) Дата начала отсчета срока действия патента:
26.03.2007
(46) Опубликовано: 10.08.2008
(56) Список документов, цитированных в отчете о поиске:
RU 2264631 С1, 31.03.2004. RU 2264630 C1, 31.03.2004. RU 2242014 С1, 31.12.2003. RU 2004109897 A, 10.10.2005. RU 2004106034 A, 10.08.2005.
Адрес для переписки:
634050, г.Томск, пр. Ленина, 30, Томский политехнический университет, отдел интеллектуальной и промышленной собственности
|
(72) Автор(ы):
Гольдштейн Ефрем Иосифович (RU), Кац Илья Маркович (RU)
(73) Патентообладатель(и):
Государственное образовательное учреждение высшего профессионального образования Томский политехнический университет (RU)
|
(54) СПОСОБ ОПРЕДЕЛЕНИЯ УГЛА СДВИГА ФАЗ МЕЖДУ ДВУМЯ СИГНАЛАМИ
(57) Реферат:
Способ заключается в измерении мгновенных значений двух сигналов, оцифровке их для одних и тех же моментов времени. Сохранении каждого цифрового отсчета как текущего и предыдущего. Определении разности и суммы каждой пары текущего и предыдущего значений. Перемножении разности и суммы, суммировании произведения. Определении реактивной квазимощности . Перемножении текущих отсчетов сигналов и определении их активной квазимощности , а угол сдвига фаз a,b между сигналами a(tj) и b(tj) определяется по соотношению:
В начальный момент времени формируют опорный сигнал s(tj) вида s(tj)=Amsin(tj), где Аm – амплитуда опорного сигнала, – круговая частота опорного сигнала. Производится определение углов сдвига фаз a,s и b,s между первым и опорным сигналами, вторым и опорным сигналами, а значения угла фазового сдвига между сигналами a(t) и b(t) определяется как разность значений a,s и b,s. Технический результат заключается в упрощении, повышении точности определения сдвига фаз между двумя сигналами синусоидальной и несинусоидальной формы. 1 ил.
Изобретение относится к радиоизмерительной технике и может быть использовано при разработке и построении цифровых фазометров.
Известен способ определения угла сдвига фаз между несинусоидальными сигналами одинаковой частоты, в котором с помощью алгоритма быстрого преобразования Фурье раскладывают исследуемые сигналы в спектр, по условию максимума амплитуды определяют главные гармоники и находят разность фаз между ними (Заявка 2004106034/28, МПК 7, G01R 25/00, опубл. 10.08.2005).
К недостаткам данного способа можно отнести сложность необходимых расчетов и неизбежную погрешность дискретного преобразования Фурье, ограничивающую точность определения угла сдвига фаз.
Известен способ определения угла сдвига фаз между двумя несинусоидальными сигналами (Патент РФ №2264631, МПК 7 G01R 25/00, опубл. 20.11.2005), выбранный в качестве прототипа, в котором измеряют мгновенные значения двух сигналов a(tj) и b(tj), оцифровывают их для одних и тех же моментов времени
где N – число разбиений на периоде Т.
Сохраняют каждый цифровой отсчет как текущий и предыдущий, далее определяют разность и сумму каждой пары текущего и предыдущего значений, перемножают разность и сумму, затем суммируют произведения, далее определяют реактивную квазимощность , затем перемножают текущие отсчеты сигналов и определяют их активную квазимощность , а угол сдвига фаз a,b между сигналами a(tj) и b(tj) определяют по соотношению:
Основным недостатком данного способа является низкая точность определения угла сдвига фаз между несинусоидальными сигналами.
Задачей изобретения является создание простого и точного способа определения угла сдвига фаз между двумя сигналами как синусоидальной, так и несинусоидальной формы.
Это достигается тем, что в способе определения угла сдвига фаз между двумя сигналами, включающем, так же как и в прототипе, измерение мгновенных значений этих сигналов, оцифровку их для одних и тех же моментов времени
где N – число разбиений на периоде Т;
сохранение каждого цифрового отсчета как текущий и предыдущий, определение разности и суммы каждой пары текущего и предыдущего значений, перемножение разности и суммы, суммирование произведения, определение реактивной квазимощности , перемножение текущих отсчетов сигналов и определение их активной квазимощности , определение угла сдвига фаз a,b между сигналами a(tj) и b(tj) по соотношению:
Согласно изобретению в начальный момент времени формируют опорный сигнал s(tj) вида
где Аm – амплитуда опорного сигнала, – круговая частота опорного сигнала.
По соотношению (4) одновременно производят определение углов сдвига фаз a,s и b,s между первым и опорным сигналами, вторым и опорным сигналами, а значения угла сдвига фаз между сигналами a(t) и b(t) определяют как разность значений a,s и b,s.
Авторами экспериментально установлено, что определение угла сдвига фаз по выражению (2) справедливо только в случае, если сигналы синусоидальной формы. Если же сигналы несинусоидальной формы, то точность вычислений по соотношению (2) достаточно мала.
Также было экспериментально установлено, что использование опорного сигнала синусоидальной формы s(tj) при определении углов сдвига фаз между сигналами несинусоидальной формы по соотношению (2) приводит к увеличению точности вычислений. При этом наблюдается вычисление углов сдвига фаз между первыми гармониками сигналов, что согласно [Меерсон A.M. Радиоизмерительная техника. – Л.: Энергия, 1978. С.247] является углом сдвига фаз между двумя сигналами несинусоидальной формы.
Таким образом, предлагаемый способ позволяет производить определения угла сдвига фаз между двумя сигналами как синусоидальной, так и несинусоидальной формы по массивам мгновенных значений токов и напряжений путем их измерения в течение одного периода.
На чертеже приведена функциональная блок-схема способа.
Предложенный способ диагностики может быть реализован, например, с помощью функциональной блок-схемы, которая представлена на чертеже. Она содержит: аналогово-цифровой преобразователь 1(АЦП); блок кольцевой памяти 2(КП); программатор вычисления углов фазовых сдвигов между сигналами 3(П); генератор опорных сигналов 4(ГОС). К входам аналогово-цифрового преобразователя 1(АЦП) присоединены трансформаторы тока и напряжения (не показаны на чертеже). Выходы преобразователя 1(АЦП) связаны с входами блока кольцевой памяти 2(КП) и с входами программатора 3(П). Выходы блока кольцевой памяти 2(КП) соединены с входами программатора 3(П).
Выходы программатора 3(П) связаны с входами сегментных индикаторов (не показаны на чертеже). Выходы генератора опорных сигналов 4(ГОС) связаны с входами программатора 3(П).
В качестве аналогово-цифрового преобразователя 1(АЦП) может быть выбран, например, аналого-цифровой преобразователь серии МАХ186 (12 бит). Блок кольцевой памяти 2(КП) может быть реализован на внешней перезаписываемой памяти данных Amtel AT25L256 (32 кБайта). Программатор 3 (П) может быть выполнен на микроконтроллере серии 51 производителя Atmel AT89S53. Генератор опорных сигналов 4 (ГОС) может быть выполнен на базе микросхемы МАХ038.
При поступлении аналоговых сигналов a(t) и b(t) с трансформаторов тока и напряжения (не показаны на чертеже) на вход аналого-цифрового преобразователя 1(АЦП) происходит их преобразование в цифровой сигнал при числе N дискретных значений выходного сигнала аналого-цифрового преобразователя на периоде входного сигнала. Оцифрованные сигналы a(tj) и b(tj), измеренные в момент времени tj, поступают одновременно на вход блока кольцевой памяти 2 (КП) и на вход программатора 3 (П). В блоке кольцевой памяти 2 (КП) запоминают текущие значения сигналов a(tj) и b(tj) для использования их в последующий момент времени. С выходов блока кольцевой памяти 2 (КП) сигналы a(tj-1) и b(tj-1), соответствующие моменту времени tj-1, поступают на вход программатора 3 (П). Одновременно с этим на соответствующие входы программатора 3 (П) поступают сигналы a(tj) и b(tj) с выходов аналогово-цифрового преобразователя 1 (АЦП). Одновременно с этим на вход программатора с выходов генератора опорного сигнала 4 (ГОС) поступает сигнал s(tj).
В программаторе 3 (П) происходит вычисления угла сдвига фаз между сигналом a(tj) и s(tj) по следующему выражению:
Одновременно с этим в программаторе 3(П) происходит вычисление угла сдвига фаз между сигналом b(tj) и s(tj) по следующему выражению:
Далее происходит вычисление угла сдвига фаз a,b между сигналами a(tj) и b(tj)
С выходов программатора 3 (П) сигнал поступает на соответствующие входы сегментных индикаторов (не показаны на чертеже).
В качестве примера, возьмем два тестовых сигнала несинусоидальной формы вида:
где =314-1,
которые подаются на входы аналогово-цифрового преобразователя 1 (АЦП). При этом число дискретных значений выходного сигнала аналого-цифрового преобразователя на периоде входного сигнала равно 32.
В качестве опорного сигнала возьмем сигнал вида
который задается с помощью генератора опорного сигнала 4(ГОС) и подается на вход программатора 3(П), в котором одновременно по соотношениям (6) и (7) происходит вычисления значений a,s и b,s.
Затем вычисляют значения угла сдвига фаз a,b между сигналами a(tj) и b(tj)
Относительную погрешность вычисляют по формуле [Бронштейн И.Н., Семендяев К.А. Справочник по математике для инженеров и учащихся Втузов. – М.: Наука, 1980. -976 с.]
где – значение угла сдвига фаз между первыми гармониками.
Таким образом, получен простой и точный способ определения сдвига фаз между двумя сигналами.
Формула изобретения
Способ определения угла сдвига фаз между двумя сигналами a(t) и b(t), заключающийся в том, что измеряют мгновенные значения этих сигналов, оцифровывают их для одних и тех же моментов времени
tj=t1, t2…tN,
где N – число разбиений на периоде Т,
сохраняют каждый цифровой отсчет как текущий и предыдущий, определяют разность и сумму каждой пары текущего и предыдущего значений, перемножают разность и сумму, суммируют произведения, далее определяют реактивную квазимощность , перемножают текущие отсчеты сигналов и определяют их активную квазимощность , а угол сдвига фаз a,b, между сигналами a(tj) и b(tj) определяют по соотношению
,
отличающийся тем, что в начальный момент времени формируют опорный сигнал s(tj) вида
s(tj)=Amsin(tj),
где Аm – амплитуда опорного сигнала;
– круговая частота опорного сигнала,
по вышеприведенному соотношению одновременно производят определение углов сдвига фаз a,s и b,s между первым и опорным сигналами, вторым и опорным сигналами, а значения угла фазового сдвига между сигналами a(t) и b(t) определяют как разность значений a,s и b,s.
РИСУНКИ
|