信号类型（雷达）——脉冲雷达（四）

已于 2023-10-31 11:00:43 修改
阅读量990 收藏 21
点赞数 3
分类专栏： 信号处理 # 雷达 代码实现 文章标签： 信号处理 学习 matlab
于 2023-09-30 22:54:31 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/yjh_2019/article/details/133442895
版权
信号处理 同时被 3 个专栏收录
41 篇文章 155 订阅
订阅专栏
代码实现
33 篇文章 26 订阅
订阅专栏
雷达
11 篇文章 19 订阅
订阅专栏

系列文章目录

《信号类型（雷达+通信）》

《信号类型（雷达）——雷达波形认识（一）》

《信号类型（雷达）——连续波雷达（二）》

《信号类型（雷达）——脉冲雷达（三）》

       为实现各种辐射源信号的分选识别以及定位，首先需要对辐射源信号特性有一个充分了解，根据各类辐射源特点从众多信号冲分离出来，甚至去除辐射源调制信息的影响，实现信号长时间相干定位处理，提高定位精度。本文主要记录目前了解到的雷达脉冲信号时域、频域、时频域以及模糊函数特性，作为后续辐射源分选识别以及定位研究的笔记，当然，若对相关行业的人有所启发，也欢迎能在评论区提出建议，或自身研究心得。

一、调频类脉冲雷达

1.1、常规脉冲雷达

       矩形脉冲包络，没有相位调制信息。下图分别展示了这类脉冲波形的时域、频谱、时频图以及模糊函数特性。

​​

时域实部和虚部                                   时域模值和相位 

​​

频谱图                                                                时频图

​​

模糊函数                                                模糊函数（俯视图）

​​

模糊函数（0频）                                       模糊函数（0时延）

clear;close all;clc;

%% 
% 矩形脉冲波形（非调制）

%------------------------------波形分析------------------------------------
%% 波形参数
waveform = phased.RectangularWaveform;
waveform.SampleRate = 5e6;                             % 采样率
waveform.DurationSpecification = 'Pulse width';           % 脉冲持续时间设置方式，'Pulse width' or 'Duty cycle'
waveform.PulseWidth = 20e-6;                            % 脉宽
waveform.PRF = 5000;                                      % PRF间隔
xt = step(waveform);

%% 信号实部/虚部
figure;
plot(real(xt));
hold on;
plot(imag(xt));
title('信号实部/虚部');
legend('实部','虚部');

%% 信号模值/相位
figure;
plot(abs(xt));
hold on;
plot(angle(xt));
title('信号模值/相位');
legend('模值','相位');

%% 信号频谱
f=linspace(-waveform.SampleRate/2,waveform.SampleRate/2,length(xt));
figure;
plot(f,abs(fftshift(fft(xt))));
xlabel('采样率（Hz）');
title('信号频谱');

%% 信号时频图
window = blackman(128);
noverlap = 32;
nfft = 512;
[s,f,t] = spectrogram(xt,window,noverlap,nfft,waveform.SampleRate,'centered'); 
figure;
imagesc(t,f,abs(s));
xlabel('时间（s）');
ylabel('频率（Hz）');
title('信号时频图');

%% 模糊函数
[afmag,delay,doppler] = ambgfun(xt,waveform.SampleRate,waveform.PRF); 
[ambrect1,delayrect] = ambgfun(xt,waveform.SampleRate,...,
    waveform.PRF,'Cut','Doppler','CutValue',0);
[ambrect2,dopplerrect] = ambgfun(xt,waveform.SampleRate,...,
    waveform.PRF,'Cut','Delay','CutValue',0);

figure;
mesh(delay,doppler,abs(afmag));
xlabel('时间（s）');
ylabel('频率（Hz）');
title('信号模糊函数');

figure;
plot(delayrect*10^6,ambrect1);
xlabel('时延 (\mu sec)');
title('0Hz多普勒频点处的自相关函数');

figure;
plot(dopplerrect/10^6,ambrect2);
xlabel('多普勒频率 (MHz)');
title('0s时延处的自相关函数');

1.2、线性调频脉冲雷达

       矩形脉冲包络，相位存在线性调频调制。下图分别展示了这类脉冲波形的时域、频谱、时频图以及模糊函数特性。

 ​​

时域实部和虚部                                   时域模值和相位 

​

频谱图                                                                时频图

​​

模糊函数                                                模糊函数（俯视图）

​​

模糊函数（0频）                                       模糊函数（0时延）

