СПОСОБ ИЗМЕРЕНИЯ ВЕКТОРА ГАРМОНИЧЕСКОГО СИГНАЛА Российский патент 2009 года по МПК G01R17/00 

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

Известен способ измерения параметров многоэлементных двухполюсников путем уравновешивания разветвленной мостовой измерительной цепи при одновременном действии нескольких тестовых воздействий с разными частотами, разделяемых с помощью аналоговых фильтров [Шеремет Л.П. Принципы построения мостовых измерительных цепей для одновременного уравновешивания на нескольких частотах. //Проблемы технической электродинамики, вып.54, Киев: Наукова думка, 1975. - С.14-19].

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

Известен также способ измерения основной (первой) гармоники периодического сигнала, принятый автором за прототип, согласно которому измерение синусной и косинусной составляющих этой гармоники, т.е. ее проекций на два ортогональных, т.е. сдвинутых относительно друг друга по фазе на 90° опорных сигнала, осуществляют в реальном масштабе времени путем неравномерной дискретизации данного сигнала на множестве моментов времени, являющихся линейными комбинациями полупериодов его гармоник с номерами, равными простым числам: 1, 2, 3, 5, 7, 11,…, и суммирования множества полученных на этом множестве времени дискрет [Пат. Р.246351, Польша, МКИ G01R, Sposob i uklad do pomiaru wektora harmonicznej podstawowej przebiegu odksztalconego / Sawicki J. - Warszawa: Biuletyn Urzedu Patentowego. - 1985. - Nr 19 (307). - P.39].

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

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

Поставленная цель достигается тем, что в способе измерения вектора гармонического сигнала, действующего вместе с М-1 другими гармоническими сигналами и имеющего, как и они, известный период T_j и неизвестные амплитуду A_j и начальный фазовый сдвиг φ_0j, при котором проекции и измеряемого сигнала

S_j(t)=A_jsin(2πt/T_j+φ_0j)

на два ортогональных вектора опорных сигналов, связанные с A_j и φ_0j, например, соотношениями

и ,

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

с помощью мгновенных импульсов, действующих в моменты времени, образующие множества {t_i ^'} и {t_i ^''=t_i ^'+T_j/4}_j, причем {t_i ^'}_j формируют так, чтоб

и

где K_j - коэффициент, зависящий от процедуры формирования

{t_i'}_j, a N - число t_i ^', в качестве измеряемого выбирают любой из совместно действующих сигналов, частота следования которого 1/T_j не кратна частотам других сигналов, а множество {t_i ^'}_j моментов времени действия мгновенных импульсов выборки дискрет суммарного сигнала σ(t) формируют с помощью пошаговой процедуры из начального множества, состоящего из одного произвольного начального момента t₀, путем получения на первом шаге дополнительного множества, состоящего также из одного момента, посредством сдвига начального множества на нечетное число полупериодов первого сигнала, затем получения на втором шаге дополнительного множества, состоящего уже из двух моментов, посредством сдвига двух моментов суммарного множества, полученного на первом шаге, на нечетное число полупериодов второго сигнала и продолжают эту процедуру с образованием на каждом m-м шаге дополнительного множества посредством сдвига суммарного множества, сформированного на предыдущем шаге, на нечетное число n_m полупериодов очередного подавляемого сигнала, и удвоением числа моментов выборки дискрет сигнала σ(t), до тех пор, пока число шагов не станет равным М-1, а амплитуду A_j и фазовый сдвиг φ_0j сигнала S_j(t) определяют по соотношениям

при

где t_i ^'=t_i+t₀, t₁=0, t₂=t₁+n₁T₁/2,

t₃=t₁+n₂T₂/2, t₄=t₂+n₂T₂/2, t₅=t₁+n₃T₃/2, t₆=t₂+n₃T₃/2, t₇=t₃+n₃T₃/2, t₈=t₄+n₃T₃/2,……, m≠j.

Изобретение поясняется графическим материалом: на фиг.1 отражена процедура формирования множества моментов времени {t_i'}_j выборки дискрет суммарного сигнала

а фиг.2 поясняет методику получения математического выражения для K_j.

Сущность изобретения состоит в том, что предлагаемая процедура формирования множества моментов времени {t_i'}_j позволяет точно и быстро (в реальном масштабе времени) измерить проекции на ортогональные векторы опорных сигналов любого из входящих в σ(t) сигналов, в частности, проекции p^' _j и р^'' _j сигнала S_j(t), если его частота не кратна частотам других сигналов, инвариантно по отношению к остальным, т.е. исключить влияние всех этих сигналов на точность измерения p^' _j и р^'' _j, а значит и A_j и

φ_oj.

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

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

Поясним это математически. Пусть имеется исходное множество {t_i} моментов времени выборки дискрет (подавляемого) гармонического сигнала, сформированное на (m-1)-м шаге. И пусть образовано дополнительное множество {t_i*} путем сдвига исходного на интервал времени, равный n_mT_m/2, где n_m - нечетное число. Покажем, что сумма дискрет гармонического сигнала с периодом Т_m, полученных на суммарном множестве моментов времени, состоящем из исходного и дополнительного множеств, равна нулю. Для этого пронумеруем в произвольном порядке моменты времени t_i исходного множества. Тогда моментами времени дополнительного множества будут t_i*=t_i+n_mT_m/2. Но сдвиг на интервал времени n_mТ_m/2, если n_m - нечетное число, эквивалентен сдвигу по фазе на π, а sin(φ+π)=-sinφ и значит суммы всех пар дискрет подавляемого сигнала с одинаковыми номерами моментов времени t_i и t_i* будут равны нулю, что и требовалось доказать.

Заметим, что равенство сдвига нечетному числу полупериодов, а не одному, нужна для исключения возможного совпадения моментов времени выборки дискрет при сдвигах множеств {t_i} и {t_i*}, например, в случае Т₃=Т₂+Т₁ и ему подобных исключениях.

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

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

Для этого рассмотрим множество, сформированное на очередном шаге, в частности, на следующем шаге по отношению к описанному выше. Оно также состоит из двух одинаковых множеств, основного (исходного) и дополнительного, сдвинутых во времени на n_mТ_m+1/2, и благодаря этому позволяет подавить сигнал с периодом

T_m+1. Но из приведенного выше доказательства следовало, что первое из этих множеств (исходное) позволяет подавить сигнал с периодом Т_m. Сдвиг же этого множества на произвольный отрезок времени Δt не влияет на его свойство подавления этого сигнала, так как sin(φ+2πΔt/T_m+π)=-sin(φ+2πΔt/T_m). Таким образом, дополнительное множество, как и исходное, тоже обладает свойством подавления сигнала с периодом

Т_m, а значит и суммарное множество будет обладать этим свойством. Отсюда получаем, что суммарное множество будет обладать свойством подавления сигналов с периодами и Т_m, и T_m+1.

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

с периодом Т₁, второе - два с периодами

Т₁ и Т₂, m-е - m с периодами Т₁, Т₂,…, Т_m, а окончательное - (М-1) сигналов соответственно с периодами Т₁, Т₂,…, Т_m, Т_M, где m≠j, a j - номер измеряемого сигнала, что и требовалось доказать.

Проиллюстрируем эту процедуру графически с помощью фиг.1, соответствующей случаю М=5, где на верхней строке показано исходное (начальное) множество, состоящее из одного изображенного стрелкой импульса выборки дискреты, инициирующего формирование всего множества дискрет и соответствующего произвольному начальному моменту времени t₀. Таким образом, отсчет времени t осуществляется относительно t₀, т.е. t_i'=t_i+t₀. На второй строке показано первое множество, состоящее из двух импульсов выборки дискрет, соответствующих множеству из двух моментов времени, отстоящих друг от друга на интервал τ₁=n₁Т₁/2, и значит позволяющее подавить сигнал с периодом Т₁.

На третьей строке показано второе множество, состоящее из двух подмножеств, а именно изображенного на второй строке и сдвинутого относительно него на τ₂=n₂T₂/2 того же самого множества. Согласно вышеизложенному это множество содержит четыре дискреты, соответствующие четырем моментам времени, т.е. вдвое большему числу моментов, чем у предыдущего, и позволяет подавить сигналы с периодами Т₁ и

Т₂. Подтвердим это численно подстановкой имеющих место значений моментов времени в математические выражения сигналов S₁(t)=A₁sin(2πt/T₁+φ₀₁) и S₂(t)=A₂sin(2πt/T₂+φ₀₂) и суммированием полученных дискрет, приняв для упрощения выкладок n_i=1

На четвертой строке показано третье множество, состоящее из двух изображенных на третьей строке множеств, сдвинутых относительно друг друга на τ₃=n₃Т₃/2, и позволяющее подавить сигналы с периодами Т₁, T₂ и T₃. Убедиться в этом легко путем подстановки моментов времени t_i в выражения для сигналов S₁(t), S₂(t) и S₃(t) аналогично тому, как это было сделано выше для третьей строки.

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

τ₄=n₄Т₄/2, и позволяющее подавить четыре сигнала с периодами Т₁, T₂, Т₃ и Т₄. Пары формируемых на очередных шагах множеств отмечены штриховыми линиями над стрелками.

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

при будет составлять пару его проекций (p^' _j и p^'' _j) на оси ортогональной системы координат, так как в силу периодичности синусоиды значения входящих в данные суммы дискрет приводятся к одному периоду S_j(t), a Δt_j=Т_j/4 означает фазовый сдвиг сигналов на 90°. Найдем теперь формулу для коэффициента K_j.

Чтобы найти математическое выражение для К_j, воспользуемся векторной моделью гармонических сигналов, а для упрощения выкладок преобразуем форму представления сигнала S_j(t)=A_jsin(2πt/T_j+φ_0j) к удобному для анализа виду с учетом инвариантности K_j по отношению к амплитуде A_j и фазовому сдвигу φ_0j сигнала S_j(t), положив A_j=1 и φ_0j=0

S_j(t)⇒s(t)=sin2πt/T_j=sinω_jt,

а поскольку процедура формирования {t_i}_j одинакова для всех сигналов, опустим индекса и запишем {t_i}_j в виде {t_i},

где

а также положим

K_j=К, T_j=Т и ω_j=ω.

Кроме этого учтем, что, как отмечалось выше, благодаря периодичности синусоиды множество {t_i} преобразуется в , где - дробная часть частного от деления t_i на T=2π/ω, т.е. моменты времени выборки дискрет сигнала S_j(t) приводятся к одному периоду синусоиды S_j(t). В результате множество дискрет {s_i=s(t_i)} преобразуется в множество дискрет , приведенных к одному периоду синусоиды s(t), а множество фазовых углов {φ(t_i)=ωt_i}, соответствующих дискретам S_j(t_i), преобразуется в множество , приведенных к четырем квадрантам декартовой системы координат. При этом каждому из фазовых углов соответствует дискрета и наоборот.

Воспользуемся теперь формулой Эйлера, имеющей в данном случае вид: e^jωt=cosωt+jsinωt, согласно которой комплексное число выступает как единичный вектор в полярной системе координат, а дискрета - как проекция единичного вектора на мнимую ось комплексной плоскости .

Связанные с применением этой формулы построения приведены на фиг.2, на которой показано множество единичных векторов , где ω=2π/T_j, с действительными и мнимыми составляющими (проекциями на оси координат) s*_ix, s*_iy и их сумма , полученная по правилам векторной алгебры. Поворот на любой угол не приводит к изменению его модуля

где N=2^M-1, т.е. (в отличие от проекций

S_Σx и S_Σy) он является инвариантом по отношению к

Далее, поскольку математическое выражение для получено для случая единичного вектора, т.е. для синусоиды с единичной амплитудой (A_j=1), то выражение

где N=2^M-1 выступает в качестве нормирующего множителя.

Учитывая произведенные выше упрощения обозначений величин, а именно ω_j=ω, получим

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

Изобретение относится к области электроизмерительной техники. Способ может быть применен в средствах измерений пассивных и активных комплексных величин переменного тока, например в мостах переменного тока для измерения параметров многоэлементных двухполюсников, путем уравновешивания разветвленной мостовой измерительной цепи одновременно при нескольких тестовых воздействиях с разными частотами, а также при измерениях в случае действия нескольких гармонических помех. Способ состоит в том, что проекции и измеряемого сигнала S_j(t)=A_jsin(2πt/T_j+φ_0j) на два ортогональных опорных сигнала получают путем сложения дискрет суммы совместно действующих сигналов σ(t), моменты выборки которых составляют соответственно множества {t_i'}_j и {T_i ^''=t_i ^'+T_j/4}, которые формируют с помощью пошаговой процедуры из начального множества, состоящего из одного произвольного начального момента времени t₀, путем получения на первом шаге дополнительного множества, состоящего тоже из одного момента, посредством сдвига начального множества на нечетное число полупериодов первого сигнала, затем получения на втором шаге дополнительного множества, состоящего уже из двух моментов, посредством сдвига двух моментов суммарного множества, полученного на первом шаге, на нечетное число полупериодов второго сигнала и продолжают эту процедуру с образованием дополнительного множества на каждом шаге посредством сдвига суммарного множества, сформированного на предыдущем шаге, на нечетное число полупериодов очередного подавляемого сигнала и с удвоением на каждом шаге числа моментов времени выборки дискрет суммы совместно действующих сигналов до тех пор, пока число шагов не станет равным М-1, а амплитуду A_j и фазовый сдвиг φ_0j сигнала S_j(t) определяют по соотношениям при где t_i ^'=t_i+t₀, t₁=0, t₂=t₁+n₁T₁/2, t₃=t₁+n₂T₂/2, t₄=t₂+n₂T₂/2, t₅=t₁+n₃T₃/2, t₆=t₂+n₃T₃/2, t₇=t₃+n₃T₃/2, t₈=t₄+n₃T₃/2,……, m≠j. Технический результат - обеспечение возможности точного измерения в реальном масштабе времени любого из совместно действующих гармонических сигналов с известными периодами, а также возможность одновременного измерения нескольких сигналов при условии, что частота измеряемого сигнала не кратна частотам действующих совместно с ним сигналов. 2 ил.

Способ измерения вектора гармонического сигнала, действующего вместе с М-1 другими гармоническими сигналами и имеющего, как и они, известный период Т_j и неизвестные амплитуду A_j, и начальный фазовый сдвиг φ_0j, при котором проекции и измеряемого сигнала
S_j(t)=A_jsin(2πt/T_j+φ_0j) на два ортогональных вектора опорных сигналов, связанные с A_j и φ_0j, например, соотношениями и измеряют путем выборки и суммирования дискрет суммарного сигнала с помощью мгновенных импульсов, действующих в моменты времени, образующие множества {t_i ^'} и {t_j ^''=t_i ^'+T_j/4}_j, причем {t_j ^'}_j, формируют так, чтоб и
где K_j - коэффициент, зависящий от процедуры формирования {t_i ^'}, а N - число t_i ^', отличающийся тем, что в качестве измеряемого выбирают любой из совместно действующих сигналов, частота следования которого 1/T_j не кратна частотам других сигналов, а множество {t_i ^'} моментов времени действия мгновенных импульсов выборки дискрет суммарного сигнала σ(t) формируют с помощью пошаговой процедуры из начального множества, состоящего из одного произвольного начального момента t₀, путем получения на первом шаге дополнительного множества, состоящего также из одного момента, посредством сдвига начального множества на нечетное число полупериодов первого сигнала, затем получения на втором шаге дополнительного множества, состоящего уже из двух моментов, посредством сдвига двух моментов суммарного множества, полученного на первом шаге, на нечетное число полупериодов второго сигнала и продолжают эту процедуру с образованием на каждом m-ом шаге дополнительного множества посредством сдвига суммарного множества, сформированного на предыдущем шаге, на нечетное число n_m полупериодов очередного подавляемого сигнала и с удвоением числа моментов выборки дискрет сигнала σ(_t) до тех пор, пока число шагов не станет равным М-1, а амплитуду A_j и фазовый сдвиг φ_0j сигнала S_j(t) определяют по соотношениям:
при
где t_i ^'=t_i+t₀, t₁=0, t₂=t₁+n₁T₁/2,
t₃=t₁+n₂T₂/2, t₄=t₂+n₂T₂/2, t₅=t₁+n₃T₃/2,
t₆=t₂+n₃T₃/2, t₇=t₃+n₃T₃/2, t₈=t₄+n₃T₃/2,……,
m≠j.