扩频通信变换域抗干扰技术matlab实现

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