%------------------------------波形分析------------------------------------
%% 波形参数
waveform = phased.LinearFMWaveform;
waveform.SampleRate = 5e6;                             % 采样率
waveform.DurationSpecification = 'Pulse width';           % 脉冲持续时间设置方式，'Pulse width' or 'Duty cycle'
waveform.PulseWidth = 20e-6;                            % 脉宽
waveform.PRF = 5000;                                      % PRF间隔
waveform.SweepBandwidth = 1e6;                            % 频扫带宽
waveform.SweepDirection = 'Up';                           % 'Up' or 'Down'
waveform.SweepInterval = 'Symmetric';                     % 'Positive' or 'Symmetric'
% waveform.Envelope = 'Gaussian';                           % 'Rectangular' or 'Gaussian'.
xt = step(waveform);

%% 信号实部/虚部
figure;
plot(real(xt));
hold on;
plot(imag(xt));
title('信号实部/虚部');
legend('实部','虚部');

%% 信号模值/相位
figure;
plot(abs(xt));
hold on;
plot(angle(xt));
title('信号模值/相位');
legend('模值','相位');

%% 信号频谱
f=linspace(-waveform.SampleRate/2,waveform.SampleRate/2,length(xt));
figure;
plot(f,abs(fftshift(fft(xt))));
xlabel('采样率（Hz）');
title('信号频谱');

%% 信号时频图
window = blackman(16);
noverlap = 10;
nfft = 512;
[s,f,t] = spectrogram(xt,window,noverlap,nfft,waveform.SampleRate,'centered'); 
figure;
imagesc(t,f,abs(s));
xlabel('时间（s）');
ylabel('频率（Hz）');
title('信号时频图');

%% 模糊函数
[afmag,delay,doppler] = ambgfun(xt,waveform.SampleRate,waveform.PRF); 
[ambLFM1,delayLFM] = ambgfun(xt,waveform.SampleRate,...,
    waveform.PRF,'Cut','Doppler','CutValue',0);
[ambLFM2,dopplerLFM] = ambgfun(xt,waveform.SampleRate,...,
    waveform.PRF,'Cut','Delay','CutValue',0);

figure;
mesh(delay,doppler,abs(afmag));
xlabel('时间（s）');
ylabel('频率（Hz）');
title('信号模糊函数');

figure;
plot(delayLFM*10^6,ambLFM1);
xlabel('时延 (\mu sec)');
title('0Hz多普勒频点处的自相关函数');

figure;
plot(dopplerLFM/10^6,ambLFM2);
xlabel('多普勒频率 (MHz)');
title('0s时延处的自相关函数');

1.3、步进频雷达

       矩形脉冲包络，脉间频率逐渐递增。下图分别展示了这类脉冲波形的时域、频谱、时频图以及模糊函数特性。

 ​​

时域实部和虚部                                   时域模值和相位 

​

频谱图                                                                时频图

​​

模糊函数                                                模糊函数（俯视图）

​​

模糊函数（0频）                                       模糊函数（0时延）

%------------------------------波形分析------------------------------------
%% 波形参数
waveform = phased.SteppedFMWaveform;
waveform.SampleRate = 5e6;                             % 采样率
waveform.DurationSpecification = 'Pulse width';           % 脉冲持续时间设置方式，'Pulse width' or 'Duty cycle'
waveform.PulseWidth = 20e-6;                            % 脉宽
waveform.PRF = 5000;                                      % PRF间隔
waveform.FrequencyStep = 500e3;                           % 步进间隔
waveform.NumSteps = 3;                                    % 步进数量
waveform.OutputFormat = 'Pulses';                         % 输出信号的格式为'Pulses' or 'Samples'
waveform.NumPulses = 3;
xt = step(waveform);

%% 信号实部/虚部
figure;
plot(real(xt));
hold on;
plot(imag(xt));
title('信号实部/虚部');
legend('实部','虚部');

%% 信号模值/相位
figure;
plot(abs(xt));
hold on;
plot(angle(xt));
title('信号模值/相位');
legend('模值','相位');

%% 信号频谱
f=linspace(-waveform.SampleRate/2,waveform.SampleRate/2,length(xt));
figure;
plot(f,abs(fftshift(fft(xt))));
xlabel('采样率（Hz）');
title('信号频谱');

%% 信号时频图
window = blackman(16);
noverlap = 10;
nfft = 512;
[s,f,t] = spectrogram(xt,window,noverlap,nfft,waveform.SampleRate,'centered'); 
figure;
imagesc(t,f,abs(s));
xlabel('时间（s）');
ylabel('频率（Hz）');
title('信号时频图');

