RU2017162C1 - Способ измерения частоты сигнала и устройство для его осуществления - Google Patents

Способ измерения частоты сигнала и устройство для его осуществления Download PDF

Info

Publication number
RU2017162C1
RU2017162C1 SU4928486A RU2017162C1 RU 2017162 C1 RU2017162 C1 RU 2017162C1 SU 4928486 A SU4928486 A SU 4928486A RU 2017162 C1 RU2017162 C1 RU 2017162C1
Authority
RU
Russia
Prior art keywords
input
signal
frequency
samples
output
Prior art date
Application number
Other languages
English (en)
Inventor
В.Г. Ципоренко
Е.Д. Ципоренко
Original Assignee
Житомирский научно-исследовательский институт радиосистем
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Житомирский научно-исследовательский институт радиосистем filed Critical Житомирский научно-исследовательский институт радиосистем
Priority to SU4928486 priority Critical patent/RU2017162C1/ru
Application granted granted Critical
Publication of RU2017162C1 publication Critical patent/RU2017162C1/ru

Links

Images

Landscapes

  • Measuring Frequencies, Analyzing Spectra (AREA)

Abstract

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

Description

Изобретение относится к радиотехнике и контрольно-измерительной технике и может быть использовано в радиотехнических и информационно-вычислительных системах, функционирующих в реальном масштабе времени.
Известен способ определения частоты сигнала по положению максимума амплитуды на оси частот в измеряемом спектре сигнала (Мартынов В.А., Селихов Ю. И. Панорамные приемники и анализаторы спектра. Под ред. Г.Д.Заварина, 2-е изд. перераб. и доп. М.: Советское радио, 1980, с.352).
Наиболее близким к предлагаемому способу является способ определения частоты, заключающийся в получении дискретных значений входного сигнала, определении номера максимальной спектральной составляющей, полученного с помощью дискретного преобразования Фурье (ДПФ), и определении частоты сигнала как частоты частотной составляющей с максимальной амплитудой [1].
Недостатком этих способов является недостаточная точность определения частоты при ограниченных времени анализа и аппаратурных затратах.
Наиболее близким к предлагаемому устройству является устройство измерения частоты сигнала, состоящее из аналого-цифрового преобразователя (АЦП), анализатора спектра на основе ДПФ и решающего устройства. Оценка частоты сигнала определяется по номеру частотного канала с максимальной амплитудой выходного напряжения [2].
Недостатком данного устройства является недостаточная точность определения частоты при ограниченном времени анализа и больших аппаратурных затратах.
Целью изобретения является повышение точности определения частоты сигнала и уменьшение аппаратурных затрат при неизменном времени анализа.
Цель достигается тем, что в способе измерения частоты сигнала, основанном на аналого-цифровом преобразовании сигнала и подсчете количества мгновенных отсчетов, дополнительно производят формирование вспомогательного сигнала и подсчет количества его отсчетов, начиная с момента его формирования, прекращение подсчета количества отсчетов вспомогательного сигнала при появлении нулевого его значения, формирование значения частоты сигнала путем уменьшения значения частоты дискретизации сигнала пропорционально количеству отсчетов вспомогательного сигнала, накопленному между началом его формирования и появлением нулевого отсчета. При этом вспомогательный сигнал формируется путем суммирования значений мгновенных отсчетов анализируемого сигнала.
Причем в устройство, реализующее способ измерения частоты, содержащее аналого-цифровой преобразователь, вход которого соединен с входной шиной, введены накапливающий сумматор, узел сравнения, триггер, счетчик, делитель, ключ, генератор импульсов, выход которого подключен к второму входу аналого-цифрового преобразователя и к первому входу ключа, второй вход ключа соединен с выходом триггера, вход ключа подсоединен к вторым входам накапливающего сумматора, узла сравнения и счетчика, выход аналого-цифрового преобразователя подключен к первому входу накапливающего сумматора и к первому входу триггера, выход накапливающего сумматора соединен с первым входом узла сравнения, выход которого подключен к третьему входу накапливающего сумматора, второму входу триггера, к первому входу счетчика и второму входу делителя, вход делителя соединен с выходом счетчика.
Именно использование формирования вспомогательного сигнала с определенными свойствами и определение количества его отсчетов между началом формирования и появлением нулевого значения с последующим формированием значения частоты сигнала позволили достичь поставленной цели изобретения в предложенном способе. Также именно введение накапливающего сумматора, узла сравнения делителя, триггера, счетчика, генератора импульсов, ключа с соответствующими связями в устройство, реализующее способ, позволили достичь поставленной цели изобретения. Это позволяет сделать вывод, что заявляемые изобретения связаны между собой единым изобретательским замыслом.
Сравнение заявляемых технических решений с прототипами позволило установить соответствие их критерию "новизна". При изучении других известных технических решений в данной области техники установлено, что признаки, отличающие заявляемые изобретения от прототипов, не выявлены в них, и поэтому эти признаки обеспечивают заявляемому техническому решению соответствие критерию "существенные отличия".
На фиг.1 приведена структурная схема устройства, реализующего заявляемый способ; на фиг.2 приведены временные диаграммы работы устройства.
Пусть необходимо определить частоту fs гармонического или квазигармонического узкополосного сигнала, описываемого соотношением Sвх(t) = Assin( ωst+φн). Частота сигнала fs = ω/2π лежит в пределах полосы анализа Δ fa = { 0,fmax } и кратна частоте дискретизации Fд. Для определения частоты fs входного сигнала производят его следующую обработку.
Производим аналого-цифровое преобразование сигнала Sвх(t) с частотой дискретизации Fд≥2fa. В результате анализируемый сигнал Sвх(t) преобразуется в последовательность цифровых отсчетов Sвх(K):
Sвх(n)=Assin(ωsΔtnн)=Assin(ωsTдn+φn),(1) где n = 0,1,2...;
Тд = 1/Fд - период дискретизации.
Затем сформируем вспомогательный сигнал S(m) путем суммирования цифровых отсчетов последовательности Sвх(K), описываемый следующим образом:
S(m) =
Figure 00000001
Assin(ωsTдK+φн), (2) где φн≠2πTдl/Ts, l - целое число;
m - текущее значение количества просуммированных отсчетов.
Одновременно с началом формирования вспомогательного сигнала S(m) производим подсчет количества m суммируемых отсчетов Sвх(n).
Во вспомогательном сигнале S(m) контролируем появление нулевого значения с начала его формирования. При выявлении этого нулевого значения S(mo) = 0 накопление отсчетов прекращаем и фиксируем соответствующее количество mo накопленных отсчетов.
Частоту fs сигнала Sвх(t) формируем следующим образом:
fs=
Figure 00000002
(3)
Покажем, что частота сигнала fs и момент появления (mo +1)Tднулевого значения S(mo) имеют однозначную взаимосвязь. Для этого преобразуем уравнение (2) следующим образом:
S(m) =
Figure 00000003
Ascosφнsin(ωsTдn)+
Figure 00000004
Assinφncos(ωsTдn) = Sc(m)+Ss(m), где Sc(m) =
Figure 00000005
Assinφнcos(ωsTдn);
Ss(m) =
Figure 00000006
Ascosφнsin(ωsTдn).
Нулевое значение вспомогательного сигнала S(m) появится в случае:
S (mo ) = 0 при
Figure 00000007
Figure 00000008
Figure 00000009
(5)
Составляющие Sc(m) и Ss(m) можно описать суммами рядов с тригонометрическими членами, сходящимися по признаку Дирихле. При этом суммы синусного и косинусного рядов равны (4):
Zc(X) =
Figure 00000010
cospx =
Figure 00000011

