Способ определения угла сдвига фаз между синусоидальными сигналами (варианты)



Способ определения угла сдвига фаз между синусоидальными сигналами (варианты)
Способ определения угла сдвига фаз между синусоидальными сигналами (варианты)
Способ определения угла сдвига фаз между синусоидальными сигналами (варианты)
Способ определения угла сдвига фаз между синусоидальными сигналами (варианты)
Способ определения угла сдвига фаз между синусоидальными сигналами (варианты)
Способ определения угла сдвига фаз между синусоидальными сигналами (варианты)
Способ определения угла сдвига фаз между синусоидальными сигналами (варианты)
Способ определения угла сдвига фаз между синусоидальными сигналами (варианты)
Способ определения угла сдвига фаз между синусоидальными сигналами (варианты)
Способ определения угла сдвига фаз между синусоидальными сигналами (варианты)
Способ определения угла сдвига фаз между синусоидальными сигналами (варианты)
Способ определения угла сдвига фаз между синусоидальными сигналами (варианты)
Способ определения угла сдвига фаз между синусоидальными сигналами (варианты)
Способ определения угла сдвига фаз между синусоидальными сигналами (варианты)
Способ определения угла сдвига фаз между синусоидальными сигналами (варианты)
Способ определения угла сдвига фаз между синусоидальными сигналами (варианты)
Способ определения угла сдвига фаз между синусоидальными сигналами (варианты)
Способ определения угла сдвига фаз между синусоидальными сигналами (варианты)
Способ определения угла сдвига фаз между синусоидальными сигналами (варианты)
Способ определения угла сдвига фаз между синусоидальными сигналами (варианты)
Способ определения угла сдвига фаз между синусоидальными сигналами (варианты)
Способ определения угла сдвига фаз между синусоидальными сигналами (варианты)
Способ определения угла сдвига фаз между синусоидальными сигналами (варианты)

 


Владельцы патента RU 2563556:

Федеральное государственное бюджетное образовательное учреждение высшего образования "Иркутский национальный исследовательский технический университет" (ФГБОУ ВО "ИРНИТУ") (RU)

Изобретение относится к электротехнике, в частности к электрооборудованию, установленному на электрических станциях и подстанциях в системах производства, передачи и потребления электроэнергии, и может быть использовано во всех электроустановках, использующих цифровую обработку данных. Способ определения угла сдвига фаз между двумя синусоидальными сигналами путем измерения N/2 раз в течение полупериода Т/2 и в каждый текущий момент времени tj, j=1, 2, …, N/2 мгновенного значения одного из двух синусоидальных сигналов a(tj), изменяющегося во времени t по следующей зависимости: a(t)=A m sin(ωt). При наступлении момента выполнения условия, при котором мгновенное значение a(tj)=0, осуществления измерения и фиксации мгновенного значения другого синусоидального сигнала - b(tj)|а=0 той же частоты, изменяющегося во времени t по следующей зависимости: b(t)=B m sin(ωt+φ). Определяют значение угла сдвига фаз φ: где φ - угол сдвига фаз между двумя синусоидальными сигналами a(t) и b(t); b(tj)|а=0 - значение синусоидального сигнала b(t) в течение одного полупериода Т/2 в момент времени tj, когда значение синусоидального сигнала a(t) равно нулю, единицы измерения сигнала b(t); Вm - амплитудное значение синусоидального сигнала b(t), единицы измерения сигнала b(t), взятое со знаком плюс, если выполняется условие где b(tj-1) - предыдущее мгновенное значение синусоидального сигнала b(t), измеренное в момент времени tj-1, и со знаком минус, если Технический результат заключается в повышении быстродействия и точности определения сдвига фаз. 4 н.п. ф-лы, 3 табл.

 

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

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

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

Заявляемое изобретение относится к приоритетному направлению развития науки и технологий «Технологии создания энергосберегающих систем транспортировки, распределения и потребления тепла и электроэнергии» [Алфавитно-предметный указатель к Международной патентной классификации по приоритетным направлениям развития науки и технологий / Ю.Г. Смирнов, Е.В. Скиданова, С.А. Краснов. - М.: ПАТЕНТ, 2008. - с.97], так как позволяет определить угол сдвига фаз между двумя синусоидальными сигналами, необходимые для мгновенного и суммарного учета электроэнергии, для регулирования режимов работы компенсаторов реактивной энергии в линиях электропередачи, для управления электроэнергетической системой.

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

К прямому способу, основанному на непосредственном измерении угла сдвига фаз между двумя синусоидальными сигналами по наблюдениям мгновенных значений тока и напряжения можно отнести, например, способ, рассмотренный в следующем источнике [Овчаренко Н.И. Цифровые аппаратные и программные элементы микропроцессорной релейной защиты и автоматики энергосистем. М.: 2006, стр. 61-62, 67].

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

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

