Заградительный передатчик помех

Поддержка моделирования заградительного передатчика помех

barrageJammer объектные модели широкополосный передатчик помех. Выход barrageJammer комплексная последовательность белого Гауссова шума. Модифицируемые свойства заградительного передатчика помех:

  • ERP — Эффективная излучаемая мощность в ваттах

  • SamplesPerFrameSource — Источник количества отсчетов на систему координат

  • SamplesPerFrame — Количество отсчетов на систему координат

  • SeedSource — Источник seed для генератора случайных чисел

  • Seed — Отберите для генератора случайных чисел

Действительные и мнимые части комплексной последовательности белого Гауссова шума у каждого есть отклонение, равное 1/2 эффективная излучаемая мощность в ваттах. Обозначьте эффективную излучаемую мощность в ваттах P. Заградительный передатчик помех выход:

w[n]=P2x[n]+jP2y[n]

В этом уравнении, x[n] и y[n] некоррелированые последовательности нулевых средних Гауссовых случайных переменных с модульным отклонением.

Заградительный передатчик помех модели Выход

Этот пример исследует статистические свойства заградительного передатчика помех выход и как они относятся к эффективной излучаемой мощности (ERP). Создайте заградительный передатчик помех с помощью эффективной излучаемой мощности 5 000 ватт. Сгенерируйте выход на 500 выборках на систему координат. Затем вызовите step функционируйте однажды, чтобы сгенерировать одну систему координат комплексных данных. Используя histogram функция, покажите распределение заградительных выходных значений передатчика помех. BarrageJammer Системный объект использует генератор случайных чисел. В этом примере начальное значение генератора случайных чисел фиксируется в иллюстративных целях и может быть удалено.

rng default
jammer = barrageJammer('ERP',5000,...
    'SamplesPerFrame',500);
y = jammer();
subplot(2,1,1)
histogram(real(y))
title('Histogram of Real Part')
subplot(2,1,2)
histogram(imag(y))
title('Histogram of Imaginary Part')
xlabel('Watts')

Figure contains 2 axes objects. Axes object 1 with title Histogram of Real Part contains an object of type histogram. Axes object 2 with title Histogram of Imaginary Part contains an object of type histogram.

Средние значения действительных и мнимых частей

mean(real(y))
ans = -1.0961
mean(imag(y))
ans = -2.1671

которые являются эффективно нулевыми. Стандартные отклонения действительных и мнимых частей

std(real(y))
ans = 50.1950
std(imag(y))
ans = 49.7448

которые соглашаются с ожидаемым значением ERP/2.

Эффект модели заградительного передатчика помех на целевом эхе

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

antenna = phased.ULA(4); 
Fs = 1e6;
fc = 1e9;
rng('default')
waveform = phased.RectangularWaveform('PulseWidth',100e-6,...
    'PRF',1e3,'NumPulses',5,'SampleRate',Fs);
transmitter = phased.Transmitter('PeakPower',1e4,'Gain',20,...
    'InUseOutputPort',true);
radiator = phased.Radiator('Sensor',antenna,'OperatingFrequency',fc);
jammer = barrageJammer('ERP',1000,...
    'SamplesPerFrame',waveform.NumPulses*waveform.SampleRate/waveform.PRF);
target = phased.RadarTarget('Model','Nonfluctuating',...
    'MeanRCS',1,'OperatingFrequency',fc);
targetchannel = phased.FreeSpace('TwoWayPropagation',true,...
    'SampleRate',Fs,'OperatingFrequency', fc);
jammerchannel = phased.FreeSpace('TwoWayPropagation',false,...
    'SampleRate',Fs,'OperatingFrequency', fc);
collector = phased.Collector('Sensor',antenna,...
    'OperatingFrequency',fc);
amplifier = phased.ReceiverPreamp('EnableInputPort',true);

Примите, что массив, цель и передатчик помех являются стационарными. Массив расположен в глобальном начале координат, (0,0,0). Цель расположена в (1000,500,0), и передатчик помех расположен в (2000 2000 100). Определите направления от массива до цели и передатчика помех.

targetloc = [1000 ; 500; 0];
jammerloc = [2000; 2000; 100];
[~,tgtang] = rangeangle(targetloc);
[~,jamang] = rangeangle(jammerloc);

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

wav = waveform();
% Transmit waveform
[wav,txstatus] = transmitter(wav);
% Radiate pulse toward the target
wav = radiator(wav,tgtang);
% Propagate pulse toward the target
wav = targetchannel(wav,[0;0;0],targetloc,[0;0;0],[0;0;0]);
% Reflect it off the target
wav = target(wav);
% Collect the echo
wav = collector(wav,tgtang);

Сгенерируйте создающий затор сигнал

jamsig = jammer();
% Propagate the jamming signal to the array
jamsig = jammerchannel(jamsig,jammerloc,[0;0;0],[0;0;0],[0;0;0]);
% Collect the jamming signal
jamsig = collector(jamsig,jamang);

% Receive target echo alone and target echo + jamming signal
pulsewave = amplifier(wav,~txstatus);
pulsewave_jamsig = amplifier(wav + jamsig,~txstatus);

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

subplot(2,1,1)
t = unigrid(0,1/Fs,size(pulsewave,1)*1/Fs,'[)');
plot(t*1000,abs(pulsewave(:,1)))
title('Magnitudes of Pulse Waveform Without Jamming--Element 1')
ylabel('Magnitude')
subplot(2,1,2)
plot(t*1000,abs(pulsewave_jamsig(:,1)))
title('Magnitudes of Pulse Waveform with Jamming--Element 1')
xlabel('millisec')
ylabel('Magnitude')

Figure contains 2 axes objects. Axes object 1 with title Magnitudes of Pulse Waveform Without Jamming--Element 1 contains an object of type line. Axes object 2 with title Magnitudes of Pulse Waveform with Jamming--Element 1 contains an object of type line.