1、伪随机序列由mgen 函数产生,共有1024个码元。对已得到的基带信号进行扩频调制,直接把基带信号S与产生的伪随机序列相乘,得到扩频信号。然后对已作扩频处理的信号作BPSK 载波调制,得到发射信号。发射信号通过存在高斯白噪声的信道,到达接到端,接收端首先对信号进行解调,将信号解调到中频段,然后将解调的信号与伪随机信号码作乘法运算,对解调信号进行解扩。解扩信号中由于存在高斯白噪声,所有要对噪声进行滤波处理,由于信号为低频基带信号,所以用由汉明窗设计成的FIR低通数字滤波器对解扩信号作低通滤波,产生的信号在经过检验判决,最后还原为原始信号。误码率的计算是由一个比较程序完成的,对判决得到的信号与原始
2、信号进 行位对位的比较,计算判决得到的信号中相比原始信号不同的码元的个数,然后将错误码元个数除以总的码元个数,得到在此信噪比下的误码率。3.系统参数信息码元数 8位二进制码调制方式 BPSK信息码速率 128bps扩频码速率 8bps采样率 128Hz信道参数 高斯白噪声信道,信噪比为自定义值4.系统结构 扩频通信系统的结构框图如图1所示:信息码产生恢复信道滤波扩频解扩白噪声扩频码产生图1 扩频通信系统的系统结构 原始信息码是8位双极性二进制信号,然后再通过时域上的抽样,最后变成信息码。扩频码是由函数mgen产生的十六位伪随机信号,然后对随机信号进行抽样,然后再复制加长,最后产生扩频码。mge
3、n函数的程序代码如下:. functionout=mgen(g,state,N)gen=dec2bin(g)-48;%二进制数转换成相应的数字向量M=length(gen);curState=dec2bin(state,M-1)-48;for k=1:N out(k)=curState(M-1); a=rem(sum(gen(2:end).*curState),2); curState=a curState(1:M-2);end 信息码与扩频码进行乘法运算,完成扩频过程。对扩频信号进行载波调制,应用BPSK(二进制相位键控)把扩频信号调制到发射频率上,形成发射信号。发射信号加高斯白噪声,通过函
4、数 awgn(signal,snr,measured)对扩频信号加入噪声干扰。这个函数的好处是可以通过对snr 的改变,实现信号的信噪比的变化。 加入噪声后的信号进入接收端,成为接受信号,对接受信号进行BPSK解调处理,然后对调后的信号进行解扩。解扩的方法是用与发射端进行扩频处理时使用的伪随机码完全一样的伪随机序列与解调的信号相乘,将信号解扩。 解扩后得到的信号由于仍存在噪声的干扰,因此进行低通滤波,滤波器选用比较常用的FIR 数字滤波器。 滤波后的信号进行检验判决,最后恢复原始信号。最后,通过对恢复信号码与原始信号码逐位的比较,得到误码的个数,最后经过计算,得出误码率。5系统工作原理(1)程
5、序代码流程图(2)扩频通信系统的误码率仿真程序如下:code_length=8;x=1 1 -1 1 -1 1 -1 -1; %产生信号码for i=1:8 %对信号进行抽样,抽样率为128 s(1+(i-1)*128:i*128)=x(i);end n=1:1:1024;subplot(3,2,1) plot(n,s);axis(1 1024 -1.5 1.5) ylabel(信息码);figure(1)PN_length=16;Pn=sign(mgen(19,8,PN_length)-0.5); %产生伪随机码,为双极型码 PN_length %对伪随机码进行抽样编码,抽样率为8 PN2(
6、1+(i-1)*8:i*8)=Pn(i);PN=PN2 PN2 PN2 PN2 PN2 PN2 PN2 PN2;subplot(3,2,2)plot(n,PN);axis(1 1024 -1.5 1.5)扩频码figure(1);k_signal=s.*PN; %对已编码的信号进行扩频 subplot(3,2,3)plot(n,k_signal);扩频信号)bpsk_m=zeros(1,1024); %对扩频信号作BPSK调制 j=1;k=1;1024 if(j=2) k=k+1;end; % 0.05*pi 为初始相位,可以任意改变,这里设为0.05*pi bpsk_m(1,i)=k_sig
7、nal(1,k)*sin(2*pi*0.5*i+0.05*pi)+k_signal(1,k)*cos(2*pi*0.5*i+0.05*pi);j=j+1;subplot(3,2,4)plot(n,bpsk_m);BPSK调制信号y=awgn(bpsk_m,-4,measured %加入白噪声,此处的信噪比为-4dBsubplot(3,2,5)plot(n,y);axis(1 1024 -2 2)加入高斯白噪声后的调制信号bpsk_m1=y.*sin(2*pi*0.5*n); %解调bpsk_m2=y.*cos(2*pi*0.5*n);signal_jk1=bpsk_m1.*PN; %解扩sig
8、nal_jk2=bpsk_m2.*PN;signal_j=signal_jk1+signal_jk2;subplot(3,2,6)plot(n,signal_j);解扩的信号N=1024;h0=zeros(1,N);f=0 0.1 0.1 1; w0 = 1 1 0 0;b=fir2(40,f,w0); %由汉明窗构成的40阶低通滤波器h,w=freqz(b,1,N/2);h0(1,1:N/2)=abs(hN/2h0(1,N-i+1)=h0(1,i);tempx=fft(signal_jk1);tempx=tempx.*h0; %低通滤波tempx=ifft(tempx);real_x=rea
9、l(tempx);tempx=h0.*fft(signal_jk2);real_x1=real(tempx);figure(2)subplot(2,2,1)plot(n,real_x1+real_x);axis(1 N -1.5 1.5);滤波后的信号In=real_x1+real_x; % 取了两路之和N % 滤波后整形if(In(1,i)0) % 判决,得到解调结果In(1,i)=1;elseIn(1,i)=-1;subplot(2,2,2)plot(n,In);判决后的最终信号%计算误码率r=0;for j=1:N %计算错误码元的个数if(s(j)=In(j)r=r+1;Pe=r/N; %计算误码率,误码率=错误码元的个数/所有信号码元个数Pe6.计算机仿真结果在信噪比为-4dB的情况下,仿真结果如图:(2)信号的恢复,恢复信号如图: 有前后两幅仿真结果可以看出,信号经过扩频-调制-加噪声-解调-解扩-滤波-判决一系列过程之后,最终恢复为原来的信号。程序运行的误码率为0.88%,这个结果可以说明在信噪比较低的环境中,扩频通信系统也具有良好的性能。
copyright@ 2008-2022 冰豆网网站版权所有
经营许可证编号:鄂ICP备2022015515号-1