К косвенному способу, основанному на вычислениях активной и реактивной мощностей по результатам измерений мгновенных значений тока и напряжения и последующем определении угла сдвига фаз между двумя синусоидальными сигналами, можно отнести, например, способ, рассмотренный в следующем источнике [Патент РФ №2264631, МПК G01R 25/00, «Способ определения сдвига фаз между двумя синусоидальными сигналами», Гольдштейн Е.И., Сулайманов А.О., Бацева Н.Л., Панкратов А.В., опубликованному 20.11.2005].

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

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

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

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

- сократить временные затраты на определения угла сдвига фаз между двумя синусоидальными сигналами в эксплуатации;

- использовать значения угла сдвига фаз между двумя синусоидальными сигналами для текущего контроля устойчивости электроэнергетической системы;

- использовать значения угла сдвига фаз между двумя синусоидальными сигналами для мгновенного и суммарного учета электроэнергии;

- управлять текущим состоянием электроэнергетического комплекса.

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

- возможность постоянного мониторинга процесса изменения сдвига фаз электрических сигналов во времени для всех трех фаз передачи электроэнергии;

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

- повышение точности вычисления значения угла сдвига фаз между двумя синусоидальными сигналами и, как следствие, повышение качества управления электроэнергетическим объектом;

- возможность контроля соотношения активной и реактивной мощности в загрузке линий электропередачи, с целью принятия или непринятия решения об установке на ней компенсаторов реактивной мощности;

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

Технический результат достигается тем, что в первом варианте формулы изобретения раскрыта техническая сущность способа определения угла сдвига фаз φ между двумя синусоидальными сигналами путем измерения N/2 раз в течение полупериода T/2 и в каждый текущий момент времени tj, j=1, 2, …, N/2 мгновенного значения одного из двух синусоидальных сигналов a(tj), изменяющегося во времени t по следующей зависимости:

где a(t) - синусоидальный сигнал, единицы измерения сигнала a(t);

Am - амплитудное значение синусоидального сигнала a(t), единицы измерения сигнала a(t);

ω - круговая частота синусоидального сигнала, рад/с;

при наступлении момента выполнения условия, при котором мгновенное значение a(tj)=0, осуществления измерения и фиксации мгновенного значения другого синусоидального сигнала - b(tj)|a=0 той же частоты, изменяющегося во времени t по следующей зависимости:

где b(t) - синусоидальный сигнал, единицы измерения сигнала b(t);

Вm - амплитудное значение синусоидального сигнала b(t), единицы измерения сигнала b(t);

ω - круговая частота синусоидального сигнала, рад/с;

и определения значения угла сдвига фаз φ по следующему математическому выражению

где φ - угол сдвига фаз между двумя синусоидальными сигналами a(t) и b(t);

b(tj)|а=0, - значение синусоидального сигнала b(t) в течение одного полупериода Т/2 в момент времени tj, когда значение синусоидального сигнала a(t) равно нулю, единицы измерения сигнала b(t);

Вm - амплитудное значение синусоидального сигнала b(t), единицы измерения сигнала b(t), взятое со знаком плюс, если выполняется условие

где b(tj-1) - предыдущее мгновенное значение синусоидального сигнала b(t), измеренное в момент времени tj-1,

и со знаком минус, если

Технический результат достигается тем, что во втором варианте формулы изобретения раскрыта техническая сущность способа определения угла сдвига фаз φ между двумя синусоидальными сигналами путем измерения N/2 раз в течение полупериода Т/2 и в каждый текущий момент времени tj, j=1, 2, …, N/2 мгновенного значения одного из двух синусоидальных сигналов b(tj), изменяющегося во времени t по следующей зависимости:

где b(t) - синусоидальный сигнал, единицы измерения сигнала b(t);

Вm - амплитудное значение синусоидального сигнала b(t), единицы измерения сигнала b(t)\

ω - круговая частота синусоидального сигнала, рад/с;

φ - угол сдвига сигнала b(t), рад,

постоянной при каждом измерении проверки выполнения следующего условия для этого сигнала:

где b(tj-1) предыдущее мгновенное значение синусоидального сигнала, то есть измеренное в момент времени tj-1,

при наступлении указанного условия фиксируют текущие и предыдущие значения сигналов a(t) и b(t) и осуществляют вычисление значения - а|b≈0 сигнала a(tj) той же частоты, изменяющегося во времени t по следующей зависимости:

где a(t) - синусоидальный сигнал, единицы измерения сигнала a(t);

Am - амплитудное значение синусоидального сигнала a(t), единицы измерения сигнала a(t);

ω - круговая частота синусоидального сигнала, рад/с;

согласно математическому выражению

где запись abs(Z) обозначает абсолютное, то есть без учета знака, значение некоторой величины Z;

b(tj) - текущее мгновенное значение синусоидального сигнала b(t) в момент времени tj, единицы измерения сигнала b(t);

b(tj-1) - предыдущее мгновенное значение синусоидального сигнала b(t), измеренное в момент времени tj-1, единицы измерения сигнала b(t);

a(tj) - мгновенное значение синусоидального сигнала a(t) в момент времени tj, соответствующее измерению b(tj) синусоидального сигнала b(t), единицы измерения сигнала a(t);

