Способ обнаружения подводного источника широкополосного шума

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

Известен способ обнаружения шумящих в море объектов (Патент РФ №2300118, МПК G01S 3/80, G01S 15/04, опубликован 27.05.2007 г. бюл. №15), включающий прием первичного поля шумоизлучения объектов многоканальной антенной, при котором осуществляют частотно-временную обработку принятых шумовых сигналов для каждого канала наблюдения в горизонтальной плоскости, квадрируют, усредняют по времени, центрируют и нормируют сигналы к помехе, осуществляют наблюдение на каждом цикле обзора отметок принятых шумовых сигналов всей совокупности пространственных каналов в горизонтальной плоскости и принимают решение об обнаружении путем сравнения с пороговым значением отношения сигнал-помеха.

В известном способе введены новые операции, а именно:

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

во второй последовательности осуществляют сопровождение информационных параметров шумового сигнала по пеленгу, дисперсии пеленга, по скорости и ускорению изменения пеленга и отбраковку пеленгов локальных ложных максимумов шумовых сигналов, для чего выполняют следующие операции: фиксируют пеленг пространственного канала, в котором наблюдается каждый локальный максимум сигнала; вычисляют по заданному закону аппроксимации уточненное значение пеленга шумового сигнала по данным нескольких откликов пространственных каналов в окрестности данного локального максимума сигнала, образующих сигнальный фрагмент; составляют из совокупности оценки пеленга и величины изменения пеленга (ВИП) вектор параметров движения локальных максимумов сигнала; вычисляют матрицу взаимно корреляционных функций для вектора параметров движения локальных максимумов сигнала; вычисляют скорость изменения ВИПа и пеленга за время между предыдущим и текущим циклами обзора; в результате определяют прогнозные оценки пеленга и ВИПа локальных максимумов сигнала за время между предыдущим и текущим циклами обзора для заданного времени накопления; определяют дисперсию прогнозной оценки пеленга для заданного времени накопления и вычисляют ширину строба по пеленгу, в пределах которого осуществляется наблюдение каждого сигнала; вычисляют плотность вероятности смещения измеренного пеленга для заданного времени накопления; по результатам выполнения двух последовательностей операций вычисляют обобщенный вес локальных максимумов сигнала; сравнивают обобщенный вес локальных максимумов сигнала с порогом обнаружения сигнала, который соответствует пороговому отношению сигнал-помеха.

Недостатком известного способа является малая помехоустойчивость и малая дальность действия приемной системы при работе на низких частотах, когда размер приемной системы соизмерим с длиной волны, и при работе в мелком море, когда алгоритмы формирования пространственной направленности становятся неэффективными из-за дисперсионных искажений сигналов, а также возможность реализации способа только для сигналов фиксированного частотного диапазона, в котором апертура антенны обеспечивает необходимую пространственную направленность.

Известен также способ обнаружения шумящих объектов в мелком и глубоком море (патент РФ №2653189, МПК G01S 3/80, G01S 15/04, опубликован 07.05.2018 г. бюл. №13), включающий прием шумового сигнала комбинированным приемником, содержащим приемник звукового давления и трехкомпонентный приемник вектора колебательной скорости, при котором выполняют частотно-временной обработку в заданном фиксированном частотном диапазоне принятого шумового сигнала; вычисляют в каждом частотном канале, сформированном в результате частотно-временной обработки принятого шумового сигнала, текущие значения комплексных амплитуд звукового давления, трех компонент вектора колебательной скорости, трех компонент вещественной составляющей вектора интенсивности и трех компонент мнимой составляющей вектора интенсивности в локальной системе координат комбинированного приемника для суммарного процесса сигнал плюс помеха и для помехи отдельно; формируют в каждом частотном канале усредненные за время T₁ значения трех компонент вещественной составляющей вектора интенсивности, трех компонент мнимой составляющей вектора интенсивности и квадрата звукового давления для суммарного процесса сигнал плюс помеха и для помехи отдельно; формируют в каждом частотном канале усредненные за время T₂=10T₁ комплексные амплитуды нулевой и первой гармоник вторичного спектра для трех компонент вещественной составляющей вектора интенсивности, трех компонент мнимой составляющей вектора интенсивности и квадрата звукового давления для суммарного процесса сигнал плюс помеха и для помехи отдельно; нормируют все 21 информативный параметр, вычисленных для суммарного процесса сигнал плюс помеха, на соответствующие значения информативных параметров, вычисленные для помехи; вычисляют максимальное отношение сигнал/помеха для одного из 21 информативных параметров и принимают решения об обнаружении путем его сравнения с пороговым значением отношения сигнал/помеха.