%% 模糊函数
[afmag,delay,doppler] = ambgfun(xt,waveform.SampleRate,waveform.PRF); 
[ambsfm1,delaysfm] = ambgfun(xt,waveform.SampleRate,...,
    waveform.PRF,'Cut','Doppler','CutValue',0);
[ambsfm2,dopplersfm] = ambgfun(xt,waveform.SampleRate,...,
    waveform.PRF,'Cut','Delay','CutValue',0);

figure;
mesh(delay,doppler,abs(afmag));
xlabel('时间（s）');
ylabel('频率（Hz）');
title('信号模糊函数');

figure;
plot(delaysfm*10^6,ambsfm1);
xlabel('时延 (\mu sec)');
title('0Hz多普勒频点处的自相关函数');

figure;
plot(dopplersfm/10^6,ambsfm2);
xlabel('多普勒频率 (MHz)');
title('0s时延处的自相关函数');

二、相位编码雷达

2.1、Barker编码

        矩形脉冲包络，相位存在Barker编码，已知的barker码包括：2位（++或+-），3位（++-），4位（+++-或++-+），5位（+++-+），7位（+++--+-），11位（+++---+--+-），13位（+++++--++-+-+）。下图分别展示了13位Barker相位编码脉冲波形的时域、频谱、时频图以及模糊函数特性。

 ​​

时域实部和虚部                                   时域模值和相位 

​

频谱图                                                                时频图

​​

模糊函数                                                模糊函数（俯视图）

​​

模糊函数（0频）                                       模糊函数（0时延）

%------------------------------波形分析------------------------------------
%% 波形参数
waveform = phased.PhaseCodedWaveform;
waveform.SampleRate = 5e6;                              % 采样率
waveform.Code = 'Barker';                                 % 编码方式
waveform.ChipWidth = 10/waveform.SampleRate;              % 码片时间宽度
waveform.NumChips = 13;                                   % 2, 3, 4, 5, 7, 11, or 13
waveform.PRF = 5000;                                      % PRF间隔
xt = step(waveform);

%% 信号实部/虚部
figure;
plot(real(xt));
hold on;
plot(imag(xt));
title('信号实部/虚部');
legend('实部','虚部');

%% 信号模值/相位
figure;
plot(abs(xt));
hold on;
plot(angle(xt));
title('信号模值/相位');
legend('模值','相位');

%% 信号频谱
f=linspace(-waveform.SampleRate/2,waveform.SampleRate/2,length(xt));
figure;
plot(f,abs(fftshift(fft(xt))));
xlabel('采样率（Hz）');
title('信号频谱');

%% 信号时频图
window = blackman(16);
noverlap = 10;
nfft = 512;
[s,f,t] = spectrogram(xt,window,noverlap,nfft,waveform.SampleRate,'centered'); 
figure;
imagesc(t,f,abs(s));
xlabel('时间（s）');
ylabel('频率（Hz）');
title('信号时频图');

%% 模糊函数
[afmag,delay,doppler] = ambgfun(xt,waveform.SampleRate,waveform.PRF); 
[ambbk1,delaybk] = ambgfun(xt,waveform.SampleRate,...,
    waveform.PRF,'Cut','Doppler','CutValue',0);
[ambbkt2,dopplerbk] = ambgfun(xt,waveform.SampleRate,...,
    waveform.PRF,'Cut','Delay','CutValue',0);

figure;
mesh(delay,doppler,abs(afmag));
xlabel('时间（s）');
ylabel('频率（Hz）');
title('信号模糊函数');

figure;
plot(delaybk*10^6,ambbk1);
xlabel('时延 (\mu sec)');
title('0Hz多普勒频点处的自相关函数');

figure;
plot(dopplerbk/10^6,ambbkt2);
xlabel('多普勒频率 (MHz)');
title('0s时延处的自相关函数');

2.2、Frank编码

        矩形脉冲包络，相位存在Frank编码，是LFM脉冲压缩波形的阶跃近似，需要是整数的平方。下图分别展示了36位Frank相位编码脉冲波形的时域、频谱、时频图以及模糊函数特性。

 ​

时域实部和虚部                                   时域模值和相位 

​

频谱图                                                                时频图

​​

模糊函数                                                模糊函数（俯视图）

​​​​​​​​​​​​​​​​​​​​​​​​​​​​​​​

模糊函数（0频）                                       模糊函数（0时延）