a(tj-1) - предыдущее мгновенное значение синусоидального сигнала a(t) в момент времени tj-1, соответствующее измерению b(tj-1) синусоидального сигнала b(t), единицы измерения сигнала a(t),

а определение значения угла сдвига фаз φ между сигналами a(t) и b(t) осуществляют по следующему математическому выражению

где Аm - амплитудное значение синусоидального сигнала a(t), единицы измерения сигнала, взятое со знаком минус, если выполняется условие

где b(tj-1) - предыдущее мгновенное значение синусоидального сигнала b(t), измеренное в момент времени tj-1,

и со знаком плюс, если

Технический результат достигается тем, что в третьем варианте формулы изобретения раскрыта техническая сущность способа определения угла сдвига фаз φ между двумя синусоидальными сигналами путем измерения N раз в течение периода Т и в каждый текущий момент времени tj, j=1, 2, …, N мгновенного значения одного из двух синусоидальных сигналов a(tj), изменяющегося во времени t по следующей зависимости:

где a(t) - синусоидальный сигнал, единицы измерения сигнала a(t);

Am - амплитудное значение синусоидального сигнала a(t), единицы измерения сигнала a(t);

ω - круговая частота синусоидального сигнала, рад/с;

φа - угол сдвига сигнала a(t), рад,

постоянной проверки выполнения условия для синусоидального сигнала a(t)

где a(tj-1) - предыдущее мгновенное значение синусоидального сигнала a(t), измеренное в момент времени tj-1,

при наступлении указанного условия фиксируют текущие и предыдущие значения сигналов a(t) и b(t) и осуществляют вычисление значения - b|a≈0 сигнала b(t) той же частоты, изменяющегося во времени t по следующей зависимости:

где b(t) - синусоидальный сигнал, единицы измерения сигнала b(t);

Вm - амплитудное значение синусоидального сигнала b(t), единицы измерения сигнала b(t);

ω - круговая частота синусоидального сигнала, рад/с;

φb - угол сдвига сигнала b(t), рад,

согласно математическому выражению

где запись abs(Z) обозначает абсолютное, то есть без учета знака, значение некоторой величины Z;

a(tj) - мгновенное значение синусоидального сигнала a(t) в момент времени tj, единицы измерения сигнала a(t);

a(tj-1) - предыдущее мгновенное значение синусоидального сигнала a(t) в момент времени tj-1, единицы измерения сигнала a(t);

b(tj) - мгновенное значение синусоидального сигнала b(t) в момент времени tj, соответствующему измерению a(tj) синусоидального сигнала a(t), единицы измерения сигнала b(t);

b(tj-1) - предыдущее мгновенное значение синусоидального сигнала b(t) в момент времени tj-1, соответствующему измерению a(tj-1) синусоидального сигнала a(t), единицы измерения сигнала b(t),

а определение значения угла сдвига фаз φ между сигналами a(t) и b(t) осуществляют по следующему математическому выражению

где φа и φb - углы сдвига фаз сигналов a(t) и b(t), соответственно, а знак ± означает, что суммирование осуществляется по правилам электротехники, то есть из индуктивных углов вычитаются емкостные, если углы известны; b|а≈0, - значение синусоидального сигнала b(t) в момент времени tj, когда значение синусоидального сигнала a(t) изменило знак, то есть почти равно нулю, единицы измерения сигнала b(t);

Вm - амплитудное значение синусоидального сигнала b(t), единицы измерения сигнала b(t), взятое со знаком плюс, если выполняется условие

где b(tj-1) - предыдущее мгновенное значение синусоидального сигнала b(t), измеренное в момент времени tj-1,

и со знаком минус, если

Технический результат достигается тем, что в четвертом варианте формулы изобретения раскрыта техническая сущность способа определения угла сдвига фаз φ между двумя синусоидальными сигналами путем измерения N раз в течение периода T и в каждый текущий момент времени tj, j=1, 2, …, N мгновенного значения одного из двух синусоидальных сигналов b(tj), изменяющегося во времени t по следующей зависимости:

где b(t) - синусоидальный сигнал, единицы измерения сигнала b(t);

Вm - амплитудное значение синусоидального сигнала b(t), единицы измерения сигнала b(t);

ω - круговая частота синусоидального сигнала, рад/с;

φb - угол сдвига сигнала b(t), рад,

постоянной при каждом измерении проверки выполнения следующего условия для этого сигнала:

где b(tj-1) - предыдущее мгновенное значение синусоидального сигнала b(t), измеренное в момент времени tj-1,

при наступлении указанного условия фиксируют текущие и предыдущие значения сигналов a(t) и b(t) и осуществляют вычисление значения - а|b≈0 сигнала a(tj) той же частоты, изменяющегося во времени t по следующей зависимости:

где a(t) - синусоидальный сигнал, единицы измерения сигнала a(t);

Am - амплитудное значение синусоидального сигнала a(t), единицы измерения сигнала a(t);

ω - круговая частота синусоидального сигнала, рад/с;

φа - угол сдвига сигнала a(t), рад,

согласно математическому выражению

