Swerling 1 по сравнению с Swerling 2 модели

Для Swerling 1 и целевых моделей Swerling 2, общий ЭПР является результатом многих независимых маленьких рассеивателей приблизительно равного отдельного ЭПР. Общий ЭПР может меняться в зависимости от каждого импульса в скане (Swerling 2) или может быть постоянным по полному скану, состоящему из нескольких импульсов (Swerling 1). В любом случае статистические данные выполняют функцию плотности вероятности в квадрате хи с двумя степенями свободы.

Время задержки и радиолокационный обзор

Для простоты начните с вращающегося радара, имеющего время вращения 5 секунд, соответствуя уровню вращения или частоте развертки 72 степеней/секунда.

Trot = 5.0;
rotrate = 360/Trot;

Радар имеет основную ширину луча на уровне половинной мощности (HPBW) 3,0 градусов. В течение времени, когда цель освещается основным лучом, радарные импульсы поражают цель и отражаются назад к радару. Период времени, во время которого освещается цель, называется временем задержки. Этот период времени также называется сканом. Пример обработает 3 скана цели.

HPBW = 3.0;
Tdwell = HPBW/rotrate;
Nscan = 3;

Количество импульсов, которые прибывают в цель во время времени задержки, зависит от импульсной частоты повторения (PRF). PRF является инверсией импульсного интервала повторения (PRI). Примите, что 5 000 импульсов передаются в секунду.

prf = 5000.0;
pri = 1/prf;

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

Np = floor(Tdwell*prf);

Настройте радарную Swerling 1 модель

Вы создаете цель Swerling 1 путем соответствующего использования step метод RadarTarget Система object™. Чтобы произвести модель Swerling 1, установите Model свойство phased.RadarTarget Системный объект к любому 'Swerling1' или 'Swerling2'. Оба эквивалентны. Затем в первом вызове step метод в начале скана, набор updatercs аргумент к true. Установите updatercs к false для остающихся вызовов step во время скана. Это означает, что эффективная площадь рассеивания только обновляется в начале скана и остается постоянной для остатка от скана.

Установите целевую модель на 'Swerling1'.

tgtmodel = 'Swerling1';

Настройте радарные компоненты Системного объекта модели

Настройте исходящую антенну. Примите, что рабочая частота антенны составляет 1 ГГц.

fc = 1e9;
antenna = phased.IsotropicAntennaElement('BackBaffled',true);
radiator = phased.Radiator('OperatingFrequency',fc, ...
    'Sensor',antenna);

Задайте местоположение стационарной антенны.

radarplatform = phased.Platform('InitialPosition',[0;0;0]);

Задайте местоположение стационарной цели.

targetplatform = phased.Platform('InitialPosition',[2000; 0; 0]);

Переданный сигнал является линейной формой волны FM. Передайте один импульс на вызов step метод.

waveform = phased.LinearFMWaveform('PulseWidth',50e-6,...
    'OutputFormat','Pulses','NumPulses',1);

Настройте усилитель передачи.

transmitter = phased.Transmitter('PeakPower',1000.0,'Gain',40);

Настройте среду распространения, чтобы быть свободным пространством.

channel = phased.FreeSpace('OperatingFrequency',fc,'TwoWayPropagation',true);

Задайте радарную цель, чтобы иметь средний ЭПР 1 m2 и быть типа модели Swerling 1 или 2. Можно использовать Swerling 1 или 2 взаимозаменяемо.

target = phased.RadarTarget('MeanRCS',1,'OperatingFrequency',fc,...
    'Model',tgtmodel);

Настройте радарный коллектор.

collector = phased.Collector('OperatingFrequency',1e9,...
    'Sensor',antenna);

Задайте согласованный фильтр, чтобы обработать входящий сигнал.

wav = step(waveform);
filter = phased.MatchedFilter('Coefficients',getMatchedFilter(waveform));

Обработка цикла для 3 сканов цели Swerling 1

Обеспечьте память для радара возвращают амплитуды.

z = zeros(Nscan,Np);
tp = zeros(Nscan,Np);

Введите цикл. Установите updatercs к true только для первого импульса скана.

for m = 1:Nscan
    t0 = (m-1)*Trot;
    t = t0;
    for k = 1:Np

        if k == 1
            updatercs = true;
        else
            updatercs = false;
        end
        t = t + pri;
        txwav = transmitter(wav);

Найдите радар и предназначайтесь для положений

        [xradar,vradar] = radarplatform(t);
        [xtgt,vtgt] = targetplatform(t);

Излучите форму волны, чтобы предназначаться

        [~,ang] = rangeangle(xtgt,xradar);
        radwav = radiator(txwav,ang);

Распространите форму волны к и от цели

        propwav = channel(radwav,xradar,xtgt,vradar,vtgt);

Отразите форму волны от цели. Установите updatercs флаг.

        reflwav = target(propwav,updatercs);

Соберите принятую форму волны

        collwav = collector(reflwav,ang);

Примените согласованный фильтр к входящему сигналу

        y = filter(collwav);
        z(m,k) = max(abs(y));
        tp(m,k) = t;

    end
end

Постройте импульсные амплитуды

Постройте амплитуды импульсов для скана в зависимости от времени.

plot(tp(:),z(:),'.')
xlabel('Time (sec)')
ylabel('Pulse Amplitudes')

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