至此，扩频系统仿真完成，比较最终判决信号与原发送信号完全吻合。
- 21 -
（2）信道中加入窄带强干扰
加入窄带干扰的matlab系统流程图
100/7Hz 二进制比特信息 100Hz 7位双极性m序列 100Hz 扩频序列 2000Hz 载波cos4000πt BPSK调制信号 高斯白噪声 2040~2050Hz 窄带强干扰 恢复载波cos4000πt 凯萨尔滤波器低通滤波 100Hz 7位双极性m序列 采样、判决 15、在信道中加入2040~2050Hz窄带强干扰并乘以恢复载波
产生方法为：生成单位冲激信号，其频谱为白色。之后，用2040~2050Hz的带通滤波器对其进行滤波。之后，与恢复载波相乘。
figure(15) subplot(2,1,1) fd=200000;
Wp1=2*2040/fd; 40~2050Hz带通滤波器 Wp2=2*2050/fd; Wc1=2*2030/fd; Wc2=2*2060/fd; Ap=1; As=100;
W1=(Wp1+Wc1)/2; W2=(Wp2+Wc2)/2;
wdth=min((Wp1-Wc1),(Wc2-Wp2)); Nd=ceil(11*pi/wdth)+1; bd=fir1(Nd,[W1 W2]); zd(1)=1; %产生单位冲激信号 for i=2:1:350000 zd(i)=0; end
ds=abs(freqz(bd,1,400000,fd))';
ybz=fft(zd,N)*100000; magz=abs(ybz);
dz=ds.*magz; %冲激信号频谱乘以带通滤波器 dsz=maga+dz; %magrelz=magrel;
plot(freb,dz(1:N/2)*2/N,freb,maga(1:N/2)*2/N); xlabel('Hz');
axis([1700,2300,0,0.6]);
title('经信道后，加上一窄带强干扰的频谱'); subplot(2,1,2)
rez=real(ifft(dz,N)); ts=(1:N)/N*4;
yzz=rez.*cos(2*pi*2000*ts); %信号乘以恢复载波 yz=fft(yzz,N); magyz=abs(yz);
renz=real(ifft(maga,N));
ynzz=renz.*cos(2*pi*2000*ts);
- 23 -
ynz=fft(ynzz,N); magynz=abs(ynz);
plot(freb,magyz(1:N/2)*2/N,freb,magynz(1:N/2)*2/N); axis([0,5000,0,0.2]);
title('加窄带干扰的信号乘以恢复载波后频谱');
如图，绿色部分为原扩频信号，蓝色部分为新加上窄带干扰频谱。从图上可以看出，窄带干扰非常强。
16、对加窄带干扰的信号进行低通滤波并解扩
figure(16) subplot(3,1,1)
magyzl=bs.*magyz; %用凯赛尔窗低通滤波器滤波 magynzl=bs.*magynz;
plot(freb,magyzl(1:N/2)*2/N,freb,magynzl(1:N/2)*2/N); axis([0,200,0,0.2]); xlabel('Hz');
title('信号与窄带干扰经过凯萨尔窗低通滤波后频谱');
- 24 -
subplot(3,1,2)
yrnzl=real(ifft(bs.*yre,400000)); %对滤波后带干扰信号进行ifft变换 yrzl=real(ifft(magynzl,400000)); tm=(1:N)/N*4;
yrnzlj=jj.*yrnzl(1:350000); yrzlj=jj.*yrzl(1:350000);
plot(ts(1:350000),yrnzlj+yrzlj); xlabel('t');
axis([0,4,-1.5,1.5]);
title('加窄带干扰后解扩信号波形');
subplot(3,1,3) %分别画出信号与干扰经低通滤波后频谱 yzj=fft(yrzlj,N); magzj=abs(yzj); ynzj=fft(yrnzlj,N); magnzj=abs(ynzj);
plot(freb,magzj(1:N/2)*2/N,freb,magnzj(1:N/2)*2/N); axis([0,500,0,0.2]);
title('信号与窄带干扰解扩后频谱'); xlabel('Hz');
能够看出，由于窄带干扰与扩频序列不相关，经解扩后频谱幅度明显降低，对信号无太大影响，信号的时域波形只是增加了幅度很小的干扰，对采样判决没有影响。说明扩频系统能够很好抑制窄带强干扰。
- 25 -