%——————————————-------------------跳頻信號（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）;???%根據信噪比重新調整信