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

5、观察调制信号频谱，并经awgn信道加高斯白噪和窄带强干扰

figure(5)

subplot(3,1,1) N=400000;

yb=fft(s_bpsk,N); mag=abs(yb);

fb=(1:N/2)*100000/N; plot(fb,mag(1:N/2)*2/N); axis([1000,3000,0,0.25]);

title('码分复用扩频后调制信号频谱'); xlabel('f/Hz'); subplot(3,1,2)

s_bpska=awgn(s_bpsk,3,'measured');%经过awgn信号，信噪比3dB plot(ts,s_bpsk,ts,s_bpska);

title('码分复用扩频后经加噪过信道后的信号与原信号时域波形对比'); xlabel('t/s');

axis([0.0775,0.0825,-2.4,2.4]); subplot(3,1,3)

yba=fft(s_bpska,N); maga=abs(yba);

fb=(1:N/2)*100000/N; fd=200000;

Wp1=2*2040/fd;%%%%%%%%%%带通滤波器 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)*40000; magz=abs(ybz); dz=ds.*magz;

dsz=maga+dz;%将窄带干扰叠加到信号上

- 41 -

plot(fb,dsz(1:N/2)*2/N); axis([1000,3000,0,0.4]); xlabel('f/Hz');

title('码分复用扩频后经信道调制加窄带强干扰信号频谱');

- 42 -

6、接收机信号乘以恢复载波，观察时域和频域

figure(6)

subplot(2,1,1)

rez=real(ifft(dz,N));

re=(s_bpska+rez).*cos(2*pi*fs*ts); plot(ts,re);

axis([0.065,0.095,-2.8,2.8]); xlabel('t');

title('扩频后接收信号乘以恢复载波'); xlabel('t/s'); subplot(2,1,2) yre=fft(re,N); magre=abs(yre);

plot(fb,magre(1:N/2)*2/N);

title('码分复用扩频后乘以恢复载波后信号频谱'); axis([0,5000,0,0.16]); xlabel('f/Hz');

- 43 -

7、信号经凯萨尔窗低通滤波器

figure(7)

subplot(2,1,1)

magrel=bs.*magre;

plot(fb,magrel(1:N/2)*2/N);

title('码分复用信号经过凯萨尔窗函数低通滤波'); axis([0,200,0,0.25]); xlabel('f/Hz'); subplot(2,1,2) tm=(1:N)/N*4;

yrel=real(ifft(bs.*yre,400000)); plot(tm,yrel); xlabel('t/s');

title('经过凯萨尔窗函数滤波后时域波形');

可以看出，信号经过滤波后的时域波形前后部分畸变较大，说明用改善的walsh码进行码分复用时，不同用户间仍然后一定的码间干扰。

- 44 -