RU2444740C1 - Способ определения местоположения и мощностей источников излучения - Google Patents

Способ определения местоположения и мощностей источников излучения Download PDF

Info

Publication number
RU2444740C1
RU2444740C1 RU2011102806/07A RU2011102806A RU2444740C1 RU 2444740 C1 RU2444740 C1 RU 2444740C1 RU 2011102806/07 A RU2011102806/07 A RU 2011102806/07A RU 2011102806 A RU2011102806 A RU 2011102806A RU 2444740 C1 RU2444740 C1 RU 2444740C1
Authority
RU
Russia
Prior art keywords
power
vector
radiation
sources
axis
Prior art date
Application number
RU2011102806/07A
Other languages
English (en)
Inventor
Марина Витальевна Самойленко (RU)
Марина Витальевна Самойленко
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 RU2011102806/07A priority Critical patent/RU2444740C1/ru
Application granted granted Critical
Publication of RU2444740C1 publication Critical patent/RU2444740C1/ru

Links

Landscapes

  • Radar Systems Or Details Thereof (AREA)

Abstract

Способ определения местоположения и мощностей источников излучения позволяет по измеренной мощности на выходе приемной антенны однопозиционной локационной станции определять местоположение и мощности источников излучения. Достигаемый технический результат - расширение области применения способа определения местоположения источников на пассивные однопозиционные локационные системы, упрощение измерений и повышение информативности за счет определения мощностей источников. Способ заключается в разбиении контролируемой области пространства на элементы разрешения по местоположению, определении коэффициентов усиления, создаваемых приемной антенной для каждого элемента разрешения при выбранных априори направлениях оси диаграммы направленности (ДН) антенны, формировании матрицы усилений ! ! где wi(ej) - коэффициент усиления по мощности, создаваемый антенной при приеме излучения из i-го элемента разрешения при направлении оси ее ДН ej, измерении мощности на выходе антенны при выбранных априори направлениях оси ДН, оценивании вектора сигналов исходя из уравнения измерений p=Wf+n, где р - вектор измерений мощности, n - вектор ошибок измерений мощности на выходе приемной антенны, f - вектор мощностей источников, число компонент которого равно числу элементов разрешения, причем i-я компонента равна нулю, если в i-м элементе разрешения нет источника излучения, и равна мощности источника излучения, если он в i-м элементе разрешения есть, определении местоположения и мощностей источников излучения по оценке вектора мощностей источников. 1 ил.

Description

