СПОСОБ АДАПТИВНОГО ВЫБОРА ОПТИМАЛЬНОГО ПАРАМЕТРА АЛГОРИТМА КОРРЕКЦИИ НА ОСНОВЕ РАЗБРОСА ФАЗ КОРРЕКТИРУЕМОГО СИГНАЛА Российский патент 2019 года по МПК H03H21/00 

Изобретение относится к технике связи и может быть использовано в системах передачи данных с адаптивной коррекцией сигналов для выбора параметра алгоритма коррекции.

Во многих системах передачи данных для компенсации искажений, внесенных каналом связи, применяют алгоритмы адаптивной коррекции сигналов. Для этого в передаваемый сигнал осуществляют периодические вставки известного на приемной стороне тестового сигнала. Такой подход используется, например, в стандарте авиационной передачи данных ARINC-635 [1].

Известно большое количество различных методов, алгоритмов и их модификаций, используемых для коррекции сигналов, например, метод наименьших квадратов или алгоритм LMS [2], алгоритм RLS [3], метод регуляризации Тихонова [4]. Во всех этих и многих других алгоритмах для обеспечения устойчивости и сходимости вводят некоторый параметр, в частности, параметр регуляризации (алгоритм RLS, метод регуляризации Тихонова), размер шага сходимости (алгоритм LMS). Выбор того или иного параметра оказывает существенное влияние не только на устойчивость решения, но и на вероятность ошибки на бит после демодуляции откорректированного сигнала, т.е. на помехоустойчивость.

Как известно задача адаптивной коррекции сводится к решению двух уравнений, которые можно записать в следующей форме:

где K(t) - передаваемый тестовый сигнал, u(t) - принимаемый тестовый сигнал, K_m(t) - передаваемый информационный сигнал, u_m(t) - принимаемый информационный сигнал, h(t) - импульсная характеристика канала, * - оператор свертки.

Из уравнения (1) получают приближенную импульсную характеристику канала в общем случае в виде:

а результат коррекции в этом случае можно записать в виде:

где R{} - некоторый алгоритм расчета, α₁, α₂ - параметры, используемые для устойчивости алгоритма. Отметим, что в большинстве практических случаев допустимо принять:

тогда вместо (5) запишем:

Известны различные способы выбора оптимального значения этого параметра.

Известен способ невязки описанный, в [5]. Этот способ часто применяют для выбора параметра регуляризации в методе регуляризации Тихонова. В условиях (3) (6) способ невязки заключается в том, что, используя алгоритм R{} на основе поступившего тестового сигнала u₀(t) рассчитывают импульсную характеристику канала h(t, α) и корректирующего фильтра h_кop(t, α), с помощью которой, используя некоторый алгоритм R{} корректируют поступивший информационный сигнал u_m(t), в результате чего получают откорректированный информационный сигнал K_m(t, α), после чего определяют значение ошибки е, в качестве которой служит разница среднеквадратичного уклонения откорректированного информационного сигнала K_m(t, α), свернутого с рассчитанной импульсной характеристикой канала h(t, α), от принятого информационного сигнала u_m(t) и дисперсии шумовой составляющей Δu, т.е. после чего определяют зависимость значения ошибки е от параметра α, путем изменения значения этого параметра, в результате чего получают оптимальное значение параметра α_opt, обеспечивающего минимальное значение ошибки е.

Недостатком указанного способа является необходимость знания определенной априорной информации, а именно дисперсии шумовой составляющей, оценка которой является отдельной достаточно сложной задачей и имеет определенную погрешность. Кроме того при вычислении значения ошибки е вносится дополнительная погрешность при свертке приближенных (рассчитанных) значений K_m(t, α) и h(t, α).

Известен способ адаптивного выбора оптимального параметра алгоритма коррекции сигналов описанный в патенте РФ №2653485 от 08.05.2018 и принятый за прототип. Способ адаптивного выбора оптимального параметра алгоритма коррекции сигналов заключается в том, что, используя алгоритм R{} на основе поступившего тестового сигнала u₀(t) рассчитывают импульсную характеристику корректирующего фильтра h_коp(t, α), с помощью импульсной характеристики корректирующего фильтра h_кop(t, α), используя алгоритм R{} корректируют поступившие тестовые сигналы u₁(t)…u_n(t), задержанные на интервал, равный длине информационного сигнала L_И, в результате чего получают откорректированные тестовые сигналы K₁(t, α)…K_n(t, α), определяют значения ошибки е₁…e_n, в качестве которой служит среднеквадратичное уклонение откорректированного тестового сигнала K₁(t, α)…K_n(t, α) от образцового тестового сигнала после чего определяют зависимость значения ошибки е₁…e_n от параметра а, путем изменения значения этого параметра, в результате чего получают массив значений параметров α₁…α_n, обеспечивающие соответствующее минимальное значение ошибки е₁…e_n для каждого откорректированного тестового сигнала K₁(t, α)…K_n(t, α), после чего из массива α₁…a_n осуществляют окончательный выбор оптимального значения параметра α_opt, в качестве которого, в зависимости от конкретного алгоритма R{} и диапазона значений параметров α₁…α_n, берут среднее арифметическое значение или медианное значение из массива α₁…α_n.