Zs(X) =
Figure 00000012
sinpx =
Figure 00000013
(6)
Уравнения (4) с учетом уравнений (6) при x = ωs Tд преобразуются к виду:
Figure 00000014
Figure 00000015
Figure 00000016
(7)
Совместим решением системы (7) для случая (5) является значение mo= KFд/fs-1. При этом для Kmin = 1 выполняется равенство:
mo=
Figure 00000017
- 1 (8)
Так как значение Fд известно заранее, то из уравнения (8) можно определить искомое значение частоты входного сигнала:
fs=
Figure 00000018
(9)
Таким образом, по значению mo номера отсчета формируемого вспомогательного сигнала S(m), при котором последний принимает нулевое значение, можно однозначно судить о частоте анализируемого сигнала.
Определим погрешность δ fs измерения частоты сигнала предлагаемым способом. Время анализа сигнала равно времени накопления его отсчетов:
Та = Тд˙ Ns, (10) где Та - время анализа;
Ns - количество накапливаемых отсчетов.
Тогда достижимая погрешность измерения частоты сигнала Sвх(t) не превысит половины формируемого дискрета по частоте:
δfs
Figure 00000019
_
Figure 00000020
=
Figure 00000021
(11)
По способу-прототипу при неизменном времени анализа (см.10) погрешность измерения частоты составит (1):
δf * s =
Figure 00000022
=
Figure 00000023
(12)
Из сравнения равенств (11) и (12) следует, что при использовании предлагаемого способа точность измерения частоты может быть повышена в М = 1/(Ns-1) раза по сравнению с прототипом без увеличения времени анализа.
Обычно количество обрабатываемых отсчетов в процедуре ДПФ лежит в пределах (28-212). Поэтому ожидаемый выигрыш по точности при использовании предлагаемого способа составит соответственно 2-3 порядка.
Важным параметром аппаратуры измерения параметров сигналов является время реакции, равное интервалу времени между началом анализа и получением искомого значения частоты. При использовании предлагаемого способа вспомогательный сигнал S(m) может формироваться в темпе поступления отсчетов последовательности Sвх(n) и время реакции Тр будет равно времени анализа:
Tp = Ta (13)
Для повышения быстродействия в способе-прототипе при Ns = 2lможно использовать процедуру быстрого преобразования Фурье (БПФ). В этом случае время реакции практически будет не менее:
Тр * ≥ Та + ТБПФ (14) где ТБПФ - время обработки накопленного массива отсчетов по алгоритму БПФ.
При режиме обработки в реальном масштабе времени Та = ТБПФ и, следовательно, Тр * ≥ 2Та. Таким образом, время реакции Тр в предлагаемом способе меньше времени реакции Тр * в прототипе
Тр < Тр* (15)
Таким образом, при использовании предлагаемого способа обеспечивается возможность существенно сократить время реакции измерения частоты сигнала при сохранении неизменным времени анализа.
Для реакции предлагаемого способа необходимы простые операции, такие как аналого-цифровое преобразование, накопление и сравнение.
Для реализации известного способа помимо аналого-цифрового преобразования и накопления необходимы процедуры ДПФ и поиска максимума с вычислением значения частоты, которые намного сложнее процедуры,реализующей предлагаемый способ, как по количеству операций, так и по сложности их алгоритмических взаимосвязей.
Устройство, реализующее способ, содержит (фиг.1) генератор 1 импульсов, АЦП 20,ключ 3, накапливающий сумматор 4, счетчик 5 импульсов, узел 6 сравнения, делитель 7, триггер 8.
Устройство работает следующим образом (фиг.2). Входной анализируемый сигнал Sвх(t) поступает на первый (сигнальный) вход АЦП 2. Запуска АЦП 2 производится тактовыми импульсами с генератора 1 импульсов, поступающими на его второй вход (вход запуска). Период следования импульсов равен периоду дискретизации Тд входного сигнала Sвх(t). В исходном состоянии (после включения питания) триггер 8 находится в нулевом состоянии и его выходной сигнал, поступая на второй вход ключа 3, препятствует прохождению через него тактовых импульсов с генератора 1. Накапливающий сумматор 4, счетчик 5 импульсов и делитель 7 также находится в нулевом состоянии. АЦП 2 производит преобразование входного сигнала Sвх(t) в эквивалентную последовательность цифровых отсчетов Sвх(n), которая поступает на накапливающий сумматор 4. При этом знаковый разряд цифровых отсчетов с АЦП 2 соединен с первым синхровходом триггера 8. Триггер 8 является Т-триггером с двумя входами: первым - динамическим синхровходом (С-вход) и вторым - установочным в нулевое состояние (Р-вход). В момент перехода значения входного сигнала Sвх(t) через ноль его знак изменяется с отрицательного на положительный. При этом знаковый разряд соответствующих цифровых отсчетов Sвх(n) изменяет свое значение с нулевого на единичное (при использовании прямого двоичного кодирования) (фиг.2, т.А). В этом случае на первом (С-входе) входе триггера 8 формируется положительный фронт (переход из логического нуля в единицу), по которому триггер 8 переходит в единичное состояние (на его выходе формируется сигнал с уровнем логической единицы). Единичный сигнал с выхода триггера 8, поступая на второй вход ключа 1, разрешает прохождение через него тактовых импульсов с генератора 4 импульсов. Таким образом обеспечивается синхронизация начала измерения частоты с моментом появления сигнала. Цифровые отсчеты с АЦП 2 по синхроимпульсам с ключа 3 принимаются накапливающим сумматором 4, на выходе которого формируется вспомогательный сигнал S(m) (см. (2)). Одновременно с формированием значений вспомогательного сигнала S(m) счетчик 5 производит подсчет текущего количества m отсчетов, накопленных в накапливающем сумматоре 4. Значения отсчетов вспомогательного сигнала S(m) поступают на узел 6 сравнения, который сравнивает их значения с порогом Пo. Работа узла 6 сравнения стробируется импульсами запуска с генератора 1 импульсов. При появлении значения формируемого сигнала S(mo), равного нулю или меньше по модулю порога По, на выходе узла 6 сравнения формируется импульс, который поступает на второй (R-вход) вход триггера 8 и переводит его в нулевое состояние, тем самым закрывая ключ 1 и препятствуя дальнейшему накоплению отсчетов с АЦП 2 в накапливающем сумматоре 4 и подсчету их количества в счетчике 5. Одновременно с этим полученное значение количества накопленных отсчетов (mo + 1) из счетчика 5 фиксируется (записывается) в делителе 7, который формирует на выходе измеренное значение частоты fs входного сигнала Sвх(t) (см.(3)). Накапливающий сумматор 4 и счетчик 5 обнуляются, переходя в исходное состояние. Очередной цикл измерения частоты входного сигнала повторится аналогично описанному выше порядку. Таким образом устройство производит измерение частоты входного сигнала.
Время анализа предлагаемого устройства с учетом (9) и (10) равно:
Та = (mo + 1)Тд при max{m} = Ns (15)
Время реакции предлагаемого устройства равно интервалу времени между моментами появления перехода значения входного сигнала Sвх(t) через ноль и выдачи результата измерения:
Tp=Ta467,(16) где τ4 - время задержки сумматора 4;
τ6 - время задержки узла 6 сравнения;
τ7- время задержки делителя 7.
Для реального масштаба времени измерений при большом значении mo(более 10) практически выполняется равенство:
Тр ≈ Та (17)
Достижимая погрешность δf измерения частоты сигнала не превышает:
δf
Figure 00000024
(18)
Для устройства-прототипа время реакции равно:
Figure 00000025
+ Tf , (19) где ТДПФ - время вычисления ДПФ накопленного массива из отсчетов;
Tf - время вычисления частоты сигнала.
Для режима реального масштаба времени при Та = ТДПФ = Тf
Тр * ≥ Та(2-3) (20)
При этом для устройства-прототипа достижимая погрешность измерения частоты сигнала не превысит:
δf * s
Figure 00000026
(21)
Из сравнения параметров предлагаемого устройства и прототипа видно, что при одинаковом времени анализа Та точность измерения частоты сигнала и время реакции у предлагаемого устройства существенно лучше. По аппаратурным затратам предлагаемое устройство требует существенно меньших аппаратурных затрат по сравнению с прототипом (примерно один-два модуля), а все его узлы (АЦП, счетчик, сумматор, делитель, триггер и ключ) являются типовыми для аппаратуры цифровой обработки сигналов.
Устройство-прототип измерения частоты сигнала состоит из АЦП, прецессора ДПФ, устройства поиска максимума и вычислителя частоты. Процессор ДПФ является наиболее сложным вычислительным устройством. Для реального масштаба времени при Ns ≥ 1024 его объем составляет не менее целого блока с количеством микросхем до нескольких тысяч. Устройства поиска максимума и вычисления частоты также довольно сложные. Следовательно, аппаратурные затраты предлагаемого устройства существенно меньше аналогичных затрат устройства-прототипа.
Таким образом, предлагаемое устройство измерения частоты сигнала при неизменном времени анализа обеспечивает более высокую точность измерения частоты при существенном сокращении аппаратурных затрат и времени реакции.