Изобретение относится к области антенных измерений и может быть использовано для высокоточного определения местоположения и мощностей источников излучения однопозиционной активной или пассивной локационной системой.
Под местоположением источников излучения понимается пространственное положение источников в заданной системе координат. Определение местоположения источников излучения в контролируемой области пространства является основной задачей радиолокации. Определение мощностей источников излучения, также входящее в задачу заявляемого способа, весьма полезно в ряде приложений, например, в задаче идентификации.
Под источниками излучения будем понимать как источники собственного излучения, так и источники переизлучения (отраженного объектом излучения). Под приемной антенной будем понимать как антенну, работающую только на прием сигналов, так и приемо-передающую антенну в режиме приема сигналов.
Известны способы определения местоположения источников излучения, основанные на обработке сигналов, принимаемых направленной антенной или антеннами. Местоположение источника определяется как точка пересечения линий его положения или линии и поверхности положения, определенных с помощью этих антенн. При этом используют ту или иную комбинацию решения двух независимых задач - определения направления на источник и определение его дальности.
Так, в пеленгационном способе для определения местоположения источника определяют из двух разнесенных в пространстве точек направления на этот источник и местоположение источника получают как точку пересечения в пространстве двух прямых, соответствующих найденным направлениям.
Недостатком этого способа является невозможность во многих случаях использовать двухпозиционную систему, включающую две пространственно разнесенные станции пеленгации, причем разнесенные на достаточно большое расстояние для обеспечения приемлемой точности. По этой причине пеленгационный способ определения местоположения используется, как правило, только в системах навигации.
Наиболее близким к заявляемому способу является дальномерно-пеленгационный способ определения местоположения источников излучения (прототип) [1]: он применим в однопозиционных локационных системах, т.е. позволяет определять местоположение источников из одной точки, с использованием одной антенны. Согласно [1], прототип является единственным известным способом, позволяющим однозначно определять местоположение источника из одной точки, т.е. однопозиционной локационной станцией. Таким же свойством обладает и заявляемый способ, однако основаны заявляемый способ и прототип на разных принципах.
Способ-прототип заключается в раздельном определении направления на источник излучения и его наклонной дальности. В результате решения первой задачи определяется линия положения - прямая в пространстве, указывающая направление на источник. В результате решения второй задачи определяется поверхность положения - сфера с центром в точке измерений, радиус которой равен наклонной дальности источника. Местоположение источника определяют как точку пересечения прямой и сферы.
При этом, если для определения направлений на источники достаточно измерений мощности на выходе приемной антенны, то для определения дальностей источников необходимы специальные измерительные схемы, которые, в зависимости от принципа действия, измеряют либо запаздывание фазы монохроматического излучения при использовании фазового метода дальнометрии [2], либо частоту биений при частотном методе [3], либо время запаздывания зондирующего импульса в импульсной дальнометрии [4].
Способ-прототип обладает следующими недостатками.
1) Он может быть реализован только в активной локационной системе, работающей по принципу излучения зондирующего сигнала и последующего приема этого сигнала, отраженного объектом наблюдения (источником отраженного излучения). Необходимость активной системы обусловлена тем, что все известные способы определения дальности основаны на сопоставлении отраженного сигнала с излученным. В то же время, в ряде случаев желательно применять пассивные локационные станции, работающие только на прием излучения источника.
2) Для реализации прототипа необходимы достаточно сложные измерения запаздывания фазы, частоты биений или запаздывания зондирующего импульса и соответствующая измерительная аппаратура.
3) Для однозначного определения дальности источника излучения по задержке зондирующего сигнала необходимо предварительное обнаружение источника с грубой оценкой его местоположения. Это требует соответствующих аппаратурных и временных затрат.
4) Прототип, как и другие методы определения местоположения источников излучения, не включает оценивание мощностей их излучения, которые представляют собой важную информацию для ряда прикладных задач.
Технической задачей данного изобретения является расширение области применения способа определения местоположения источников излучения на пассивные однопозиционные локационные системы, упрощение измерений, сокращение измерительной аппаратуры, уменьшение времени оценивания, а также получение информации о мощностях источников излучения.
Поставленная задача достигается тем, что контролируемую область пространства разбивают на малые объемы - элементы разрешения по местоположению, нумеруют их и фиксируют местоположение каждого элемента разрешения, например, с помощью координатных векторов r1, r2, …, rN, где N - число элементов разрешения в контролируемой области пространства, определяют коэффициенты усиления по мощности, создаваемые приемной антенной для каждого элемента разрешения при множестве априорно выбранных направлений оси ее диаграммы направленности (ДН) е1, е2, …, eK, где ej - направляющий вектор оси ДН, К - число выбранных направлений оси ДН, формируют из полученных коэффициентов усиления матрицу усилений
Figure 00000001
где wi(ej) - коэффициент усиления по мощности, создаваемый антенной при приеме излучения из i-го элемента разрешения при направлении оси ее ДН ej, устанавливают априорно выбранные направления оси ДН приемной антенны е1, е2, …, eK и измеряют при каждом направлении мощность на выходе приемной антенны р1, р2, …pK, находят оценку
Figure 00000002
вектора мощностей источников f исходя из уравнения измерений р=Wf+n, где р=[р1р2…рК]T - вектор измерений мощности на выходе приемной антенны, n - вектор ошибок измерений, индекс Т обозначает транспонирование, f - вектор, число компонент которого равно числу элементов разрешения N, причем i-я компонента равна нулю, если в i-м элементе разрешения нет источника излучения, и равна мощности источника излучения, если он в i-м элементе разрешения есть, определяют местоположение и мощности источников излучения по оценке вектора мощностей источников
Figure 00000002
, в котором значение каждой компоненты есть оценка мощности источника излучения, находящегося в элементе разрешения, номер которого равен номеру этой компоненты.
Особенностью заявляемого способа является возможность по одним лишь измерениям мощности на выходе приемной антенны пассивной или активной однопозиционной локационной станции определить местоположение источников излучения (или переизлучения) с точностью размера элемента разрешения, а также оценить мощности источников.
Обоснование способа.
Разобьем контролируемую область пространства на маленькие объемы - элементы разрешения по местоположению. Пронумеруем их и зафиксируем местоположение каждого элемента разрешения, например, с помощью координатных векторов их центров r1, r2, …, rN, где N - число элементов разрешения в контролируемой области пространства. Качественно такое разбиение представлено на фиг.1, где 1 - контролируемая область пространства, в которой показаны некоторые из элементов разрешения, на которые она разбита. Местоположения элементов разрешения задаются координатными векторами 2 с началом в точке базирования приемной антенны 3.
Определим коэффициенты усиления по мощности, создаваемые приемной антенной для каждого элемента разрешения. Это можно сделать экспериментально или по известной ДН. Определим такие коэффициенты усиления при множестве выбранных априори направлений оси ДН е1, е2, …, eK, где ei - направляющий вектор оси ДН, К - число выбранных направлений. При этом ДН должна быть достаточно широкоугольной, чтобы, по возможности, охватывать большую часть контролируемой области пространства.
Обозначим коэффициенты усиления, создаваемые антенной при j-м направлении оси ДН w1(ej), w2(ej), …wN(ej) и составим из полученных коэффициентов усиления вектор усилений
Figure 00000003
где индекс Т обозначает транспонирование.
Векторы усилений (1), полученные для всех К выбранных направлений оси ДН w(e1), w(e2), …w(eK), представляют собой априорную информацию, необходимую для реализации заявляемого способа. Эта информация определяется направленными свойствами антенны и выбранными направлениями оси ДН, которые будут затем устанавливаться при проведении измерений. Заметим, что вся априорная информация может быть получена заранее, до начала измерений, при практически неограниченном времени на ее получение.
Найдем выражение, определяющее выходную мощность приемной антенны при направлении оси ее ДН ej. Если в контролируемой области пространства находится только один источник излучения, расположенный в i-м элементе разрешения и мощность этого источника fi, то мощность на выходе приемной антенны определится выражением
p(ej)=fiwi(ej).
В общем случае источников в контролируемой области пространства может быть несколько. Для определения в общем случае выходной мощности антенны составим вектор мощностей источников
Figure 00000004
в котором число компонент равно числу элементов разрешения по местоположению в контролируемой области пространства. Значения компонент вектора (2) определим следующим образом: если в i-м элементе разрешения источника нет, то fi=0. Если же источник в i-м элементе разрешения есть, то компонента fi равна мощности этого источника.
С учетом значений компонент векторов (1) и (2), в общем случае произвольного числа источников, выходная мощность приемной антенны при направлении оси ее ДН ej определится выражением
Figure 00000005
Перейдем теперь к процессу измерений. Будем устанавливать ДН антенны в ранее выбранных направлениях е1, е2, …, eK и при каждом таком направлении оси ДН измерять мощность на выходе приемной антенны. Эти измерения, согласно (3), связаны с вектором мощностей источников f соотношениями
Figure 00000006
где введены обозначения wi=w(ei) и pi=p(ei).
В системе уравнений (4) известны измеренные мощности на выходе антенны p1, p2, …, pK, а также векторы усилений w1, w2, …,wK, которые представляют собой ранее определенную априорную информацию. Неизвестным является вектор мощностей источников f. Задача определения местоположения и мощностей источников излучения состоит в нахождении этого вектора. Действительно, каждая его компонента дает информацию о мощности источника, расположенного в элементе разрешения, номер которого равен номеру компоненты (если источника в этом элементе разрешения нет, то соответствующая компонента вектора f равна нулю).
Таким образом, задача свелась к нахождению вектора мощностей источников из системы уравнений (4). Для ее решения введем в рассмотрение вектор измерений
Figure 00000007
и матрицу усилений
Figure 00000008
С учетом (5) и (6) запишем систему уравнений (4) в виде векторно-матричного уравнения
Figure 00000009
С учетом ошибок измерений выходной мощности приемной антенны, уравнение (7) приобретает вид уравнения измерений
Figure 00000010
где n - вектор ошибок измерений мощности на выходе приемной антенны.
Уравнение (8) позволяет оценить вектор мощностей источников f. Это можно сделать разными известными методами. Так, метод винеровского оценивания [5] дает искомую оценку в виде выражения
Figure 00000011
где Rff и Rnn - ковариационные матрицы, соответственно, сигналов и ошибок измерений.
При отсутствии данных о вероятностных характеристиках сигналов и ошибок измерений, более грубую оценку можно найти из уравнения (7) методом псевдообращения [6]:
Figure 00000012
где индекс + обозначает псевдообратную матрицу.
Полученная оценка вектора мощностей источников
Figure 00000013
представляет собой оценку распределения источников излучения по элементам разрешения: каждая отличная от нуля компонента есть оценка мощности источника, расположенного в элементе разрешения с номером этой компоненты. Полученное распределение позволяет определить положение источников в контролируемой области пространства с точностью до размера элемента разрешения, который выбирается априори из соображений требуемых точности и разрешения.
Преимущества предлагаемого способа по сравнению с прототипом следующие.
1. Заявляемый способ обладает большей универсальностью, поскольку позволяет определять местоположение источников излучения не только с помощью активной (как прототип), но и с помощью пассивной однопозиционной локационной системы.
2. Существенно упрощается процедура измерений: уменьшается число измеряемых величин и, соответственно, сокращается необходимая для измерений аппаратура. Действительно, в прототипе необходимо измерять мощность на выходе приемной антенны для определения направлений на источники, а также с использованием специальных схем - задержку фазы, частоту биений или запаздывание импульса для определения дальностей источников. В заявляемом способе единственными измерениями являются измерения мощности на выходе приемной антенны.
3. В заявляемом способе отсутствует режим поиска с грубым предварительным оцениванием местоположения источников, которое необходимо в прототипе для однозначного определения дальности. Это сокращает время определения местоположения источников.
4. Поскольку в заявляемом способе используется широкоугольная ДН, появляется возможность применять антенны значительно меньших размеров, чем в прототипе, при одинаковой точности оценивания местоположения источников, что расширяет возможности применения способа. Это происходит потому, что в прототипе точность определяется шириной ДН, а в заявляемом способе - размером элемента разрешения.
5. Заявляемый способ является более информативным, чем прототип, поскольку, помимо местоположения источников излучения, определяет также их мощности.
Источники информации.
1. Сайбель А.Г. Основы радиолокации. - М., «Советское радио», 1961, с.15-17.
2. Бакулев П.А. Радиолокационные системы. - М., «Радиотехника», 2007, с.242-246.
3. Бакулев П.А. Радиолокационные системы. - М., «Радиотехника», 2007, с.246-252.
4. Бакулев П.А. Радиолокационные системы. - М., «Радиотехника», 2007, с.252-263.
5. Самойленко В.И., Пузырев В.А., Грубрин И.В. Техническая кибернетика: Учеб. пособие. - М.: Изд-во МАИ, 1994, c.130-132.
6. Гантмахер Ф.Р. Теория матриц. 4-е изд. - М.: Наука. Гл. ред. физ.-мат. лит., 1988, с.35.

