lzchirpshipei11.m

利用MATLAB作的完整实例

subplot(2,3,1);
plot(t,s1);                                      % 时域信号图
title('f 0？=0时Chirp信号');
xlabel('time[s]');
ylabel('scope');
% 信号FFT变换后分析频域特性：
S1=fft(s1);                                      % 傅立叶变换
angS1=angle(S1);                                 % 相频特性
S11=fftshift(S1);                                % 移动零点到频谱中心 
magS1=abs(S11/max(S11));                         % 幅频特性
k=(-length(S1)/2:(length(S1)-1)/2)/length(S1);   % 单位长度点数
f=fs*k;                                          % 频率
subplot(2,3,2);
plot(f,magS1);                                   % 幅频图
axis([-B B 0 1.2]);
title('幅频谱');
xlabel('frequency[Hz]');
ylabel('scope');
subplot(2,3,3);
plot(f,unwrap(angS1));                           % 相频图
grid;
title('相频谱');
xlabel('frequency[Hz]');
ylabel('angle[drgree]');
% Chirp信号通过匹配滤波器，实现脉冲压缩
R1=S11.*H1;                                      % 通过滤波器的信号
magR1=abs(R1/max(R1));                           % 幅频特性
angR1=angle(R1);                                 % 相频特性
r1=ifft(R1);                                     % 变换到时域                   
r1=ifftshift(r1);
r11=abs(r1)/max(abs(r1));                        % 归一化
r111=20*log10(r11);                              % 用对数表示
subplot(2,3,5);
plot(f,magR1);
axis([-B B 0 1.2]);
title('通过失配滤波器幅频谱');
xlabel('frequency[Hz]');
ylabel('scope');
% 用对数形式表示，并插值
% 原理：频域补零相当于时域插值
W1=S1.*H;  
W1=fftshift(W1);
[intp_sig,intp_logsig]=intp(W1);
w11=intp_sig;
w111=intp_logsig;
subplot(2,3,4);
plot(t_add0,w11); 
grid;
axis([-Tao/10-Tao/10 Tao/10 0 1]);
title('信号通过失配滤波器');
xlabel('time[s]');
ylabel('scope');
subplot(2,3,6)
plot(t_add0,w111);
grid;
axis([-Tao/10-Tao/10 Tao/10 -30 0]);
title('信号通过失配滤波器对数图');
xlabel('time[s]');
ylabel('scope[dB]');
% 加评估函数                                   
[O,width3dB,ratio,peak1,peak_1,peak2,peak_2,peak3,peak_3]=para(w11,t_add0);   % 评估参数
O                   
width3dB
ratio
peak1
peak_1
peak2
peak_2
peak3
peak_3
%-------------------------------------------------------------

figure(7);
% 7.2
% u1?=u时Chirp复信号的建模：
s2=exp(j*pi*(B+1.e+7)*(t.^2)/(Tao));            % 线性调频信号建摸
subplot(2,3,1);
plot(t,s2);                                     % 时域信号图
title('u 1？=u时Chirp信号');
xlabel('time[s]');
ylabel('scope');
% 信号FFT变换后分析频域特性：
S2=fft(s2);                                      % 傅立叶变换
angS2=angle(S2);                                 % 相频特性
S22=fftshift(S2);                                % 移动零点到频谱中心 
magS2=abs(S22/max(S22));                         % 幅频特性
k=(-length(S2)/2:(length(S2)-1)/2)/length(S2);   % 单位长度点数
f=fs*k;                                          % 频率
subplot(2,3,2);
plot(f,magS2);                                   % 幅频图
axis([-B B 0 1.2]);
title('幅频谱');
xlabel('frequency[Hz]');
ylabel('scope');
subplot(2,3,3);
plot(f,unwrap(angS2));                           % 相频图
grid;
title('相频谱');
xlabel('frequency[Hz]');
ylabel('angle[drgree]');
% Chirp信号通过匹配滤波器，实现脉冲压缩
R2=S22.*H1;                                      % 通过滤波器的信号
magR2=abs(R2/max(R2));                           % 幅频特性
angR2=angle(R2);                                 % 相频特性
r2=ifft(R2);                                     % 变换到时域                   
r2=ifftshift(r2);
r21=abs(r2)/max(abs(r2));                        % 归一化
r211=20*log10(r21);                              % 用对数表示
subplot(2,3,5);
plot(f,magR2);
axis([-B B 0 1.2]);
title('通过失配滤波器幅频谱');
xlabel('frequency[Hz]');
ylabel('scope');    
% 用对数形式表示，并插值
% 原理：频域补零相当于时域插值
W2=S2.*H;  
W2=fftshift(W2);
[intp_sig,intp_logsig]=intp(W2);
w21=intp_sig;
w211=intp_logsig;
subplot(2,3,4);
plot(t_add0,w21); 
grid;
axis([-Tao/4 Tao/4 0 1]);
title('信号通过失配滤波器');
xlabel('time[s]');
ylabel('scope');
subplot(2,3,6)
plot(t_add0,w211);
grid;
axis([-Tao/4 Tao/4 -30 0]);
title('信号通过失配滤波器对数图');
xlabel('time[s]');
ylabel('scope[dB]');
% 加评估函数                                   
[O,width3dB,ratio,peak1,peak_1,peak2,peak_2,peak3,peak_3]=para(w21,t_add0);   % 评估参数
O                   
width3dB
ratio
peak1
peak_1
peak2
peak_2
peak3
peak_3
%-------------------------------------------------------------

