RU2660320C2 - Способ дискретизации и восстановления непрерывного сигнала - Google Patents
Способ дискретизации и восстановления непрерывного сигнала Download PDFInfo
- Publication number
- RU2660320C2 RU2660320C2 RU2015129169A RU2015129169A RU2660320C2 RU 2660320 C2 RU2660320 C2 RU 2660320C2 RU 2015129169 A RU2015129169 A RU 2015129169A RU 2015129169 A RU2015129169 A RU 2015129169A RU 2660320 C2 RU2660320 C2 RU 2660320C2
- Authority
- RU
- Russia
- Prior art keywords
- samples
- sampling
- spectrum
- formula
- continuous signal
- Prior art date
Links
Images
Classifications
-
- G—PHYSICS
- G11—INFORMATION STORAGE
- G11B—INFORMATION STORAGE BASED ON RELATIVE MOVEMENT BETWEEN RECORD CARRIER AND TRANSDUCER
- G11B20/00—Signal processing not specific to the method of recording or reproducing; Circuits therefor
Landscapes
- Engineering & Computer Science (AREA)
- Signal Processing (AREA)
- Data Exchanges In Wide-Area Networks (AREA)
- Mobile Radio Communication Systems (AREA)
- Digital Transmission Methods That Use Modulated Carrier Waves (AREA)
Abstract
Изобретение относится к области измерительной технике и предназначено для использования в системах контроля окружающей среды и технологических процессов. При дискретизации отсчеты берут пакетами по m отсчетов с задержками от первого , последовательность из N пакетов отсчетов регистрируют или передают в канал связи. При восстановлении из последовательности с пакетами отсчетов формируют m последовательностей одиночных отсчетов {х(k)}, задержанных на интервал времени τν(τ1=0) от первой, дискретным преобразованием Фурье, вычисляют спектр с алиасингом каждой из последовательностей {xν(k)} по формуле где для каждой из частот неизвестного спектра формируют суммы взвешенных задержками гармоник частот наложения по уравнению где из которого формируют систему линейных алгебраических уравнений относительно неизвестных составляющих гармоник (-m≤i≤m), решением системы определяют неизвестные, обратным преобразованием Фурье восстанавливают непрерывный сигнал x(t) по формуле Технический результат заключается в расширении спектра сигнала. 6 ил.
Description
Изобретение относится к области измерительной техники и предназначено для использования в системах контроля окружающей среды и технологических процессов.
Известен способ дискретизации и восстановления непрерывных сигналов с ограниченным спектром, вытекающий из теоремы Котельникова [Котельников В.Л. О пропускной способности эфира и проволоки в электросвязи. Материалы к 1-му Всесоюзному съезду, ВЭК, 1933, Харкевич А.А. Спектры и анализ. Издание четвертое. Изд. ФМЛ. - М., 1962. - 236 с.].
По этому способу цифровые отсчеты непрерывного сигнала берут с частотой, в два раза превышающей верхнюю граничную частоту в спектре сигнала, а восстанавливают по интерполяционной формуле (ряду Котельникова), содержащей бесконечную сумму взвешенных по функциям sin x/x отсчетов.
Этот способ принят в качестве прототипа. К недостаткам прототипа можно отнести следующее.
Сигналы с ограниченным спектром используются только в технических системах. Реальные непрерывные сигналы от естественных источников окружающей среды обычно имеют бесконечный, хотя и затухающий, спектр. Если допустим определенный уровень погрешности при восстановлении сигнала, возможно ограничение бесконечного спектра некоторой частотой среза. Если это ограничение до половины частоты дискретизации не производится низкочастотной фильтрацией перед дискретизацией, то возникает известное явление наложения частот (апиасинг), приводящее к не устраняемой составляющей погрешности дискретизации при восстановлении непрерывного сигнала любыми временными интерполяторами, в том числе и рядом Котельникова [Ж. Макс. Методы и техника обработки сигналов при физических измерениях: В 2-х т. Пер. с франц. - М.: Мир, 1983. - Т.1. - 312 с.].
В системах циклического контроля окружающей среды (например, зондированием прибором толщи океана, полигонной съемкой региона последовательным обходом станций и т.д.) возможности низкочастотной фильтрации процессов и полей обычно ограничены или она в принципе не возможна. Одновременно ограничена сверху и частота временной дискретизации экономическими затратами или возможностями практической peaлизации. В этих случаях, когда возможности низкочастотной фильтрации и повышения частоты дискретизации исчерпаны, при априорно известных типовых бесконечных спектрах непрерывных сигналов можно оценить возникающую погрешность дискретизации от алиасинга [Гайский В.А., Гайский П.В. Погрешности дискретизации случайных многомерных полей со степенными спектрами. Морской гидрофизический журнал, №6, 1994]. Если эта погрешность оказывается выше допустимого уровня, то нет известных способов ее уменьшения.
Предлагаемый способ решает эту задачу.
Целью предлагаемого изобретения является снижение погрешности дискретизации, обусловленной алиасингом, путем определения параметров гармоник частот наложения и восстановления непрерывного сигнала по спектру в расширенной полосе частот.
Эта цель достигается тем, что при дискретизации сигнала x(t), максимальная значимая частота в спектре которого в m раз превышает половину частоты дискретизации цифровые отсчеты сигнала берут пакетами по m отсчетов с задержками τν от первого и интервалом времени Т между пакетами, последовательности отсчетов регистрируют или передают в канал связи, при восстановлении непрерывного сигнала из последовательности с N пакетами отсчетов формируют m последовательностей {xν(k)} с N одиночными отсчетами, (m-1) которых задержаны на интервал от первой дискретным преобразованием Фурье вычисляют спектр с алиасингом каждой из последовательностей {xν(k)} по формуле
для каждой из частот неизвестного спектра сигнала формируют суммы взвешенных задержками гармоник частот наложения по уравнению
при из последнего уравнения формируют систему с неизвестными при -m<i<m, решением системы уравнений находят неизвестные составляющие расширенного спектра исходного процесса восстанавливают непрерывный сигнал x(t) обратным преобразованием Фурье по формуле
Изобретение поясняется рисунками.
На фиг. 1 показаны схемы равномерной (а) и равномерной пакетной (б) дискретизации непрерывного сигнала мгновенными отсчетами, на фиг. 2 представлена схема алиасинга для положительных (ω) и отрицательных (-ω) частот, фиг. 3 иллюстрирует схему формирования рядов задержанных одиночных отсчетов, на фиг. 4 показаны, для примера, графики спектров и функций модуляции для пакета из 5 отсчетов с кратными задержками, на фиг. 5 показан график вертикального зондирования верхнего слоя моря непрерывного (1) и с остановками на некоторых горизонтах (2), на фиг. 6 представлен график выполнения циклического обхода станций гидрологического полигона с одним (1) и несколькими (2) зондированиями на станциях.
Суть предлагаемого способа дискретизации мгновенными отсчетами и восстановления непрерывного сигнала по этим отсчетам состоит в следующем.
Предполагается, что непрерывный процесс x(t) на отрезке времени TN с достаточной точностью представляется ограниченным частотой Ωc спектром
После дискретизации сигнала x(t) рядом {х(k)} (фиг. 1а) из N отсчетов с интервалом Т его спектр определяется дискретным преобразованием Фурье [Сергиенко А.Б. Цифровая обработка сигналов: Учебник для вузов. 2-е изд. - Спб.: Питер, 2006. - 751 с.: ил.]
Если частота дискретизации то алиасинга нет, дискретный спектр не содержит частот наложения и непрерывный сигнал x(t) восстанавливается по интерполяционной формуле Котельникова или тригонометрическим рядом Фурье.
Если с то имеет место алиасинг и вычисляемые в рабочей полосе частот спектральные составляющие содержат суммы гармоник частот наложения в соответствии с известной формулой
Эта формула наглядно иллюстрируется «пилой» алиасинга на фиг. 2, где горизонтальная ось соответствует шкале частот ω или Ω, а вертикальная ось направлением вниз соответствует шкале частот рабочей полосы
Для каждой частоты ω в рабочей полосе частоты гармоник наложения определяется точками пересечения горизонтальной прямой, проходящей через эту частоту на вертикальной оси с «пилой» алиасинга так, как это показано на фиг. 2.
Для нахождения спектральных составляющих необходимо придать им вес и сформировать дополнительные уравнения, это обеспечивается задержками.
на τν, то при дискретизации берут пакет из т отсчетов с задержками от первого (фиг. 1б). При восстановлении формируют m отдельных ν, k последовательностей {xν(k)} с единичными отсчетами, интервалом дискретизации Т и сдвигами (фиг. 3).
Для каждой ν-й последовательности {xν(k)} вычисляют дискретный спектр
и записывают уравнение с наложением частот, аналогичное выражению (3),
При этом наложение частот происходит с модуляцией функцией Gν(ω), графики для вещественных и для мнимых составляющих которой, показаны на фиг. 4, для примера, из 5 отсчетов с кратными задержками.
Для каждой ω формируют из выражения (5) при систему линейных алгебраических уравнений с неизвестными спектральными компонентами исходного непрерывного процесса x(t), расширенные матрицы которых имеют вид
для ω
для ω
Решением систем (6) и (7) находят компоненты расширенного спектра исходного сигнала, которых стало mN вместо N, что эквивалентно увеличению частоты дискретизации в m раз.
Используя обратное преобразование Фурье, восстанавливают непрерывный сигнал x(t) по конечному ряду со сколь угодно высоким временным разрешением по формуле
содержащей конечное число членов, в отличие от формулы Котельникова.
Таким образом, для реализации предложенного способа дискретизации и восстановления непрерывного сигнала необходимо обеспечить дискретизацию пакетами. При этом общее число отсчетов равно числу отсчетов при увеличении частоты дискретизации в m раз.
Рассмотрим примеры использования способа.
В многоканальных системах аналого-цифрового преобразования с одним аналого-цифровым преобразователем и коммутатором на N каналов время цикла равно (период дискретизации одного сигнала) (фиг. 1а)
где τk - время коммутации одного канала;
τn - время одного аналого-цифрового преобразования. Если τn<τk, τn=ατk, α<1, что характерно при использовании квазиследящих или сигма-дельта АЦП, то для повышения быстродействия целесообразно использовать пакетную дискретизацию (фиг. 1б). В этом случае при m отсчетах в пакете эквивалентная частота дискретизации возрастает в m раз, но и период опроса возрастает и будет равен
Выгрыш в увеличении эквивалентной частоты дискретизации составит
Следовательно, при реальных значениях времен преобразования и коммутации предложенный способ пакетной дискретизации и восстановления дает выигрыш в эквивалентной частоте дискретизации.
В океанографии контроль морской среды осуществляется путем выполнения с помощью экспедиционного судна повторяющегося полигона станций с вертикальным зондированием верхнего слоя моря и измерениями по траектории зонда. График такого зондирования показан на фиг. 5(1). Если выполняется одно зондирование на станции, то частота временной дискретизации каждой точки среды на фиксированной глубине составит где T1 - периодичность съемок, которая в лучшем случае не будет меньше времени выполнения полигона. Это время фактически составляет время сканирования пространства, поскольку измерения выполняются параллельно и, в принципе, время одного измерения может быть существенно меньше времени перемещения между двумя точками пространства на расстояние требуемого шага пространственной дискретизации.
Если озаботиться при этом эквивалентной частотой временной дискретизации в конкретных точках пространства (например, на «стандартных» горизонтах гидрологической станции), то это можно выполнить с использованием предложенного способа дискретизации и восстановления непрерывного сигнала путем реализации пакета отсчетов в этих точках при остановке зонда на стандартных горизонтах (фиг. 5) или выполнением нескольких зондирований на каждой станции полигона (фиг. 6).
На фиг. 5 показаны графики многократного зондирования верхнего слоя моря непрерывного (1) и с остановками на некоторых горизонтах (2). В первом случае измерение в каждой точке осуществляется при движении вниз параллельно с движением зонда и затрат времени на него нет. Время между отсчетами равно τ3 и частота дискретизации равна Для выбранного числа горизонтов g (например, при принятых в гидрологии стандартных глубин g=16) возможно увеличение эквивалентной частоты дискретизации, если осуществить остановку зонда на этих горизонтах для взятия пакета из m измерений за время При этом период дискретизации будет равен а эквивалентная частота дискретизации будет равна при α<<1. Выигрыш в эквивалентной частоте дискретизации составит
Таким образом, практически всегда возможно использовать предлагаемый способ дискретизации и восстановления непрерывного сигнала на выбранных горизонтах океанографической станции для расширения полосы рабочих частот.
Второй пример использования предложенного способа при океанографической съемке путем циклического обхода станций полигона с непрерывным зондированием на станциях одним (1) и несколькими зондированиями (2) показан на фиг. 6.
При этом период дискретизации с одним зондированием составит с m зондированиями - , где - время перехода между станциями; обычно α<1.
Выигрыш в эквивалентной частоте дискретизации составит
Эта формула аналогична формуле (9), поэтому выводы (10) для нее справедливы.
Claims (7)
- Способ дискретизации и восстановления непрерывного сигнала x(t) с использованием его цифровых отсчетов через равные интервалы времени Т, отличающийся тем, что при дискретизации отсчеты берут пакетами по m отсчетов с задержками от первого, последовательность из N пакетов отсчетов регистрируют или передают в канал связи, при восстановлении из последовательности с пакетами отсчетов формируют m последовательностей одиночных отсчетов {х(k)}, задержанных на интервал времени τν(τ1=0) от первой, дискретным преобразованием Фурье, вычисляют спектр с алиасингом каждой из последовательностей {xν(k)} по формуле
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
RU2015129169A RU2660320C2 (ru) | 2015-07-16 | 2015-07-16 | Способ дискретизации и восстановления непрерывного сигнала |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
RU2015129169A RU2660320C2 (ru) | 2015-07-16 | 2015-07-16 | Способ дискретизации и восстановления непрерывного сигнала |
Publications (3)
Publication Number | Publication Date |
---|---|
RU2015129169A RU2015129169A (ru) | 2017-01-25 |
RU2015129169A3 RU2015129169A3 (ru) | 2018-03-13 |
RU2660320C2 true RU2660320C2 (ru) | 2018-07-05 |
Family
ID=58451009
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
RU2015129169A RU2660320C2 (ru) | 2015-07-16 | 2015-07-16 | Способ дискретизации и восстановления непрерывного сигнала |
Country Status (1)
Country | Link |
---|---|
RU (1) | RU2660320C2 (ru) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4121295A (en) * | 1977-04-07 | 1978-10-17 | Wittronics, Inc. | Integer weighted impulse equivalent coded signal processing apparatus |
SU1288726A2 (ru) * | 1984-06-26 | 1987-02-07 | Ленинградский Электротехнический Институт Связи Им.Проф.М.А.Бонч-Бруевича | Устройство дл восстановлени непрерывных функций по дискретным отсчетам |
RU2051414C1 (ru) * | 1991-04-09 | 1995-12-27 | Алексей Веденеевич Заворин | Способ формирования сигналов произвольной формы |
RU2259591C2 (ru) * | 2003-05-05 | 2005-08-27 | Государственное образовательное учреждение высшего профессионального образования "Оренбургский государственный университет" | Способ восстановления исходных сигналов |
-
2015
- 2015-07-16 RU RU2015129169A patent/RU2660320C2/ru not_active IP Right Cessation
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4121295A (en) * | 1977-04-07 | 1978-10-17 | Wittronics, Inc. | Integer weighted impulse equivalent coded signal processing apparatus |
SU1288726A2 (ru) * | 1984-06-26 | 1987-02-07 | Ленинградский Электротехнический Институт Связи Им.Проф.М.А.Бонч-Бруевича | Устройство дл восстановлени непрерывных функций по дискретным отсчетам |
RU2051414C1 (ru) * | 1991-04-09 | 1995-12-27 | Алексей Веденеевич Заворин | Способ формирования сигналов произвольной формы |
RU2259591C2 (ru) * | 2003-05-05 | 2005-08-27 | Государственное образовательное учреждение высшего профессионального образования "Оренбургский государственный университет" | Способ восстановления исходных сигналов |
Also Published As
Publication number | Publication date |
---|---|
RU2015129169A3 (ru) | 2018-03-13 |
RU2015129169A (ru) | 2017-01-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Kânoğlu et al. | Initial value problem solution of nonlinear shallow water-wave equations | |
Bahreinimotlagh et al. | Application of shallow-water acoustic tomography to measure flow direction and river discharge | |
JP5007391B2 (ja) | 津波波源推定方法及び津波波高予測方法並びにその関連技術 | |
Taniguchi et al. | Measuring the Kuroshio Current with ocean acoustic tomography | |
JP5507903B2 (ja) | 震度推定方法及び装置 | |
CN112985639A (zh) | 基于去啁啾和时域子啁啾脉冲提取的分布式光纤传感方法 | |
RU170068U1 (ru) | Адаптивное устройство для подавления помех | |
Chi et al. | Utilizing orthogonal coprime signals for improving broadband acoustic Doppler current profilers | |
RU2549207C2 (ru) | Устройство обнаружения шумовых гидроакустических сигналов на основе квадратурного приемника | |
CN103353610A (zh) | 基于广义s变换子波能量衰减特征检测烃类的方法及装置 | |
RU2660320C2 (ru) | Способ дискретизации и восстановления непрерывного сигнала | |
Pawłowski | Reconstruction of input signal of sensor with frequency output | |
RU2460093C1 (ru) | Способ измерения дистанции гидролокатором | |
Peng et al. | Application of a Threshold‐Based Earthquake Early Warning Method to the M w 6.6 Lushan Earthquake, Sichuan, China | |
RU2524673C1 (ru) | Способ измерения вариаций фазового сдвига центральной частоты сигналов с минимальной частотной манипуляцией и устройство для его осуществления | |
RU2346300C1 (ru) | Способ обнаружения возможности наступления катастрофических явлений | |
CN101533045A (zh) | 一种中子脉冲序列的频谱分析方法 | |
CN107870034B (zh) | 一种基于相位差的水声声速测量方法 | |
RU2482517C1 (ru) | Трассопоисковый приемник | |
RU2593622C1 (ru) | Способ измерения радиальной скорости объекта по его шумоизлучению | |
Kundu et al. | Artificial neural network based estimation of moment magnitude with relevance to Earthquake Early Warning | |
RU2550757C1 (ru) | Устройство обнаружения шумовых гидроакустических сигналов на основе квадратурного приемника | |
Kandula et al. | Field testing of indirect displacement estimation using accelerometers | |
RU2650747C1 (ru) | Способ и устройство определения места прохождения трубопровода | |
Saulig et al. | Nonstationary signals information content estimation based on the local Rényi entropy in the time-frequency domain |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
MM4A | The patent is invalid due to non-payment of fees |
Effective date: 20200717 |