Claims (1)

  1. Способ определения местоположения и мощностей источников излучения, отличающийся тем, что контролируемую область пространства разбивают на малые объемы - элементы разрешения по местоположению, нумеруют их и фиксируют местоположение каждого элемента разрешения, например, с помощью координатных векторов r1, r2, …, rN, где N - число элементов разрешения в контролируемой области пространства, определяют коэффициенты усиления по мощности, создаваемые приемной антенной для каждого элемента разрешения при множестве априорно выбранных направлений оси ее диаграммы направленности (ДН) e1, е2, …, еK, где ej - направляющий вектор оси ДН, K - число выбранных направлений оси ДН, формируют из полученных коэффициентов усиления матрицу усилений
    Figure 00000014

    где wi(ej) - коэффициент усиления по мощности, создаваемый антенной при приеме излучения из i-го элемента разрешения при направлении оси ее ДН ej, устанавливают априорно выбранные направления оси ДН приемной антенны e1, е2, …, еK и измеряют при каждом направлении мощность на выходе приемной антенны p1, р2, … рK, находят оценку
    Figure 00000015
    вектора мощностей источников f исходя из уравнения измерений p=Wf+n, где p=[p1p2…pK]T - вектор измерений мощности на выходе приемной антенны, n - вектор ошибок измерений, индекс Т обозначает транспонирование, f - вектор, число компонент которого равно числу элементов разрешения N, причем i-я компонента равна нулю, если в i-м элементе разрешения нет источника излучения, и равна мощности источника излучения, если он в i-м элементе разрешения есть, определяют местоположение и мощности источников излучения по оценке вектора мощностей источников
    Figure 00000016
    в котором значение каждой компоненты есть оценка мощности источника излучения, находящегося в элементе разрешения, номер которого равен номеру этой компоненты.
