СПОСОБ ОПРЕДЕЛЕНИЯ СДВИГА ФАЗ МЕЖДУ ДВУМЯ СИНУСОИДАЛЬНЫМИ СИГНАЛАМИ Российский патент 2005 года по МПК G01R25/00 

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

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

Известен способ определения сдвига фаз двух синусоидальных сигналов [Патент РФ №2039360, МПК 6 G 01 R 25/00, опубл. 1995.07.09], заключающийся в том, что измеренные мгновенные значения отфильтровывают от постоянной составляющей сигналов X(t) и Y(t), имеющих период колебаний Т, измеряют два мгновенных значения одного из сигналов, принятого за измерительный, в моменты времени которые выбирают на интервале полуволны другого сигнала, принятого за опорный, а значение разности фаз определяют по формуле F₀=m(g+πn).

Недостатком известного способа являются многоэтапность и сложность его реализации.

Известен способ определения разности фаз двух синусоидальных сигналов [А.С. №1503025, МПК 4 G 01 R 25/00, опубл. 1987.04.27], выбранный в качестве прототипа, заключающийся в том, что измеряют мгновенные значения синусоидальных сигналов, синусоидальные сигналы отфильтровывают от постоянной составляющей, сдвигают каждый из них по фазе на угол в сторону опережения без изменения амплитуды, а разность фаз определяют по формуле

где U₁₁, U₁₂ - мгновенные значения соответственно первого и сдвинутого по отношению к нему на угол сигналов, измеренных в один момент времени;

U₂₁, U₂₂ - мгновенные значения соответственно второго и сдвинутого по отношению к нему на угол сигналов, измеренных одновременно со значениями U₁₁ и U₁₂.

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

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

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

t_j=t₁, t₂, ..., t_N,

где N - число разбиений на периоде Т,

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

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

Известна теорема Телледжена о квазимощности между любыми двумя синусоидальными сигналами, которые необязательно существуют в цепи одно и то же время и необязательно связаны одной зоной цепи [П.Пенфилд и др. Энергетическая теория электрических цепей / П.Пенфилд, Р.Спенс, С.Дюинкер. - М.: Энергия, 1974. - 152 с.] Также известно, что реактивная мощность может быть определена двояко:

1) по общеизвестной формуле

2) по методике О.А.Маевского [Маевский О.А. Энергетические показатели вентильных преобразователей. - М.: Энергия, 1978. - 320 с.]

где F_BAX - площадь вольт-амперной характеристики для исследуемых сигналов, найденная в данном случае по формуле для определения площади многоугольника, заданного координатами концов отрезков [Бронштейн И.Н., Семендяев К.А. Справочник по математике для инженеров и учащихся Втузов. - М.: Наука, 1980. - 976 с.]

Подставив формулу (4) в формулу (3) получим формулу для вычисления реактивной квазимощности

Экспериментально установлено, что выражения (2) и (5) вполне справедливы для вычисления реактивной квазимощности , тогда, приравняв правые части формул (2) и (5), можно найти sinϕ_ab

Описанный способ определения сдвига фаз обладает рядом преимуществ, в частности нет необходимости отфильтровывать сигналы от постоянной составляющей и сдвигать сигналы по фазе на угол 90°. Относительная погрешность при определении сдвига фаз в среднем составляет 0,77%.

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

На фиг.2 приведена схема, для которой измерены мгновенные значения сигналов a(t_j)=i₁(t_j), b(t_j)=i₂(t_j).

На фиг.3 приведена векторная диаграмма для наглядности проверки сдвига фаз между сигналами.

В табл.1 приведены цифровые отсчеты мгновенных значений сигналов a(t_j)=i₁(t_j), b(t_j)=i₂(t_j).

В табл.2 приведены результаты проверки работоспособности предлагаемого способа определения сдвига фаз.

Способ может быть осуществлен с помощью устройства (фиг.1). Оно включает в себя первое устройство выборки-хранения 1 (УВХ 1), второе устройство выборки-хранения 2 (УВХ 2), третье устройство выборки-хранения 3 (УВХ 3), четвертое устройство выборки-хранения 4 (УВХ 4), инвертор 5 (Инвертор), первый сумматор 6 (Сумматор 1), второй сумматор 7 (Сумматор 2), перемножитель 8 (Перемножитель), интегратор 9 (Интегратор), перемножитель-делитель 10 (Перемножитель-делитель), тактовый генератор 11 (ТГ), первый выпрямитель 12 (Выпрямитель 1), второй выпрямитель 13 (Выпрямитель 2), первый фильтр низких частот 14 (ФНЧ 1), второй фильтр низких частот 15 (ФНЧ 2).