%------------------------------波形分析------------------------------------
%% 波形参数
waveform = phased.PhaseCodedWaveform;
waveform.SampleRate = 5e6;                             % 采样率
waveform.Code = 'Frank';                                  % 编码方式
waveform.ChipWidth = 10/waveform.SampleRate;              % 码片时间宽度
waveform.NumChips = 36;                                   % 整数的平方
waveform.PRF = 5000;                                      % PRF间隔
xt = step(waveform);

%% 信号实部/虚部
figure;
plot(real(xt));
hold on;
plot(imag(xt));
title('信号实部/虚部');
legend('实部','虚部');

%% 信号模值/相位
figure;
plot(abs(xt));
hold on;
plot(angle(xt));
title('信号模值/相位');
legend('模值','相位');

figure;
plot(phase(xt));

%% 信号频谱
f=linspace(-waveform.SampleRate/2,waveform.SampleRate/2,length(xt));
figure;
plot(f,abs(fftshift(fft(xt))));
xlabel('采样率（Hz）');
title('信号频谱');

%% 信号时频图
window = blackman(16);
noverlap = 10;
nfft = 512;
[s,f,t] = spectrogram(xt,window,noverlap,nfft,waveform.SampleRate,'centered'); 
figure;
imagesc(t,f,abs(s));
xlabel('时间（s）');
ylabel('频率（Hz）');
title('信号时频图');

%% 模糊函数
[afmag,delay,doppler] = ambgfun(xt,waveform.SampleRate,waveform.PRF); 
[ambfk1,delayfk] = ambgfun(xt,waveform.SampleRate,...,
    waveform.PRF,'Cut','Doppler','CutValue',0);
[ambfk2,dopplerfk] = ambgfun(xt,waveform.SampleRate,...,
    waveform.PRF,'Cut','Delay','CutValue',0);

figure;
mesh(delay,doppler,abs(afmag));
xlabel('时间（s）');
ylabel('频率（Hz）');
title('信号模糊函数');

figure;
plot(delayfk*10^6,ambfk1);
xlabel('时延 (\mu sec)');
title('0Hz多普勒频点处的自相关函数');

figure;
plot(dopplerfk/10^6,ambfk2);
xlabel('多普勒频率 (MHz)');
title('0s时延处的自相关函数');

2.3、P1编码

        P1、P2、P3和P4多相位编码 。这些编码是LFM脉冲压缩波形的阶跃近似 ，其距离副瓣低，且具有与LFM编码一样的多普勒容忍性，需要是整数的平方。下图分别展示了36位P1相位编码脉冲波形的时域、频谱、时频图以及模糊函数特性。

 ​​​​​​​​​​​​​​​

时域实部和虚部                                   时域模值和相位 

​

频谱图                                                                时频图

​​

模糊函数                                                模糊函数（俯视图）

​​​​​​​​​​​​​​​​​​​​​​​

模糊函数（0频）                                       模糊函数（0时延）

waveform = phased.PhaseCodedWaveform;
waveform.SampleRate = 5e6;                             % 采样率
waveform.Code = 'P1';                                     % 编码方式
waveform.ChipWidth = 10/waveform.SampleRate;              % 码片时间宽度
waveform.NumChips = 36;                                   % 整数的平方
waveform.PRF = 5000;                                      % PRF间隔
xt = step(waveform);

%% 信号实部/虚部
figure;
plot(real(xt));
hold on;
plot(imag(xt));
title('信号实部/虚部');
legend('实部','虚部');

%% 信号模值/相位
figure;
plot(abs(xt));
hold on;
plot(angle(xt));
title('信号模值/相位');
legend('模值','相位');

figure;
plot(phase(xt));

%% 信号频谱
f=linspace(-waveform.SampleRate/2,waveform.SampleRate/2,length(xt));
figure;
plot(f,abs(fftshift(fft(xt))));
xlabel('采样率（Hz）');
title('信号频谱');

%% 信号时频图
window = blackman(16);
noverlap = 10;
nfft = 512;
[s,f,t] = spectrogram(xt,window,noverlap,nfft,waveform.SampleRate,'centered'); 
figure;
imagesc(t,f,abs(s));
xlabel('时间（s）');
ylabel('频率（Hz）');
title('信号时频图');

%% 模糊函数
[afmag,delay,doppler] = ambgfun(xt,waveform.SampleRate,waveform.PRF); 
[ambP11,delayP1] = ambgfun(xt,waveform.SampleRate,...,
    waveform.PRF,'Cut','Doppler','CutValue',0);
