Документация

Импортируйте данные о ландшафте

Импортируйте данные о ландшафте DTED-формата для области вокруг Булдера, Колорадо, США. Файл ландшафта был загружен с "SRTM Пусто Заполненный" набор данных, доступный от Геологической службы США (USGS). Файл является форматом уровня 1 DTED и имеет разрешение выборки приблизительно 90 метров. Один файл DTED задает область, которая охватывает 1 степень и в области широты и в области долготы.

dtedfile = "n39_w106_3arc_v2.dt1";
attribution = "SRTM 3 arc-second resolution. Data available from the U.S. Geological Survey.";
addCustomTerrain("southboulder",dtedfile, ...
    "Attribution",attribution)

Открытое Средство просмотра Сайта с помощью импортированного ландшафта. Визуализация со спутниковым формированием изображений карты с высоким разрешением требует Интернет-соединения.

viewer = siteviewer("Terrain","southboulder");

Покажите радарные сайты кандидата

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

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

names = "Radar site" + (1:5);
rdrlats = [39.6055 39.6481 39.7015 39.7469 39.8856];
rdrlons = [-105.1602 -105.1378 -105.1772 -105.2000 -105.2181];

% Create transmitter sites associated with radars
rdrtxs = txsite("Name",names, ...
    "AntennaHeight",10, ...
    "Latitude",rdrlats, ...
    "Longitude",rdrlons);

% Create receiver sites associated with radars
rdrrxs = rxsite("Name",names, ...
    "AntennaHeight",10, ...
    "Latitude",rdrlats, ...
    "Longitude",rdrlons);

% Show just the radar transmitter sites
show(rdrtxs);

Масштабируйте и вращайте карту, чтобы просмотреть 3-D ландшафт вокруг радарных сайтов кандидата. Выберите сайт, чтобы просмотреть местоположение, высоту антенны и наземное повышение.

Разработайте моностатическую радиолокационную систему

Разработайте основную моностатическую импульсную радиолокационную систему, чтобы обнаружить неколеблющиеся цели с радарным сечением (RCS) на 0,1 квадратных метра на расстоянии до 35 000 метров от радара с разрешением области значений 5 метров. Желаемый индекс производительности является вероятностью обнаружения (Фунт) 0,9 и вероятность ложного предупреждения (PFA) ниже 1e-6. Радары приняты, чтобы быть поворотными и поддержать то же усиление антенны во всех направлениях, где усиление антенны соответствует очень направленной антенной решетке.

pd = 0.9;         % Probability of detection
pfa = 1e-6;       % Probability of false alarm
maxrange = 35000; % Maximum unambiguous range (m)
rangeres = 5;     % Required range resolution (m)
tgtrcs = .1;      % Required target radar cross section (m^2)

Используйте импульсное интегрирование, чтобы уменьшать необходимый ОСШ в радарном получателе. Используйте 10 импульсов и вычислите ОСШ, требуемый обнаружить цель.

numpulses = 10;
snrthreshold = albersheim(pd, pfa, numpulses); % Unit: dB
disp(snrthreshold);

    4.9904

Задайте радарную частоту центра и усиление антенны, приняв очень направленную антенную решетку.

fc = 10e9;    % Transmitter frequency: 10 GHz
antgain = 38; % Antenna gain: 38 dB
c = physconst('LightSpeed');
lambda = c/fc;

Вычислите необходимую пиковую импульсную степень (Ватты) радарного передатчика с помощью основного уравнения радиолокации.

pulsebw = c/(2*rangeres);
pulsewidth = 1/pulsebw;
Ptx = radareqpow(lambda,maxrange,snrthreshold,pulsewidth,...
    'RCS',tgtrcs,'Gain',antgain);
disp(Ptx)

   3.1521e+05

Задайте целевые положения

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

% Define region of interest
latlims = [39.5 40];
lonlims = [-105.6 -105.1];

% Define grid of target locations in region of interest
tgtlatv = linspace(latlims(1),latlims(2),50);
tgtlonv = linspace(lonlims(1),lonlims(2),50);
[tgtlons,tgtlats] = meshgrid(tgtlonv,tgtlatv);
tgtlons = tgtlons(:);
tgtlats = tgtlats(:);

Вычислите минимум, максимум, и имейте в виду наземное повышение для целевых местоположений.

% Create temporary array of sites corresponding to target locations and query terrain
Z = elevation(txsite("Latitude",tgtlats,"Longitude",tgtlons));
[Zmin, Zmax] = bounds(Z);
Zmean = mean(Z);

disp("Ground elevation (meters): Min    Max    Mean" + newline + ...
    "                           " + round(Zmin) + "   " + round(Zmax) + "   " + round(Zmean))

Ground elevation (meters): Min    Max    Mean
                           1253   3953   2373

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

% Target altitude above ground level (m)
tgtalt = 500;

Покажите видимую область твердой зеленой зоной на карте.

viewer.Name = "Radar Coverage Region of Interest";
helperContourMap(rdrtxs, tgtlats, tgtlons)

Вычислите ОСШ для целевых положений с ландшафтом

Основное уравнение радиолокации включает потерю пути к свободному пространству и имеет параметр за дополнительные потери. Используйте модель распространения ландшафта, чтобы предсказать дополнительную потерю пути по ландшафту. Используйте Ландшафт Интегрированная Грубая Земля Model™ (TIREM™) от Науки Alion, если это доступно, или иначе используйте Лонгли-Райса (иначе ITM) модель. TIREM™ поддерживает частоты до 1 000 ГГц, тогда как Лонгли-Райс - допустимых до 20 ГГц. Вычислите общую дополнительную потерю включая распространение от радара до цели и затем обратно с цели на получатель.

