资源描述：
《实验九-Rake接收机仿真.docx》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、实验九Rake接收机仿真一、实验目的1.了解Rake接收机原理。2.比较3种合并算法性能。二、实验内容3.编写MATLAB程序，实现Rake接收机。4.修改SNR观察误码率变化。三、实验代码clearall;Numusers=1;Nc=16;ISI_Length=1;EbN0db=[0:1:30];Tlen=8000;Bit_Error_Number1=0;Bit_Error_Number2=0;Bit_Error_Number3=0;power_unitary_factor1=sqrt(6/9);power_unitary_factor2=sqrt(2/9);power_unitary_fa
2、ctor3=sqrt(1/9);s_initial=randsrc(1,Tlen);wal2=[11;1-1];wal4=[wal2wal2;wal2wal2*(-1)];wal8=[wal4wal4;wal4wal4*(-1)];wal16=[wal8wal8;wal8wal8*(-1)];s_spread=zeros(Numusers,Tlen*Nc);ray1=zeros(Numusers,2*Tlen*Nc);ray2=zeros(Numusers,2*Tlen*Nc);ray3=zeros(Numusers,2*Tlen*Nc);4fori=1:Numusersx0=s_initia
3、l(i,:).'*wal16(8,:);x1=x0.';s_spread(i,:)=(x1(:)).';endray1(1:2:2*Tlen*Nc-1)=s_spread(1:Tlen*Nc);ray1(2:2:2*Tlen*Nc)=ray1(1:2:2*Tlen*Nc-1);ray2(ISI_Length+1:2*Tlen*Nc)=ray1(1:2*Tlen*Nc-ISI_Length);ray2(2*ISI_Length+1:2*Tlen*Nc)=ray1(1:2*Tlen*Nc-2*ISI_Length);fornEN=1:length(EbN0db)en=10^(EbN0db(nEN)
4、/10);sigma=sqrt(32/(2*en));demp=power_unitary_factor1*ray1+power_unitary_factor2*ray2+power_unitary_factor3*ray3+(rand(1,2*Tlen*Nc)+randn(1,2*Tlen*Nc)*i)*sigma;dt=reshape(demp,32,Tlen)';wal16_d(1:2:31)=wal16(8,1:16);wal16_d(2:2:32)=wal16(8,1:16);rdata1=dt*wal16_d(1,:).';wal16_delay1(1,2:32)=wal16_d(
5、1,1:31);rdata2=dt*wal16_delay1(1,:).';wal16_delay2(1,3:32)=wal16_d(1,1:30);wal16_delay2(1,1:2)=wal16_d(1,31:32);rdata3=dt*wal16_delay2(1,:).';p1=rdata1'*rdata1;p2=rdata2'*rdata2;p3=rdata3'*rdata3;p=p1+p2+p3;u1=p1/p;u2=p2/p;u3=p3/p;rd_m1=real(rdata1*u1+rdata2*u2+rdata3*u3);rd_m2=(real(rdata1+rdata2+r
6、data3))/3;u=[u1,u2,u3];maxu=max(u);if(maxu==u1)rd_m3=real(rdata2);elseif(maxu==u2)rd_m3=real(rdata2);elserd_m3=real(rdata3);4endendr_Data1=sign(rd_m1)';r_Data2=sign(rd_m2)';r_Data3=sign(rd_m3)';Bit_Error_Number1=length(find(r_Data1(1:Tlen)~=s_initial(1:Tlen)));Bit_Error_Rate1(nEN)=Bit_Error_Number1/
7、Tlen;Bit_Error_Number2=length(find(r_Data2(1:Tlen)~=s_initial(1:Tlen)));Bit_Error_Rate2(nEN)=Bit_Error_Number2/Tlen;Bit_Error_Number3=length(find(r_Data3(1:Tlen)~=s_initial(1:Tlen)));Bit_Error_Rate3(n
