线性调频信号产生及压缩.m

雷达系统中经常学习的东西！！！！ 也是初学者最基本的东西

disp('脉冲压缩仿真程序')
M=[];A=[];B=[];B1=[];W=[];B2=[];W1=[];
b=input('input 调频带宽 b=')
T=input('input 脉冲宽度 T=')
fc=input('input 采样频率 fc=')
Ts=T*fc;
if Ts<=50
      w=128;
elseif Ts<=100
      w=256;
elseif Ts<=200
      w=512;
elseif Ts<=500
      w=1024;
elseif Ts<=1000
      w=2048; 
elseif Ts<=2000
       w=4096;
end
K=-b/T;
N=T*fc;
for n=-(N-1)/2:1:(N-1)/2
       if n<(N-1)/2                                                %??????????????????
Sn=10*(cos(K*pi*n*n/fc/fc)+j*sin(K*pi*n*n/fc/fc)+(randn+j*randn));%Sn1=10*(cos(K*pi*n1^2/fc^2+fd*n1/fc)+j*sin(K*pi*n1^2/fc^2+fd*n1/fc))+(randn+j*randn);
hn=cos(K*pi*n*n/fc/fc)-j*sin(K*pi*n*n/fc/fc);
Wn=0.54+0.46*cos(2.0*pi*n*K/b/fc);
Wb=0.42+0.5*cos(2.0*pi*n*K/b/fc)+0.08*cos(4.0*pi*n*K/b/fc);
Hn=hn*Wn;
Qn=hn*Wb;
      else
Sn=(randn+j*randn);
Wn=0;
hn=0;
Hn=0;
Wb=0;
Qn=0;
     end
W=[W,Wn];
A=[A,Sn];
B=[B,Hn];
B1=[B1,Qn];
B2=[B2,hn];
W1=[W1,Wb];
end
C=fft(A,w);

D=[fft(B,w)];
D1=[fft(B1,w)];
D2=[fft(B2,w)];
E=D2.*C;
E0=D.*C;
E1=D1.*C;
E2=[ifft(E,w)];
E3=[ifft(E0,w)];
E4=[ifft(E1,w)];
E5=abs(E2);
E6=max(E5);
E7=20*log10(E5/E6);    %Divide by zero
E8=abs(E3);
E9=max(E8);
E10=20*log10(E8/E9);   %Divide by zero
E11=abs(E4);
E12=max(E11);
E13=20*log10(E11/E12); %Divide by zero
subplot(3,2,1);
plot(E7)
grid
title('T微秒信号压缩结果')
xlabel('样本个数N')
Ylabel('主副瓣比db')
subplot(3,2,2);
plot(E5)
grid
title('T微秒信号压缩结果')
xlabel('样本个数N')
Ylabel('加权函数值')
subplot(3,2,3);
plot(E10)
grid
title('T微秒信号压缩结果(海明窗)')
xlabel('样本个数N')
Ylabel('主副瓣比db')
subplot(3,2,4);
plot(E8)
grid
title('T微秒信号压缩结果(海明窗)')
xlabel('样本个数N')
Ylabel('主副瓣比db')
subplot(3,2,5);
plot(E13)
grid
title('T微秒信号压缩结果(布莱克曼窗)')
xlabel('样本个数N')
Ylabel('主副瓣比db')
subplot(3,2,6);
plot(E11)
grid
title('T微秒信号压缩结果(布莱克曼窗)')
xlabel('样本个数N')