Claims (2)

1. Способ измерения частоты сигнала, заключающийся в аналого-цифровом преобразовании сигнала и подсчете количества мгновенных отсчетов, отличающийся тем, что, с целью повышения точности и упрощения при неизменном времени анализа, дополнительно производят формирование вспомогательного сигнала и подсчет количества его отсчетов, начиная с момента начала измерения, сравнение значений отсчетов вспомогательного сигнала с нулем, прекращение подсчета количества отсчетов вспомогательного сигнала при появлении нулевого его значения, формирование значения частоты сигнала путем уменьшения значения частоты дискретизации сигнала пропорционально количеству отсчетов вспомогательного сигнала, накопленных между началом его формирования и появлением нулевого отсчета, при этом вспомогательный сигнал формируется путем суммирования значения мгновенных отсчетов анализируемого сигнала.
2. Устройство измерения частоты сигнала, содержащее аналого-цифровой преобразователь, вход которого соединен с входной шиной, отличающееся тем, что, с целью повышения точности измерения частоты и сокращения аппаратурных затрат при неизменном времени анализа, в него дополнительно введены накапливающий сумматор, узел сравнения, триггер, счетчик, делитель, ключ, генератор импульсов, выход которого подключен к второму входу аналого-цифрового преобразователя и к первому входу ключа, второй вход ключа соединен с выходом триггера, выход ключа подсоединен к вторым входам накапливающего сумматора, узла сравнения и счетчика, выход аналого-цифрового преобразователя подключен к первому входу накапливающего сумматора и к первому входу триггера, выход накапливающего сумматора соединен с первым входом узла сравнения, выход которого подключен к третьему входу накапливающего сумматора, второму входу триггера, к первому входу счетчика и второму входу делителя, вход делителя соединен с выходом счетчика.
SU4928486 1991-04-22 1991-04-22 Способ измерения частоты сигнала и устройство для его осуществления RU2017162C1 (ru)