Отличительными существенными признаками этого известного способа являются следующие операции: используют в качестве приемной системы комбинированный приемник, содержащий, кроме приемника звукового давления, трехкомпонентный приемник вектора колебательной скорости; формируют методами частотно-временной обработки сигналов набор частотных каналов в заданном фиксированном частотном диапазоне в векторных каналах комбинированного приемника; - вычисляют в каждом частотном канале, сформированном в результате частотно-временной обработки принятых шумовых сигналов, текущие значения комплексных амплитуд трех компонент вектора колебательной скорости, текущие значения амплитуд трех компонент вещественной составляющей вектора интенсивности, трех компонент мнимой составляющей вектора интенсивности для суммарного процесса сигнал плюс помеха;- усредняют за заранее определенный временной интервал T₁ значения трех компонент вещественной составляющей вектора интенсивности, трех компонент мнимой составляющей вектора интенсивности и квадрата звукового давления для суммарного процесса сигнал плюс помеха; выделяют из текущих значений суммарного случайного процесса сигнал плюс помеха текущие значения помехи; вычисляют в каждом частотном канале текущие значения комплексных амплитуд трех компонент вектора колебательной скорости, текущие значения амплитуд трех компонент вещественной составляющей вектора интенсивности, трех компонент мнимой составляющей вектора интенсивности для помехи; усредняют за заранее определенный временной интервал T₁ значения трех компонент вещественной составляющей вектора интенсивности, трех компонент мнимой составляющей вектора интенсивности и квадрата звукового давления для помехи; вычисляют в каждом частотном канале за заранее определенный временной интервал T₂=10T₁ текущие значения комплексных амплитуд нулевой и первой гармоник вторичного спектра для трех компонент вещественной составляющей вектора интенсивности, трех компонент мнимой составляющей вектора интенсивности и квадрата звукового давления для суммарного процесса сигнал плюс помеха; вычисляют в каждом частотном канале за заранее определенный временной интервал Т₂=10T₁ текущие значения комплексных амплитуд нулевой и первой гармоник вторичного спектра для трех компонент вещественной составляющей вектора интенсивности, трех компонент мнимой составляющей вектора интенсивности и квадрата звукового давления для помехи; нормируют квадрат звукового давления и компоненты комплексного вектора интенсивности, усредненные за время Т₁ вычисленные для суммарного процесса сигнал плюс помеха, на соответствующие значения квадрата звукового давления и компоненты комплексного вектора интенсивности, усредненные за время Т₁ вычисленные для помехи; нормируют вычисленные за время Т₂=10Т₁ текущие значения комплексных амплитуд нулевой и первой гармоник вторичного спектра для трех компонент вещественной составляющей вектора интенсивности, трех компонент мнимой составляющей вектора интенсивности и квадрата звукового давления для суммарного процесса сигнал плюс помеха на соответствующие текущие значения комплексных амплитуд нулевой и первой гармоник вторичного спектра для трех компонент вещественной составляющей вектора интенсивности, трех компонент мнимой составляющей вектора интенсивности и квадрата звукового давления для помехи; вычисляют максимальное отношение сигнал/помеха из набора 21 информативных параметров, 7 информативных параметров для усредненных за время T₁ нормированных на помеху значений комплексного вектора интенсивности и квадрата звукового давления и 14 информативных параметров для усредненных за время T₂=10T₁ нормированных на помеху значений комплексных амплитуд нулевой и первой гармоник вторичного спектра для комплексного вектора интенсивности и квадрата звукового давления; принимают в качестве модельной статистики поля помехи в гидрофонном канале и в каналах вектора колебательной скорости гауссову статистику; принимают в качестве модельной статистики поля помехи в каналах вектора интенсивности Лапласову статистику; вычисляют на основе принятых статистик аналитическую зависимость вероятности правильного обнаружения при заданной вероятности ложной тревоги от порогового отношения сигнал-помеха по методу максимального правдоподобия; принимают решение об обнаружении путем сравнения с пороговым значением отношения сигнал/помеха максимального отношения сигнал/помеха, вычисленного из набора 21 информативных параметров.

Известный способ является наиболее близким к заявленному изобретению и принят за прототип.

Для данного известного способа обнаружения источника шумового сигнала характерна низкая помехоустойчивость и высокая чувствительность к ложным срабатываниям из-за отсутствия достоверных данных о характеристиках помехи, на основе которых формируются пороговые уровни в устройстве обнаружения, а также применения пороговых энергетических механизмов обнаружения, основанных на вычислении величины не всей энергии сигнала, принимаемого от источника, а лишь отдельных наиболее сильных ее спектральных составляющих.