[ambP12,dopplerP1] = ambgfun(xt,waveform.SampleRate,...,
    waveform.PRF,'Cut','Delay','CutValue',0);

figure;
mesh(delay,doppler,abs(afmag));
xlabel('时间（s）');
ylabel('频率（Hz）');
title('信号模糊函数');

figure;
plot(delayP1*10^6,ambP11);
xlabel('时延 (\mu sec)');
title('0Hz多普勒频点处的自相关函数');

figure;
plot(dopplerP1/10^6,ambP12);
xlabel('多普勒频率 (MHz)');
title('0s时延处的自相关函数');

2.4、P2编码

        需要是偶数的平方，下图分别展示了36位P2相位编码脉冲波形的时域、频谱、时频图以及模糊函数特性。

 ​​​​​​​​​​​​​​​

时域实部和虚部                                   时域模值和相位 

​

频谱图                                                                时频图

​​

模糊函数                                                模糊函数（俯视图）

​​​​​​​​​​​​​​​​​​​​​​​​​​​​​​

模糊函数（0频）                                       模糊函数（0时延）

%------------------------------波形分析------------------------------------
%% 波形参数
waveform = phased.PhaseCodedWaveform;
waveform.SampleRate = 5e6;                             % 采样率
waveform.Code = 'P2';                                     % 编码方式
waveform.ChipWidth = 10/waveform.SampleRate;              % 码片时间宽度
waveform.NumChips = 36;                                   % 偶数的平方
waveform.PRF = 5000;                                      % PRF间隔
xt = step(waveform);

%% 信号实部/虚部
figure;
plot(real(xt));
hold on;
plot(imag(xt));
title('信号实部/虚部');
legend('实部','虚部');

%% 信号模值/相位
figure;
plot(abs(xt));
hold on;
plot(angle(xt));
title('信号模值/相位');
legend('模值','相位');

figure;
plot(phase(xt));

%% 信号频谱
f=linspace(-waveform.SampleRate/2,waveform.SampleRate/2,length(xt));
figure;
plot(f,abs(fftshift(fft(xt))));
xlabel('采样率（Hz）');
title('信号频谱');

%% 信号时频图
window = blackman(16);
noverlap = 10;
nfft = 512;
[s,f,t] = spectrogram(xt,window,noverlap,nfft,waveform.SampleRate,'centered'); 
figure;
imagesc(t,f,abs(s));
xlabel('时间（s）');
ylabel('频率（Hz）');
title('信号时频图');

%% 模糊函数
[afmag,delay,doppler] = ambgfun(xt,waveform.SampleRate,waveform.PRF); 
[ambP21,delayP2] = ambgfun(xt,waveform.SampleRate,...,
    waveform.PRF,'Cut','Doppler','CutValue',0);
[ambP22,dopplerP2] = ambgfun(xt,waveform.SampleRate,...,
    waveform.PRF,'Cut','Delay','CutValue',0);

figure;
mesh(delay,doppler,abs(afmag));
xlabel('时间（s）');
ylabel('频率（Hz）');
title('信号模糊函数');

figure;
plot(delayP2*10^6,ambP21);
xlabel('时延 (\mu sec)');
title('0Hz多普勒频点处的自相关函数');

figure;
plot(dopplerP2/10^6,ambP22);
xlabel('多普勒频率 (MHz)');
title('0s时延处的自相关函数');

2.5、P3编码

        需要是整数的平方，下图分别展示了36位P3相位编码脉冲波形的时域、频谱、时频图以及模糊函数特性。

 ​​​​​​​​​​​​​​​

时域实部和虚部                                   时域相位 

​

频谱图                                                                时频图

​​

模糊函数                                                模糊函数（俯视图）

​​​​​​​​​​​​​​​​​​​​​​​​​​​​​​

模糊函数（0频）                                       模糊函数（0时延）

%------------------------------波形分析------------------------------------
%% 波形参数
waveform = phased.PhaseCodedWaveform;
waveform.SampleRate = 5e6;                             % 采样率
waveform.Code = 'P3';                                     % 编码方式
waveform.ChipWidth = 10/waveform.SampleRate;              % 码片时间宽度
waveform.NumChips = 36;                                   % 正整数
waveform.PRF = 5000;                                      % PRF间隔
xt = step(waveform);

%% 信号实部/虚部
figure;
plot(real(xt));
hold on;
plot(imag(xt));
title('信号实部/虚部');
legend('实部','虚部');

%% 信号模值/相位
figure;
plot(abs(xt));
hold on;
plot(angle(xt));
title('信号模值/相位');
legend('模值','相位');

