matlab开发-数字脉冲压缩MATLAB的实现

matlab开发-数字脉冲压缩MATLAB的实现。本项目的主要目标是设计脉冲压缩滤波器（匹配滤波器）。

tic
clc;
clf;
BW=5*10^6;%bandwidth
T=51.2*10^-6;%pulse width
fc=0; %zero centered frequency
n=fix（2*BW*T+1）;% sampling points
sampling_interval = T / n;%sampling interval
freqlimit = 0.5/ sampling_interval;
alpha=BW/（T）;%chirp rate
t=linspace（-T/2T/2n）;%time interval
ichannel1=cos（2*pi*（（fc）*t+alpha*t.^2））;%real part
qchannel2=sin（2*pi*（（fc）*t+alpha*t.^2））;%imaginary part
LFMD=0.0892*（ichannel1+（1i.*qchannel2））;%Complex base band time domain signal
LFMD（2:1:10）=0;
LFMD（504:1:511）=0;
freq = linspace（-freqlimit/2freqlimit/2n）;
FLFMD=fftshift（fft（LFMDn））;%Frequency domain of base band signal
[lfmfreiflfm]=lffm1（BWnT）;%Function calling for input Lfm signal 
figure（1）
subplot（311）
plot（freiabs（flfm））;
plot（freqabs（（FLFMD）））;
xlabel（‘Frequency domain‘）
ylabel（‘Amplitude‘）
title（‘Input signal of pulse compression filter in Frequency domain‘）;
grid;
%plot（t（（LFMD））‘k‘）
% plot（t（lfm）‘k‘）;
rep=fliplr（LFMD）;
%  disp（length（rep））;
j=0;
y1=zeros（110001）;%10Khzs Sampling domain   
for i=2246:1:2758
    j=j+1;
    y1（i）=LFMD（j）;
end
%     j=0;
%      for i=3976:1:4488
%     j=j+1;
%     y1（i）=LFMD（j）;
%      end
 
%    for i=9246:1:9758
%     j=j+1;
%     y1（i）=LFMD（j）;
%   end

Q=length（lfm）;%length of lfm signal h（k）
P=length（y1）;%length of base band signal domain X（k） 
N=P+Q-1;
N=2.^nextpow2（N）;% FFTbuffer
x1=zeros（1N）;%creating a  empty 16k Buffer
% creating stored replica buffer of 16k
for k1=1:1:N
    if k1<=Q
        h1（k1）=rep（k1）;
    else
        h1（k1）=0;
    end
    
end
hn=fftshift（fft（h1N））; %taking FFt of 16k H（k）buffer
hnc=（conj（hn））;%
%creating the 16k buffer for X（k）signal
n3=0;
for k2=1:1:N
     if k2<=N-P+1
         x1（k2）=0;
     else
         n3=n3+1;
        x1（k2）=y1（n3）;
     end
end

  x1=awgn（x10.1）;
 subplot（211）;
 tt=linspace（01*10^-6N）;
 plot（ttx1）;
 xlabel（‘Time domain‘）
 ylabel（‘Amplitude‘）
title（‘Noise along with the signal ‘）;
%  grid;
 xf=（fftshift（fft（x1N）））; %fft of the 16k X（k） buffer
%computing the product of X（k）*h（k）
zn=xf.*hnc;
%IFFT of the Multiplier output buffer
zk=（ifft（znN））;
zk=abs（（zk））;
zk1=zk（1:1:8192）;
var1=var（zk1）;
zk2=zk（9590:1:16000）;
var2=var（zk2）;
n12=var1-var2;
t4=linspace（0NN）;
subplot（212）
plot（t4（zk））;
xlabel（‘Time domain‘）
ylabel（‘Amplitude‘）
title（‘ Mainlobe width analysis for the pulse compression filter‘）;
grid;
toc;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        2478  2010-05-18 07:24  dd1.asv
     文件        2441  2010-05-14 00:52  dd1.m
     文件        2319  2010-05-13 23:57  direct.m
     文件         620  2010-05-13 22:18  lffm1.m
     文件       42555  2010-06-20 16:51  llll.mdl
     文件        1883  2010-06-22 00:39  matched2.m
     文件        1890  2010-05-14 01:15  matched21.m
     文件        4741  2010-04-05 15:24  ovsav1.m
     文件        2500  2010-04-13 13:02  ovsav2.m
     文件        2050  2010-04-05 15:24  ovsav3.m
     文件        1836  2010-04-05 14:49  ovsav4.m
     文件        2150  2010-05-14 01:12  ovsav5 （2）.m
     文件        2033  2010-05-19 06:58  ovsav5.asv
     文件        2150  2010-05-14 01:12  ovsav5.m
     文件        2124  2010-05-14 01:16  ovsav51 （2）.m
     文件        2124  2010-05-14 01:16  ovsav51.m
     文件       22546  2010-06-20 17:27  sineee.mdl
     文件        1306  2014-02-12 13:15  license.txt