Создайте библиотеку спецификаций импульсного сжатия
The pulseCompressionLibrary System 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 формой волны. В согласованных фильтрах используются окно Ханна и окно Тейлора, соответственно.Чтобы использовать функцию объекта, задайте системный объект в качестве первого входного параметра. Например, чтобы освободить системные ресурсы системного объекта с именем 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);
Создайте библиотеку обработки сжатия с помощью p ulseCompressionLibrary Системный 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, поражающего три цели на областях значений 2000, 4000 и 5500 метров. Принимая максимальную дальность действия радара 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');