Priority Applications (1)

Application Number Priority Date Filing Date Title
SU4928486 RU2017162C1 (ru) 1991-04-22 1991-04-22 Способ измерения частоты сигнала и устройство для его осуществления

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
SU4928486 RU2017162C1 (ru) 1991-04-22 1991-04-22 Способ измерения частоты сигнала и устройство для его осуществления

Publications (1)

Publication Number Publication Date
RU2017162C1 true RU2017162C1 (ru) 1994-07-30

Family

ID=21570322

Family Applications (1)

Application Number Title Priority Date Filing Date
SU4928486 RU2017162C1 (ru) 1991-04-22 1991-04-22 Способ измерения частоты сигнала и устройство для его осуществления

Country Status (1)

Country Link
RU (1) RU2017162C1 (ru)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2478213C2 (ru) * 2011-10-17 2013-03-27 Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Санкт-Петербургский государственный политехнический университет" (ФГБОУ ВПО "СПбГПУ") Способ измерения частоты сигнала
RU2791820C1 (ru) * 2022-05-20 2023-03-13 АО "Автограф" Способ определения несущей частоты дискретного сигнала

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
1. Финкельштейн М.И. Основы радиолокации. М.: Радио и связь, 1983, с.324-344. *
2. Кузьмин С.З. Основы проектирования систем цифровой обработки радиолокационной информации. М.: Радио и связь, 1986, с.103. *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2478213C2 (ru) * 2011-10-17 2013-03-27 Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Санкт-Петербургский государственный политехнический университет" (ФГБОУ ВПО "СПбГПУ") Способ измерения частоты сигнала
RU2791820C1 (ru) * 2022-05-20 2023-03-13 АО "Автограф" Способ определения несущей частоты дискретного сигнала

