Способ совместного измерения частоты, амплитуды, фазы и начальной фазы гармонического сигнала
Изобретение относится к измерительной технике и может быть использовано для одновременного измерения частоты, амплитуды, фазы и начальной фазы непрерывного или импульсного гармонического сигнала по одному и тому же минимальному набору исходных данных. Способ включает в себя дискретизацию аналогового сигнала, представление его фрагмента тройкой цифровых кодов S1, S2, S3, сформированных в моменты времени t1, t2, t3 и используемых для вычисления частоты сигнала f по формуле f = 1 2 π τ arccos S 1 + S 3 2 S 2 , где τ - интервал дискретизации. При этом фрагмент сигнала и соответствующую ему тройку кодов выбирают так, чтобы код S2 не равнялся нулю. При изменении частоты те же коды используют для вычисления амплитуды U, фазы φ и начальной фазы сигнала φ0 в соответствии с выражениями:
U = 2 | S 2 | ( S 2 2 − S 1 S 3 ) [ 4 S 2 2 − ( S 1 + S 3 ) 2 ] , φ = arccos ( S 2 U ) и φ 0 = φ − 2 π f t 2 ,
где фаза φ соответствует моменту времени t2.
Техническим результатом изобретения является ускорение измерений.
Реферат
Предлагаемое изобретение относится к измерительной технике и может быть использовано для одновременного измерения частоты, амплитуды, фазы и начальной фазы непрерывного или импульсного гармонического сигнала по одному и тому же минимальному набору исходных данных.
Известен способ определения частоты сигнала (патент РФ на изобретение №2117306), заключающийся в том, что сигнал дискретизируют, вычисляют его спектр, определяют номер максимальной спектральной составляющей, измеряют ее амплитуду, а также номер и амплитуду большей из смежных с ней и используют эти исходные данные в формуле вычисления частоты.
В этом способе-аналоге, так же как и в предлагаемом изобретении, сигнал дискретизируют по времени и представляют дискретные значения сигнала цифровыми кодами. Кроме того, используемые для вычисления частоты исходные данные несут в себе информацию и о частоте, и об амплитуде сигнала.
Недостаток способа состоит в том, что в нем используется не непосредственно сигнал, а его спектр (т.е. выполняется дополнительное преобразование - Фурье-преобразование сигнала). Исходные данные (в данном случае амплитуды спектральных составляющих) используются в формуле вычисления частоты, и, следовательно, погрешности измерения амплитуд вносят дополнительный вклад в погрешности вычисления частоты. Кроме метрологических факторов, определяющих погрешности измерения амплитуд (таких, например, как разрядность цифрового кода), возможны еще и объективные факторы, связанные, например, с тем, что найденная максимальная спектральная составляющая может неточно соответствовать частоте сигнала, и потому ее измеренное значение будет занижено относительно истинного.
Известен также способ измерения амплитуды гармонических сигналов (патент РФ на изобретение №2060475), заключающийся в том, что в спектре сигнала регистрируют гармонику с максимальной амплитудой, определяют ее частоту и используют эту частоту в формуле вычисления амплитуды гармонических колебаний.
В этом способе-аналоге, как и в предлагаемом изобретении, для измерения амплитуды сигнала и его частоты используются одни и те же исходные данные (в этом аналоге это спектральная составляющая с максимальной амплитудой).
Недостаток этого способа-аналога, так же как и предыдущего, состоит в использовании для измерений не сигнала, а его Фурье-преобразования и использовании частоты спектральной составляющей с максимальной амплитудой в формуле вычисления амплитуды. Из-за этого погрешности измерения частоты, связанные, например, с тем, что частота спектральной составляющей с максимальной амплитудой может не совпадать с истинной частотой сигнала, вносят дополнительный вклад в погрешности вычисления амплитуды.
Наиболее близким к заявляемому способу по технической сути представляется способ, реализованный в "Устройстве определения частоты переменного напряжения" (Авторское свидетельство СССР №1185260. Опубл. 30.10.88. Бюл. №40) и включающий в себя дискретизацию аналогового сигнала, представление его фрагмента тройкой цифровых кодов S1, S2, S3, сформированных в моменты времени t1, t2, t3 соответственно и используемых для вычисления частоты сигнала f по формуле f = 1 2 π τ arccos S 1 + S 3 2 S 2 , где τ - интервал дискретизации; при этом фрагмент сигнала и соответствующую ему тройку кодов выбирают так, чтобы код S2 не равнялся нулю.
Признаки, совпадающие с признаками предлагаемого изобретения: дискретизация аналогового сигнала, представление его фрагмента тройкой цифровых кодов S1, S2, S3, сформированных в моменты времени t1, t2, t3 соответственно; использование кодов для вычисления частоты сигнала f по формуле f = 1 2 π τ arccos S 1 + S 3 2 S 2 , где τ - интервал дискретизации; выбор такого фрагмента сигнала и соответствующей ему тройки кодов, при котором код S2 не равнялся бы нулю.
Недостаток этого способа, выбранного в качестве прототипа, состоит в недоиспользовании информации о сигнале, содержащейся в сформированных цифровых кодах S1, S2, S3. Между тем, зная коды S1, S2 и S3, можно не только измерять частоту, но и дополнительно к найденной частоте измерять амплитуду, фазу и начальную фазу непрерывного или импульсного гармонического сигнала. Кроме того, измерение амплитуды, фазы и начальной фазы можно выполнять одновременно (параллельно) с измерением частоты. Это позволит ускорить измерения.
Задачей, на решение которой направлено предлагаемое изобретение, является ускорение измерений за счет одновременного использования сформированных цифровых кодов S1, S2, S3 и для измерения частоты и для измерения амплитуды, фазы и начальной фазы анализируемого непрерывного или импульсного гармонического сигнала.
Этот технический результат достигается тем, что имеющиеся цифровые коды S1, S2, S3 одновременно с измерением частоты используют для вычисления амплитуды U, фазы φ и начальной фазы сигнала φ0 в соответствии с выражениями: U = 2 | S 2 | ( S 2 2 − S 1 S 3 [ 4 S 2 2 − ( S 1 + S 3 ) 2 ] , ϕ = arccos ( S 2 U ) и φ0=φ-2πft2, где фаза φ соответствует моменту времени t2.
Для достижения технического результата в способе совместного измерения частоты, амплитуды, фазы и начальной фазы гармонического сигнала, включающем дискретизацию аналогового сигнала, представление его фрагмента тройкой цифровых кодов S1, S2, S3, сформированных в моменты времени t1, t2, t3 соответственно и используемых для вычисления частоты сигнала f по формуле f = 1 2 π τ arccos S 1 + S 3 2 S 2 , где τ - интервал дискретизации; при этом фрагмент сигнала и соответствующую ему тройку кодов выбирают так, чтобы код S2 не равнялся нулю, одновременно с измерением частоты, те же коды используют для вычисления амплитуды U, фазы φ и начальной фазы сигнала φ0 в соответствии с выражениями:
U = 2 | S 2 | ( S 2 2 − S 1 S 3 ) [ 4 S 2 2 − ( S 1 + S 3 ) 2 ] , ϕ = arccos ( S 2 U ) и φ0=φ-2πft2, где фаза φ соответствует моменту времени t2.
Сравнение предлагаемого способа с прототипом показывает, что он содержит новые признаки, т.е. соответствует критерию новизны. Из сравнения с аналогами следует, что заявляемый способ соответствует критерию «существенные отличия», так как в аналогах не обнаружены новые заявляемые признаки.
Для доказательства существования причинно-следственной связи между заявляемыми признаками и достигаемым техническим результатом рассмотрим сущность предлагаемого способа совместного измерения частоты, амплитуды, фазы и начальной фазы гармонического сигнала и сопоставим его со способом-прототипом и способами-аналогами.
Гармонический сигнал
S ( t ) = U cos ( 2 π f t + ϕ ) ( 1 )
считается полностью известным, если известны амплитуда U, частота f и начальная фаза φ сигнала. Этот сигнал считается полностью известным потому, что, зная U, f и φ, можно вычислить (и это прямая задача) значение сигнала в любой момент времени t.
С другой стороны, между сигналом S(t) и его параметрами U, f и φ существует взаимно-однозначное соответствие. Поэтому, если сигнал S(t) известен (например, задан последовательностью дискретных отсчетов, полученных с интервалом времени τ), то параметры U, f и φ могут быть определены. Для их определения (обратная задача) достаточно составить и решить относительно U, f и φ три уравнения вида:
S ( t j ) = U cos ( 2 π f t j + ϕ ) ( 2 )
с тремя неизвестными. Покажем это.
Обозначим:
t j = τ j ; S ( t j ) = S j ; 2 π f t j = 2 π f τ j = α j ; г д е α = 2 π f τ . ( 3 )
С учетом этих обозначений уравнение (2) примет вид:
S j = U cos ( α j + ϕ ) ( 4 )
Составим 3 уравнения вида (4) для -1≤j≤1:
S − 1 = U cos ( − α + ϕ ) ( 5 )
S 0 = U cos ϕ ( 6 )
S 1 = U cos ( α + ϕ ) ( 7 )
где S-1, S0, S1 - отсчеты сигнала.
Раскроем уравнения (5) и (7):
S − 1 = U cos ( − α ) cos ϕ − U sin ( − α ) sin ϕ , ( 8 )
S 1 = U cos α cos ϕ − U sin α sin ϕ . ( 9 )
Учитывая четность функции cos(x) и нечетность функции sin(x), сложим уравнения (8) и (9):
S − 1 + S 1 = 2 U cos α cos ϕ ( 10 )
Подставив в это уравнение вместо Ucosφ величину S0 из уравнения (6) и решив полученное уравнение относительно cosα, получим:
cos α = ( S − 1 + S 1 ) / ( 2 S 0 ) . ( 11 )
Угол α и частота f определятся теперь по формулам:
α = arccos ( ( S − 1 + S 1 ) / ( 2 S 0 ) ) , ( 12 )
f = α / ( 2 π τ ) = 1 2 π τ arccos S − 1 + S 1 2 S 0 ( 13 )
U и φ могут быть найдены из уравнений (5), (6), (7) и (9) методом подстановки. Найдем амплитуду U. Для этого из уравнения (6) получим cosφ:
cos ϕ = S 0 / U ( 14 )
и подставим его в (9). Имея в виду, что sin2x=1-cos2x, возведем уравнение (9) в квадрат и после преобразований получим:
U = S 0 2 + S 1 2 − 2 S 0 S 1 cos α sin α ( 15 )
или с учетом (11) и формулы sin2x=1-cos2x:
U = 2 | S 0 | ( S 0 2 + S − 1 S 1 ) [ 4 S 0 2 − ( S − 1 + S 1 ) 2 ] ( 16 )
Теперь найдем начальную фазу φ из уравнения (14):
ϕ = arccos ( S 0 U ) ( 17 )
Фаза, вычисленная по формуле (17), соответствует моменту получения отсчета S0. Зная положение отсчета S0 на оси времени, можно вычислить набег фазы сигнала относительно нулевого момента времени. С его помощью можно пересчитать φ к нулевому моменту времени, т.е. вычислить начальную фазу сигнала φ0.
Учитывая, что от замены индексов суть вышеизложенного не изменится, заменим в формулах (13), (16) и (17): S-1 на S1; S0 на S2 и S1 на S3 и получим формулы, приведенные в формуле предлагаемого изобретения. Кроме того, учтем, что набег фазы, о котором говорилось в предыдущем абзаце, может быть в контексте формулы изобретения представлен выражением 2πft2.
Внедрение заявляемого способа позволит по трем отсчетам S1, S2, S3 сигнала без дополнительных затрат, преобразований и действий одновременно вычислить и частоту и остальные параметры непрерывного или импульсного гармонического сигнала (амплитуду, фазу и начальную фазу). Одновременное выполнение указанных измерений позволит уменьшить общее время измерений, т.е. ускорить измерения.
Способ совместного измерения частоты, амплитуды, фазы и начальной фазы гармонического сигнала, включающий дискретизацию аналогового сигнала, представление его фрагмента тройкой цифровых кодов S1, S2, S3, сформированных в моменты времени t1, t2, t3 и используемых для вычисления частоты сигнала f по формуле f = 1 2 π τ arccos S 1 + S 3 2 S 2 , где τ - интервал дискретизации; при этом фрагмент сигнала и соответствующую ему тройку кодов выбирают так, чтобы код S2 не равнялся нулю, отличающийся тем, что одновременно с измерением частоты те же коды используют для вычисления амплитуды U, фазы φ и начальной фазы сигнала φ0 в соответствии с выражениями: U = 2 | S 2 | ( S 2 2 − S 1 S 3 ) [ 4 S 2 2 − ( S 1 + S 3 ) 2 ] , ϕ = arccos ( S 2 U ) и φ0=φ-2πft2,где фаза φ соответствует моменту времени t2.