Входные шины устройства подключены к входам устройств выборки-хранения: первого 1 (УВХ 1) и второго 2 (УВХ 2), выходы которых - к входам третьего 3 (УВХ 3), четвертого 4 (УВХ 4) устройств выборки-хранения. Первое устройство выборки-хранения 1 (УВХ 1) подключено к входам первого сумматора 6 (Сумматор 1). Второе устройство выборки-хранения 2 (УВХ 2) подсоединено к входу второго сумматора 7 (Сумматор 2). Вход третьего устройства выборки-хранения 3 (УВХ 3) подключен к входу инвертора 5 (Инвертор), выход которого соединен с входом первого сумматора 6 (Сумматор 1). Выход четвертого устройства выборки-хранения 4 (УВХ 4) подключен к входу второго сумматора 7 (Сумматор 2). Выходы первого 6 (Сумматор 1) и второго 7 (Сумматор 2) сумматоров связаны с входами перемножителя 8 (Перемножитель), выход которого соединен с входом интегратора 9 (Интегратор). Выход интегратора 9 (Интегратор) связан с входом перемножителя-делителя 10 (Перемножитель-делитель). Выходы тактового генератора 11 (ТГ) связаны с входами управления первого 1 (УВХ 1), второго 2 (УВХ 2), третьего 3 (УВХ 3) и четвертого 4 (УВХ 4) устройств выборки-хранения. Входы первого 12 (Выпрямитель 1) и второго 13 (Выпрямитель 2) выпрямителей подключены к входным шинам, а их выходы к входам первого 14 (ФНЧ 1) и второго 15 (ФНЧ 2) фильтров низких частот. Выходы первого 14 (ФНЧ 1) и второго 15 (ФНЧ 2) фильтров низких частот подключены к множительным входам перемножителя-делителя 10 (Перемножитель-делитель). Выход перемножителя-делителя 10 (Перемножитель-делитель) соединен с входом сегментного индикатора для вывода значения сдвига фаз.

Первое 1 (УВХ 1), второе 2 (УВХ 2), третье 3 (УВХ 3) и четвертое 4 (УВХ 4) устройства выборки-хранения могут быть реализованы на микросхемах 1100СК2. Инвертор 5 (Инвертор) может быть реализован на микросхеме 140УД17А. Первый 6 (Сумматор 1) второй 7 (Сумматор 2), сумматоры могут быть реализованы на операционных усилителях 140УД17А. В качестве перемножителя 8 (Перемножитель) и перемножителя-делителя 10 (Перемножитель-делитель) может быть использована микросхема 525ПС3. Интегратор 9 (Интегратор) может быть реализован на операционном усилителе 140УД17А. Тактовый генератор 11 (ТГ) может быть реализован на микроконтроллере АТ80С2051. Выпрямители 12 (Выпрямитель 1) и 13 (Выпрямитель 2), а также фильтры низких частот 14 (ФНЧ 1) и 15 (ФНЧ 2) могут быть выполнены операционных усилителях 140УД17А.

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

На вход первого устройства выборки-хранения 1 (УВХ 1) поступает сигнал, пропорциональный первому одночастотному синусоидальному сигналу, например a(t_j)=i₁(t_j)=10,9329sin(ωt_j-30°), а на вход второго устройства выборки-хранения 2 (УВХ 2) сигнал, пропорциональный второму одночастотному синусоидальному сигналу, например b(t_j)=i₂(t_j)=9,1926sin(ωt_j+40°),

где

t_j=t₁, t₂, ..., t_N,

- число разбиений на периоде Т,

Δt=1·10^-4 - дискретность массивов значений сигналов,

