Изобретение относится к области дискретного спектрального анализа, к области систем обработки информации и измерительной техники, и может быть использовано для доплеровской фильтрации (выделения) лучевой структуры ионосферных сигналов. Способ включает прием ионосферного сигнала с помощью приемника, усиление и преобразование по частоте, преобразование аналогового сигнала на выходе приемника в цифровую форму с помощью аналого-цифрового преобразователя. При этом снижают промежуточную частоту сигнала до нулевого значения с помощью преобразования Фурье. Формируют элементы корреляционной матрицы
A
m
,
p
=
Y
^
n
+
m
−
1
Y
^
n
+
p
−
1
∗
¯
(черта сверху означает суммирование по индексу n) и правого столбца
b
^
p
=
Y
^
n
+
p
−
1
e
−
i
ω
t
n
¯
векторного уравнения
A
^
c
¯
=
b
¯
. Определяют элементы
c
^
m
вектора неопределенных коэффициентов
c
¯
,
решая в вычислительном устройстве векторное уравнение
A
^
c
¯
=
b
¯
. Формируют частотную зависимость функционала правдоподобия
Δ
(
ω
)
=
1
1
−
∑
m
=
1
M
c
^
m
Y
^
n
+
m
−
1
∗
перебирая частоты с заданным шагом в заданном интервале частот. Оценивают частоты доплеровских составляющих ионосферного сигнала и их достоверность по максимумам частотной зависимости функционала правдоподобия Δ(ω). Формируют элементы матрицы
A
^
1
m
,
p
=
e
i
(
ω
m
−
ω
p
)
t
n
¯
и правого столбца
b
^
1
m
=
Y
^
n
e
−
i
ω
m
t
n
¯
векторного уравнения
A
^
1
U
¯
=
b
¯
1
.
Оценивают комплексные амплитуды доплеровских составляющих ионосферного сигнала
U
^
1
÷
U
^
M
,
решая в вычислительном устройстве векторное уравнение
A
^
1
U
¯
=
b
¯
1
.
Технический результат заключается в
повышении точности и достоверности оценок доплеровского спектра многолучевых ионосферных сигналов, и в расширении возможностей спектрального анализа в область малых интервалов обработки сигналов, где критерий Рэлея не выполняется. 6 ил.