RU2645445C2 - Дробно-дифференцирующий электрический фильтр порядка 1/2 по методу Прони - Google Patents
Дробно-дифференцирующий электрический фильтр порядка 1/2 по методу Прони Download PDFInfo
- Publication number
- RU2645445C2 RU2645445C2 RU2015106763A RU2015106763A RU2645445C2 RU 2645445 C2 RU2645445 C2 RU 2645445C2 RU 2015106763 A RU2015106763 A RU 2015106763A RU 2015106763 A RU2015106763 A RU 2015106763A RU 2645445 C2 RU2645445 C2 RU 2645445C2
- Authority
- RU
- Russia
- Prior art keywords
- inverting amplifier
- fractional
- inputs
- filter
- repeaters
- Prior art date
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
- G01R23/16—Spectrum analysis; Fourier analysis
- G01R23/165—Spectrum analysis; Fourier analysis using filters
Landscapes
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- General Physics & Mathematics (AREA)
- Networks Using Active Elements (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
Изобретение относится к области радиотехники и радиолокации и может быть использовано для оперативного контроля средней частоты по критерию центра тяжести энергетического спектра широкополосных доплеровских радиосигналов во временной области без спектральной обработки. Изобретение представляет собой дробно-дифференцирующий электрический фильтр порядка 1/2 по методу Прони, состоящий из М интегрирующих цепей, повторителей, инвертирующего усилителя и инвертирующего усилителя-сумматора. Причем входы интегрирующих цепей с соответствующими постоянными времени и инвертирующего усилителя с заданным коэффициентом усиления объединены и подключены к входу дробно-дифференцирующего фильтра, выходы интегрирующих цепей подключены к входам повторителей, а выходы повторителей и инвертирующего усилителя с заданным коэффициентом усиления подключены к входам инвертирующего усилителя-сумматора, выход которого подключен к выходу дробно-дифференцирующего фильтра. Технический результат заключается в повышении точности и скорости измерения средней частоты. 8 ил.
Description
Изобретение относится к области радиотехники и радиолокации и может быть использовано для оперативного контроля средней частоты по критерию центра тяжести энергетического спектра широкополосных доплеровских радиосигналов во временной области без спектральной обработки.
Из существующего уровня техники известен метод дискретного счета [1], используемый для измерения частоты F электронно-счетными частотомерами. Принцип работы метода дискретного счета основан на подсчете числа N пересечений нулевого уровня сигналом x(t) за временной интервал Т и формировании отношения FT=N/T (квазичастота). Метод дискретного счета прост и хорошо работает в реальном масштабе времени при измерении частоты монохроматических сигналов.
В случае широкополосных радиосигналов (например, доплеровских сигналов в радиолокации) [2] используется понятие «средней частоты» - центра тяжести физического спектра сигнала:
где - 
спектральная плотность сигнала; ω=2πF - круговая частота. Однако, при измерении частоты широкополосных радиосигналов, носящих случайный характер, метод дискретного счета дает завышенную оценку [3] в сравнении с (1), что является его существенным недостатком.
Известен способ оценки средней частоты ω0 [2], в котором средняя частота рассчитывается путем непосредственного вычисления по соотношению (1); при этом необходим предварительный расчет спектра сигнала (спектральный анализ) с помощью алгоритмов быстрого преобразования Фурье. Однако спектральный анализ требует большого объема вычислений и значительно увеличивает время получения результата, что в ряде случаев (например, в ближней радиолокации) неприемлемо.
Наиболее близким к заявляемому техническому решению является способ оценки средней частоты широкополосных доплеровских сигналов [4], используемый в качестве прототипа и основанный на вычислении центра тяжести спектра во временной области по алгоритму:
схематически изображенному на фиг. 1, где D1/2 - оператор дробного дифференцирования, представляющий собой некоторый линейный стационарный оператор во временной области, что дает возможность реализовать его в виде электрического фильтра и получать оценку средней частоты в реальном масштабе времени без спектральной обработки. Частотная характеристика этого фильтра должна описываться выражением
Недостатком данного способа является то, что частотная характеристика фильтра вида (3) абсолютно неинтегрируема и реализовать ее в виде электрической схемы (четырехполюсника) технически невозможно.
В основу изобретения положена задача, заключающаяся в разработке дробно-дифференцирующего электрического фильтра порядка 1/2 по методу Прони, лишенного вышеуказанных недостатков и в котором обеспечивается возможность реализации алгоритма дробного дифференцирования доплеровских радиосигналов с помощью типовых RC схем в реальном масштабе времени. Это, в свою очередь, позволяет реализовать схему измерения центра тяжести спектра сигнала, позволяющую производить оценку средней частоты энергетического спектра широкополосных радиосигналов во временной области без спектральной обработки, тем самым повысить точность и скорость измерения средней частоты.
Указанный технический результат достигается тем, что в состав дробно-дифференцирующего электрического фильтра порядка 1/2 по методу Прони включены М интегрирующих цепей, повторителей, инвертирующий усилитель и инвертирующий усилитель-сумматор, причем входы интегрирующих цепей с соответствующими постоянными времени и инвертирующего усилителя с заданным коэффициентом усиления объединены и подключены ко входу дробно-дифференцирующего фильтра, выходы интегрирующих цепей подключены ко входам повторителей, а выходы повторителей и инвертирующего усилителя с заданным коэффициентом усиления подключены ко входам инвертирующего усилителя-сумматора, выход которого подключен к выходу дробно-дифференцирующего фильтра.
Сущность изобретения поясняется чертежами, на которых представлено;
На фиг. 1 - структурная схема измерения центра тяжести спектра сигнала x(t) по методу дробного дифференцирования.
На фиг. 2 - импульсная характеристика дробно-дифференцирующего фильтра порядка 1/2.
На фиг. 3 - структурная схема дробно-дифференцирующего фильтра порядка 1/2 по методу Прони.
На фиг. 4-7 - основные электрические звенья, используемые для реализации дробно-дифференцирующего фильтра: интегрирующая цепь (фиг. 4); повторитель (фиг. 5); инвертирующий усилитель(фиг. 6); инвертирующий усилитель-сумматор (фиг. 7).
На фиг. 8 - схема дробно-диффференцирующего электрического фильтра порядка 1/2 по методу Прони.
Известно [5], что собственная функция оператора D1/2 (импульсная характеристика дробно-дифференцирующего фильтра) описывается соотношением
где δ(t) - дельта-функция Дирака, σ(t) - функция включения Хевисайда. Функция h1/2(t) представлена на фиг. 2 и удовлетворяет условию
В описываемом устройстве эту характеристику предлагается формировать конечным числом звеньев с экспоненциальными импульсными характеристиками. Для этой цели величина ε в соотношении (4) полагается достаточно малой, но конечной ε<<1/ωmax, где ωmax - максимальная частота спектра сигнала x(t), и используется разложение h1/2(t) для t>0 в ряд Прони [6] по экспонентам:
где М - прядок разложения, определяемый требуемой точностью аппроксимации. В дальнейшем фильтр, импульсная характеристика которого при t>0 описывается суммой экспонент, полученных разложением в ряд Прони (6), будем называть «фильтром Прони».
Структура дробно-дифференцирующего фильтра Прони представлена на фиг. 3 и содержит М динамических элементов с импульсными характеристиками hm(t)=σ(t)ехр(λmt) соответственно, М+1 усилительных элементов с коэффициентами усиления h0,А1,А2,…Ам, сигналы с выходов которых суммируются. Число динамических элементов М определяется требуемой точностью аппроксимации импульсной характеристики h1/2(t) при t>0.
Коэффициенты Am и параметры λm определяются разложением в ряд Прони [6]. Значение коэффициента h0 определяется из условия (5):
Импульсная характеристика фильтра аппроксимируется суммой экспонент, каждое из слагаемых которой реализуется в виде электрического звена первого порядка.
Для схемной реализации дробно-дифференцирующего фильтра предлагается использовать наиболее распространенные электрические звенья [7], обладающие передаточными характеристиками:
- интегрирующая цепь K(р)=(1+pR1C)-1
- повторитель K(р)=1
- инвертирующий усилитель K(р)=-R/R1.
- инвертирующий усилитель-сумматор Km(р)=-R/R2,m, m=l, M,
где М - число входов сумматора.
Осуществление изобретения
Дробно-дифференцирующий фильтр порядка 1/2 по методу Прони состоит из М электрических звеньев: интегрирующих цепей, повторителей, инвертирующего усилителя и инвертирующего усилителя-сумматора, причем входы интегрирующих цепей с параметрами R1,mC=-1/λm и инвертирующего усилителя с коэффициентом усиления h0 объединены и подключены ко входу фильтра, выходы интегрирующих цепей подключены ко входам повторителей, а выходы повторителей и инвертирующего усилителя с коэффициентом усиления h0 подключены ко входам инвертирующего усилителя-сумматора, выход которого подключен к выходу фильтра.
На фиг. 4-7 изображены электрические схемы звеньев на базе операционных усилителей (ОУ), реализующих структуру дробно-дифференцирующего фильтра, представленную на фиг. 3. Соответствующие импульсные характеристики звеньев имеют вид:
- интегрирующая цепь (фиг. 4) 
;
- повторитель на ОУ (фиг. 5) 
;
- инвертирующий усилитель на ОУ (фиг. 6) 
;
- инвертирующий усилитель-сумматор на ОУ (фиг. 7) 
;
m=1, M.
На фиг. 8 представлена схема дробно-дифференцирующего фильтра порядка 1/2 по методу Прони.
Схема состоит из М параллельно включенных ветвей (m=1,М), на входы которых подается сигнал U1(t), а выходы суммируются на входе инвертирующего усилителя-сумматора 10, выход которого является выходом схемы.
Каждая из ветвей представляет собой последовательное включение интегрирующей RC цепи 8 и повторителя 9 на ОУ, служащего для развязки выхода интегрирующей цепи от остальной части схемы. Параметры m-й RC цепи определяются из условия R1.mC=-1/λm.
«Нулевая» ветвь (m=0) состоит из инвертирующего усилителя (фиг. 6), в схеме которого R=R1.
Работает устройство следующим образом.
Входной сигнал U1(t) поступает одновременно на интегрирующие цепи 8 и инвертирующий усилитель 7. Выходные сигналы интегрирующих цепей 8 поступают на входы повторителей напряжения 9, служащих для развязки. Выходные напряжения операционных усилителей, включая «нулевую» ветвь 7 (см. фиг. 8), подаются на инвертирующий усилитель-сумматор 10, где суммируются с коэффициентами Am. Величина сопротивлений R2.m для m=l, M на схеме фиг. 8 выбирается из условия
Поскольку все значения Am, λm в разложении (6) получаются отрицательными, значения сопротивлений R2.m>0. Величина сопротивления R2.0 в «нулевой» ветви 7 на схеме фиг. 8 выбирается из условия
Таким образом, электрическая схема, изображенная на фиг. 8, реализует импульсную характеристику
аппроксимирующую характеристику (4) дробно-дифференцирующего фильтра порядка 1/2.
Источники информации
1. Кушнир Ф.В. Радиотехнические измерения: Учебник для техникумов связи. -М.: Связь, 1980. - 176 с.
2. Гоноровский И.С. Радиотехнические цепи и сигналы. - М.: Радио и связь, 1986. - 512 с.
3. Баскаков С.И. Радиотехнические цепи и сигналы. - М.: Высшая школа, 1988. - 448 с.
4. Захарченко В.Д. Способ оценки средней частоты широкополосных доплеровских сигналов. - Патент на изобретение (РФ) №2114440 от 12.05.1995. Опубл. 27.06.1998. //Изобретения. Заявки и патенты. 1998, №18(11), с. 344. (прототип).
5. Zakharchenko V.D., Kovalenko I.G. On protecting the planet against cosmic attack: ultrafast real-time estimate of the asteroid's radial velocity //Acta Astronautica. - 2014. - V. 98C. - P. 158-162; DOI:10.1016/j.actaastro. 2014.02.002
6. Марпл. - мл. С.Л. Цифровой спектральный анализ и его приложения: Пер. с англ. - М: Мир, 1990. - 584 с.
7. Титце У., Шенк К. Полупроводниковая схемотехника: Справочное руководство. Пер. с нем. - М: Мир, 1983. - 512 с.
Claims (1)
- Дробно-дифференцирующий электрический фильтр порядка 1/2 по методу Прони, состоящий из М интегрирующих цепей, повторителей, инвертирующего усилителя и инвертирующего усилителя-сумматора, отличающийся тем, что входы интегрирующих цепей с соответствующими постоянными времени и инвертирующего усилителя с заданным коэффициентом усиления объединены и подключены к входу дробно-дифференцирующего фильтра, выходы интегрирующих цепей подключены к входам повторителей, а выходы повторителей и инвертирующего усилителя с заданным коэффициентом усиления подключены к входам инвертирующего усилителя-сумматора, выход которого подключен к выходу дробно-дифференцирующего фильтра.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| RU2015106763A RU2645445C2 (ru) | 2015-02-26 | 2015-02-26 | Дробно-дифференцирующий электрический фильтр порядка 1/2 по методу Прони |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| RU2015106763A RU2645445C2 (ru) | 2015-02-26 | 2015-02-26 | Дробно-дифференцирующий электрический фильтр порядка 1/2 по методу Прони |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| RU2015106763A RU2015106763A (ru) | 2016-09-20 |
| RU2645445C2 true RU2645445C2 (ru) | 2018-02-21 |
Family
ID=56891839
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| RU2015106763A RU2645445C2 (ru) | 2015-02-26 | 2015-02-26 | Дробно-дифференцирующий электрический фильтр порядка 1/2 по методу Прони |
Country Status (1)
| Country | Link |
|---|---|
| RU (1) | RU2645445C2 (ru) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| SU720695A1 (ru) * | 1977-12-02 | 1980-03-05 | Предприятие П/Я Р-6155 | Дробный фильтр нижних частот |
| RU2114440C1 (ru) * | 1995-05-12 | 1998-06-27 | Волгоградский государственный университет | Способ оценки средней частоты широкополосных доплеровских сигналов |
| EP0971516A2 (en) * | 1998-07-06 | 2000-01-12 | Metalink Limited | Channel precoder using an infinite impulse response filter |
| US8014541B1 (en) * | 2004-10-08 | 2011-09-06 | Kind of Loud Technologies, LLC. | Method and system for audio filtering |
-
2015
- 2015-02-26 RU RU2015106763A patent/RU2645445C2/ru not_active IP Right Cessation
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| SU720695A1 (ru) * | 1977-12-02 | 1980-03-05 | Предприятие П/Я Р-6155 | Дробный фильтр нижних частот |
| RU2114440C1 (ru) * | 1995-05-12 | 1998-06-27 | Волгоградский государственный университет | Способ оценки средней частоты широкополосных доплеровских сигналов |
| EP0971516A2 (en) * | 1998-07-06 | 2000-01-12 | Metalink Limited | Channel precoder using an infinite impulse response filter |
| US8014541B1 (en) * | 2004-10-08 | 2011-09-06 | Kind of Loud Technologies, LLC. | Method and system for audio filtering |
Also Published As
| Publication number | Publication date |
|---|---|
| RU2015106763A (ru) | 2016-09-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP2797226B1 (en) | Signal playback device and signal playback method | |
| He et al. | Parallel calibration method for phased array with harmonic characteristic analysis | |
| RU170068U1 (ru) | Адаптивное устройство для подавления помех | |
| RU2645445C2 (ru) | Дробно-дифференцирующий электрический фильтр порядка 1/2 по методу Прони | |
| Ding et al. | Parameter identification of multi-input, single-output systems based on FIR models and least squares principle | |
| Bai | Frequency domain identification of Hammerstein models | |
| Kochańska | Adaptive identification of time-varying impulse response of underwater acoustic communication channel | |
| JP5797135B2 (ja) | 濾波装置および濾波方法 | |
| US7561639B2 (en) | Method and device for estimating channel properties of a transmission channel | |
| Punchalard | Steady-state analysis of a complex adaptive notch filter using modified gradient algorithm | |
| RU2280878C2 (ru) | Способ статистической оценки искажений, вносимых четырехполюсником, и двухканальный статистический анализатор (варианты) | |
| JP2013246107A (ja) | 信号処理装置および信号処理方法 | |
| RU2541919C1 (ru) | Способ повышения точности аппроксимации при выделении полезного сигнала в условиях априорной неопределенности и устройство, его реализующее | |
| Goos et al. | Continuous time frequency domain LPV state space identification via periodic time-varying input-output modeling | |
| Thakral | Design of low power current mode square-root circuit | |
| Shali et al. | Approximate solutions of the nonlinear stochastic differential equations | |
| RU2552117C2 (ru) | Устройство для оценки частотных искажений случайных сигналов (варианты) | |
| RU2321004C2 (ru) | Устройство для статистической оценки нелинейных искажений | |
| RU167937U1 (ru) | Мост для измерения параметров многоэлементного произвольного двухполюсника | |
| Tanc | Sparsity regularized recursive total least-squares | |
| Weatherspoon et al. | Accurate and efficient modeling of FET cold noise sources using ANNs | |
| RU2321005C2 (ru) | Устройство для статистической оценки нелинейных искажений | |
| Gower et al. | Analysis of sensor impulse response effects on Cramèr–Rao lower bounds for signal parameter estimators | |
| RU164241U1 (ru) | Устройство управления коэффициентом усиления передающего устройства рлс на основе восстановления параметров усиливаемого сигнала | |
| Oke | On Adomian Decomposition Method for solving general wave equations on transmission lines |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| MM4A | The patent is invalid due to non-payment of fees |
Effective date: 20180312 |