Создайте библиотеку импульсных технических требований сжатия
pulseCompressionLibrary Система object™ создает импульсную библиотеку сжатия. Библиотека содержит наборы параметров, которые описывают импульсные операции сжатия, выполняемые на полученных сигналах сгенерировать их ответ области значений. Можно пользоваться этой библиотекой, чтобы выполнить обработка фрагмента или согласованная фильтрация. Этот объект может обработать формы волны, созданные pulseWaveformLibrary объект.
Сделать импульсную библиотеку сжатия
Создайте pulseCompressionLibrary объект и набор его свойства.
Вызовите объект с аргументами, как будто это была функция.
Чтобы узнать больше, как Системные объекты работают, смотрите то, Что Системные объекты?
complib = pulseCompressionLibrary()complib, со значениями свойств по умолчанию.
complib = pulseCompressionLibrary(Name,Value)Name установите на заданный Value. Можно задать дополнительные аргументы пары "имя-значение" в любом порядке как (Name1, Value1..., NameN, ValueN). Заключите каждое имя свойства в одинарные кавычки.
complib = pulseCompressionLibrary('SampleRate',1e9,'WaveformSpecification',{{'Rectangular','PRF',1e4,'PulseWidth',100e-6},{'SteppedFM','PRF',1e4}},'ProcessingSpecification',{{'MatchedFilter','SpectrumWindow','Hann'},{'MatchedFilter','SpectrumWindow','Taylor'}}) создает библиотеку с двумя согласованными фильтрами. Каждый является соответствующим к прямоугольному радиоимпульсу и другому к ступенчатой форме волны FM. Согласованные фильтры используют окно Hann и окно Тейлора, соответственно.Чтобы использовать объектную функцию, задайте Системный объект как первый входной параметр. Например, чтобы выпустить системные ресурсы Системного объекта под названием obj, используйте этот синтаксис:
release(obj)
Создайте прямоугольный радиоимпульс и линейную форму волны FM. Используйте методы обработки в импульсной библиотеке сжатия к процессу области значений формы волны. Используйте согласованную фильтрацию для прямоугольного радиоимпульса и обработку фрагмента для линейной формы волны FM.
Создайте две формы волны с помощью pulseWaveformLibrary Система object™. Частота дискретизации составляет 1 МГц, и импульсная частота повторения для обеих форм волны составляет 1 кГц. Ширина импульса является также тем же самым в 50 микросекунд.
fs = 1.0e6;
prf = 1e3;
pw = 50e-6;
waveform1 = {'Rectangular','PRF',prf,'PulseWidth',pw};
waveform2 = {'LinearFM','PRF',prf,'PulseWidth',pw,...
'SweepBandwidth',1e5,'SweepDirection','Up',...
'SweepInterval', 'Positive'};
pulselib = pulseWaveformLibrary('WaveformSpecification',...
{waveform1,waveform2},'SampleRate',fs);Получите формы волны для обработки импульсной библиотекой сжатия.
rectwav = pulselib(1); lfmwav = pulselib(2);
Создайте библиотеку обработки сжатия с помощью pulseCompressionLibrary Система object™ с двумя техническими требованиями обработки. Первой спецификацией обработки является согласованная фильтрация, и второй является обработка фрагмента.
mf = getMatchedFilter(pulselib,1);
procspec1 = {'MatchedFilter','Coefficients',mf};
procspec2 = {'StretchProcessor','ReferenceRange',5000,...
'RangeSpan',200,'RangeWindow','Hamming'};
comprlib = pulseCompressionLibrary( ...,
'WaveformSpecification',{waveform1,waveform2}, ...
'ProcessingSpecification',{procspec1,procspec2}, ...
'SampleRate',fs,'PropagationSpeed',physconst('Lightspeed'));Процесс обе формы волны.
rect_out = comprlib(rectwav,1); lfm_out = comprlib(lfmwav,2); nsamp = fs/prf; t = [0:(nsamp-1)]/fs; plot(t*1000,real(rect_out)) hold on plot(t*1000,real(lfm_out)) hold off title('Pulse Compression Output') xlabel('Time (millsec)') ylabel('Amplitude')
Постройте ответ области значений сигнала LFM достижение трех целевых показателей в областях значений 2 000, 4000, и 5 500 метров. Принятие максимальной области значений радара составляет 10 км, определите импульсный интервал повторения из максимальной области значений.
% Create the pulse waveform. rmax = 10.0e3; c = physconst('Lightspeed'); pri = 2*rmax/c; fs = 1e6; pri = ceil(pri*fs)/fs; prf = 1/pri; nsamp = pri*fs; rxdata = zeros(nsamp,1); t1 = 2*2000/c; t2 = 2*4000/c; t3 = 2*5500/c; idx1 = floor(t1*fs); idx2 = floor(t2*fs); idx3 = floor(t3*fs); lfm = phased.LinearFMWaveform('PulseWidth',10/fs,'PRF',prf, ... 'SweepBandwidth',(30*fs)/40); w = lfm(); %% % Imbed the waveform part of the pulse into the received signal. x = w(1:11); rxdata(idx1:idx1+10) = x; rxdata(idx2:idx2+10) = x; rxdata(idx3:idx3+10) = x; %% % Create the pulse waveform library. w1 = {'LinearFM','PulseWidth',10/fs,'PRF',prf,... 'SweepBandwidth',(30*fs)/40}; wavlib = pulseWaveformLibrary('SampleRate',fs,'WaveformSpecification',{w1}); wav = wavlib(1); %% % Generate the range response signal. p1 = {'MatchedFilter','Coefficients',getMatchedFilter(wavlib,1),'SpectrumWindow','None'}; idx = 1; complib = pulseCompressionLibrary( ... 'WaveformSpecification',{w1}, ... 'ProcessingSpecification',{p1}, ... 'SampleRate',fs, ... 'PropagationSpeed',c); y = complib(rxdata,1); %% % Plot range response of processed data plotResponse(complib,rxdata,idx,'Unit','mag');