где запись abs(Z) обозначает абсолютное, то есть без учета знака, значение некоторой величины Z;

b(tj) - текущее мгновенное значение синусоидального сигнала b(t) в момент времени tj, единицы измерения сигнала b(t);

b(tj-1) - предыдущее мгновенное значение синусоидального сигнала b(t), измеренное в момент времени tj-1, единицы измерения сигнала b(t);

a(tj) - мгновенное значение синусоидального сигнала a(t) в момент времени tj, соответствующее измерению b(tj) синусоидального сигнала b(t), единицы измерения сигнала a(t);

a(tj-1) - предыдущее мгновенное значение синусоидального сигнала a(t) в момент времени tj-1, соответствующее измерению b(tj-1) синусоидального сигнала b(t), единицы измерения сигнала a(t),

а определение значения угла сдвига фаз φ между сигналами a(t) и b(t) осуществляют по следующему математическому выражению

где Am - амплитудное значение синусоидального сигнала a(t), единицы измерения сигнала, взятое со знаком минус, если выполняется условие

где a(tj-1) - предыдущее мгновенное значение синусоидального сигнала a(t), измеренное в момент времени tj-1,

и со знаком плюс, если

Предлагаемый способ определения угла сдвига фаз φ между изменяющимися во времени t синусоидальными сигналами a(t) и b(t) основан на измерении оцифрованных мгновенных значений этих сигналов a(tj), b(tj), для одних и тех же моментов времени tj, где j=1, 2, … N, N - число разбиений на периоде Т, Δt=T/N - шаг дискретизации сигнала по времени.

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

В реальных измерениях значений токов, напряжений и других параметров трехфазной цепи с помощью цифровых приборов напряжение одной из фаз (обычно фазы А) берется в качестве опорного, то есть его фаза равна нулю, а значения остальных углов сдвига фаз токов и напряжений рассчитывается относительно этого опорного напряжения. Такой подход для определения угла сдвига фаз φ между изменяющимися во времени t сигналами a(t)=Amsin(ωt) - опорным синусоидальным сигналом и b(t)=Bm·sin(ωt+φ) - другим синусоидальным сигналом одной и той же частоты, где Аm, Вm - амплитудные значения сигналов a(t) и b(t), соответственно, ω=2πf рад/с - круговая частота, а f=50 Гц - промышленная частота гармонического сигнала, позволяет получить наименьшую погрешность определения угла сдвига фаз в момент времени tj, когда опорный сигнал a(t) пересекает ось абсцисс. Это обусловлено тем, что при цифровой обработке значение N всегда кратно двум, обычно это цифра 2 в некоторой степени или сумма таких чисел, а процесс измерения синхронизирован с началом периода опорного сигнала. Поэтому точка пересечения опорным сигналом оси абсцисс всегда приходится точно на момент измерения, и значение другого, а также любого последующего синусоидального сигнала трехфазной цепи, будет точно измерено для последующих вычислений.

Итак, в момент времени tj, когда один сигнал a(t) пересекает ось абсцисс, то есть равен нулю, осуществляют измерение и фиксацию значения другого сигнала b(tj)|а=0. Учитывая установленный факт, что отношение полученного мгновенного значения сигнала к его амплитудному значению равно синусу угла сдвига фаз, то есть

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

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

Есть же использовать точку пересечения другим сигналом оси абсцисс в момент времени, близкий к tj, то она находится между двумя некоторыми значениями сигнала b(t), одно из которых b(tj), например, положительное и полученное в момент времени tj, и другое - предыдущее значение b(tj-1), которое отрицательное и полученное в предыдущий момент времени tj-1. Знаки поменялись, так как этот сигнал пересек ось абсцисс. И определять значение опорного синусоидального сигнала a(t) в примерной точке пересечения вторым сигналом оси абсцисс - а|b≈0, приходится из следующей пропорции с использованием двух значений второго сигнала b(tj), b(tj-1) и соответствующих значений первого сигнала a(tj) и a(tj-1)

Такой подход является источником дополнительной погрешности, однако ничего более точного пока предложить не удается, но и в этом случае погрешность предлагаемого изобретения остается меньше, чем при определении угла сдвига фаз другим способом, например, способами, указанными в ссылках выше по тексту. Значение опорного синусоидального сигнала a(t) в примерной точке пересечения вторым сигналом оси абсцисс а|b≈0, полученное из указанной пропорции, имеет вид

где запись abs(Z) обозначает абсолютное, то есть без учета знака, значение некоторой величины Z.

Используя факт, что отношение полученного мгновенного значения опорного сигнала а|b≈0 к его амплитудному значению Аm равно синусу угла сдвига фаз, то есть

найдем угол сдвига фаз

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

В варианте три и в варианте четыре формулы предлагаемого изобретения рассмотрены случаи определения угла сдвига фаз без использования опорного сигнала. То есть и a(t)=Amsin(ωt+φ1) - один изменяющийся во времени t синусоидальный сигнал, и b(t)=Bm·sin(ωt+φ2) - другой изменяющийся во времени t синусоидальный сигнал имеют свои углы сдвига фаз, φ1 и φ2, соответственно, которые, как правило, неизвестны. Иначе искомый угол сдвига фаз φ равен