Similar Documents

Publication Publication Date Title
EP0710846B1 (en) Apparatus for generating a signal representative of the total harmonic distortion in waveforms of an AC electrical system
Wagdy Effect of ADC quantization errors on some periodic signal measurements
Friedlander Parametric signal analysis using the polynomial phase transform
Glisson et al. The digital computation of discrete spectra using the fast Fourier transform
RU2017162C1 (ru) Способ измерения частоты сигнала и устройство для его осуществления
US4296374A (en) Wideband digital spectrometer
JPS5819068B2 (ja) デンシシキデンリヨクリヨウケイ
Sienkowski et al. Simple, fast and accurate four-point estimators of sinusoidal signal frequency
Cennamo et al. Dynamic testing and diagnostics of digitizing signal analyzers
US4181949A (en) Method of and apparatus for phase-sensitive detection
US4823293A (en) Correlation time-difference detector
RU2591742C1 (ru) Способ измерения частоты гармонического сигнала и устройство для его осуществления
CN104950168A (zh) 一种基于二次平均的低信噪比正弦信号高精度频率测量方法
RU2096788C1 (ru) Статистический анализатор качества параметров электрической энергии
Weiss et al. Confidence on the modified Allan variance and the time variance
JPS5526413A (en) Waveform analyzing unit
SU656018A1 (ru) Устройство дл измерени длительности импульсов со случайным периодом следовани
RU2256928C2 (ru) Способ измерения нестабильности частоты и устройство для его осуществления
SU1027734A1 (ru) Устройство дл определени двумерной плотности веро тности случайного процесса
Pawłowski Software-based method of increasing the effective resolution of a measurement chain for a transducer with a pulse frequency output
SU690403A1 (ru) Частотометр
SU842694A1 (ru) Цифровой измеритель временных интер-ВАлОВ C ВЕСОВыМ уСРЕдНЕНиЕМ
SU953573A1 (ru) Способ стробоскопического преобразовани периодических электрических сигналов
SU1013905A1 (ru) Устройство дл определени центра т жести повтор ющихс импульсов
SU813290A1 (ru) Устройство дл измерени центральнойчАСТОТы СпЕКТРА СигНАлА