figure;
plot(phase(xt));

%% 信号频谱
f=linspace(-waveform.SampleRate/2,waveform.SampleRate/2,length(xt));
figure;
plot(f,abs(fftshift(fft(xt))));
xlabel('采样率（Hz）');
title('信号频谱');

%% 信号时频图
window = blackman(16);
noverlap = 10;
nfft = 512;
[s,f,t] = spectrogram(xt,window,noverlap,nfft,waveform.SampleRate,'centered'); 
figure;
imagesc(t,f,abs(s));
xlabel('时间（s）');
ylabel('频率（Hz）');
title('信号时频图');

%% 模糊函数
[afmag,delay,doppler] = ambgfun(xt,waveform.SampleRate,waveform.PRF); 
[ambP31,delayP3] = ambgfun(xt,waveform.SampleRate,...,
    waveform.PRF,'Cut','Doppler','CutValue',0);
[ambP32,dopplerP3] = ambgfun(xt,waveform.SampleRate,...,
    waveform.PRF,'Cut','Delay','CutValue',0);

figure;
mesh(delay,doppler,abs(afmag));
xlabel('时间（s）');
ylabel('频率（Hz）');
title('信号模糊函数');

figure;
plot(delayP3*10^6,ambP31);
xlabel('时延 (\mu sec)');
title('0Hz多普勒频点处的自相关函数');

figure;
plot(dopplerP3/10^6,ambP32);
xlabel('多普勒频率 (MHz)');
title('0s时延处的自相关函数');

2.6、P4编码

        需要是整数的平方，下图分别展示了36位P4相位编码脉冲波形的时域、频谱、时频图以及模糊函数特性。

 ​​​​​​​​​​​​​​​

时域实部和虚部                                   时域相位 

​

频谱图                                                                时频图

​​

模糊函数                                                模糊函数（俯视图）

​​​​​​​​​​​​​​​​​​​​​​​​​​​​​​

模糊函数（0频）                                       模糊函数（0时延）

%------------------------------波形分析------------------------------------
%% 波形参数
waveform = phased.PhaseCodedWaveform;
waveform.SampleRate = 5e6;                             % 采样率
waveform.Code = 'P4';                                     % 编码方式
waveform.ChipWidth = 10/waveform.SampleRate;              % 码片时间宽度
waveform.NumChips = 36;                                   % 正整数
waveform.PRF = 5000;                                      % PRF间隔
xt = step(waveform);

%% 信号实部/虚部
figure;
plot(real(xt));
hold on;
plot(imag(xt));
title('信号实部/虚部');
legend('实部','虚部');

%% 信号模值/相位
figure;
plot(abs(xt));
hold on;
plot(angle(xt));
title('信号模值/相位');
legend('模值','相位');

figure;
plot(phase(xt));

%% 信号频谱
f=linspace(-waveform.SampleRate/2,waveform.SampleRate/2,length(xt));
figure;
plot(f,abs(fftshift(fft(xt))));
xlabel('采样率（Hz）');
title('信号频谱');

%% 信号时频图
window = blackman(16);
noverlap = 10;
nfft = 512;
[s,f,t] = spectrogram(xt,window,noverlap,nfft,waveform.SampleRate,'centered'); 
figure;
imagesc(t,f,abs(s));
xlabel('时间（s）');
ylabel('频率（Hz）');
title('信号时频图');

%% 模糊函数
[afmag,delay,doppler] = ambgfun(xt,waveform.SampleRate,waveform.PRF); 
[ambP41,delayP4] = ambgfun(xt,waveform.SampleRate,...,
    waveform.PRF,'Cut','Doppler','CutValue',0);
[ambP42,dopplerP4] = ambgfun(xt,waveform.SampleRate,...,
    waveform.PRF,'Cut','Delay','CutValue',0);

figure;
mesh(delay,doppler,abs(afmag));
xlabel('时间（s）');
ylabel('频率（Hz）');
title('信号模糊函数');

figure;
plot(delayP4*10^6,ambP41);
xlabel('时延 (\mu sec)');
title('0Hz多普勒频点处的自相关函数');

figure;
plot(dopplerP4/10^6,ambP42);
xlabel('多普勒频率 (MHz)');
title('0s时延处的自相关函数');

总结

        本文简单介绍常见脉冲雷达信号的时域、频域、时频域以及模糊函数特性，其中相位编码雷达信号主要有phased工具箱产生，可能和定义中的P1,P2,P3,P4码有所区别，有问题欢迎评论区留言。转载请附链接【杨（_> <_)】的博客_CSDN博客-信号处理,SAR,代码实现领域博主