где φ1 и φ2 - углы сдвига фаз сигналов a(t) и b(t), соответственно, а знак ± означает, что суммирование осуществляется по правилам электротехники, то есть из индуктивных углов вычитаются емкостные.

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

Замечание, касающееся уменьшения погрешности вычислений и повышения точности измерения угла сдвига фаз между двумя синусоидальными сигналами. Усреднение вычисленного значения некоторого параметра Rk, k=1, …, n, по выражению

где суммирование ведется по k=1, …, n,

R - среднее значение некоторого параметра, например, угла сдвига фаз между двумя синусоидальными сигналами,

n - количество вычислений этого параметра, по которым производится усреднение,

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

Пример использования способа

Для сравнения относительных погрешностей, предлагаемых вариантом три и вариантом четыре формулы изобретения, в качестве примера взяты массивы значений сигналов из патента РФ №2264631, МПК G01R 25/00, «Способ определения сдвига фаз между двумя синусоидальными сигналами», Гольдштейн Е.И., Сулайманов А.О., Бацева Н.Л., Панкратов А.В., опубликованного 20.11.2005 (далее по тексту - пример), которые представлены в табл. 1. Как и в примере, представлен только один период, включающий N=200, то есть двести измерений в течение времени Т=0,02 секунды, при частоте f=50 Гц.

В табл. 1 представлены из примера результаты измерения a(tj), b(tj) двух синусоидальных сигналов a(t), b(t) со следующими параметрами

где φ1=-30° - угол сдвига первого синусоидального сигнала относительно опорного, как правило, φ1 - неизвестен;

φ2=40° - угол сдвига второго синусоидального сигнала относительно опорного, как правило, φ2 - неизвестен,

которые занимают, соответственно, вторую и третью колонки таблицы. В первой колонке располагается время с шагом дискретизации Δt=T/N=0,0001 с.

Для сравнения относительных погрешностей, предлагаемых вариантом один и вариантом два формулы изобретения, приведены данные, аналогичные вышеприведенному примеру. Также представлен только один период, включающий N=200, то есть двести измерений в течение времени Т=0,02 секунды, при частоте f=50 Гц, которые приведены в табл. 2.

В табл. 2 представлены результаты измерения a(tj), b(tj) двух синусоидальных сигналов a(t), b(t) со следующими параметрами

где φ - неизвестный угол сдвига второго синусоидального сигнала относительно первого, в табл. 2 данные для b(t) получены с φ=40°, которые занимают, соответственно, вторую и третью колонки таблицы. В первой колонке располагается время с шагом дискретизации Δt=T/N=0,0001 с.

Для экономии места обе табл. 1 и табл. 2, каждая длиной двести строк и состоящая из трех колонок, «разорваны» пополам, затем сложены рядом и представлены в виде таблицы длиной сто строк из шести колонок.

Результаты расчетов сведены в табл. 3 для углов сдвига фаз между двумя синусоидальными сигналами, взятых с шагом 10°. Чистая активная нагрузка (угол сдвига фаз равен нулю) и чистая реактивная нагрузка (угол сдвига фаз равен 90°) не рассматривались, так как это идеальные случаи, не типичные в эксплуатации. Для каждого угла сдвига фаз между двумя синусоидальными сигналами вначале вычислялись таблицы, аналогичные табл. 1 и табл. 2, а затем по выражениям, приведенным в вариантах 1-4 предлагаемой формулы изобретения, определялся расчетный угол сдвига фаз.

Относительную погрешность измерения - δ в процентах определяли согласно [РМГ 29-99 Метрология. Основные термины и определения. М. 2001] умножением на 100% отношения abs(Δx) - абсолютной погрешности измерения к x - действительному значению измеряемой величины, то есть из выражения

δ=abs(Δx)/x 100%,

где abs(Δx) - абсолютная погрешность измерений;

x - действительное значение измеряемой величины.

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

Используя данные табл. 2, найдем угол сдвига фаз согласно первому варианту формулы предлагаемого изобретения. В табл. 2 есть три точки, в которых a(tj)=0, это при tj=0, tj=0,01 и tj=0,02. Во всех трех точках значение второго синусоидального сигнала равно b(tj)|а=0=5,908889381. В точке tj=0,01 значение второго синусоидального сигнала отрицательное, но в этот момент выполняется условие b(tj)<b(tj-1), то есть -5,9088894<-5,684781, и знак у амплитудного значения Вm тоже должен быть отрицательным, согласно варианту один формулы предлагаемого изобретения, и в последующих вычислениях эти знаки сократятся, так как берется их отношение. В результате получим

Переводя вычисленное значение из радиан в градусы, согласно выражению 180°·0,698132/n, получим 40°, что равно точному значению угла сдвига.

