直接序列扩频系统MATLAB仿真（BPSK调制）

13、比较解扩前后信号功率谱密度

figure(13) subplot(2,1,1) yjb=fft(yrel,N);

prelb=yjb.*conj(yjb)/N; plot(freb,prelb(1:N/2)*2/N); axis([0,200,0,0.01]);

title('解扩前信号功率谱');xlabel('Hz'); subplot(2,1,2) yj=fft(yrej,N);

prel=yj.*conj(yj)/N;

plot(freb,prel(1:N/2)*2/N); axis([0,200,0,0.01]);

title('解扩后信号功率谱'); xlabel('Hz');

如图，解扩后信号的频谱被压缩，功率幅度增加，符合理论分析结果。

- 21 -

14、对解扩信号进行采样、判决

figure(14) subplot(2,1,1) for i=1:1:350

ij=i*1000-500; ss(i)=yrej(ij); end

stem(ss);

title('解扩信号采样'); subplot(2,1,2) for i=1:1:350 %判决信号算法 if ss(i)>0.2 ss(i)=1; elseif ss(i)<-0.2 ss(i)=-1; else ss(i)=0; end end

for i=1:1:50 ij=7*i-6; if ss(ij)==0

ss(ij)=ss(ij+4); end end

for i=1:1:348 if ss(i)==0

ss(i)=ss(i+2); end end

for i=1:1:50

S(i)=ss(i*7-3); if S(i)==0 S(i)=S(i)+1; end

S(i)=(1-S(i))/2; end stem(S);

title('判决后的最终信号');

- 22 -

至此，扩频系统仿真完成，比较最终判决信号与原发送信号完全吻合。

- 23 -

（2）信道中加入窄带强干扰

加入窄带干扰的matlab系统流程图

100/7Hz 二进制比特信息 100Hz 7位双极性m序列 100Hz 扩频序列 2000Hz 载波cos4000πt BPSK调制信号 高斯白噪声 2040~2050Hz 窄带强干扰 恢复载波cos4000πt 凯萨尔滤波器低通滤波 100Hz 7位双极性m序列 采样、判决 - 24 -