Массивы значений сигналов представлены в табл.1. Значения сигналов записываются в блоки выборки-хранения 1 (УВХ 1) и 2 (УВХ 2) и хранятся там, как текущие, затем с выхода устройства выборки-хранения 1 (УВХ 1) сигнал i₁(t_j) поступает в устройство выборки-хранения 3 (УВХ 3) и становится предыдущим значением, а с выхода устройства выборки-хранения 2 (УВХ 2), значение сигнала i₂(t_j) поступает в устройство выборки-хранения 4 (УВХ 4) и становится предыдущим значением. С выхода устройства выборки-хранения 3 (УВХ 3) предыдущее значение сигнала i₁(t_j) поступает в инвертор 5 (Инвертор). С помощью инвертора 5 (Инвертор) отрицательное значение предыдущего сигнала i₁(t_j) преобразуется в положительное. С выхода инвертора 5 (Инвертор) значение сигнала i₁(t_j) поступает на вход сумматора 6 (Сумматор 1). В то же время с выхода устройства выборки-хранения 1 (УВХ 1) текущее значение сигнала i₁(t_j) поступает на вход сумматора 6 (Сумматор 1). С помощью сумматора 6 (Сумматор 1) определяют разность текущего и предыдущего значений сигнала i₁(t_j). Одновременно с описанным выше процессом, с выхода устройства выборки-хранения 4 (УВХ 4) предыдущее значение сигнала i₂(t_j) поступает на вход сумматора 7 (Сумматор 2), а с выхода устройства выборки-хранения 2 (УВХ 2), текущее значение сигнала i₂(t_j) поступает на вход сумматора 7 (Сумматор 2). С помощью сумматора 7 (Сумматор 2) определяют сумму текущего и предыдущего значений сигнала i₂(t_j). С выхода сумматора 6 (Сумматор 1) разность текущего и предыдущего значений сигнала i₁(t_j) поступает на вход перемножителя 8 (Перемножитель), а с выхода сумматора 7 (Сумматор 2) сумма текущего и предыдущего значений сигнала i₂(t_j) поступает на вход перемножителя 8 (Перемножитель). С помощью перемножителя 8 (Перемножитель) значения разности и суммы сигналов перемножают и подают на вход интегратора 9 (Интегратор). С помощью интегратора 9 (Интегратор) суммируют произведения разности и суммы сигналов и определяют значение реактивной квазимощности . В данном случае (формула 5). С выхода интегратора 9 (Интегратор) значение реактивной квазимощности поступает на вход перемножителя-делителя 10 (Перемножитель-делитель). Параллельно с вышеописанным процессом сигналы i₁(t_j) и i₂(t_j) поступают на входы выпрямителей 12 (Выпрямитель 1) и 13 (Выпрямитель 2) соответственно. С выходов выпрямителей выпрямленные сигналы и подаются на входы фильтров низких частот 14 (ФНЧ 1) и 15 (ФНЧ 2). С помощью фильтра низких частот 14 (ФНЧ 1) определяют действующее значение сигнала . В данном случае I₁=7,7404. С помощью фильтра низких частот 15 (ФНЧ 2) определяют действующее значение сигнала . В данном случае I₂=6,5136. С выходов фильтров низких частот 14 (ФНЧ 1) и 15 (ФНЧ 2) действующие значения сигналов поступают на входы перемножителя-делителя 10 (Перемножитель-делитель),с помощью которого определяют сдвиг фазы между сигналами (формула 6). В данном случае

Если , то

Для наглядности проверки использовали векторную диаграмму, приведенную на фиг.3 и формулу

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

где

является приближенным значением числа

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

Табл.1Время t, ca(t_j)=i₁(t_j), Ab(t_j)=i₂(t_j), A1230-5,466455,9088890,0001-5,166356,1271660,0002-4,861156,3393970,0003-4,551166,5453710,0004-4,236676,7448850,0005-3,9186,9377430,0006-3,595477,1237550,0007-3,269387,3027360,0008-2,940077,474510,0009-2,607867,6389080,001-2,273087,7957670,0011-1,936057,9449330,0012-1,597118,0862570,0013-1,25668,2196020,0014-0,914848,3448350,0015-0,572188,4618330,0016-0,228968,570480,00170,1144878,6706690,00180,4578228,7623010,00190,8007068,8452850,0021,1427998,919540,00211,4837658,9849930,00221,8232669,0415790,00232,1609689,0892420,00242,4965379,1279350,00252,8296439,1576190,00263,1599569,1782670,00273,487159,1898560,00283,8109039,1923760,00294,1308969,1858240,0034,4468119,1702070,00314,7583389,145540,00325,0651699,1118480,00335,3670029,0691630,00345,6635379,0175280,00355,9544848,9569940,00366,2395548,8876210,00376,5184678,8094760,00386,7909478,7226380,00397,0567248,627191

Продолжение табл.1230,0047,3155388,523230,00417,5671328,4108580,00427,8112588,2901860,00438,0476768,1613320,00448,2761518,0244230,00458,4964597,8795960,00468,7083827,7269930,00478,9117117,5667630,00489,1062457,3990670,00499,2917927,2240680,0059,4681697,041940,00519,6352036,8528630,00529,7927276,6570220,00539,9405876,4546120,005410,078646,2458320,005510,206746,0308880,005610,324775,8099930,005710,432615,5833630,005810,530165,3512240,005910,617315,1138030,00610,693994,8713360,006110,760114,6240610,006210,815624,3722230,006310,860454,116070,006410,894563,8558550,006510,917923,5918350,006610,93053,324270,006710,93233,0534240,006810,923312,7795650,006910,903542,5029630,00710,873012,2238910,007110,831751,9426240,007210,77981,659440,007310,717211,3746190,007410,644041,088440,007510,560370,8011880,007610,466280,5131450,007710,361860,2245950,007810,24721-0,064180,007910,12245-0,352880,0089,987701-0,641240,00819,843095-0,92897