Недостатком прототипа является то, что в каналах с быстрыми замираниями импульсная характеристика канала для тестовых сигналов пришедших много «раньше» или много «позже» обрабатываемого в данный момент тесового сигнала может значительно отличаться от текущей импульсной характеристики канала, а, следовательно, и полученная текущая импульсная характеристика корректирующего фильтра не будет являться оптимальной для более «ранних» и «поздних» тестовых сигналов. Это приведет к сильной погрешности в оценках значений ошибки е₁…e_n, соответствующих более «ранним» и «поздним» тестовым сигналам, и соответствующей погрешности в определении оптимального значения алгоритма коррекции. Чтобы избежать данного эффекта, можно уменьшить число анализируемых тестовых сигналов при вычислении значений ошибок е₁…e_n, однако при малых отношениях сигнал/шум объема поучаемой статистики е₁…e_n может быть недостаточно для получения достаточно точной оценки оптимального значения алгоритма коррекции.

Целью изобретения является выбор оптимального параметра алгоритма коррекции сигналов на основе анализа разброса фаз корректируемого информационного сигнала без знания априорной информации о мощности шума, отношении сигнал/шум и корректируемого информационного сигнала.

Поставленная цель достигается тем, что способ адаптивного выбора оптимального параметра алгоритма коррекции на основе разброса фаз корректируемого сигнала заключается в том, что, используя алгоритм R{} на основе поступившего тестового сигнала u₀(t) рассчитывают импульсную характеристику корректирующего фильтра h_кop(t, α), при этом с помощью импульсной характеристики корректирующего фильтра h_кop(t, α), используя алгоритм R{} корректируют информационный сигнал u_m{t) с фазовой манипуляцией, переносящий N информационных символов, в результате чего получают откорректированный информационный сигнал K_m(t, α), определяют абсолютное значение фазы и амплитуды каждого символа откорректированного сигнала, получая значения фаз Ф₁…Ф_N и значения амплитуд A₁…A_N, преобразовывают полученные значения фаз, таким образом, чтобы все преобразованные фазы всех символов располагались в секторе соответствующего символу «0», получая значения преобразованных фаз после чего вычисляют экспоненты A₁⋅exp(i⋅Ф₁)…A_N⋅exp(i⋅Ф_N), суммируют полученные значения экспонент, вычисляют модуль полученной суммы |Q|, после чего определяют зависимость |Q| от параметра α, путем изменения значения этого параметра, в результате чего получают оптимальное значение параметра а_опт, соответствующее максимуму |Q|.

На фиг. 1 представлена структурная схема способа адаптивного выбора оптимального параметра алгоритма коррекции на основе разброса фаз корректируемого сигнала. Она содержит:

1 - первая линия задержки;

2 - блок расчета импульсной характеристики;

3 - корректирующий фильтр;

4 - вторая линия задержки;

5(1), 5(2)…5(N) - фазовый детектор;

6 - преобразователь фазы;

7(1), 7(2)…7(N) - блок вычисления экспоненты;

8 - сумматор;

9 - решающее устройство.

Работа способа осуществляется следующим образом. На вход первой линии задержки 1 поступает сигнал, содержащий периодически повторяющиеся тестовые и информационные сигналы. Структура такого сигнала представлена на фиг. 2. С первого выхода первой линии задержки 1 на первый вход блока расчета импульсной характеристики 2 подают поступивший тестовый сигнал u₀(t) и рассчитывают импульсную характеристику корректирующего фильтра h_коp(t, α). Полученную импульсную характеристику корректирующего фильтра h_коp(t, α) с выхода блока расчета импульсной характеристики 2 подают на второй вход корректирующего фильтра 3. Со второго выхода первой линии задержки на первый вход корректирующего фильтра 3 подают поступивший информационный сигнал u_m(t) с фазовой манипуляцией, переносящий N информационных символов. В результате на выходе корректирующего фильтра 3 получают откорректированный информационный сигнал K_m(t, α), который подают на вход второй линии задержки 4. При этом в блоке расчета импульсной характеристики 2 и корректирующем фильтре 3 используют один и тот же алгоритм обозначенный ранее как R{}.

С 1, 2…N выходов второй линии задержки 4 на вход соответствующего фазового детектора 5(1), 5(2)…5(N) подают сигналы соответствующие информационным символам 1, 2…N. В фазовом детектора 5(1), 5(2)…5(N) получают значения фаз Ф₁…Ф_N и значения амплитуд A₁…A_N, соответствующего символа. При этом каждое значение фазы представляет собой смесь следующих компонент:

где k - номер символа, ϕ_k - значение фазы символа (0, π - для двухпозиционной фазовой манипуляции, - для четырехпозиционной фазовой манипуляции и т.д.), ψ - постоянное смещение фазы, ξ_k - погрешность фазы, связанная с шумом и неоптимальным выбором параметра α.

