跳频信号时频分析源代码

%——————————————-------------------跳频信号（Alpha稳定分布）参数估计，对于alpha>1无需低阶，直接stft


clearclcclose all





%------------------------------------产生跳频信号
fs=4*10^3;   %采样频率
fk = [1.1 1.3 1.6 1.0 1.7 1.5 1.2 1.4]*10^3;   %跳频图案集
th=50*10^（-3）;   %跳频周期，每跳驻留时间

ze = zeros（10）;   %0表示信号延迟为零
x = []; sig = ze;
t = 0:1/fs:th-1/fs;
for i = fk;
    x = cos（2*pi*i*t）;
    sig = [sig x];   %sig为输出跳频信号
end

N = length（sig）;
   
%----------------------------------加入Alpha稳定分布噪声
         U=unifrnd（-pi/2pi/211600）;
        W=exprnd（111600）;

       alpha=0.8;

        X=（sin（alpha*U）./cos（U）.^（1/alpha））.*（cos（U-alpha*U）./W）.^（（1-alpha）/alpha）; %产生Alpha稳定分布

       m=２;
        
        sig_temp=sig.*（10.^m/20）./std（sig）;   %根据信噪比重新调整信