Кроме того, для данного способа характерна высокая вероятность ложных срабатываний, которая формируется из-за малого объема объективных данных об источнике, к которым относятся данные о спектральных характеристиках источника, его угловом положении, текущих координатных оценках и их изменении во времени.

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

Для решения поставленной задачи в способе обнаружения подводного источника широкополосного шума, включающем прием шумового сигнала комбинированным приемником, содержащим приемник звукового давления и трехкомпонентный приемник вектора колебательной скорости, при котором выполняют частотно-временную обработку принятого сигнала в наблюдаемом диапазоне частот, вычисляют в каждом частотном канале, сформированном в результате частотно-временной обработки принятого шумового сигнала, текущие значения комплексных амплитуд звукового давления и трех компонент вектора колебательной скорости, а также для текущего временного фрагмента принятого сигнала заданной длительности вычисляют в каждом частотном канале средние величины компонент вектора плотности потока энергии и принимают решение об обнаружении источника, прием шумового сигнала выполняют двумя комбинированными приемниками, выставленными в акватории в двух точках с известными координатами, приемники снабжены средствами контроля углового положения их осей, а также устройствами связи, по которым от каждого приемника результаты синхронизированного вычисления для текущего временного фрагмента принятого сигнала средних величин горизонтальных компонент вектора плотности потока энергии в каждом сформированном частотном канале и текущие данные положения осей приемников относительно севера передают в пункт обработки данных, в котором для первого и второго приемников в каждом сформированном частотном канале вычисляют величину горизонтального потока энергии (интенсивность потока) Е(ϖ) и углы β(ϖ) прихода потока по формулам:

где I_x, I_y - средние величины горизонтальных компонент вектора плотности потока энергии;

- знак комплексного сопряжения;

, ϖ_н, ϖ_в - нижняя и верхняя частоты наблюдаемого диапазона;

для каждого приемника суммируют энергии потоков частотных каналов, углы прихода потоков энергии которых группируются в выделенном угловом секторе шириной Δβ, при этом формируют равные горизонтальные угловые секторы шириной Δβ, охватывающие весь горизонт наблюдения 0-360 град, и строят угловое распределение величины потока энергии по всем сформированным угловым секторам Е₁(Δβ), Е₂(Δβ); в каждом приемнике определяют среднюю секторную величину потока энергии по всем сформированным угловым секторам наблюдения, назначают ее пороговой величиной и выделяют из сформированных угловых секторов доминирующие секторы, в которых величина потока энергии превышает пороговую; для каждого приемника устанавливают информативные частотные каналы, которые формируют энергию в доминирующих угловых секторах, фиксируют спектральные энергетические портреты принятых в доминирующих угловых секторах сигналов по величине потоков энергии информативных частотных каналов, формируют оценки направлений на источник от каждого приемника β₁(ϖ), β₂(ϖ) в информативных частотных каналах, совпадающих в обоих приемниках; в информативных частотных каналах, совпадающих в обоих приемниках, рассчитывают k оценок текущих координат источника шума x_k(t₁, ϖ), y_k(t₁, ϖ) по сформированным парам пеленгов на источник шума:

где k - число совпадающих частотных каналов в массивах ϖ₁ и ϖ₂,

x_a, y_a - координаты второго приемника относительно первого,

пеленги γ₁(ϖ)=δ₁-δ₁(ϖ) и γ₂(ϖ)=δ₂-β₂(ϖ) определены через углы прихода потоков энергии β₁(ϖ), β₂(ϖ) и значения положения осей комбинированных приемников относительно севера δ₁, δ₂; выполняют описанный выше цикл обработки текущего временного фрагмента принятого сигнала и для последующих временных фрагментов той же длительности, но опережающих предыдущие на фиксированный интервал времени Δt=t₂-t₁, формируют новые спектральные энергетические портреты принятых сигналов и оценки координат источника x_k(t₂, ϖ), y_k(t₂, ϖ), x_k(t₃, ϖ), y_k(t₃, ϖ)…; по результатам обработки временных фрагментов принятого сигнала проверяют временную устойчивость спектрального портрета сигнала, принятого каждым приемником, сравнивают спектральные энергетические портреты сигналов, принятых первым E₁(t₁, ϖ), E₁(t₂, ϖ), E₁(t₃, ϖ)… и вторым E₂(t₁, ϖ), E₂(t₂, ϖ), E₂(t₃, ϖ)… приемниками, проверяют концентрацию в локальных областях средних значений координатных оценок x(t), y(t) и дисперсии их распределения за время наблюдения, а также контролируют динамику изменения средних значений координатных оценок подбором условий