Продолжение табл. 11230,00829,688775-1,215780,00839,524893-1,501390,00849,351612-1,785520,00859,169101-2,067890,00868,977542-2,348210,00878,777123-2,626220,00888,568042-2,901640,00898,350506-3,174190,0098,124728-3,443610,00917,890932-3,709630,00927,649349-3,971990,00937,400217-4,230430,00947,143782-4,48470,00956,880297-4,734540,00966,610022-4,979710,00976,333223-5,219960,00986,050174-5,455060,00995,761155-5,684780,015,46645-5,908890,01015,16635-6,127170,01024,861152-6,33940,01034,551156-6,545370,01044,236669-6,744890,01053,918001-6,937740,01063,595466-7,123750,01073,269383-7,302740,01082,940073-7,474510,01092,607862-7,638910,0112,273078-7,795770,01111,93605-7,944930,01121,597111-8,086260,01131,256596-8,21960,01140,914841-8,344840,01150,572184-8,461830,01160,228961-8,570480,0117-0,11449-8,670670,0118-0,45782-8,76230,0119-0,80071-8,845290,012-1,1428-8,919540,0121-1,48376-8,984990,0122-1,82327-9,041580,0123-2,16097-9,089240,0124-2,49654-9,127930,0125-2,82964-9,15762

Продолжение табл.11230,0126-3,15996-9,178270,0127-3,48715-9,189860,0128-3,8109-9,192380,0129-4,1309-9,185820,013-4,44681-9,170210,0131-4,75834-9,145540,0132-5,06517-9,111850,0133-5,367-9,069160,0134-5,66354-9,017530,0135-5,95448-8,956990,0136-6,23955-8,887620,0137-6,51847-8,809480,0138-6,79095-8,722640,0139-7,05672-8,627190,014-7,31554-8,523230,0141-7,56713-8,410860,0142-7,81126-8,290190,0143-8,04768-8,161330,0144-8,27615-8,024420,0145-8,49646-7,87960,0146-8,70838-7,726990,0147-8,91171-7,566760,0148-9,10624-7,399070,0149-9,29179-7,224070,0159,46817-7,041940,0151-9,6352-6,852860,0152-9,79273-6,657020,0153-9,94059-6,454610,0154-10,0786-6,245830,0155-10,2067-6,030890,0156-10,3248-5,809990,0157-10,4326-5,583360,0158-10,5302-5,351220,0159-10,6173-5,11380,016-10,694-4,871340,0161-10,7601-4,624060,0162-10,8156-4,372220,0163-10,8604-4,116070,0164-10,8946-3,855860,0165-10,9179-2,591830,0166-10,9305-3,324270,0167-10,9323-3,053420,0168-10,9233-2,779570,0169-10,9035-2,50296

Окончание табл.1230,017-10,873-2,223890,0171-10,8317-1,942620,0172-10,7798-1,659440,0173-10,7172-1,374620,0174-10,644-1,088440,0175-10,5604-0,801190,0176-10,4663-0,513140,0177-10,3619-0,22460,0178-10,24720,0641760,0179-10,12250,3528840,018-9,98770,6412430,0181-9,84310,928970,0182-9,688781,215780,0183-9,524891,501390,0184-9,351611,7855190,0185-9,16912,0678850,0186-8,977542,3482110,0187-8,777122,6262190,0188-8,568042,9016350,0189-8,350513,1741880,019-8,124733,4436090,0191-7,890933,709630,0192-7,649353,9719910,0193-7,400224,2304330,0194-7,143784,4846990,0195-6,88034,7345390,0196-6,610024,9797070,0197-6,333225,219960,0198-6,050175,4550620,0199-5,761165,6847810,02-5,466455,908889

Табл.2СигналыПроизведение действующих значений I₁I₂Квазимощность Угол, град, Проверка угла, градПогрешность, %a(t_j)=i₁(t_j)=10,9329sin(ωt_j-30°);50,417747,212669,461700,77b(t_j)=i₂(t_j)=9,1926sin(ωt_j+40°)

Способ может быть использован для определения фазового угла между синусоидальными сигналами одинаковой частоты в однофазной цепи переменного тока при диагностике работоспособности электротехнических и электромеханических систем и устройств. Сигналы оцифровывают для одних и тех же моментов времени. Используют сумму текущего и предыдущего отсчетов первого сигнала и разность текущего и предыдущего отсчетов второго сигнала для определения реактивной квазимощности на интервале времени, равном периоду сигналов. Определяют действующие значения сигналов. Фазовый угол между сигналами определяют по отношению реактивной квазимощности к произведению действующих значений сигналов. По сравнению с аналогичными решениями данный способ не требует фильтрации постоянной составляющей сигналов и сдвига сигналов на угол 90° и позволяет достичь средней погрешности определения сдвига фаз 0,77%. 3 ил., 2 табл.

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

t_j=t₁, t₂, ..., t_N,

где N - число разбиений на периоде Т,

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

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