实验二 RAKE接收机仿真.doc

实验二 RAKE接收机仿真
一、实验目的
（1）了解RAKE接收机的原理；
（2）分析比较3种不同合并算法的性能。
二、实验内容
（1）编写MATLAB程序，实现RAKE接收机；
（2）修改信噪比，观察3种合并算法 误码率。
三、实验原理
见教材P212-213。
四、实验步骤
（1）画出仿真程序流程图；
（2）运行MATLAB开发环境，编写程序；
（3）运行程序，观察实验结果；
（4）分别修改扩频因子、信噪比、数据长度、功率因子等参数，观察并比较实验结果；
（5）画出RAKE接收机的性能曲线并分析原因。
五、思考题
（1）本仿真采用三径支路，实际系统中径数要多得多，是否径数越多越好？（2）三种合并算法哪个最佳？
参考代码如下：
Numusers=1;
Nc=16; %扩频因子
ISI_Length=1; %每径延时为ISI_Length/2
EbN0db=[0:2:10]; %信噪比，单位dB
Tlen=8000; %数据长度
%误比特率的初始值
Bit_Error_Number1=0;
Bit_Error_Number2=0;
Bit_Error_Number3=0;
%每径功率因子
power_unitary_factor1=sqrt(6/9);
power_unitary_factor2=sqrt(3/9);
power_unitary_factor3=sqrt(1/9);
s_initial=randsrc(1,Tlen); %数据源
%产生walsh矩阵
Wal2=[1 1;1 -1];
Wal4=[Wal2 Wal2;Wal2 Wal2*(-1)];
Wal8=[Wal4 Wal4;Wal4 Wal4*(-1)];
Wal16=[Wal8 Wal8;Wal8 Wal8*(-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);
for i=1:Numusers
x0=s_initial(i,:).’*Wal16(8,:);
x1=x0.’;
s_Spread(i,:)=(x1(:)).’;
end
%将每个扩频后的输出重复为两次，有利于下面的延迟（延迟半个码元）
ray1(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);
ray3(2*ISI_Length+1:2*Tlen*Nc)=ray1(1:2*Tlen*Nc-2*ISI_Length);
for nEN=1:length(EbN0db)
en=10^(EbN0db(nEN)/10); %将Eb/N0的dB值转化成十进制数值
sigma=sqrt((32/(2*en)));
%接收到的信号demp
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)’;
%将walsh码重复为两次
Wal16_d(1:2:31)=Wal16(8,1:16);
Wal16_d(2:2:32)=Wal16(8,1:16);
%解扩后rdata1为第一径输出
rdata1=dt*Wal16_d(1,:).’;
%将walsh码延迟半个码片
Wal16_delay1(1,2:32)=Wal16_d(1,1:31);
%解扩后rdata2为第二径输出
rdata2=dt*Wal16_delay1(1,:).’;
%将walsh码延迟一个码片
Wal16_delay2(1,3:32)=Wal16_d(1,1:30);
Wal16_delay2(1,1:2)=Wal16_d(1,31:32);
%解扩后rdata3为第三径输出
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+rdata3))/3;
%选择式合并
u=[u1,u2,u3];
maxu=max(u);
if(maxu==u1)
rd_m3=real(rdata1);
else if(maxu==u2)
rd_m3=real(rdata2);
else rd_m3=real(rdata3);
end
end
%三种方法判决输出
r_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/(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(nEN)=Bit_Error_Number3/(Tlen);
end
semilogy(EbN0db,Bit_Error_Rate1,’*-’);hold on;
semilogy(EbN0db,Bit_Error_Rate2,’o-’);hold on;
semilogy(EbN0db,Bit_Error_Rate3,’+-’);
grid on;
legend(’最大比合并’,’等增益合并’,’选择式合并’);
xlabel(’信噪比(dB)’);
ylabel(’误比特率’);
title(’3种主要分集合并方式性能比较’);
EbN0db=[0:2:20]; %信噪比，单位dB