где V_x, V_y - параметр траектории движения шумящего объекта, который подбирается с учетом априорных сведений о предполагаемом источнике; анализируют полученные результаты обработки и принимают решение об обнаружении источника широкополосного шума по совокупности данных о временной устойчивости спектральных энергетических портретов, полученных приемниками в течение времени наблюдения, и степени их идентичности в каждом приемнике, концентрации координатных оценок на различных частотах в локальной области, размер которой существенно меньше ее дальности от установленных приемников, а также соответствия временной изменчивости средних значений координатных оценок предполагаемой динамике движения источника, и устанавливают местоположение источника по текущим значениям координатных оценок, а по выбранному параметру траектории движения объекта V_x, V_y вычисляют скорость и направление движения источника.

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

Новизна заявленного способа обнаружения подводного источника широкополосного шума заключается в том, что для расширения объема данных об источнике при решении задачи обнаружения предлагается использование двух разнесенных приемников с известной базой, что позволяет получать оценки координат источника на всех выделенных информативных частотах наблюдаемого диапазона, обнаружение источника основывается на предположении, что его спектральные характеристики слабо меняются в процессе обнаружения, и учета того обстоятельства, что поток энергии от источника концентрируется в определенном угловом секторе, значение которого в течение процесса обнаружения может меняться незначительно и соответствовать предполагаемой динамике движения источника, а координатные оценки, рассчитанные для каждой информативной частоте, группируются в локальной области малых размеров, решение об обнаружении источника широкополосного шума принимается на основе увеличенного времени накопления данных и энергетического накопления по всему информативному частотному диапазону, что позволяет повысить помехоустойчивость приемной системы, увеличение вероятности правильного обнаружения обеспечивается последовательной обработкой ряда временных фрагментов принятого сигнала и совместным анализом ее результатов на временную стабильность спектральных энергетических характеристик источника, оценок его местоположения и характеристик движения.

Структурная схема, поясняющая заявленный способ обнаружения подводного источника широкополосного шума, приведена на чертеже, где обозначены следующие элементы:

- I и II - первая и вторая идентичные приемные системы, содержащие каждая:

1 - комбинированный приемник;

2 - блок расчета спектральных характеристик принятого сигнала с усреднением (быстрое преобразование Фурье, например);

3 - измеритель углового положения осей приемника (в простейшем случае компас);

4 - блок передачи системы связи;

- III - пункт обработки данных, содержащий:

5 - блок приема системы связи;

6 - блок расчета величины и углов прихода потока энергии в каждом частотном канале для каждого приемника;

7 - блок формирования углового распределения потоков энергии, выделения угловых секторов с высокой величиной энергии потока, фиксации состава и уровня спектральных компонент, формирующих угловые секторы с высокой величиной энергии потока для каждого приемника;

8 - блок расчета текущих координатных оценок положения подводного источника широкополосного шума;

9 - блок накопления результатов обработки для ряда временных фрагментов принятого сигнала;

10 - блок анализа и принятия решений об обнаружении, оценки спектрального состава шумов источника, его местоположения и параметров движения.

Способ обнаружения подводного источника широкополосного шума осуществляется следующим образом.

В акватории в двух точках с известными координатами выставляют две приемные системы I и II, которые системой связи связаны с пунктом обработки данных III. В каждой точке акватории акустический сигнал от источника, излучающего неизвестный широкополосный шум, синхронизировано принимают комбинированными приемниками 1, с выхода которых принятый сигнал звукового давления акустического поля и трех ортогональных компонент вектора колебательной скорости частиц водной среды (компоненты X и Y составляют горизонтальную плоскость) поступает в блок расчета спектральных характеристик 2.

В блоках 2 приемных систем формируют для принятого сигнала методами частотно-временной обработки сигналов набор частотных каналов ϖ, охватывающий весь наблюдаемый частотный диапазон, (с шириной полосы каналов Δf=1 Гц, например), определяют для принятого сигнала методами частотно-временной обработки в каждом частотном канале текущие комплексные значения амплитуд звукового давления и трех компонент вектора колебательной скорости, по которым вычисляют текущую величину трех компонент вектора плотности потока акустической энергии, вычисляют в каждом частотном канале для текущего временного фрагмента принятого сигнала длительностью Т (Т=10 сек, например) по текущим величинам средние величины трех компонент вектора плотности потока акустической энергии, которые наряду с текущими данными блока 3 об угловой ориентации оси векторного приемника относительно севера, передают блоком 4 в пункт обработки данных III, где принимаются блоком 5.