figure(8);
% 7.3
% r(t)?=0时Chirp复信号的建模：
s3=exp(j*pi*B*(t.^2)/Tao+j*1.e+19*t.^3);        % 线性调频信号建摸
subplot(2,3,1);
plot(t,s3);                                     % 时域信号图
title('r(t)？=0时Chirp信号');
xlabel('time[s]');
ylabel('scope');
% 信号FFT变换后分析频域特性：
S3=fft(s3);                                      % 傅立叶变换
angS3=angle(S3);                                 % 相频特性
S33=fftshift(S3);                                % 移动零点到频谱中心 
magS3=abs(S33/max(S33));                         % 幅频特性
k=(-length(S3)/2:(length(S3)-1)/2)/length(S3);   % 单位长度点数
f=fs*k;                                          % 频率
subplot(2,3,2);
plot(f,magS3);                                   % 幅频图
axis([-B B 0 1.2]);
title('幅频谱');
xlabel('frequency[Hz]');
ylabel('scope');
subplot(2,3,3); 
plot(f,unwrap(angS3));                           % 相频图
grid;
title('相频谱');
xlabel('frequency[Hz]');
ylabel('angle[drgree]');
% Chirp信号通过匹配滤波器，实现脉冲压缩
R3=S33.*H1;                                      % 通过滤波器的信号
magR3=abs(R3/max(R3));                           % 幅频特性
angR3=angle(R3);                                 % 相频特性
r3=ifft(R3);                                     % 变换到时域                   
r3=ifftshift(r3);
r31=abs(r3)/max(abs(r3));                        % 归一化
r311=20*log10(r31);                              % 用对数表示
subplot(2,3,5);
plot(f,magR3);
axis([-B B 0 1.2]);
title('通过失配滤波器幅频谱');
xlabel('frequency[Hz]');
ylabel('scope');
% 用对数形式表示，并插值
% 原理：频域补零相当于时域插值
W3=S3.*H;  
W3=fftshift(W3);
[intp_sig,intp_logsig]=intp(W3);
w31=intp_sig;
w311=intp_logsig;
subplot(2,3,4);
plot(t_add0,w31); 
grid;
axis([-Tao/10 Tao/10 0 1]);
title('信号通过失配滤波器');
xlabel('time[s]');
ylabel('scope');
subplot(2,3,6)
plot(t_add0,w311);
grid;
axis([-Tao/10 Tao/10 -30 0]);
title('信号通过失配滤波器对数图');
xlabel('time[s]');
ylabel('scope[dB]');
% 加评估函数                                   
[O,width3dB,ratio,peak1,peak_1,peak2,peak_2,peak3,peak_3]=para(w31,t_add0);   % 评估参数
O                   
width3dB
ratio
peak1
peak_1
peak2
peak_2
peak3
peak_3
%-------------------------------------------------------------


figure(9);
% 7.3
% 高频干扰Chirp复信号的建模：
s4=exp(j*pi*B*(t.^2)/Tao+j*sin(2.e+10*t));        % 线性调频信号建摸
subplot(2,3,1);
plot(t,s4);                                     % 时域信号图
title('高频干扰时Chirp信号');
xlabel('time[s]');
ylabel('scope');
% 信号FFT变换后分析频域特性：
S4=fft(s4);                                      % 傅立叶变换
angS4=angle(S4);                                 % 相频特性
S43=fftshift(S4);                                % 移动零点到频谱中心 
magS4=abs(S43/max(S43));                         % 幅频特性
k=(-length(S4)/2:(length(S4)-1)/2)/length(S4);   % 单位长度点数
f=fs*k;                                          % 频率
subplot(2,3,2);
plot(f,magS4);                                   % 幅频图
axis([-B B 0 1.2]);
title('幅频谱');
xlabel('frequency[Hz]');
ylabel('scope');
subplot(2,3,3); 
plot(f,unwrap(angS4));                           % 相频图
grid;
title('相频谱');
xlabel('frequency[Hz]');
ylabel('angle[drgree]');
% Chirp信号通过匹配滤波器，实现脉冲压缩
R4=S43.*H1;                                      % 通过滤波器的信号
magR4=abs(R4/max(R4));                           % 幅频特性
angR4=angle(R4);                                 % 相频特性
r4=ifft(R4);                                     % 变换到时域                   
r4=ifftshift(r4);
r41=abs(r4)/max(abs(r4));                        % 归一化
r411=20*log10(r41);                              % 用对数表示
subplot(2,3,5);
plot(f,magR4);
axis([-B B 0 1.2]);
title('通过失配滤波器幅频谱');
xlabel('frequency[Hz]');
ylabel('scope');
% 用对数形式表示，并插值
% 原理：频域补零相当于时域插值
W4=S4.*H;  
W4=fftshift(W4);
[intp_sig,intp_logsig]=intp(W4);
w41=intp_sig;
w411=intp_logsig;
subplot(2,3,4);
plot(t_add0,w41); 
grid;
%axis([-Tao/10 Tao/10 0 1]);
title('信号通过失配滤波器');
xlabel('time[s]');
ylabel('scope');
subplot(2,3,6)
plot(t_add0,w411);
grid;
%axis([-Tao/10 Tao/10 -30 0]);
title('信号通过失配滤波器对数图');
xlabel('time[s]');
ylabel('scope[dB]');
% 加评估函数                                   
[O,width3dB,ratio,peak1,peak_1,peak2,peak_2,peak3,peak_3]=para(w41,t_add0);   % 评估参数
O                   
width3dB
ratio
peak1
peak_1
peak2
peak_2
peak3
peak_3
%-------------------------------------------------------------