barrageJammer
объектные модели широкополосный передатчик помех. Выход barrageJammer
комплексная последовательность белого Гауссова шума. Модифицируемые свойства заградительного передатчика помех:
ERP
— Эффективная излучаемая мощность в ваттах
SamplesPerFrameSource
— Источник количества отсчетов на систему координат
SamplesPerFrame
— Количество отсчетов на систему координат
SeedSource
— Источник seed для генератора случайных чисел
Seed
— Отберите для генератора случайных чисел
Действительные и мнимые части комплексной последовательности белого Гауссова шума у каждого есть отклонение, равное 1/2 эффективная излучаемая мощность в ваттах. Обозначьте эффективную излучаемую мощность в ваттах P. Заградительный передатчик помех выход:
В этом уравнении, 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')
Средние значения действительных и мнимых частей
mean(real(y))
ans = -1.0961
mean(imag(y))
ans = -2.1671
которые являются эффективно нулевыми. Стандартные отклонения действительных и мнимых частей
std(real(y))
ans = 50.1950
std(imag(y))
ans = 49.7448
которые соглашаются с ожидаемым значением .
Этот пример демонстрирует, как симулировать эффект заградительного передатчика помех на целевом эхе. Во-первых, создайте требуемые объекты. Вам нужны массив, передатчик, излучатель, цель, передатчик помех, коллектор и приемник. Кроме того, необходимо задать два пути к распространению: один от массива до цели и назад, и другой путь от передатчика помех до массива.
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')