RU2011102806/07A 2011-01-26 2011-01-26 Способ определения местоположения и мощностей источников излучения RU2444740C1 (ru)

Priority Applications (1)

Application Number Priority Date Filing Date Title
RU2011102806/07A RU2444740C1 (ru) 2011-01-26 2011-01-26 Способ определения местоположения и мощностей источников излучения

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
RU2011102806/07A RU2444740C1 (ru) 2011-01-26 2011-01-26 Способ определения местоположения и мощностей источников излучения

Publications (1)

Publication Number Publication Date
RU2444740C1 true RU2444740C1 (ru) 2012-03-10

Family

ID=46029146

Family Applications (1)

Application Number Title Priority Date Filing Date
RU2011102806/07A RU2444740C1 (ru) 2011-01-26 2011-01-26 Способ определения местоположения и мощностей источников излучения

Country Status (1)

Country Link
RU (1) RU2444740C1 (ru)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2499273C1 (ru) * 2012-06-18 2013-11-20 Марина Витальевна Самойленко Способ определения местоположений и мощностей источников излучения однопозиционной локационной станцией
CN113534077A (zh) * 2021-07-16 2021-10-22 北京华力创通科技股份有限公司 一种雷达辐射源威力反演方法、装置及电子设备
RU2796396C1 (ru) * 2022-12-08 2023-05-23 Александр Сергеевич Дмитриев Устройство долговременного непрерывного измерения плотности мощности неионизирующего электромагнитного излучения и плотности суммарной энергии
WO2024123213A1 (ru) * 2022-12-08 2024-06-13 Александр Сергеевич Дмитриев Устройство для измерения плотности мощности и суммарной энергии неионизирующего электромагнитного излучения

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0477798A3 (en) * 1990-09-27 1993-04-07 Alliant Techsystems Inc. Apparatus and method for determining the position and attitude of an object
RU2251803C1 (ru) * 2004-07-20 2005-05-10 Панов Владимир Петрович Способ определения информативных параметров и характеристик радиосигналов передатчиков
JP3802422B2 (ja) * 2000-04-24 2006-07-26 ロッキード・マーティン・コーポレイション 受動可干渉位置決めシステムおよびそのシステムの使用方法
CA2739723A1 (en) * 2008-10-06 2010-04-15 Tektrap Systems, Inc. Method and device for tracing objects and detecting change in configuration of objects
RU2392634C1 (ru) * 2009-05-18 2010-06-20 Марина Витальевна Самойленко Способ определения направлений на источники излучения и углового разрешения источников
RU2405166C2 (ru) * 2009-02-17 2010-11-27 Закрытое акционерное общество "ИРКОС" Способ определения местоположения передатчика переносным пеленгатором

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0477798A3 (en) * 1990-09-27 1993-04-07 Alliant Techsystems Inc. Apparatus and method for determining the position and attitude of an object
JP3802422B2 (ja) * 2000-04-24 2006-07-26 ロッキード・マーティン・コーポレイション 受動可干渉位置決めシステムおよびそのシステムの使用方法
RU2251803C1 (ru) * 2004-07-20 2005-05-10 Панов Владимир Петрович Способ определения информативных параметров и характеристик радиосигналов передатчиков
CA2739723A1 (en) * 2008-10-06 2010-04-15 Tektrap Systems, Inc. Method and device for tracing objects and detecting change in configuration of objects
RU2405166C2 (ru) * 2009-02-17 2010-11-27 Закрытое акционерное общество "ИРКОС" Способ определения местоположения передатчика переносным пеленгатором
RU2392634C1 (ru) * 2009-05-18 2010-06-20 Марина Витальевна Самойленко Способ определения направлений на источники излучения и углового разрешения источников

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
САЙБЕЛЬ Ф.Г. Основы радиолокации. - М.: Советское радио, 1961. С.15-17. *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2499273C1 (ru) * 2012-06-18 2013-11-20 Марина Витальевна Самойленко Способ определения местоположений и мощностей источников излучения однопозиционной локационной станцией
CN113534077A (zh) * 2021-07-16 2021-10-22 北京华力创通科技股份有限公司 一种雷达辐射源威力反演方法、装置及电子设备
CN113534077B (zh) * 2021-07-16 2024-03-12 北京华力创通科技股份有限公司 一种雷达辐射源威力反演方法、装置及电子设备
RU2796396C1 (ru) * 2022-12-08 2023-05-23 Александр Сергеевич Дмитриев Устройство долговременного непрерывного измерения плотности мощности неионизирующего электромагнитного излучения и плотности суммарной энергии
WO2024123213A1 (ru) * 2022-12-08 2024-06-13 Александр Сергеевич Дмитриев Устройство для измерения плотности мощности и суммарной энергии неионизирующего электромагнитного излучения
RU225252U1 (ru) * 2023-12-19 2024-04-16 Публичное акционерное общество "Россети Центр" (ПАО "Россети Центр") Устройство для фиксации электромагнитного излучения