Аналогично, из табл. 1 видно, что, например, в точках tj-1=0,0016 и tj=0,0017 первый синусоидальный сигнал пересекает ось абсцисс, то есть его значения в этих точках поменяли знак: a(tj-1)=-0,2289614, a(tj)=0,11448697. Значения второго синусоидального сигнала в этих точках соответственно равны: b(tj-1)=8,57047978, b(tj)=8,67066861. Подставляя полученные значения гармонических сигналов в математические выражения варианта три формулы предлагаемого изобретения, также вычислим почти точные значения угла сдвига фаз.

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

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

1. Способ определения угла сдвига фаз между двумя синусоидальными сигналами путем измерения N/2 раз в течение полупериода Т/2 и в каждый текущий момент времени tj, j=1, 2, …, N/2 мгновенного значения одного из двух синусоидальных сигналов a(tj), изменяющегося во времени t по следующей зависимости:

где a(t) - синусоидальный сигнал, единицы измерения сигнала a(t);
Am - амплитудное значение синусоидального сигнала a(t), единицы измерения сигнала a(t);
ω - круговая частота синусоидального сигнала, рад/с;
при наступлении момента выполнения условия, при котором мгновенное значение a(tj)=0, осуществления измерения и фиксации мгновенного значения другого синусоидального сигнала - b(tj)|а=0 той же частоты, изменяющегося во времени t по следующей зависимости:

где b(t) - синусоидальный сигнал, единицы измерения сигнала b(t);
Вm - амплитудное значение синусоидального сигнала b(t), единицы измерения сигнала b(t);
ω - круговая частота синусоидального сигнала, рад/с;
и определения значения угла сдвига фаз φ по следующему математическому выражению

где φ - угол сдвига фаз между двумя синусоидальными сигналами a(t) и b(t);
b(tj)|а=0 - значение синусоидального сигнала b(t) в течение одного полупериода Т/2 в момент времени tj, когда значение синусоидального сигнала a(t) равно нулю, единицы измерения сигнала b(t);
Вm - амплитудное значение синусоидального сигнала b(t), единицы измерения сигнала b(t), взятое со знаком плюс, если выполняется условие

где b(tj-1) - предыдущее мгновенное значение синусоидального сигнала b(t), измеренное в момент времени tj-1,
и со знаком минус, если

2. Способ определения угла сдвига фаз между двумя синусоидальными сигналами путем измерения N/2 раз в течение полупериода Т/2 и в каждый текущий момент времени tj, j=1, 2, …, N/2 мгновенного значения одного из двух синусоидальных сигналов b(tj), изменяющегося во времени t по следующей зависимости:

где b(t) - синусоидальный сигнал, единицы измерения сигнала b(t);
Вm - амплитудное значение синусоидального сигнала b(t), единицы измерения сигнала b(t);
ω - круговая частота синусоидального сигнала, рад/с;
φ - угол сдвига сигнала b(t), рад,
постоянной при каждом измерении проверки выполнения следующего условия для этого сигнала:

где b(tj-1) - предыдущее мгновенное значение синусоидального сигнала, то есть измеренное в момент времени tj-1,
при наступлении указанного условия фиксируют текущие и предыдущие значения сигналов a(t) и b(t) и осуществляют вычисление значения - а|b≈0 сигнала a(tj) той же частоты, изменяющегося во времени t по следующей зависимости:

где a(t) - синусоидальный сигнал, единицы измерения сигнала a(t);
Am - амплитудное значение синусоидального сигнала a(t), единицы измерения сигнала a(t);
ω - круговая частота синусоидального сигнала, рад/с;
согласно математическому выражению

где запись abs(Z) обозначает абсолютное, то есть без учета знака, значение некоторой величины Z;
b(tj) - текущее мгновенное значение синусоидального сигнала b(t) в момент времени tj, единицы измерения сигнала b(t);
b(tj-1) - предыдущее мгновенное значение синусоидального сигнала b(t), измеренное в момент времени tj-1, единицы измерения сигнала b(t);
a(tj) - мгновенное значение синусоидального сигнала a(t) в момент времени tj, соответствующее измерению b(tj) синусоидального сигнала b(t), единицы измерения сигнала a(t);
a(tj-1) - предыдущее мгновенное значение синусоидального сигнала a(t) в момент времени tj-1, соответствующее измерению b(tj-1) синусоидального сигнала b(t), единицы измерения сигнала a(t),
а определение значения угла сдвига фаз φ между сигналами a(t) и b(t) осуществляют по следующему математическому выражению

где Аm - амплитудное значение синусоидального сигнала a(t), единицы измерения сигнала, взятое со знаком минус, если выполняется условие

где b(tj-1) - предыдущее мгновенное значение синусоидального сигнала b(t), измеренное в момент времени tj-1,
и со знаком плюс, если

3. Способ определения угла сдвига фаз между двумя синусоидальными сигналами путем измерения N раз в течение периода Т и в каждый текущий момент времени tj, j=1, 2, …, N мгновенного значения одного из двух синусоидальных сигналов a(tj), изменяющегося во времени t по следующей зависимости:

