猝发信号的matlab实现

BPSK猝发信号及其短时傅里叶变换和能量谱密度：

代码实现：

close all;

clc;clear

%% 纯BPSK信号

snr = 15;

N = 512;

sig=anabpsk(N,50,0.25);

sig2=anabpsk(N,50,0.25);

y = awgn(sig,snr);

y2 = awgn(sig2,snr);

%% 噪声

figure

n01 = randn(1,1024);

n02 = randn(1,512);

n03 = randn(1,512);

%% 猝发BPAK信号

sig_contac = [n01 y’ n02 y2’ n03]’;

[S,F,T,P] = spectrogram(sig_contac(900:1100),kaiser(128,18),120,128,1E3,‘yaxis’);%选取 时间轴上900:1100的信号画能量谱密度作存在性检测

figure(1)

plot(F,S);

title(‘BPSK调制猝发信号短时傅里叶变换’)

axis([0 500 -50 50])

figure(2)

plot(F,P);

title(‘BPSK调制猝发信号能量谱密度’)

[tfr,t,f] = tfrstft(sig_contac,1:length(sig_contac),1024);

figure(3)

imagesc(t,f,abs(tfr));

title(‘BPSK调制猝发信号时频图’)