Создание библиотеки спецификаций импульсного сжатия
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'}}) создает библиотеку с двумя подходящими фильтрами. Один соответствует прямоугольной форме сигнала, а другой - ступенчатой ЧМ форме сигнала. Соответствующие фильтры используют окно Ханна и окно Тейлора соответственно.Чтобы использовать функцию объекта, укажите объект System в качестве первого входного аргумента. Например, для освобождения системных ресурсов объекта System с именем obj, используйте следующий синтаксис:
release(obj)
Создайте прямоугольный и линейный ЧМ-сигналы. Используйте методы обработки в библиотеке сжатия импульсов для обработки диапазона сигналов. Используйте согласованную фильтрацию для прямоугольной формы сигнала и обработку растяжения для линейной ЧМ формы сигнала.
Создайте две формы сигнала с помощью 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, поражающего три цели на дальностях 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');