Полученные значения фаз Ф₁…Ф_N и значения амплитуд A₁…A_N с выхода каждого фазового детектора 5(1), 5(2)…5(N) подают на 1, 2…N входы преобразователя фазы 6, где преобразовывают полученные значения фаз Ф₁…Ф_N, таким образом, чтобы все преобразованные фазы всех символов располагались в секторе соответствующего символу «0», в результате получая значения преобразованных фаз . Например для двухпозиционной фазовой манипуляции вычисляют:

в результате чего все значения преобразованных фаз , т.е. располагаются в секторе соответствующего символу «0».

С 1, 2…N выходов преобразователя фазы 6 на вход соответствующего блока вычисления экспоненты 7(1), 7(2)…7(N) подают значения преобразованных фаз и значения амплитуд A₁…A_N. В каждом блоке вычисления экспоненты 7(1), 7(2)…7(N) вычисляют экспоненты и передают их на 1, 2…N входы сумматора 8, где суммируют полученные значения экспонент и вычисляют модуль суммы |Q|.

Полученное на выходе сумматора 8 значение |Q| подают на вход решающего устройства 9, в котором определяют зависимость |Q| от параметра а, путем изменения значения этого параметра на первом выходе решающего устройства 9, и передаче его на второй вход блока расчета импульсной характеристики 2 и третий вход корректирующего фильтра 3. В результате на втором выходе решающего устройства 9 получают оптимальное значение параметра α_onm, соответствующее максимуму |Q|.

Техническим результатом является выбор оптимального параметра алгоритма коррекции сигналов на основе анализа разброса фаз корректируемого информационного сигнала, не требующий знания априорной информации о мощности шума, отношении сигнал/шум и корректируемого информационного сигнала.

Список источников

1. ARINC Characteristic 635-4. HF Data Link Protocol. - Dec, 2003.

2. Джиган В.И. Адаптивная фильтрация сигналов: теория и алгоритмы. М.: Техносфера, 2013. - 528 с.

3. Sayed А.Н. Adaptive filters. - New Jersey: Hoboken: John Wiley & Sons, Inc., 2008. - 786 c.

4. Тихонов A.H., Арсенин В.Я. Методы решения некорректных задач / Учебное пособие для вузов. - Изд. 3-е испр. - М.: Наука, 1986. - 288 с.

5. Верлань А.Ф., Сизиков B.C. Методы решеня интегральных уравнений с программами для ЭВМ. Справочное пособие. - Киев: Наукова думка, 1978. - 292 с.

Изобретение относится к технике связи и может быть использовано в системах передачи данных с адаптивной коррекцией сигналов для выбора параметра алгоритма коррекции. Техническим результатом является выбор оптимального параметра алгоритма коррекции сигналов на основе анализа разброса фаз корректируемого информационного сигнала, без знания априорной информации о мощности шума, отношении сигнал/шум и корректируемого информационного сигнала. Используя некоторый алгоритм R{}, на основе поступившего тестового сигнала u₀(t) рассчитывают импульсную характеристику корректирующего фильтра h_кор(t, α), с помощью которой, используя алгоритм R{}, корректируют поступивший информационный сигнал u_m(t) с фазовой манипуляцией, переносящий N информационных символов, в результате чего получают откорректированный информационный сигнал K_m(t, α). Определяют абсолютное значение фазы и амплитуды каждого символа откорректированного сигнала, получая значения фаз Ф₁…Ф_N и значения амплитуд A₁…A_N. Преобразовывают полученные значения фаз таким образом, чтобы все преобразованные фазы всех символов располагались в секторе, соответствующем символу «0», получая значения преобразованных фаз после чего вычисляют экспоненты A₁·exp(i·Ф₁)… A₁·exp(i·Ф_N). Суммируют полученные значения экспонент. Вычисляют модуль полученной суммы |Q|, после чего определяют зависимость |Q| от параметра α путем изменения значения этого параметра, в результате чего получают оптимальное значение параметра α_опт, соответствующее максимуму |Q|. 2 ил.

Способ адаптивного выбора оптимального параметра алгоритма коррекции на основе разброса фаз корректируемого сигнала, заключающийся в том, что, используя алгоритм R{ }, на основе поступившего тестового сигнала u₀(t) рассчитывают импульсную характеристику корректирующего фильтра h_кор(t,α), отличающийся тем, что с помощью импульсной характеристики корректирующего фильтра h_кор(t,α), используя алгоритм R{}, корректируют информационный сигнал u_m(t) с фазовой манипуляцией, переносящий N информационных символов, в результате чего получают откорректированный информационный сигнал K_m(t,α), определяют абсолютное значение фазы и амплитуды каждого символа откорректированного сигнала, получая значения фаз Ф₁…Ф_N и значения амплитуд A₁…A_N, преобразовывают полученные значения фаз таким образом, чтобы все преобразованные фазы всех символов располагались в секторе, соответствующем символу «0», получая значения преобразованных фаз, после чего вычисляют экспоненты , суммируют полученные значения экспонент, вычисляют модуль полученной суммы |Q|, после чего определяют зависимость |Q| от параметра α путем изменения значения этого параметра, в результате чего получают оптимальное значение параметра α_опт, соответствующее максимуму |Q|.