Similar Documents

Publication Publication Date Title
US7701393B2 (en) Radio frequency navigation using frequency response matching
JP6778336B2 (ja) Rf信号の到来角度決定方法及びシステム
TWI377362B (en) System and technique for calibrating radar arrays
US9672624B2 (en) Method for calibrating absolute misalignment between linear array image sensor and attitude control sensor
Norouzi et al. Joint time difference of arrival/angle of arrival position finding in passive radar
RU2444740C1 (ru) Способ определения местоположения и мощностей источников излучения
RU2593149C1 (ru) Адаптивный способ пассивной радиолокации
Sun et al. Array geometry calibration for underwater compact arrays
Bordonaro et al. Consistent linear tracker with position and range rate measurements
US10348423B1 (en) Calculating a virtual boresight vector of an antenna array
RU2499273C1 (ru) Способ определения местоположений и мощностей источников излучения однопозиционной локационной станцией
US11774574B2 (en) Methods and system for determining an angle of a detection
Yong et al. An asynchronous sensor bias estimation algorithm utilizing targets’ positions only
KR101957291B1 (ko) 전자전 지원 시스템(Warfare Support System)의 신호도래방위각 측정 장치 및 그 방법
RU2392634C1 (ru) Способ определения направлений на источники излучения и углового разрешения источников
CN113625063A (zh) 一种整机条件下的天线单脉冲性能评估方法
Le Scouarnec et al. A positioning free calibration method for mobile laser scanning applications
RU2776078C1 (ru) Способ определения координат и параметров движения источников радиоизлучений по измерениям их несущей частоты
RU2406098C1 (ru) Способ определения наклонной дальности до движущейся цели по минимальному числу пеленгов
RU2390038C2 (ru) Способ определения скорости и направления движения носителя обзорной рлс
KR20200062787A (ko) 배열 안테나 기반 위치 추정 방법 및 그 장치
RU2231806C2 (ru) Способ оценки текущих координат источника радиоизлучения
Bulychev et al. Single-Position Location of a Radiation Source with a Curvilinear Motion and Allowance for Signal Period Evolution at the Reception Point
Elsaesser Sensor data fusion using a probability density grid
RU2603971C1 (ru) Способ измерения углов в фазовых многошкальных угломерных системах и устройство, его реализующее

Legal Events

Date Code Title Description
MM4A The patent is invalid due to non-payment of fees

Effective date: 20180127