【杨（_> <_)】
关注
  • 3
    点赞
  • 踩
  • 21
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 2
    评论
专栏目录
脉冲雷达信号处理流程
wangjie36的博客
04-13 9627
1，FPGA AD采集回波：连续M个脉冲回波经过相干解调之后的基带数据形成的一个二维数据矩阵，上图只是雷达三维数据块的一个切面，其中数据的每一行代表的是对同一个距离单元的一连串脉冲测量，每一列对应的是一个脉冲回波的连续采样，也叫连续的距离门，其中每个小方块的元素都是一个复数，代表每个个距离单元的I,Q分量。 2，IQ路信号进行数字下变频处理：核心是将中频 A /D 采样信号与 DDC中的数字控制振荡器 (NCO )产生的本地数字中频载波信号进行混频 ,将中频信号下变频到基带。 3，滤波处理
脉冲多普勒（PD）雷达信号处理流程详解
jiashaocai1986的博客
12-31 2195
脉冲多普勒雷达特点简介* 脉冲串积累到底是怎么实现的* 有限长脉冲串的频谱计算过程* 速度与距离模糊的介绍与解决方法* 通过举具体的例子来说明脉冲积累与多普勒测速的实现。
mimo雷达信号处理_雷达学术入门脉冲雷达信号处理概述
weixin_39901203的博客
12-05 2233
Reviewed by ：@甜草莓 @Robert Zhou；前置知识：概率论与统计学。面向人群：本科生、研究生/信号处理博士。 编者：对于信号处理来说，雷达和通信一直是一体两面，从MIMO通信到MIMO雷达，从OFDM通信到Multicarrier雷达，很多通信和雷达领域的前沿都在事实上的相互影响。本篇文章将会介绍Fan Liu博士即将在IEEE Communication...
雷达成像技术之一---雷达波形
最新发布
Dayang88888的博客
05-14 574
雷达波形
信号类型（雷达）——脉冲雷达（三）
yjh_2019的博客
02-26 6134
雷达信号波形介绍：本文主要介绍脉冲雷达信号波形与分辨率的关系，并通过仿真对比实验说明相参的重要性。
第1讲 -- 雷达系统概述
诗小葵的博客
07-21 5031
雷达，是Radar(RadioDetectionandRanging)的音译，意思是“无线电探测和测距”，它是通过无线电/电磁波的方式获取目标的存在与否以及空间位置，因此雷达也被称为“无线电定位”。发射机向目标发射电磁波，经目标反射之后，由接收机接收，经过信号处理可获得目标至发射机的距离、（径向）速度、方位、高度等信息。能够全天时、全天候工作，穿透能力强，不受光照、雾云雨天气的干扰，因此雷达是非常重要的传感器。...............
雷达原理---脉冲雷达和连续波雷达
小席的学习日志
11-11 3万+
一、雷达分类 雷达按雷达信号的形式可分为脉冲雷达和连续波雷达。 脉冲雷达： 发射的波形为矩形脉冲，按一定的或交错的重复周期工作，是目前应用最广泛的雷达信号形式。常规脉冲雷达发射周期性的高频脉冲。 间歇式发射脉冲周期信号，并且在发射的间隙接收反射的回波信号，即收发间隔进行。 在近距离段存在探测盲区。 连续波雷达： 发射连续的正弦波，主要用来测量目标的速度。如果同时还要测量目标的距离，则需对发射的波形进行调制，如经过频率调制的调频连续波等。 发射连续波，并且发射的同时可以接收发射回来的回波信号，即收发可
雷达系统 学习笔记（四）——脉冲多普勒雷达1
热门推荐
Fying2016的博客
05-16 4万+
第五章 脉冲多普勒雷达 脉冲多普勒雷达特点：具有脉冲雷达的距离分辨力和连续波雷达的速度分辨力，有较强的抑制杂波能力，因而能在较强的杂波背景中分辨出目标回波。 5.1 脉冲多普勒雷达的基本概念 5.1.1 PD雷达定义 PD雷达是通过脉冲发射并利用多普勒效应检测目标信息的脉冲雷达。 具有如下三点特征： - 具有足够高的脉冲重复频率，以致不论杂波或所观测到的目标都没有速度模糊；...
课程设计小项目——雷达信号处理
06-29
在本“课程设计小项目——雷达信号处理”中，我们将探讨与雷达系统相关的信号处理技术，特别是使用MATLAB进行仿真和分析。MATLAB作为一种强大的数值计算和数据分析工具，被广泛应用于信号处理领域，尤其是雷达系统的...
空管二次雷达——二次雷达原理
11-01
"空管二次雷达——二次雷达原理" 本文将详细介绍空管二次雷达的原理和工作机理。 一、雷达基本概念 雷达（RADAR）的原意为无线电检测和测距（Radio Detection and Ranging），它起到对目标（空中飞行器）空中定位...
脉冲雷达仿真
04-12
仿真脉冲多卜勒雷达的信号处理。设脉冲宽度为各学生学号末两位数，单位为μs，重复周期为200μs，雷达载频为10GHz，输入噪声为高斯白噪声。目标模拟分单目标和双目标两种情况，目标回波输入信噪比可变（-35dB～10dB...
雷达信号分选源码 matlab
05-06
雷达信号通常以脉冲形式发送，每个脉冲代表一次发射事件。脉冲的特性如宽度、频率、功率等都会影响雷达的性能。 描述中提到“仿真产生的雷达信号脉冲”，这意味着代码首先会模拟雷达发射的脉冲信号。在MATLAB中，...
MATLAB——radar.zip_matlab雷达_radar signal_雷达 仿真_雷达信号_雷达信号matlab
07-13
本资源"MATLAB——radar.zip"提供了一个全面的雷达信号处理工具箱，专为从事雷达研究和开发的工程师及学者设计。这个工具箱包含了丰富的雷达信号仿真和处理功能，有助于深入理解和实践雷达系统的理论与算法。 首先...
雷达信号处理技术汇总
yjh_2019的博客
07-31 4965
雷达信号处理技术介绍：以常规雷达信号处理流程为例简单介绍相关信号处理技术，包括I/Q技术，波束形成技术，脉冲压缩技术，杂波滤除技术等
小白学雷达系列之脉冲雷达基本组成（一）
pmz13579285的博客
08-29 1801
雷达，无线电监测与测距，主要用途就是检测、跟踪与成像。R=ct02R=2ct0​​其中c为光速，上述公式的实际模型是当雷达发射的脉冲信号传播至距离R时，雷达信号（脉冲或者其他）经过时间t雷达接收到距离R出反射的雷达信号，雷达电磁波传播的距离是2R则，所以需要乘以二分之一因子。雷达的优点：采用微波波段电磁波的极低衰减特性使得雷达可以很好的“看穿”云雾，优于光学观察，另外通过电磁波的穿透性等优点，可以达到透视、超视距观察等特点。检测概率PDP_DPD​。.........
雷达信号脉冲分选
qq_41744915的博客
03-07 510
PRI变换后脉冲分选
雷达原理---脉冲法测距
小席的学习日志
11-19 2万+
文章目录一、基本原理1. 测距的本质2. 定义回波到达时间tRt_RtR​的两种方法3. 回波脉冲中心估计的原理框图：二、影响测距精度的因素1. 电波传播速度变化产生的误差2. 大气折射引起的误差3. 测读方法误差三、距离分辨力和测距范围1. 距离分辨力1.1 定义1.2 雷达距离分辨力1.3 雷达作用距离与距离分辨力之间的矛盾1.4 匹配滤波与脉冲压缩1.5 匹配滤波是如何解决距离分辨力与能量之间的矛盾问题的？2. 测距范围2.1 最小可测距范围2.2 最大无模糊测距范围2.3 最大无模糊测距范围与积累脉
matlab 雷达信号分选 聚类_进行雷达系统设计和干扰分析
weixin_39882589的博客
11-20 894
使用 Keysight SystemVue - 进行雷达系统设计和干扰分析本应用指南列出了 Keysight SystemVue 软件在进行雷达系统设计和杂波/干扰分析方面的主要特性。将要讨论的部分关键领域包括: 如何实现雷达线性调频 (Chirp) 波形; 为发射机和接收机设计射频链路; 使用快速傅立叶变换 (FFT) 卷积进行脉冲压缩分析。最后，我们在有干扰和杂波信号的环境中对雷达系统进行了测...
信号处理的MATLAB仿真.pdf
03-07
在MATLAB中，可以模拟生成不同类型的雷达信号，如线性调频连续波(LFM)、脉冲多普勒(PD)等。 2. 噪声和杂波的产生：在接收信号中，除了目标回波外，还会包含各种噪声和杂波，如环境噪声、地面散射、海面反射、气象...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 2
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

【杨（_> <_)】

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值