RU2475765C1 - Способ определения спектра электрического сигнала по измеренным выборочным значениям этого сигнала - Google Patents
Способ определения спектра электрического сигнала по измеренным выборочным значениям этого сигнала Download PDFInfo
- Publication number
- RU2475765C1 RU2475765C1 RU2011143194/28A RU2011143194A RU2475765C1 RU 2475765 C1 RU2475765 C1 RU 2475765C1 RU 2011143194/28 A RU2011143194/28 A RU 2011143194/28A RU 2011143194 A RU2011143194 A RU 2011143194A RU 2475765 C1 RU2475765 C1 RU 2475765C1
- Authority
- RU
- Russia
- Prior art keywords
- signal
- spectrum
- vector
- sample values
- analysed
- Prior art date
Links
Landscapes
- Measurement Of Resistance Or Impedance (AREA)
Abstract
Изобретение относится к способам определения спектра электрических сигналов. Анализируемый электрический сигнал дискретизируют по времени и измеряют его дискретные выборочные значения, причем до проведения измерений интервал контролируемых частот, который включает спектр анализируемого сигнала, разбивают на элементы разрешения по частоте, размер которых определяется требуемой точностью спектрального анализа, выбирают шаг дискретизации сигнала по времени исходя из удобства проведения измерений, формируют с использованием размеров элемента разрешения и шага дискретизации весовую матрицу, измеренные выборочные значения анализируемого сигнала объединяют в вектор, составляют векторно-матричное уравнение измерений и определяют из него спектр анализируемого сигнала в виде вектора, компоненты которого представляют собой оценки дискретизированных по элементам разрешения составляющих спектра анализируемого сигнала. Технический результат заключается в расширении диапазона анализируемых частот и возможности снижения требований к быстродействию используемых АЦП.
Description
Изобретение относится к области радиоэлектроники, а именно к спектральному анализу электрических сигналов.
Определение спектра электрического сигнала (спектральный анализ) является актуальной задачей в области обработки и преобразования сигналов. Одним из подходов к решению этой задачи является определение спектра анализируемого сигнала по его временной выборке. По этому принципу реализуется дискретное преобразование Фурье (ДПФ), описанное, например, в [1] - прототип заявляемого способа.
Способ-прототип заключается в том, что анализируемый сигнал f(t) дискретизируют по времени с шагом дискретизации Т, преобразуют выборочные значения сигнала в цифровую форму с помощью аналого-цифрового преобразователя (АЦП) и определяют значения спектра сигнала, исходя из преобразования Фурье по формуле
где n - номер шага дискретизации по времени, ω - круговая частота, N - длина временной выборки сигнала.
Определение спектральной составляющей по формуле (1) может выполняться для любых значений частоты ω, однако для цифровой реализации удобно определять спектр в дискретных точках оси частот, отстоящих одна от другой на постоянную величину - шаг дискретизации по частоте Ω. В этом случае выражение (1) приобретает вид
где k - номер шага дискретизации по частоте.
Способ-прототип имеет следующие недостатки. При определении спектра сигнала f(t) по формуле (1) необходимо обеспечить однозначное представление этого сигнала его временной выборкой f(nT). Согласно теореме Котельникова, для этого шаг дискретизации сигнала по времени должен удовлетворять условию , где ωmax - максимальная частота спектра сигнала. На практике, согласно [2], это условие еще ужесточается, так что шаг дискретизации должен быть . Таким образом, чем выше частота анализируемого сигнала, тем меньший шаг дискретизации Т необходимо выбирать для выполнения этого условия и, соответственно, тем более быстродействующий и дорогой требуется АЦП, используемый для представления выборочных значений сигнала в цифровой форме. Как следствие, ДПФ имеет ограничение по частоте анализируемого сигнала. Кроме этого повышение частоты анализируемых сигналов сопровождается повышением стоимости анализаторов спектра из-за высоких требований к быстродействию АЦП.
Технической задачей данного изобретения является создание способа определения спектра электрического сигнала по измеренным его выборочным значениям, полученным с произвольным шагом дискретизации по времени, не зависящим от частоты анализируемого сигнала, что позволяет расширить диапазон анализируемых частот в сторону их увеличения, а также снизить стоимость анализаторов спектра за счет применения в них менее быстродействующих, но более дешевых АЦП.
Поставленная задача достигается тем, что в способе определения спектра электрического сигнала, который заключается в дисретизации анализируемого сигнала по времени с шагом дискретизации T и измерении его дискретизированных выборочных значений путем преобразования их в цифровую форму, согласно изобретению, до проведения измерений интервал контролируемых частот (ωн, ωк), который включает спектр анализируемого сигнала, разбивают на элементы разрешения по частоте, размер которых Ω определяется требуемой точностью спектрального анализа, выбирают шаг дискретизации сигнала по времени исходя из удобства проведения измерений, формируют весовую матрицу W с компонентами где k - номер элемента разрешения по частоте, n - номер измеренного выборочного значения сигнала, k=1, 2, …, К, n=1, 2, …, N, К - число элементов разрешения по частоте в интервале контролируемых частот, N - число измеренных выборочных значений анализируемого сигнала, измеренные выборочные значения анализируемого сигнала f(T)f(2T)…f(NT) объединяют в вектор , где индекс Т обозначает транспонирование, составляют векторно-матричное уравнение измерений и определяют из него спектр анализируемого электрического сигнала в виде оценки спектрального вектора , компоненты которого представляют собой оценки днскретизированных по элементам разрешения составляющих спектра анализируемого сигнала.
Поставленная задача решается за счет того, что вместо определения составляющих спектра сигнала исходя из прямого преобразования Фурье, как это делается в прототипе согласно (1), в заявляемом способе спектральные составляющие определяют исходя из обратного преобразования Фурье как дискретные значения подынтегральный функции S(ω). При этом определяющим оказывается шаг дискретизации по частоте Ω подынтегральной функции S(ω), а моменты измерений сигнала могут выбираться произвольно, аналогично тому, как в прототипе выражение (1) может быть вычислено для произвольного значения ω, но при достаточно малом шаге дискретизации по времени Т.
Особенностью предлагаемого способа является то, что выборочные значения сигнала могут измеряться в любые удобные для измерения моменты времени, точность же спектрального анализа задается априори размером элемента разрешения по частоте, который должен быть достаточно малым.
Обоснование способа.
Запишем обратное преобразование Фурье, связывающее анализируемый сигнал f(t) с его спектром S(ω):
где t - время, ω - круговая частота.
Выберем интервал контролируемых частот (ωн, ωк), включающий в себя спектр анализируемого сигнала. Разобьем этот интервал на малые элементы разрешения, размер которых Ω определяется требуемой точностью спектрального анализа и перейдем от интеграла (2) к интегральной сумме
где k - номер элемента разрешения, К - число элементов разрешения в интервале контролируемых частот.
Введем обозначения:
Выражение (4) определяет весовые коэффициенты при дискретных составляющих спектра сигнала (5). Для всех К элементов разрешения сформируем векторы
С учетом (4)-(7) перепишем интегральную сумму в правой части (3) в векторной форме:
Спектральный вектор (7), входящий в (8) - это вектор дискретизированного по элементам разрешения спектра анализируемого сигнала. Определив этот вектор, мы определим спектр сигнала с точностью элемента разрешения Ω, выбранного априори.
Для определения спектрального вектора проведем измерения выборочных значений анализируемого сигнала в дискретные моменты времени. Для удобства практических измерений будем полагать эти моменты равноотстоящими на шаг дискретизации Т, который, в отличие от прототипа, выбирается исходя из возможностей измерительной техники и удобства измерений, в частности - исходя из ограниченного быстродействия АЦП. Спектр сигнала при этом не влияет на выбор шага дискретизации.
Проведенные измерения выборочных значений сигнала запишем аналогично (8) в виде системы линейных уравнений:
где N - число измеренных значений сигнала.
Перепишем систему уравнений (9) в векторно-матричной форме
Компоненты весовой матрицы W определяются подстановкой в (4) моментов измерений Т, 2T, …, NT и равны
где k - номер элемента разрешения по частоте, n - номер момента измерения сигнала.
Определим оценку спектрального вектора из векторно-матричного уравнения (10). Сделать это можно, например, методом псевдообращения [3]:
где индекс + обозначает операцию псевдообращения матрицы.
Оценка (11) является оценкой искомого спектрального вектора. Компоненты вектора представляют собой оценки дискретизированных по элементам разрешения составляющих спектра анализируемого сигнала, что является решением поставленной задачи. Заметим, что спектральные составляющие - компоненты вектора - определены с дискретностью Ω, заданной априори, по значениям временного сигнала, измеренным с шагом дискретизации Т, который задается независимо от частотного спектра сигнала, из соображений удобства проведения измерений. Это позволяет использовать более «редкие», чем в прототипе, измерения при одном и том же частотном спектре сигнала, что, в свою очередь, дает возможность применять менее быстродействующие, зато более дешевые АЦП. Кроме этого возможность выбирать отсчеты временного сигнала независимо от частоты самого сигнала позволяет проводить спетральный анализ сигналов на частотах, недоступных для оцифровки даже быстродействующими АЦП при шаге дискретизации, требуемой в прототипе.
Преимущества предлагаемого способа по сравнению с прототипом следующие.
1. Снижение стоимости спектрального анализа, что обусловлено возможностью использования менее быстродействующих, но более дешевых АЦП по сравнению с прототипом.
2. Возможность расширения диапазона анализируемых частот в сторону высоких частот. Это является следствием того, что в заявляемом способе снимается условие при выборе шага дискретизации временного сигнала, необходимое в прототипе.
Источники информации
1. Голд Б., Рэйдер Ч. Цифровая обработка сигналов. Пер. с англ., под ред. A.M.Трахтмана. - М., «Советское радио», 1973, стр.187-191.
2. Гоноровский И.С. Радиотехнические цепи и сигналы: Учебник для вузов. - 4-е изд., перераб. и доп. - М.: Радио и связь, 1986, стр.65.
3. Гантмахер Ф.Р. Теория матриц. 4-е изд. - М.: Наука. Гл. ред. физ.-мат. лит., 1988, с.35.
Claims (1)
- Способ определения спектра электрического сигнала по измеренным выборочным значениям этого сигнала, заключающийся в том, что анализируемый электрический сигнал дискретизируют по времени с шагом дискретизации Т, измеряют его дискретизированные выборочные значения путем преобразования их в цифровую форму, отличающийся тем, что до проведения измерений интервал контролируемых частот (ωн, ωк), который включает спектр анализируемого сигнала, разбивают на элементы разрешения по частоте, размер которых Ω определяется требуемой точностью спектрального анализа, выбирают шаг дискретизации сигнала по времени, исходя из удобства проведения измерений, формируют весовую матрицу W с компонентами , где k - номер элемента разрешения по частоте, n - номер измеренного выборочного значения сигнала, k=1, 2, …, К, n=1, 2, …, N, К - число элементов разрешения по частоте в интервале контролируемых частот, N - число измеренных выборочных значений анализируемого сигнала, измеренные выборочные значения анализируемого сигнала f(Т)f(2Т)…f(NT) объединяют в вектор где индекс Т обозначает транспонирование, составляют векторно-матричное уравнение измерений и определяют из него спектр анализируемого электрического сигнала в виде оценки спектрального вектора , компоненты которого представляют собой оценки дискретизированных по элементам разрешения составляющих спектра анализируемого сигнала.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
RU2011143194/28A RU2475765C1 (ru) | 2011-10-26 | 2011-10-26 | Способ определения спектра электрического сигнала по измеренным выборочным значениям этого сигнала |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
RU2011143194/28A RU2475765C1 (ru) | 2011-10-26 | 2011-10-26 | Способ определения спектра электрического сигнала по измеренным выборочным значениям этого сигнала |
Publications (1)
Publication Number | Publication Date |
---|---|
RU2475765C1 true RU2475765C1 (ru) | 2013-02-20 |
Family
ID=49121099
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
RU2011143194/28A RU2475765C1 (ru) | 2011-10-26 | 2011-10-26 | Способ определения спектра электрического сигнала по измеренным выборочным значениям этого сигнала |
Country Status (1)
Country | Link |
---|---|
RU (1) | RU2475765C1 (ru) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SU917119A1 (ru) * | 1979-02-28 | 1982-03-30 | Харьковский Институт Инженеров Железнодорожного Транспорта Им.С.М.Кирова | Анализатор комплексного спектра |
SU1191785A1 (ru) * | 1984-05-24 | 1985-11-15 | Ленинградский Ордена Ленина И Ордена Трудового Красного Знамени Государственный Университет Им.А.А.Жданова | Модул ционный способ спектрального анализа |
RU1799474C (ru) * | 1991-03-11 | 1993-02-28 | Нижегородский научно-исследовательский приборостроительный институт "Кварц" | Анализатор спектра |
US20090076737A1 (en) * | 2004-10-28 | 2009-03-19 | Cerno Bioscience Llc | Qualitative and quantitative mass spectral analysis |
-
2011
- 2011-10-26 RU RU2011143194/28A patent/RU2475765C1/ru not_active IP Right Cessation
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SU917119A1 (ru) * | 1979-02-28 | 1982-03-30 | Харьковский Институт Инженеров Железнодорожного Транспорта Им.С.М.Кирова | Анализатор комплексного спектра |
SU1191785A1 (ru) * | 1984-05-24 | 1985-11-15 | Ленинградский Ордена Ленина И Ордена Трудового Красного Знамени Государственный Университет Им.А.А.Жданова | Модул ционный способ спектрального анализа |
RU1799474C (ru) * | 1991-03-11 | 1993-02-28 | Нижегородский научно-исследовательский приборостроительный институт "Кварц" | Анализатор спектра |
US20090076737A1 (en) * | 2004-10-28 | 2009-03-19 | Cerno Bioscience Llc | Qualitative and quantitative mass spectral analysis |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8441379B2 (en) | Device and method for digitizing a signal | |
Zhao et al. | Minimum rate sampling and spectrum-blind reconstruction in random equivalent sampling | |
Wang et al. | A fully fledged TDC implemented in field-programmable gate arrays | |
US20080211482A1 (en) | Method and Apparatus for Measuring the Frequency of a Received Signal | |
KR20090042776A (ko) | 신호 처리 방법 및 디바이스와 컴퓨터 판독가능 매체 | |
RU2475765C1 (ru) | Способ определения спектра электрического сигнала по измеренным выборочным значениям этого сигнала | |
Ye et al. | A field-programmable-gate-array based time digitizer for the time-of-flight mass spectrometry | |
Panek et al. | Time interval measurement device based on surface acoustic wave filter excitation, providing 1ps precision and stability | |
CN106153709A (zh) | 时间间隔测量 | |
Ozawa et al. | SAR TDC architecture with self-calibration employing trigger circuit | |
Leuciuc | Sampling time calibration method for multi-channel interleaved ADCs | |
Kobayashi et al. | Fine time resolution TDC architectures-integral and delta-sigma types | |
Comberiate et al. | A bandpass sampling receiver for wide-bandwidth, spectrally-sparse waveforms for high-accuracy range measurements | |
Zhao et al. | A sparse signal reconstruction approach for sequential equivalent time sampling | |
Zhao et al. | Model-based multichannel compressive sampling with ultra-low sampling rate | |
Wang et al. | Estiamtion and correction of mismatch errors in time-interleaved ADCs | |
EP4063873A1 (en) | Signal analysis device, control circuit, and storage medium | |
RU2455653C1 (ru) | Способ спектрального анализа электрического сигнала | |
Urekar et al. | Low resolution stochastic Flash ADC for high precision energy and RMS voltage measurements for smart grid | |
Kowalewski et al. | Fast high-impedance spectroscopy method using sinc signal excitation | |
Zhou et al. | A Residual-Remainder Coupled Unlimited Sampling Framework for High Dynamic Range Signal Conversion | |
WO2022249747A1 (ja) | Adコンバータ及びこれを用いた広帯域雑音の電力測定装置、電力スペクトル測定装置 | |
Slim et al. | Digital automatic calibration method for a time-interleaved ADCs system used in time-domain EMI measurement receiver | |
Wang et al. | DACS of Microwave Limb Sounder | |
US10476483B2 (en) | Decimation filter |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
MM4A | The patent is invalid due to non-payment of fees |
Effective date: 20181027 |