В блоке 6 пункта обработки для каждого комбинированного приемника приемных систем в каждом сформированном частотном канале для текущего временного фрагмента принятого сигнала длительностью Т вычисляют энергетическую и направленную характеристики горизонтального потока энергии - величину интенсивности Е(ϖ) и углы β(ϖ) прихода потока по формулам:

где I_x, I_y - средние значения горизонтальных компонент вектора плотности потока энергии;

* - знак комплексного сопряжения;

, ϖ_н, ϖ_в - нижняя и верхняя частоты наблюдаемого диапазона;

В блоке 7 для каждого комбинированного приемника приемных систем для текущего временного фрагмента длительностью Т формируют равные горизонтальные угловые секторы шириной Δβ, охватывающие весь горизонт наблюдения 0-360 град, вычисляют для каждого сформированного углового сектора Δβ_i; величину его секторной энергии E₁(Δβ_i), Е₂(Δβ_i) - сумма величин интенсивности потоков энергии Е₁(ϖ), Е₂(ϖ) частотных каналов, имеющих углы β₁(ϖ), β₂(ϖ) прихода потоков энергии в пределах данного углового сектора, и формируют угловое распределение энергии потоков в горизонтальной плоскости, вычисляют среднюю величину секторной энергии в сформированных угловых секторах и назначают ее пороговой величиной для секторов, из сформированных угловых секторов выделяют доминирующие угловые секторы, в которых величина секторной энергии потока превышает пороговую величину для секторов, устанавливают информативные частотные каналы, которые формируют энергию в доминирующих угловых секторах, фиксируют спектральные энергетические портреты принятых в доминирующих угловых секторах сигналов по величине потоков энергии информативных частотных каналов, формируют оценки направлений на источник от каждого приемника β₁(ϖ), β₂(ϖ) в информативных частотных каналах, совпадающих в обоих приемниках.

В блоке 8 в информативных частотных каналах, совпадающих в обоих приемниках, рассчитывают k оценок текущих координат источника шума x_k(t₁, ϖ), по сформированным парам пеленгов на источник шума:

где k - число совпадающих частотных каналов в массивах ϖ₁ и ϖ₂;

x_a, y_a - координаты второго приемника относительно первого;

пеленги γ₁(ϖ)=δ₁-β₁(δϖ) и γ₂(ϖ)=δ₂-β₂(ϖ) определены через угловые приходы энергии β₁(ϖ),β₂(ϖ) и значения положения осей комбинированных приемников относительно севера δ₁, δ₂.

В блоке 9 выполняют описанный выше цикл обработки текущего временного фрагмента принятого сигнала и для последующих временных фрагментов той же длительности Т, но опережающих предыдущие на фиксированный интервал времени Δt=t₂-t₁, формируют новые оценки спектральных портретов и координат источника x_k(t₂, ϖ), y_k(t₂, ϖ), x_k(t₃, ϖ), y_k(t₃, ϖ)…, по результатам обработки временных фрагментов принятого сигнала: проверяют временную устойчивость спектрального портрета сигнала, принятого каждым приемником; сравнивают спектральные энергетические портреты сигналов, принятых первым E₁(t₁, ϖ), E₁(t₂, ϖ), E₁(t₃, ϖ)… и вторым E₂(t₁, ϖ), E₂(t₂, ϖ), E₂(t₃, ϖ)… приемниками; проверяют концентрацию в локальных областях средних значений координатных оценок x(t), y(t) и дисперсии их распределения за время наблюдения, а также контролируют динамику изменения средних значений координатных оценок подбором условий

где V_x, V_y - параметр траектории движения шумящего объекта, который подбирается с учетом априорных сведений о предполагаемом источнике;

Параметры, вычисленные в циклах обработки временных фрагментов принятого сигнала, поступают в блок 10.

В блоке 10 анализируют полученные результаты обработки и принимают решение об обнаружении источника широкополосного шума по совокупности данных о временной устойчивости спектральных энергетических портретов, полученных приемниками в течение времени наблюдения, и степени их идентичности в каждом приемнике, концентрации координатных оценок на различных частотах в локальной области, размер которой существенно меньше ее дальности до установленных приемников, а также соответствия временной изменчивости средних значений координатных оценок предполагаемой динамике движения источника, и устанавливают местоположение источника по текущим значениям координатных оценок, а по выбранному параметру траектории движения объекта V_x, V_y вычисляют скорость V_ист и направление ϕ_фист движения источника.