где a(t) - синусоидальный сигнал, единицы измерения сигнала a(t);
Am - амплитудное значение синусоидального сигнала a(t), единицы измерения сигнала a(t);
ω - круговая частота синусоидального сигнала, рад/с;
φa - угол сдвига сигнала a(t), рад,
постоянной проверки выполнения условия для синусоидального сигнала a(t)

где a(tj-1) - предыдущее мгновенное значение синусоидального сигнала a(t), измеренное в момент времени tj-1,
при наступлении указанного условия фиксируют текущие и предыдущие значения сигналов a(t) и b(t) и осуществляют вычисление значения - b|а≈0 сигнала b(t) той же частоты, изменяющегося во времени t по следующей зависимости:

где b(t) - синусоидальный сигнал, единицы измерения сигнала b(t);
Вm - амплитудное значение синусоидального сигнала b(t), единицы измерения сигнала b(t);
ω - круговая частота синусоидального сигнала, рад/с;
φb - угол сдвига сигнала b(t), рад,
согласно математическому выражению

где запись abs(Z) обозначает абсолютное, то есть без учета знака, значение некоторой величины Z;
a(tj) - мгновенное значение синусоидального сигнала a(t) в момент времени tj, единицы измерения сигнала a(t);
a(tj-1) - предыдущее мгновенное значение синусоидального сигнала a(t) в момент времени tj-1, единицы измерения сигнала a(t);
b(tj) - мгновенное значение синусоидального сигнала b(t) в момент времени tj, соответствующее измерению a(tj) синусоидального сигнала a(t), единицы измерения сигнала b(t);
b(tj-1) - предыдущее мгновенное значение синусоидального сигнала b(t) в момент времени tj-1, соответствующее измерению a(tj-1) синусоидального сигнала a(t), единицы измерения сигнала b(t),
а определение значения угла сдвига фаз φ между сигналами a(t) и b(t) осуществляют по следующему математическому выражению

где φа и φb - углы сдвига фаз сигналов a(t) и b(t), соответственно, а знак ± означает, что суммирование осуществляется по правилам электротехники, то есть из индуктивных углов вычитаются емкостные, если углы известны; b|а≈0 - значение синусоидального сигнала b(t) в момент времени tj, когда значение синусоидального сигнала a(t) изменило знак, то есть почти равно нулю, единицы измерения сигнала b(t);
Вm - амплитудное значение синусоидального сигнала b(t), единицы измерения сигнала b(t), взятое со знаком плюс, если выполняется условие

где b(tj-1) - предыдущее мгновенное значение синусоидального сигнала b(t), измеренное в момент времени tj-1,
и со знаком минус, если

4. Способ определения угла сдвига фаз между двумя синусоидальными сигналами путем измерения N раз в течение периода Т и в каждый текущий момент времени tj, j=1, 2, …, N мгновенного значения одного из двух синусоидальных сигналов b(tj), изменяющегося во времени t по следующей зависимости:

где b(t) - синусоидальный сигнал, единицы измерения сигнала b(t);
Вm - амплитудное значение синусоидального сигнала b(t), единицы измерения сигнала b(t);
ω - круговая частота синусоидального сигнала, рад/с;
φb - угол сдвига сигнала b(t), рад,
постоянной при каждом измерении проверки выполнения следующего условия для этого сигнала:

где b(tj-1) - предыдущее мгновенное значение синусоидального сигнала b(t), измеренное в момент времени tj-1,
при наступлении указанного условия фиксируют текущие и предыдущие значения сигналов a(t) и b(t) и осуществляют вычисление значения - а|b≈0 сигнала a(tj) той же частоты, изменяющегося во времени t по следующей зависимости:

где a(t) - синусоидальный сигнал, единицы измерения сигнала a(t);
Am - амплитудное значение синусоидального сигнала a(t), единицы измерения сигнала a(t);
ω - круговая частота синусоидального сигнала, рад/с;
φa - угол сдвига сигнала a(t), рад,
согласно математическому выражению

где запись abs(Z) обозначает абсолютное, то есть без учета знака, значение некоторой величины Z;
b(tj) - текущее мгновенное значение синусоидального сигнала b(t) в момент времени tj, единицы измерения сигнала b(t);
b(tj-1) - предыдущее мгновенное значение синусоидального сигнала b(t), измеренное в момент времени tj-1, единицы измерения сигнала b(t);
a(tj) - мгновенное значение синусоидального сигнала a(t) в момент времени tj, соответствующее измерению b(tj) синусоидального сигнала b(t), единицы измерения сигнала a(t);
a(tj-1) - предыдущее мгновенное значение синусоидального сигнала a(t) в момент времени tj-1, соответствующее измерению b(tj-1) синусоидального сигнала b(t), единицы измерения сигнала a(t),
а определение значения угла сдвига фаз φ между сигналами a(t) и b(t) осуществляют по следующему математическому выражению

где Am - амплитудное значение синусоидального сигнала a(t), единицы измерения сигнала, взятое со знаком минус, если выполняется условие

где a(tj-1) - предыдущее мгновенное значение синусоидального сигнала a(t), измеренное в момент времени tj-1,
и со знаком плюс, если



 