% Create a terrain propagation model, using TIREM or Longley-Rice
tiremloc = tiremSetup;
if ~isempty(tiremloc)
    pm = propagationModel('tirem');
else
    pm = propagationModel('longley-rice');
end

% Compute additional path loss due to terrain and return distances between radars and targets 
[L, ds] = helperPathlossOverTerrain(pm, rdrtxs, rdrrxs, tgtlats, tgtlons, tgtalt);

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

% Compute SNR for all radars and targets
numtgts = numel(tgtlats);
numrdrs = numel(rdrtxs);
rawsnr = zeros(numtgts,numrdrs);
for tgtind = 1:numtgts
    for rdrind = 1:numrdrs
        rawsnr(tgtind,rdrind) = radareqsnr(lambda,ds(tgtind,rdrind),Ptx,pulsewidth, ...
            'Gain',antgain,'RCS',tgtrcs,'Loss',L(tgtind,rdrind));
    end
end

Оптимизируйте радарное покрытие

Цель обнаруживается, если радарный ОСШ получателя превышает порог ОСШ, вычисленный выше. Рассмотрите все комбинации радарных сайтов и выберите три сайта, которые производят самое большое количество обнаружений. Вычислите данные об ОСШ как лучший ОСШ, доступный в получателе любого из выбранных радарных сайтов.

bestsitenums = helperOptimizeRadarSites(rawsnr, snrthreshold);
snr = max(rawsnr(:,bestsitenums),[],2);

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

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

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

% Show selected radar sites using blue markers
hide(rdrtxs(bestsitenums))
show(rdrrxs(bestsitenums))

% Plot radar coverage
viewer.Name = "Radar Coverage";
legendTitle = "SNR" + newline + "(dB)";
helperContourMap(rdrtxs, tgtlats, tgtlons, snr, snrthreshold, legendTitle)

Отличайтесь количество импульсов, чтобы объединяться

Анализ выше оптимизированной радарной степени передатчика и местоположений сайта на основе системы, которая интегрирует 10 импульсов. Теперь исследуйте влияние на радарное покрытие для различных режимов работы системы, где количество импульсов, чтобы объединяться отличается. Вычислите пороги ОСШ, требуемые обнаружить цель для переменного количества импульсов.

% Calculate SNR thresholds corresponding to different number of pulses
numpulses = 1:10;
snrthresholds = zeros(1,numel(numpulses));
for k = 1:numel(numpulses)
    snrthresholds(k) = albersheim(pd, pfa, numpulses(k));
end

% Plot SNR thresholds vs number of pulses to integrate
plot(numpulses,snrthresholds,'-*')
title("SNR at Radar Receiver Required for Detection")
xlabel("Number of pulses to integrate")
ylabel("SNR (dB)")
grid on;

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

viewer.Name = "Radar Coverage for Multiple SNR Thresholds";
helperContourMap(rdrtxs, tgtlats, tgtlons, snr, snrthresholds([1 2 5 10]), legendTitle)

Обновите целевую высоту

Обновите сценарий так, чтобы целевые положения составили 250 метров над уровнем земли вместо 500 метров над уровнем земли. Повторно выполните тот же анализ как выше, чтобы выбрать три лучших радарных сайта и визуализировать покрытие. Новая карта покрытия показывает, что сокращение видимости целей также уменьшает зону охвата.

% Target altitude above ground (m)
tgtalt = 250;
[L, ds] = helperPathlossOverTerrain(pm, rdrtxs, rdrrxs, tgtlats, tgtlons, tgtalt);

% Compute SNR for all radars and targets
numrdrs = numel(rdrtxs);
rawsnr = zeros(numtgts,numrdrs);
for tgtind = 1:numtgts
    for rdrind = 1:numrdrs
        rawsnr(tgtind,rdrind) = radareqsnr(lambda,ds(tgtind,rdrind),Ptx,pulsewidth, ...
            'Gain',antgain,'RCS',tgtrcs,'Loss',L(tgtind,rdrind));
    end
end

% Select best combination of 3 radar sites
bestsitenums = helperOptimizeRadarSites(rawsnr, snrthreshold);
snr = max(rawsnr(:,bestsitenums),[],2);

% Restore original markers
hide(rdrrxs)
show(rdrtxs)

% Show selected radar sites using blue markers
hide(rdrtxs(bestsitenums))
show(rdrrxs(bestsitenums))

% Plot radar coverage
viewer.Name = "Radar Coverage";
helperContourMap(rdrtxs, tgtlats,tgtlons, snr, snrthreshold, legendTitle)

Покажите радарную карту покрытия для нескольких порогов ОСШ.

viewer.Name = "Radar Coverage for Multiple SNR Thresholds";
helperContourMap(rdrtxs, tgtlats,tgtlons, snr, snrthresholds([1 2 5 10]), legendTitle)

Заключение

Моностатическая радиолокационная система была разработана, чтобы обнаружить неколеблющиеся цели с радарным сечением (RCS) на 0,1 квадратных метра на расстоянии до 35 000 метров. Радарные сайты были выбраны среди пяти сайтов кандидата, чтобы оптимизировать количество обнаружений по видимой области. Две целевых высоты были рассмотрены: 500 метров над уровнем земли и 250 метров над уровнем земли. Карты покрытия предлагают важность видимости угла обзора между радаром и целью в порядке достигнуть обнаружения. Второй сценарий приводит к целям, которые ближе к уровню земли и поэтому более вероятно быть блокированными от видимости угла обзора с радаром. Это видно путем вращения карты, чтобы просмотреть ландшафт, где незоны охвата обычно располагаются в теневых областях гор.

Вымойтесь заключительным Средством просмотра Сайта и удалением импортированных данных о ландшафте.

close(viewer)
removeCustomTerrain("southboulder")