Учебное пособие 996
.pdfТаблица 4.1
Предельно допустимая загруженность ν диапазона частот, при которой в нем с вероятностью 99 % наблюдается не менее
|
Lqw |
смежных свободных радиоканалов |
|
|
|||||||
|
Число |
|
Предельно допустимая загруженность ν (%) |
||||||||
|
свободных |
|
|
при общем числе каналов L |
|
||||||
|
каналов Lqw |
|
|
|
|
|
|
|
|
|
|
|
|
600 |
300 |
150 |
100 |
|
80 |
|
62 |
50 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
16 |
|
17 |
13 |
9 |
- |
|
- |
|
- |
- |
|
12 |
|
25 |
20 |
14 |
10 |
|
8 |
|
- |
- |
|
8 |
|
38 |
32 |
25 |
21 |
|
18 |
|
14 |
- |
|
7 |
|
44 |
37 |
30 |
25 |
|
22 |
|
19 |
12 |
|
6 |
|
50 |
43 |
36 |
30 |
|
28 |
|
24 |
20 |
|
5 |
|
58 |
50 |
43 |
38 |
|
35 |
|
31 |
28 |
|
4 |
|
67 |
60 |
53 |
47 |
|
45 |
|
42 |
36 |
|
3 |
|
78 |
72 |
65 |
60 |
|
58 |
|
53 |
50 |
|
2 |
|
89 |
87 |
82 |
77 |
|
76 |
|
73 |
69 |
Если же ширина радиоканала составит лишь |
B = 2 |
отсчета, то |
|||||||||
число радиоканалов в диапазоне будет равно |
L = 300 , а на тот |
же свободный промежуток шириной 24 отсчета (т.е. 12 радиоканалов) можно рассчитывать лишь при загруженности ν ≤ 20%.
При выборе ширины участка радиодиапазона, на котором предполагается производить оценку интенсивности шума, следует принимать во внимание, что в условиях сложной радиообстановки целесообразно ориентироваться на поиск относительно небольшого по протяженности свободного от сигналов участка, который обязан практически наверняка присутствовать в диапазоне частот и, вместе с тем, ширина которого достаточна для получения оценки шума с не слишком высокой дисперсией. Для массива спектра XR (n) объемом от нескольких сотен до несколь-
ких тысяч отсчетов ширина подобного участка может составлять до нескольких десятков отсчетов.
21
Информация о минимальном размере массива спектральных отсчетов XR (n) , позволяющем при известной загру-
женности гарантировать наличие в полосе частот хотя бы одного промежутка между спектрами сигналов шириной QW
отсчетов, представлена на рис. 4.2. При формировании рис. 4.2 учтено, что для повышения скорости просмотра диапазона частот при панорамном обнаружении ширину радиоканала B обычно стараются минимизировать и приоритетными для исследований являются значения ширины радиоканала из диапазона B = 4…8 отсчетов [4]. Отметим, что согласно рис. 4.2 даже при анализе данных объемом N ≤ 103 отсчетов и при заметной занятости диапазона частот можно рассчитывать на существование между спектрами сигналов хотя бы одного шумового промежутка размером QW < 20 отсчетов.
Рис. 4.2. Размер массива спектральных данных, минимально необходимый для существования в диапазоне частот
22
4.1.2. Поиск расположения шумового участка по минимальному значению сглаженного спектра
Для оценки расположения в радиодиапазоне свободного от сигналов частотного интервала будем ориентироваться на минимальное значение спектра, сглаженного прямоугольным окном. Для пояснения данного утверждения рассмотрим сглаженный по частоте энергетический спектр
|
|
1 |
|
n+int[ WW 2 ] |
|
|
|
XɶR |
(n) = |
|
∑ |
XR (i), |
(4.2) |
||
|
|||||||
|
WW |
i = n−int (WW −1) 2 |
|
||||
|
|
|
|
|
|
|
|
где WW – ширина окна сглаживания, выбираемая по правилу |
|||||||
|
WW ≤ QW (N, B,v), |
|
(4.3) |
QW ( ) – прогнозируемая ширина максимального промежутка
между спектрами обнаруживаемых сигналов. Возможный вид сглаженного спектра представлен на рис. 4.3. Каждый отсчет
XɶR (n) этого спектра представляет собой оценку математического ожидания отсчетов спектра XR (n) на интервале шири-
ной WW с центром в n -м отсчете, а на сигнальных участках превышает математическое ожидание на участках шумовых. Это значит, что при наличии в анализируемом спектре хотя бы одного свободного от сигналов участка шириной не менее WW отсчетов, начало этого (или какого-либо аналогичного из шумовых) участков может быть найдено на основе правила
|
|
|
WW −1 |
|
|||
nш |
= (argminn |
XɶR |
(n))−int |
|
. |
(4.4) |
|
2 |
|||||||
|
|
|
|
|
|
Следует, однако, иметь в виду, что участок оси частот, начало которого определяется правилом (4.4), выбирается из, возможно, многих шумовых промежутков за счет поиска минимального значения. Соответственно, среднее арифметическое отсчетов на участке, стартующем с nш , будет равно
ɶ |
ɶ |
(n) . |
(4.5) |
xmin |
= min XR |
||
|
n |
|
|
23
XR(n) |
dn1 |
|
|
|
|
dn2 |
dnM |
|
|
|
|
|
1 |
|
|
|
2 |
|
|
σξ2 N |
|
|
|
|
|
|
|
|
n |
nmin |
n1 |
n2 |
nM |
N |
1 – приближенный результат сглаживания дискретного спектра, 2 – центр участка, который будет выбран в качестве шумового
Рис. 4.3. Результат сглаживания дискретного спектра мощности наблюдаемого случайного процесса
4.2. Образцы решения задач
Задача 1. В энергетическом спектре (3.2), используемом для контроля радиообстановки в диапазоне частот, отдельные радиоканалы представлены 5 спектральными отсчетами. Дайте рекомендации по поиску участка оси частот, на котором можно осуществлять оценивание интенсивности шума, при разных соотношениях между шириной этого диапазона и его загруженностью
ν .
Решение
а) Независимо от параметров радиообстановки для поиска шумового участка можно рекомендовать соотношение (4.4), однако требует внимания выбор ширины окна сглаживания WW.
б) Рассмотрим фрагмент приведенной выше табл. 4.1:
24
Число |
Предельно допустимая загруженность ν (%) |
||||||
свободных |
|
при общем числе каналов L |
|
||||
каналов Lqw |
600 |
300 |
150 |
100 |
80 |
62 |
50 |
|
|
|
|
|
|
|
|
16 |
17 |
13 |
9 |
- |
- |
- |
- |
12 |
25 |
20 |
14 |
10 |
8 |
- |
- |
8 |
38 |
32 |
25 |
21 |
18 |
14 |
- |
7 |
44 |
37 |
30 |
25 |
22 |
19 |
12 |
6 |
50 |
43 |
36 |
30 |
28 |
24 |
20 |
|
|
|
|
|
|
|
|
Рис. 4.4 Фрагмент таблицы, содержащей значения предельно допустимой загруженности ν диапазона частот
Выделенные цветом ячейки соответствуют случаям, когда загруженность диапазона не превышала 20%. Умножив количество Lqw каналов в ожидаемом шумовом интервале и общую ши-
рину диапазона в L каналов на B = 5, получим зависимость, представленную на рис. 4.4 линией красного цвета; эта линия характеризует рекомендации для случая ν = 20%.
Результаты аналогичной операция по ячейкам табл. 4.1 с загруженностью 30% отображены на рис. 4.4 линией черного цвета, а зависимость между шириной окна сглаживания спектральных отсчетов WW и соответствующей ему минимально необходимой шириной N анализируемого диапазона частот для случая ν = 50% показаны на рис. 4.4 линией синего цвета.
Полученные зависимости наглядно показывают, что в значительно загруженных диапазонах на наличие широких промежутков между спектрами можно рассчитывать лишь для очень широких, представленных тысячами спектральных отсчетов частотных диапазонов.
Отметим, что ту же самую задачу можно было бы решать и взяв за основу рис. 4.2. К примеру, линия «QW = 30, B = 5»
25
Рис. 4.5 Минимальная ширина спектра N, необходимая для безопасного использования окна сглаживания шириной WW
на рис. 4.2 позволяет сформировать рекомендации: «ν = 30%N = 500; ν = 40% N = 1200», а линия «QW = 20, B = 5» указывает на рекомендуемые соотношения: «ν = 50% N = 600; ν = 60% N = 1500». Соответствующие метки фиолетового цвета, отображенные на рис. 4.5, показывают согласованность разных подходов к решению текущей задачи. Отличие состоит в том, что решение с опорой на табл. 4.1 позволяет получать соотношения для окон и диапазонов большой ширины, а на рис. 4.2 представлены рекомендации, полезные для грамотного выбора параметров обработки для диапазонов частот, представленных сотнями спектральных отсчетов при окнах сглаживания шириной не более QW = 30 отсчетов.
Задача 2. Контроль радиообстановки в полосе частот шириной Шf = 10 МГц базируется на обработке выборок отсчетов, дискретизируемых с частотой 25 МГц. Выборки объемом 4096 отсчетов подвергаются быстрому преобразованию Фурье (БПФ), и за счет накопления нескольких подобных спектров
26
формируется усредненный энергетический спектр (3.2). В этом спектре производится поиск участка, наиболее вероятно занятого лишь шумом, за счет использования оценки (4.4) с шириной окна сглаживания WW = 16 отсчетов. Отдельные радиоканалы, образующие анализируемую полосу частот, имеют ширину 25 кГц. Определить, при какой максимально возможной загруженности может стать значимой вероятность отсутствия в полосе частот промежутка QW между спектрами радиостанций, превышающего по ширине WW?
Решение
а) При заданной в условии частоте дискретизации FД = = 25 МГц и размерности БПФ NБПФ = 4096 отсчетов получаемый дискретный спектр согласно (3.3) будет характеризоваться интервалом между соседними спектральными отсчетами (бином)
δ f = |
FД |
= |
25000 |
≈ 6,1кГц . |
|
NБПФ |
4096 |
||||
|
|
|
Соответственно, ширина радиоканалов в отсчетах составит
B = |
ШРК |
= |
25 |
|
≈ 4 отсчета . |
|
6,1 |
||||
|
δ f |
|
б) Размерность БПФ NБПФ = 4096 отсчетов соответствует частотному интервалу, определяемому частотой дискретизации FД, а занятая сигналами полоса частот заметно у́же и оказывается представленной лишь
N = N |
|
Шполосы |
= 4096 |
10 |
|
≈ 1638отсчетами . |
|
БПФ |
F |
25 |
|||||
|
|
|
|||||
|
|
Д |
|
|
|
|
в) Если для оценки взаимосвязи между размером спектрального массива N и загруженностью ν диапазона частот использовать рис. 4.2, то анализируемый случай будет соответствовать самой левой из сиреневых линий с пометкой «QW = 16, B
=4». Прослеживая пересечение этой линии с уровнем N =
=1638 нетрудно заметить, что абсциссой этого пересечения служит загруженность ν ≈ 62%. Таким образом, при загруженности полосы частот, не превышающей 60%, вероятность отсут-
27
ствия в полосе частот промежутка QW > WW можно считать пренебрежимо малой, а вот при большей загруженности риск возникновения аномальных ошибок выявления шумовых участков спектра начнет нарастать.
4.3. Индивидуальные контрольные задания
Задача 4.1. При обеспечиваемом спектроанализатором частотном разрешении ширина радиоканалов в анализируемой полосе частот составляет B отсчетов и указана в левой колонке табл. 4.3, а объем массива, представляющего всю полосу частот указан в правой колонке этой таблицы. Поиск шумового участка для определения интенсивности шума производится по минимуму сглаженного энергетического спектра, формируемого на основе (4.2) с помощью окна шириной WW = 18 отсчетов. При какой загруженности ν в полосе частот с вероятностью не менее 99% должен обнаружиться хотя бы один промежуток между спектрами сигналов, превышающий по ширине WW?
|
|
Таблица 4.3 |
Но- |
Ширина радиоканала |
Ширина анализируемой |
мер |
B, отсчетов |
полосы частот N, отсчетов |
|
|
|
1 |
6 |
1400 |
|
|
|
2 |
4 |
1024 |
|
|
|
3 |
2 |
720 |
|
|
|
4 |
3 |
2048 |
|
|
|
5 |
7 |
1024 |
|
|
|
6 |
5 |
512 |
|
|
|
7 |
6 |
400 |
|
|
|
8 |
9 |
512 |
|
|
|
9 |
6 |
1480 |
|
|
|
10 |
4 |
256 |
|
|
|
28
Задача 4.2. Для анализа радиообстановки, параметры которой заданы в табл. 4.4, используется энергетический спектр (3.2) с интервалом между отсчетами δ f, указанным в левой колонке той же таблицы. Если для поиска шумового промежутка в полосе частот предполагается использовать оценку (4.4), полученную с помощью окна сглаживания шириной QW = 30 отсчетов, то при какой минимальной ширине полосы частот наличие промежутка, превышающего по ширине QW, можно считать практически гарантированным?
|
|
|
Таблица 4.4 |
|
Но- |
Интервал между отс- |
Ширина радио- |
Загружен- |
|
мер |
четами спектра (бин) |
канала Шрк, кГц |
ность ν , % |
|
|
δ f , кГц |
|
|
|
1 |
2,5 |
25 |
60 |
|
2 |
1,5 |
12,5 |
52 |
|
|
|
|
|
|
3 |
2 |
12,5 |
44 |
|
|
|
|
|
|
4 |
5 |
25 |
40 |
|
|
|
|
|
|
5 |
3,2 |
12,5 |
36 |
|
|
|
|
|
|
6 |
20 |
200 |
70 |
|
|
|
|
|
|
7 |
25 |
200 |
60 |
|
|
|
|
|
|
8 |
20 |
120 |
50 |
|
|
|
|
|
|
9 |
1,25 |
6,25 |
44 |
|
|
|
|
|
|
10 |
1,6 |
6,25 |
32 |
|
|
|
|
|
|
Задача 4.3. Объем N массива XR (n) отражающего радио-
обстановку в диапазоне частот задан в левой колонке табл. 4.5; ширина радиоканала B – в центральной колонке той же таблицы. Для поиска шумового участка на оси частот используется окно шириной WW, указанной в правой колонке табл. 4.5. С помощью табл. 4.1 определите, при какой загруженности диапазона ν можно рассчитывать на то, что в диапазоне частот с вероятностью не менее 99% будет наблюдаться хотя бы один шумовой интервал, превышающий по ширине WW.
29
|
|
|
Таблица 4.5 |
|
Но- |
Размер спек- |
Ширина ра- |
Ширина окна поиска |
|
мер |
трального |
диоканала B, |
шумового интервала |
|
массива N |
отсчетов |
WW, отсчетов |
||
|
||||
|
|
|
|
|
1 |
600 |
2 |
16 |
|
|
|
|
|
|
2 |
1800 |
3 |
24 |
|
|
|
|
|
|
3 |
256 |
4 |
20 |
|
|
|
|
|
|
4 |
512 |
5 |
30 |
|
|
|
|
|
|
5 |
1800 |
6 |
30 |
|
|
|
|
|
|
6 |
1200 |
2 |
24 |
|
|
|
|
|
|
7 |
1024 |
3 |
24 |
|
|
|
|
|
|
8 |
1200 |
4 |
20 |
|
|
|
|
|
|
9 |
750 |
5 |
20 |
|
|
|
|
|
|
10 |
3600 |
6 |
18 |
|
|
|
|
|
Задача 4.4. В математической среде Mathcad, MATLAB или какой-либо среде программирования сформируйте моделирующую программу для исследования степени смещения статистики xɶmin , определяемой (4.5), при условии, что ширина WW
окна сглаживания определяется в левой колонке табл. 4.6, а параметры радиообстановки определены в последующих колонках той же таблицы.
|
|
|
|
Таблица 4.6 |
|
Но- |
Окно сгла- |
Размер спек- |
Ширина ра- |
|
Загру- |
живания WW, |
трального |
диоканала B, |
|
жен- |
|
мер |
отсчетов |
массива N |
отсчетов |
|
ность ν |
|
|
||||
|
|
|
|
|
|
1 |
12 |
1024 |
8 |
|
10% |
|
|
|
|
|
|
2 |
12 |
4096 |
3 |
|
50% |
|
|
|
|
|
|
3 |
16 |
1024 |
4 |
|
30% |
|
|
|
|
|
|
4 |
16 |
4096 |
4 |
|
40% |
|
|
|
|
|
|
30