Похожие патенты:

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

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

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

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

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

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

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

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

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

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

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

Предлагаемое устройство относится к области радиоэлектроники и может быть использовано для определения несущей частоты и вида модуляции сигналов, принимаемых в заданном диапазоне частот. Технической задачей изобретения является расширение функциональных возможностей устройства путем распознавания сигналов с амплитудой и частотной манипуляцией. Устройство содержит приемную антенну 1, входную цепь 2, блок 3 поиска, гетеродин 5, смеситель 6, усилитель 7 промежуточной частоты, амплитудный детектор 8, 13, 28, 29 и 33, видеоусилитель 9, устройство 10 формирования частотной развертки, ЭЛТ 11, ключи 12, 24, 37, 40, 41 и 42, фильтры 14, 27 и 32 верхних частот, фильтры 15, 19 и 26 нижних частот, квадраторы 16 и 20, делители 17 и 22 напряжений, частотный детектор 18, блоки 23, 30, 34, 39 и 49 сравнения, фазовый детектор 25, интегратор 35, пороговый блок 36, измеритель 38 частоты, блок 43 памяти, преобразователи 44, 47 и 50 аналог-код, блок 46 клиппирования, анализаторы 21, 45 и 48 спектра, фазоинверторы 51 и 52, элементы совпадения 53, 54, 55 и 56. 2 ил.

Изобретение относится к радиотехнической области промышленности и может быть использовано при приеме нескольких совмещенных по времени разночастотных сигналов. Способ определения частоты в матричном приемнике, в котором ко входу j-й ступени приемника, имеющей Lj каналов, подключают устройство измерения частоты, измеряющее частоту сигнала в диапазоне рабочих частот j-й ступени, и сопоставляют номера сработавших индикаторов каналов ступени с измеренными значениями частоты. Устройство измерения частоты содержит усилитель-ограничитель, K каналов обработки и устройство обработки. Каждый канал содержит последовательно включенные полосовой фильтр, частотно-зависимое устройство и детектор. С выходов каналов сигнал подается на устройство обработки. Технический результат заключается в повышении вероятности однозначного определения частоты, исключении регистрации ложных значений частоты и пропуска сигналов. 2 н. и 4 з.п. ф-лы, 2 ил.

Изобретение относится к измерительной технике и может использоваться в информационно-измерительных устройствах для измерения частоты гармонических сигналов прецизионных кварцевых и квантовых стандартов частоты. Осуществляют аналого-цифровое преобразование измеряемого и опорного сигналов с интервалом временной дискретизации, определяемой частотой сигнала дискретизации, формируемого из опорного сигнала, запоминают полученные в результате аналого-цифровых преобразований цифровые выборки в следующих одна за другой тетрадах моментов времени, осуществляют преобразование цифровых выборок тетрад в значения фаз измеряемого и опорного сигналов и определяют искомую разность частот опорного и измеряемого сигналов. Устройство содержит последовательно соединенные генератор измеряемого сигнала, первый аналого-цифровой преобразователь, первое оперативное запоминающее устройство и процессор цифровой обработки сигналов, связанный шиной обмена данными с персональной вычислительной машиной. Вход синхронизации первого аналого-цифрового преобразователя соединен с выходом синтезатора частоты сигнала квантования, сигнальный вход которого соединен с выходом генератора опорного сигнала, а вход управления - с управляющим выходом процессора цифровой обработки сигналов. Устройство также содержит второе оперативное запоминающее устройство, выход которого соединен с вторым входом процессора цифровой обработки, и второй аналого-цифровой преобразователь, сигнальный вход которого соединен с выходом генератора опорного сигнала, вход синхронизации соединен с выходом синтезатора частоты сигнала квантования, а выход соединен с входом второго оперативного запоминающего устройства. Технический результат заключается в повышении точности измерения частоты гармонического сигнала при расширении диапазона частот сличаемых сигналов. 2 н.п. ф-лы, 2 ил.

Изобретение относится к области радиотехники и предназначено для использования в среднеорбитальном сегменте космической системы поиска и спасения терпящих бедствия судов, летательных аппаратов, отдельных людей или групп. Согласно способу измерения производятся с использованием всей длительности сигнала посылки радиобуя (440 мс), а не только по участку длительностью 160 мс - участку излучения чистой несущей частоты радиобуя, и соответственно всей энергии сигнала. Для этого производится модуляция принятых наземной станцией (станцией приема и обработки информации со среднеорбитальных ИСЗ систем «Глонасс», GPS и Gallileo) сигналов аварийных радиобуев достоверной цифровой информацией, заложенной в сигналы, передаваемые тем же самым аварийным радиобуем и выделенной из принятого сигнала в процессе его демодуляции и декодирования, взятой с обратным знаком (ремодуляция сигнала). Это преобразует весь принятый сигнал посылки этого радиобуя в немодулированную синусоиду, чем и обеспечивается получение минимально возможной ошибки измерения его частоты. Технический результат заявленного изобретения заключается в повышении точности измерений частоты сигналов радиобуев. 1 з.п. ф-лы, 4 ил.
Наверх