clear;
L=1000;%信號長度
N=512;%FFT的點數
Fs=1000;%采樣率
Ts=1/Fs;%采樣間隔
Tc=N*Ts;%碼元持續時間
EbN0?=?3;
snr?=?10*log10（2）?+?EbN0?-?10*log10（N）;
Wc=2*pi/Tc;%信號頻譜主瓣
t=0:Ts:Tc*L-Ts;%單音干擾的時間
jsr=68;
pn_seq?=?randi（[01]1N）;??%產生0?1隨機分布序列
pn_seq?=?-2*pn_seq?+?1;????%極化
ori_data?=?randi（[01]1L）;
data?=?2*ori_data?-?1;
data?=?repmat（dataN1）;?%信號重復N點
data?=?reshape（data1L*N）;?%重復后的信號
data_dsss?=?data.*repmat（pn_seq1L）;??%擴頻處理
????energy?=?sum（data.^2）/（N*L）;
A1=sqrt（（10^（jsr/10）*energy/3）*2）;
A2=sqrt（（10^（jsr/10）*energy/3）*2）;
A3=sqrt（（10^（jsr/10）*energy/3）*2）;
noise1=A1*cos（（Wc/2）*t）;
noise2=A2*cos（（Wc/3）*t）;
noise3=A3*cos（（Wc/4）*t）;
data_dsss?=?awgn（data_dssssnr‘measured‘）;???%加入噪聲
data_dsss?=?data_dsss?+?noise1+noise2+noise3?;
data_dsss_delay?=?[zeros（1N/2）data_dsss（1:L*N-N/2）];
window?=?boxcar（N）.‘;
snr_lose_in_window?=?10*log10（sum（window）.^2/（N*sum（window.^2）））;
%?FFT頻域干擾抑制
for?j=1:L
data_pro_1=data_dsss（1+（j-1）*N:j*N）;?%取出兩組數據
data_pro_2=data_dsss_delay（1+（j-1）*N:j*N）;
data_pro_1?=?data_pro_1?.*?window;%分別加窗
data_pro_2?=?data_pro_2?