RU2608363C1 - Способ оценки параметров модели замирания огибающей сигнала по закону накагами по информационному многочастотному сигналу - Google Patents

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

Info

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
Application number
RU2015129277A
Other languages
English (en)
Other versions
RU2015129277A (ru
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 RU2015129277A priority Critical patent/RU2608363C1/ru
Application granted granted Critical
Publication of RU2608363C1 publication Critical patent/RU2608363C1/ru
Publication of RU2015129277A publication Critical patent/RU2015129277A/ru

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. Полученное среднее значение отношения сигнал/шум с выхода блока нахождения среднего значения отношения сигнал/шум передают одновременно на вход блока вычисления коэффициента А, на вход блока вычисления коэффициента С и на третий вход блока вычисления коэффициента В. При этом в блоке вычисления коэффициента А определяют коэффициент А по формуле
Figure 00000001
, в блоке вычисления коэффициента В определяют коэффициент В по формуле
Figure 00000002
, в блоке вычисления коэффициента С определяют коэффициент С по формуле
Figure 00000003
, где коэффициенты
Figure 00000004
,
Figure 00000005
,
Figure 00000006
,
Figure 00000007
,
Figure 00000008
,
Figure 00000009
,
Figure 00000010
,
Figure 00000011
,
Figure 00000012
, являются коэффициентами полинома, определяются заранее и равны
Figure 00000013
,
Figure 00000014
,
Figure 00000015
,
Figure 00000016
,
Figure 00000017
,
Figure 00000018
,
Figure 00000019
,
Figure 00000020
,
Figure 00000021
. Далее с выхода блока вычисления коэффициента А передают значение коэффициента А на первый вход блока нахождения параметров распределения, также с выхода блока вычисления коэффициента В передают значение коэффициента В на второй вход блока нахождения параметров распределения, а с выхода блока вычисления коэффициента С передают значение коэффициента С на третий вход блока нахождения параметров распределения. В блоке нахождения параметров распределения определяют параметр распределения m по формуле
Figure 00000022
.
На приемной стороне реальных систем связи обычно присутствует блок автоматической регулировки усиления (АРУ) для приведения уровня входного сигнала к значению, обеспечивающему оптимальную работу АЦП. В таком случае, выборочная плотность амплитуды не будет являться соответствующей состоятельной оценкой истинной плотности распределения. Таким образом, недостатком прототипа является то, что он не учитывает наличие блока АРУ и получаемая данным способом оценка будет обладать большой погрешностью.
Целью изобретения является получение оценки параметров модели замирания радиоканала по закону Накагами путем анализа принимаемого информационного многочастотного сигнала.
Поставленная цель достигается тем, что способ оценки параметров модели замирания огибающей сигнала по закону Накагами по информационному многочастотному сигналу состоит в том, что на приемной стороне оцифровывают принимаемый сигнал в аналогово-цифровом преобразователе, затем передают оцифрованный сигнал с выхода аналогово-цифрового преобразователя на вход первого блока вычисления амплитуды, при этом в нем определяют значение амплитуды смеси принимаемого сигнала и шума на всех n частотах, используемых для передачи информационного сигнала на длительности элементарной посылки, а с n выходов первого блока вычисления амплитуды на первые входы n соответствующих делителей передают вычисленные значения амплитуд, также с выхода аналогово-цифрового преобразователя передают оцифрованный сигнал на вход второго блока вычисления амплитуды, в котором определяют значение амплитуды шума на n частотах, неиспользуемых для передачи информационного сигнала на длительности элементарной посылки, а с n выходов второго блока вычисления амплитуды передают вычисленные значения амплитуд на вторые входы n соответствующих делителей, а в каждом делителе осуществляют деление значения амплитуды шума на частоте, неиспользуемой для передачи информационного сигнала, полученное по второму входу на значение амплитуды смеси сигнала и шума на частоте, используемой для передачи информационного сигнала, полученное по первому входу, а результат деления передают с выходов n делителей на n соответствующих входов блока накопления, в котором накапливают выборку полученных n значений на длительности интервала анализа, равной М посылкам, получая, таким образом, выборку размером n×М значений, а с выхода блока накопления передают накопленный массив значений на вход блока вычисления параметров распределения, в котором, например, методом наискорейшего спуска, определяют параметры модели замирания огибающей сигнала по закону Накагами m и
Figure 00000023
, являющиеся координатами максимума функции правдоподобия
Figure 00000024
, где xi - это i-е значение выборки,
Figure 00000025
- плотность распределения вероятности измеряемой случайной величины.
Структурная схема предложенного способа приведена на фиг. 1.
Способ основан на следующих предположениях.
В общем случае для определения плотности распределения огибающей сигнала в канале с замираниями, когда доступными для измерения являются только значения огибающей смеси сигнал + шум, можно использовать подход, заключающийся в том, чтобы по плотности распределения огибающей смеси сигнал + шум определить параметры распределения Райса. При этом восстановить истинную плотность распределения огибающей можно, используя выборочную плотность распределения огибающей смеси сигнал + шум, получаемую посредством измерений на приемной стороне.
В данном подходе следует учитывать техническую проблему, связанную с тем, что на приемной стороне чаще всего сигнал перед обработкой проходит через устройство автоматической регулировки усиления (АРУ). Поскольку коэффициент усиления АРУ неизвестен и динамически меняется в процессе измерений, статистические характеристики выборочной плотности распределения амплитуды сигнала значительно меняются и применение указанных выше способов напрямую дает неадекватные оценки.
Избавиться от указанной трудности при приеме сигнала с использованием АРУ можно, если для оценки параметров модели канала использовать выборку случайных величин, инвариантную к значению коэффициента усиления АРУ. В качестве такой случайной величины может быть использована случайная величина ξ, определяемая как отношение огибающих Ai и Aj, измеренных на длительности одной и той же элементарной посылки на различных субчастотах с номерами i и j:
ξ=Ai/Aj.
Такой подход можно реализовать, если информационный сигнал является многочастотным, и при этом часть субчастот не используются для передачи. Тогда на входе приемника на занятых субчастотах наблюдается смесь информационного сигнала с шумом, а на свободных только шум.
Для описания плотности распределения огибающей шума на свободных субчастотах при гипотезе, что шум является гауссовским, используется плотность распределения Рэлея:
Figure 00000026
Тогда в качестве Ai можно использовать измеренную огибающую шума, а в качестве Aj - огибающую смеси сигнал + шум.
В случае постоянного уровня информационного сигнала А на соответствующих субчастотах для модели гауссовского шума функцию распределения случайной величины ξ можно найти следующим образом:
Figure 00000027
Если уровень информационного сигнала А не постоянен, а подвержен замираниям и его плотность распределения WA(x) подчиняется закону Накагами, то в этом случае функция распределения случайной величины ξ можно найти следующим образом:
Figure 00000028
При этом величина
Figure 00000029
представляет собой среднее значение отношения сигнал/помеха.
Тогда в новых обозначениях функция распределения случайной величины ξ имеет следующий вид:
Figure 00000030
Выражение для плотности при этом имеет следующий вид:
Figure 00000031
Сформировав выборку случайной величины ξ и имея аналитическое выражение для ее плотности распределения, можно воспользоваться методом максимального правдоподобия, как одним из методов оценки неизвестных параметров распределений. В данном случае неизвестными параметрами будут
Figure 00000032
и m. Тогда функция правдоподобия L, определяется выражением:
Figure 00000033
где xi - значение случайной величины ξ, n*М - объем выборки.
В этом случае координаты максимума функции правдоподобия
Figure 00000034
являются оценками искомых величин
Figure 00000035
и m.
Таким образом, приведенные аналитические выводы показывают, что с помощью предложенного способа можно определить параметры модели замирания огибающей сигнала по закону Накагами по информационному многочастотному сигналу. При этом необходимыми данными являются измеренные значения амплитуды смеси сигнала и шума на частотах, используемых для передачи информационного сигнала и значения амплитуды шума на частотах, неиспользуемых для передачи информационного сигнала.
Способ работает следующим образом.
На приемной стороне оцифровывают принимаемый сигнал в аналогово-цифровом преобразователе 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 и
Figure 00000036
, являющиеся координатами максимума функции правдоподобия
Figure 00000037
, где xi - это i-е значение выборки,
Figure 00000038
- плотность распределения вероятности измеряемой случайной величины.
Предлагаемый способ может быть использован для систем связи, использующих сигналы с ортогональным многочастотным разделением каналов связи. Применение такого способа позволяет более точно определять параметры замирающего канала связи.
Предлагаемое устройство по сравнению с прототипом обладает следующим преимуществом: обеспечивает более точное определение параметров модели замирания огибающей сигнала по закону Накагами по информационному многочастотному сигналу в случае наличия на приемной стороне блока автоматической регулировки усиления.

Claims (1)

  1. Способ оценки параметров модели замирания огибающей сигнала по закону Накагами по информационному многочастотному сигналу, заключающийся в том, что на приемной стороне оцифровывают принимаемый сигнал в аналогово-цифровом преобразователе, затем передают оцифрованный сигнал с выхода аналогово-цифрового преобразователя одновременно на вход первого блока вычисления амплитуды и на вход второго блока вычисления амплитуды, отличающийся тем, что в первом блоке вычисления амплитуды определяют значение амплитуды смеси принимаемого сигнала и шума на всех n частотах, используемых для передачи информационного сигнала, на длительности элементарной посылки, а с n выходов первого блока вычисления амплитуды на первые входы n соответствующих делителей передают вычисленные значения амплитуд, также во втором блоке вычисления амплитуды определяют значение амплитуды шума на n частотах, неиспользуемых для передачи информационного сигнала, на длительности элементарной посылки, а с n выходов второго блока вычисления амплитуды передают вычисленные значения амплитуд на вторые входы n соответствующих делителей, в каждом делителе осуществляют деление значения амплитуды шума на частоте, неиспользуемой для передачи информационного сигнала, полученное по второму входу на значение амплитуды смеси сигнала и шума на частоте, используемой для передачи информационного сигнала полученное по первому входу, а результат деления передают с выхода каждого из n делителей на n соответствующих входов блока накопления, в котором накапливают выборку полученных n значений на длительности интервала анализа, равной М посылкам, получая, таким образом, выборку размером n×М значений, а с выхода блока накопления передают накопленный массив значений на вход блока вычисления параметров распределения, в котором методом наискорейшего спуска, определяют параметры модели замирания огибающей сигнала по закону Накагами m и
    Figure 00000039
    , являющиеся координатами максимума функции правдоподобия
    Figure 00000040
    , где xi - это i-е значение выборки, а
    Figure 00000041
    - плотность распределения вероятности измеряемой случайной величины.
RU2015129277A 2015-07-16 2015-07-16 Способ оценки параметров модели замирания огибающей сигнала по закону накагами по информационному многочастотному сигналу RU2608363C1 (ru)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (2)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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