移动通信课程设计实验报告-利用matlab进行m序列直接扩频仿真(8)

else x1(i)=0; end end

t=0:N-1;

subplot(2,1,1) stem(t*0.08,x1); xlabel('t/s');

title('扩频前待发送二进制信息序列1'); x_rand=rand(1,N); for i=1:N

if x_rand(i)>=0.5 x2(i)=1;a=a+1; else x2(i)=0; end end

subplot(2,1,2) stem(t*0.08,x2);

title('扩频前待发送二进制信息序列2'); xlabel('t/s');

- 36 -

3、用两个沃尔什码分别调制两路信号

figure(3)

y1=rectpulse(x1,8);y2=rectpulse(x2,8); for i=1:1:50

M1(8*i-7)=m1(1);M2(8*i-7)=m2(1); M1(8*i-6)=m1(2);M2(8*i-6)=m2(2); M1(8*i-5)=m1(3);M2(8*i-5)=m2(3); M1(8*i-4)=m1(4);M2(8*i-4)=m2(4); M1(8*i-3)=m1(5);M2(8*i-3)=m2(5); M1(8*i-2)=m1(6);M2(8*i-2)=m2(6); M1(8*i-1)=m1(7);M2(8*i-1)=m2(7); M1(8*i)=m1(8);M2(8*i)=m2(8); end

subplot(2,1,1) tt=0:8*N-1; for i=1:400;

x11(i)=1-2*y1(i);x12(i)=1-2*y2(i);s1(i)=M1(i)*x11(i);s2(i)=M2(i)*x12(i); end

subplot(2,1,1) stem(tt/100,s1);

title('扩频后的待发送序列码1'); xlabel('t/s'); subplot(2,1,2) stem(tt/100,s2);

title('扩频后的待发送序列码2'); xlabel('t/s');

4、两路信号相加，并进行BPSK调制 figure(4)

subplot(2,1,1) s=s1+s2;

stem(tt/100,s);

title('两路信息相加,即进行walsh码复用'); xlabel('t/s'); subplot(2,1,2) fs=2000;

ts=0:0.00001:4-0.00001; s_b=rectpulse(s,1000);

s_bpsk=s_b.*cos(2*pi*fs*ts); plot(ts,s_bpsk); xlabel('t/s');

axis([0.065,0.095,-2.2,2.2])

title('walsh码分复用做BPSK变换');

对比walsh码进行复用后的信号与无码分多址系统的信号，可以看出，无码分多址系统的信息码只有两个取值-1和1，但是复用后存在三个值-1,0,1，所以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;%将窄带干扰叠加到信号上

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

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

- 40 -