RU2608363C1 - Способ оценки параметров модели замирания огибающей сигнала по закону накагами по информационному многочастотному сигналу - Google Patents
Способ оценки параметров модели замирания огибающей сигнала по закону накагами по информационному многочастотному сигналу Download PDFInfo
- Publication number
- RU2608363C1 RU2608363C1 RU2015129277A RU2015129277A RU2608363C1 RU 2608363 C1 RU2608363 C1 RU 2608363C1 RU 2015129277 A RU2015129277 A RU 2015129277A RU 2015129277 A RU2015129277 A RU 2015129277A RU 2608363 C1 RU2608363 C1 RU 2608363C1
- Authority
- RU
- Russia
- Prior art keywords
- signal
- amplitude
- input
- noise
- transmitted
- Prior art date
Links
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B17/00—Monitoring; Testing
- H04B17/30—Monitoring; Testing of propagation channels
- H04B17/391—Modelling the propagation channel
- H04B17/3911—Fading models or fading generators
Landscapes
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Monitoring And Testing Of Transmission In General (AREA)
- Mobile Radio Communication Systems (AREA)
Abstract
Изобретение относится к области электрорадиотехники и связи и может быть использовано в системах передачи данных, использующих многочастотные сигналы с ортогональным частотным разделением каналов для оценки параметров канала связи. Техническим результатом заявленного изобретения является обеспечение более точного определения параметров модели замирания огибающей сигнала по закону Накагами по информационному многочастотному сигналу в случае наличия на приемной стороне блока автоматической регулировки усиления. Кроме того, данный способ не требует наличия тестового сигнала. Технический результат достигается благодаря тому, что в предложенном способе осуществляется измерение значений амплитуды смеси сигнала и шума на частотах, используемых для передачи информационного сигнала, и значений амплитуды шума на частотах, неиспользуемых для передачи информационного сигнала, и используется аналитическое выражение для плотности случайной величины, равной отношению измеренных величин. 1 ил.
Description
Изобретение относится к области электрорадиотехники и связи и может быть использовано в системах передачи данных, использующих многочастотные сигналы с ортогональным частотным разделением каналов, для оценки параметров канала связи.
Для обеспечения стабильной работы системы передачи данных необходимо осуществлять контроль качества используемого канала связи. Критерием качества канала в цифровых системах связи является вероятность ошибки на бит, которая однозначно связана с параметрами модели замираний. Поэтому актуальной является задача определения параметров модели замирания огибающей сигнала по закону Накагами по результатам анализа информационного многочастотного сигнала.
Наиболее близким к заявленному техническому решению является способ, описанный в [Y. Chen, N.C. Beaulieu, C. Tellambura. "Novel Nakagami-m Parameter Estimator for Noisy Channel Samples", IEEE Communications Letters, vol. 9, no. 5, may 2005], который принят за прототип. Оценка параметров формируется с помощью анализа амплитуд полезного сигнала и шума на основе метода моментов.
Известный способ определения параметров распределения Накагами работает следующим образом. На приемной стороне оцифровывают принимаемый сигнал в аналогово-цифровом преобразователе (АЦП), затем передают оцифрованный сигнал с выхода АЦП на вход первого блока вычисления амплитуды, в котором определяют амплитуду принимаемого сигнала на частоте, используемой для передачи информационного сигнала на длительности элементарной посылки. Одновременно с выхода АЦП передают оцифрованный сигнал на вход второго блока вычисления амплитуды, в котором определяют амплитуду принимаемого сигнала на частоте, неиспользуемой для передачи информационного сигнала на длительности элементарной посылки, определяя, таким образом, амплитуду шума. Затем с выхода первого блока вычисления амплитуды вычисленное значение амплитуды передают на вход первого квадратора, в котором полученное значение возводят в квадрат. Далее с выхода первого квадратора полученное значение передают одновременно на вход первого блока накопления, в котором накапливают последние N значений. С выхода первого блока накопления накопленный массив значений передают на вход первого сумматора, в котором вычисляют сумму всех N значений массива. Затем с выхода первого сумматора вычисленное значение передают на вход первого делителя, в котором делят полученное значение на N, а результат деления m2 передают одновременно на первый вход блока нахождения среднего значения отношения сигнал/шум и на первый вход блока вычисления коэффициента В. При этом одновременно с выхода первого блока вычисления амплитуды вычисленное значение амплитуды передают на вход второго блока накопления, в котором накапливают последние N значений. С выхода второго блока накопления накопленный массив значений передают на вход второго сумматора, в котором вычисляют сумму всех N значений массива. Затем с выхода второго сумматора вычисленное значение передают на вход второго делителя, в котором делят полученное значение на N, а результат деления m1 передают на второй вход блока вычисления коэффициента В. При этом с выхода второго блока вычисления амплитуды вычисленное значение амплитуды шума передают на вход второго квадратора, в котором полученное значение возводят в квадрат. Далее с выхода второго квадратора полученное значение передают одновременно на вход третьего блока накопления, в котором накапливают последние N значений. С выхода третьего блока накопления накопленный массив значений передают на вход третьего сумматора, в котором вычисляют сумму всех N значений массива. Затем с выхода третьего сумматора вычисленное значение передают на вход третьего делителя, в котором делят полученное значение на N, а результат деления S передают на второй вход блока нахождения среднего значения отношения сигнал/шум. При этом в блоке нахождения среднего значения отношения сигнал/шум делят полученное значение в первого входа на значение, полученное со второго входа и вычитают из результата деления единицу, получая таким образом среднее значение отношения сигнал/шум h. Полученное среднее значение отношения сигнал/шум с выхода блока нахождения среднего значения отношения сигнал/шум передают одновременно на вход блока вычисления коэффициента А, на вход блока вычисления коэффициента С и на третий вход блока вычисления коэффициента В. При этом в блоке вычисления коэффициента А определяют коэффициент А по формуле , в блоке вычисления коэффициента В определяют коэффициент В по формуле , в блоке вычисления коэффициента С определяют коэффициент С по формуле , где коэффициенты , , , , , , , , , являются коэффициентами полинома, определяются заранее и равны , , , , , , , , . Далее с выхода блока вычисления коэффициента А передают значение коэффициента А на первый вход блока нахождения параметров распределения, также с выхода блока вычисления коэффициента В передают значение коэффициента В на второй вход блока нахождения параметров распределения, а с выхода блока вычисления коэффициента С передают значение коэффициента С на третий вход блока нахождения параметров распределения. В блоке нахождения параметров распределения определяют параметр распределения m по формуле .
На приемной стороне реальных систем связи обычно присутствует блок автоматической регулировки усиления (АРУ) для приведения уровня входного сигнала к значению, обеспечивающему оптимальную работу АЦП. В таком случае, выборочная плотность амплитуды не будет являться соответствующей состоятельной оценкой истинной плотности распределения. Таким образом, недостатком прототипа является то, что он не учитывает наличие блока АРУ и получаемая данным способом оценка будет обладать большой погрешностью.
Целью изобретения является получение оценки параметров модели замирания радиоканала по закону Накагами путем анализа принимаемого информационного многочастотного сигнала.
Поставленная цель достигается тем, что способ оценки параметров модели замирания огибающей сигнала по закону Накагами по информационному многочастотному сигналу состоит в том, что на приемной стороне оцифровывают принимаемый сигнал в аналогово-цифровом преобразователе, затем передают оцифрованный сигнал с выхода аналогово-цифрового преобразователя на вход первого блока вычисления амплитуды, при этом в нем определяют значение амплитуды смеси принимаемого сигнала и шума на всех n частотах, используемых для передачи информационного сигнала на длительности элементарной посылки, а с n выходов первого блока вычисления амплитуды на первые входы n соответствующих делителей передают вычисленные значения амплитуд, также с выхода аналогово-цифрового преобразователя передают оцифрованный сигнал на вход второго блока вычисления амплитуды, в котором определяют значение амплитуды шума на n частотах, неиспользуемых для передачи информационного сигнала на длительности элементарной посылки, а с n выходов второго блока вычисления амплитуды передают вычисленные значения амплитуд на вторые входы n соответствующих делителей, а в каждом делителе осуществляют деление значения амплитуды шума на частоте, неиспользуемой для передачи информационного сигнала, полученное по второму входу на значение амплитуды смеси сигнала и шума на частоте, используемой для передачи информационного сигнала, полученное по первому входу, а результат деления передают с выходов n делителей на n соответствующих входов блока накопления, в котором накапливают выборку полученных n значений на длительности интервала анализа, равной М посылкам, получая, таким образом, выборку размером n×М значений, а с выхода блока накопления передают накопленный массив значений на вход блока вычисления параметров распределения, в котором, например, методом наискорейшего спуска, определяют параметры модели замирания огибающей сигнала по закону Накагами m и , являющиеся координатами максимума функции правдоподобия , где xi - это i-е значение выборки, - плотность распределения вероятности измеряемой случайной величины.
Структурная схема предложенного способа приведена на фиг. 1.
Способ основан на следующих предположениях.
В общем случае для определения плотности распределения огибающей сигнала в канале с замираниями, когда доступными для измерения являются только значения огибающей смеси сигнал + шум, можно использовать подход, заключающийся в том, чтобы по плотности распределения огибающей смеси сигнал + шум определить параметры распределения Райса. При этом восстановить истинную плотность распределения огибающей можно, используя выборочную плотность распределения огибающей смеси сигнал + шум, получаемую посредством измерений на приемной стороне.
В данном подходе следует учитывать техническую проблему, связанную с тем, что на приемной стороне чаще всего сигнал перед обработкой проходит через устройство автоматической регулировки усиления (АРУ). Поскольку коэффициент усиления АРУ неизвестен и динамически меняется в процессе измерений, статистические характеристики выборочной плотности распределения амплитуды сигнала значительно меняются и применение указанных выше способов напрямую дает неадекватные оценки.
Избавиться от указанной трудности при приеме сигнала с использованием АРУ можно, если для оценки параметров модели канала использовать выборку случайных величин, инвариантную к значению коэффициента усиления АРУ. В качестве такой случайной величины может быть использована случайная величина ξ, определяемая как отношение огибающих Ai и Aj, измеренных на длительности одной и той же элементарной посылки на различных субчастотах с номерами i и j:
ξ=Ai/Aj.
Такой подход можно реализовать, если информационный сигнал является многочастотным, и при этом часть субчастот не используются для передачи. Тогда на входе приемника на занятых субчастотах наблюдается смесь информационного сигнала с шумом, а на свободных только шум.
Для описания плотности распределения огибающей шума на свободных субчастотах при гипотезе, что шум является гауссовским, используется плотность распределения Рэлея:
Тогда в качестве Ai можно использовать измеренную огибающую шума, а в качестве Aj - огибающую смеси сигнал + шум.
В случае постоянного уровня информационного сигнала А на соответствующих субчастотах для модели гауссовского шума функцию распределения случайной величины ξ можно найти следующим образом:
Если уровень информационного сигнала А не постоянен, а подвержен замираниям и его плотность распределения WA(x) подчиняется закону Накагами, то в этом случае функция распределения случайной величины ξ можно найти следующим образом:
Тогда в новых обозначениях функция распределения случайной величины ξ имеет следующий вид:
Выражение для плотности при этом имеет следующий вид:
Сформировав выборку случайной величины ξ и имея аналитическое выражение для ее плотности распределения, можно воспользоваться методом максимального правдоподобия, как одним из методов оценки неизвестных параметров распределений. В данном случае неизвестными параметрами будут и m. Тогда функция правдоподобия L, определяется выражением:
где xi - значение случайной величины ξ, n*М - объем выборки.
Таким образом, приведенные аналитические выводы показывают, что с помощью предложенного способа можно определить параметры модели замирания огибающей сигнала по закону Накагами по информационному многочастотному сигналу. При этом необходимыми данными являются измеренные значения амплитуды смеси сигнала и шума на частотах, используемых для передачи информационного сигнала и значения амплитуды шума на частотах, неиспользуемых для передачи информационного сигнала.
Способ работает следующим образом.
На приемной стороне оцифровывают принимаемый сигнал в аналогово-цифровом преобразователе 1, затем передают оцифрованный сигнал с выхода аналогово-цифрового преобразователя 1 на вход первого блока вычисления амплитуды 2, в котором определяют значение амплитуды смеси принимаемого сигнала и шума на всех n частотах, используемых для передачи информационного сигнала на длительности элементарной посылки. С n выходов первого блока вычисления амплитуды 2 на первые входы n соответствующих делителей 4(1)…4(N) передают вычисленные значения амплитуд. При этом с выхода аналогово-цифрового преобразователя 1 также передают оцифрованный сигнал на вход второго блока вычисления амплитуды 3, в котором определяют значение амплитуды шума на n частотах, неиспользуемых для передачи информационного сигнала на длительности элементарной посылки, а с n выходов второго блока вычисления амплитуды 3 передают вычисленные значения амплитуд на вторые входы n соответствующих делителей 4(1)…4(N). В каждом делителе 4(1)…4(N) осуществляют деление значения амплитуды шума на частоте, неиспользуемой для передачи информационного сигнала, полученное по второму входу на значение амплитуды смеси сигнала и шума на частоте, используемой для передачи информационного сигнала, полученное по первому входу, а результат деления передают с выходов n делителей 4(1)…4(N) на n соответствующих входов блока накопления 5, в котором накапливают выборку полученных n значений на длительности интервала анализа, равной М посылкам, получая, таким образом, выборку размером n×М значений, а с выхода блока накопления 5 передают накопленный массив значений на вход блока вычисления параметров распределения 6, в котором, например, методом наискорейшего спуска, определяют параметры модели замирания огибающей сигнала по закону Накагами m и , являющиеся координатами максимума функции правдоподобия , где xi - это i-е значение выборки, - плотность распределения вероятности измеряемой случайной величины.
Предлагаемый способ может быть использован для систем связи, использующих сигналы с ортогональным многочастотным разделением каналов связи. Применение такого способа позволяет более точно определять параметры замирающего канала связи.
Предлагаемое устройство по сравнению с прототипом обладает следующим преимуществом: обеспечивает более точное определение параметров модели замирания огибающей сигнала по закону Накагами по информационному многочастотному сигналу в случае наличия на приемной стороне блока автоматической регулировки усиления.
Claims (1)
- Способ оценки параметров модели замирания огибающей сигнала по закону Накагами по информационному многочастотному сигналу, заключающийся в том, что на приемной стороне оцифровывают принимаемый сигнал в аналогово-цифровом преобразователе, затем передают оцифрованный сигнал с выхода аналогово-цифрового преобразователя одновременно на вход первого блока вычисления амплитуды и на вход второго блока вычисления амплитуды, отличающийся тем, что в первом блоке вычисления амплитуды определяют значение амплитуды смеси принимаемого сигнала и шума на всех n частотах, используемых для передачи информационного сигнала, на длительности элементарной посылки, а с n выходов первого блока вычисления амплитуды на первые входы n соответствующих делителей передают вычисленные значения амплитуд, также во втором блоке вычисления амплитуды определяют значение амплитуды шума на n частотах, неиспользуемых для передачи информационного сигнала, на длительности элементарной посылки, а с n выходов второго блока вычисления амплитуды передают вычисленные значения амплитуд на вторые входы n соответствующих делителей, в каждом делителе осуществляют деление значения амплитуды шума на частоте, неиспользуемой для передачи информационного сигнала, полученное по второму входу на значение амплитуды смеси сигнала и шума на частоте, используемой для передачи информационного сигнала полученное по первому входу, а результат деления передают с выхода каждого из n делителей на n соответствующих входов блока накопления, в котором накапливают выборку полученных n значений на длительности интервала анализа, равной М посылкам, получая, таким образом, выборку размером n×М значений, а с выхода блока накопления передают накопленный массив значений на вход блока вычисления параметров распределения, в котором методом наискорейшего спуска, определяют параметры модели замирания огибающей сигнала по закону Накагами m и , являющиеся координатами максимума функции правдоподобия , где xi - это i-е значение выборки, а - плотность распределения вероятности измеряемой случайной величины.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
RU2015129277A RU2608363C1 (ru) | 2015-07-16 | 2015-07-16 | Способ оценки параметров модели замирания огибающей сигнала по закону накагами по информационному многочастотному сигналу |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
RU2015129277A RU2608363C1 (ru) | 2015-07-16 | 2015-07-16 | Способ оценки параметров модели замирания огибающей сигнала по закону накагами по информационному многочастотному сигналу |
Publications (2)
Publication Number | Publication Date |
---|---|
RU2608363C1 true RU2608363C1 (ru) | 2017-01-18 |
RU2015129277A RU2015129277A (ru) | 2017-01-23 |
Family
ID=58450511
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
RU2015129277A RU2608363C1 (ru) | 2015-07-16 | 2015-07-16 | Способ оценки параметров модели замирания огибающей сигнала по закону накагами по информационному многочастотному сигналу |
Country Status (1)
Country | Link |
---|---|
RU (1) | RU2608363C1 (ru) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2706939C1 (ru) * | 2019-02-05 | 2019-11-21 | Акционерное общество "Российский институт мощного радиостроения" (АО "РИМР") | Способ оценки параметров модели замираний радиоканала по закону накагами по многочастотному сигналу |
RU2742695C1 (ru) * | 2020-07-17 | 2021-02-09 | Федеральное государственное автономное образовательное учреждение высшего образования «Национальный исследовательский Томский государственный университет» | Цифровой измеритель параметров случайных процессов с распределением Накагами |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SU1396149A2 (ru) * | 1986-11-24 | 1988-05-15 | Военная академия им.Ф.Э.Дзержинского | Устройство дл определени характеристик случайного процесса |
US7046963B2 (en) * | 2002-12-06 | 2006-05-16 | Texas Instruments Incorporated | Method and apparatus of signal estimation over generalized fading channel |
-
2015
- 2015-07-16 RU RU2015129277A patent/RU2608363C1/ru not_active IP Right Cessation
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SU1396149A2 (ru) * | 1986-11-24 | 1988-05-15 | Военная академия им.Ф.Э.Дзержинского | Устройство дл определени характеристик случайного процесса |
US7046963B2 (en) * | 2002-12-06 | 2006-05-16 | Texas Instruments Incorporated | Method and apparatus of signal estimation over generalized fading channel |
Non-Patent Citations (3)
Title |
---|
CHEN Y.et al "Novel Nakagami-m Parameter Estimator for Noisy Channel Samples", IEEE Communications Letters, vol. 9, no. 5, may 2005, p.417-419. * |
CHEN Y.et al "Novel Nakagami-m Parameter Estimator for Noisy Channel Samples", IEEE Communications Letters, vol. 9, no. 5, may 2005, p.417-419. CHENG J. et al " Maximum-Likehood Based Estimation jf the Nakagami m parameter", IEEE Communications Letters, vol. 5, no. 3, march 2001, p.101-103. * |
CHENG J. et al " Maximum-Likehood Based Estimation jf the Nakagami m parameter", IEEE Communications Letters, vol. 5, no. 3, march 2001, p.101-103. * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2706939C1 (ru) * | 2019-02-05 | 2019-11-21 | Акционерное общество "Российский институт мощного радиостроения" (АО "РИМР") | Способ оценки параметров модели замираний радиоканала по закону накагами по многочастотному сигналу |
RU2742695C1 (ru) * | 2020-07-17 | 2021-02-09 | Федеральное государственное автономное образовательное учреждение высшего образования «Национальный исследовательский Томский государственный университет» | Цифровой измеритель параметров случайных процессов с распределением Накагами |
Also Published As
Publication number | Publication date |
---|---|
RU2015129277A (ru) | 2017-01-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7729414B2 (en) | Guard interval analysis method and apparatus | |
US10749778B2 (en) | Communication system determining time of arrival using matching pursuit | |
Torrieri | The radiometer and its practical implementation | |
US20110294449A1 (en) | Signal-based gain control | |
RU2608363C1 (ru) | Способ оценки параметров модели замирания огибающей сигнала по закону накагами по информационному многочастотному сигналу | |
RU2706939C1 (ru) | Способ оценки параметров модели замираний радиоканала по закону накагами по многочастотному сигналу | |
US11677425B2 (en) | Systems and methods for decentralized link performance | |
RU2549207C2 (ru) | Устройство обнаружения шумовых гидроакустических сигналов на основе квадратурного приемника | |
KR100749446B1 (ko) | 직교주파수 분할다중화 시스템에서 초기 동기 이전의 자동이득조절 방법 및 장치 | |
CN105531600B (zh) | 无线网络中用于用户速度估计的时间分析 | |
RU2559734C1 (ru) | Способ определения параметров модели замирания радиоканала по закону райса по информационному многочастотному сигналу | |
RU167430U1 (ru) | Устройство оценки вероятности ошибки на бит для сигналов с восьмипозиционной фазовой модуляцией по четырехпозиционным сигналам | |
Knox et al. | AGC-based RF fingerprints in wireless sensor networks for authentication | |
Zabini et al. | Random sampling via sensor networks: Estimation accuracy vs. energy consumption | |
CN108718223B (zh) | 一种非合作信号的盲频谱感知方法 | |
RU175190U1 (ru) | Устройство оценки вероятности ошибки на бит для сигналов с шестнадцатипозиционной фазовой модуляцией по двухпозиционным сигналам | |
JP6373809B2 (ja) | 信号情報取得システムおよび信号情報取得方法 | |
RU155554U1 (ru) | Устройство оценки вероятности ошибки на бит для сигналов с восьмипозиционной фазовой модуляцией по двухпозиционным сигналам | |
RU2526283C2 (ru) | Способ определения вероятности ошибки на бит по флуктуациям фазы информационных сигналов | |
US11330394B2 (en) | Early warning method and device based on WIFI signal strength disturbance, and storage medium | |
Miao | Research of snr estimation and prediction method used in cognitive radio | |
JP2012149951A (ja) | レーダ装置 | |
CN113447708A (zh) | 非恒包络调制信号的功率测量方法、装置及电子设备 | |
KR101488787B1 (ko) | 이동통신 단말기에서 도플러 주파수 추정 방법 및 장치 | |
RU2768238C1 (ru) | Способ двухэтапной селекции спектральных компонент радиосигналов в многоканальной аппаратуре радиомониторинга |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
MM4A | The patent is invalid due to non-payment of fees |
Effective date: 20200717 |