Rake接收机
RAKE接收机是一种用于无线通信系统的接收机。它的名字来源于“并行相关键接收机”(Parallel Receiving and Correlating Key Receiver)的缩写。RAKE接收机的主要作用是通过多个接收天线和并行处理器来改善信号接收效率。
RAKE接收机通常用于高速数据传输系统,如手机通信系统和无线局域网(WLAN)。它可以帮助消除多径散射对信号的干扰,提高信号检测的准确性。
RAKE接收机的工作原理是,每个接收天线都会收到信号的不同版本,这些版本可能因为多径散射或其他原因而有所差异。RAKE接收机通过并行处理器对这些信号版本进行相关,然后将它们合并在一起,从而获得一个更准确的信号。
使用选择性合并、最大比值合并、等增益合并算法进行Rake接收机的仿真
选择性合并(Selective Combining)是RAKE接收机中最常用的合并算法。在选择性合并中,RAKE接收机会根据信号的信噪比(SNR)来选择最优的信号版本进行合并。
在最大比值合并(Maximum Ratio Combining)中,RAKE接收机会根据信号的功率比来选择最优的信号版本进行合并。
在等增益合并(Equal Gain Combining)中,RAKE接收机将所有信号版本的功率平均分配,然后将它们合并在一起。
这三种合并算法各有优劣,适用于不同的应用场景。在实际应用中,RAKE接收机可能会根据具体情况使用多种合并算法来获得最优的信号检测效果。
clc,clear; %使用选择性合并、最大比值合并、等增益合并算法进行Rake接收机的仿真
user = 1; %用户数
Nc = 16; %扩频因子
ISI_Length = 1; %码元长度
SNR = [0:1:10]; %信噪比(dB)
Tlen = 5000; %数据长度
Bit_Error_Number1 = 0; %误比特率初始值
Bit_Error_Number2 = 0;
Bit_Error_Number3 = 0;
power_unitary_factor1 = sqrt( 5/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( user, Tlen * Nc );
ray1 = zeros( user, 2*Tlen*Nc );
ray2 = zeros( user, 2*Tlen*Nc );
ray3 = zeros( user, 2*Tlen*Nc );
for i = 1:user
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( SNR )%EbN0dB
EbNo = 10^( SNR(nEN)/10 ); %信噪比转换
pn = sqrt( 32/(2*EbNo) );
%接收信号
demp = power_unitary_factor1*ray1+...
power_unitary_factor2*ray2+...
power_unitary_factor3*ray3+...
( randn( 1,2*Tlen*Nc )+randn( 1,2*Tlen*Nc )*i )*pn;
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 = dt*wal16_d(1,:).';
%WALSH码延迟
wal16_delay1(1,2:32) = wal16_d(1,1:31);
%解扩第二径输出
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 = 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
figure(1)%绘制图像
plot(SNR,Bit_Error_Rate1,'r*-');hold on;grid on;
plot(SNR,Bit_Error_Rate2,'bo-');hold on;grid on;
plot(SNR,Bit_Error_Rate3,'g.-');grid on;
legend('最大比合并','选择式合并','等增益合并');
xlabel('信噪比'),ylabel('误比特率');
title('三种分集合并方式性能比较');
figure(2)
semilogy(SNR,Bit_Error_Rate1,'b*-');hold on;grid on;
semilogy(SNR,Bit_Error_Rate2,'go-');hold on;grid on;
semilogy(SNR,Bit_Error_Rate3,'k+-');grid on;
legend('最大比合并','选择式合并','等增益合并');
xlabel('信噪比'),ylabel('误比特率');
title('三种分集合并方式性能比较');
%采用最大比合并时误码率最低
%Rake接收机可以较好的解决多径问题
%白噪声对Rake接收机误码率影响不大,随着信噪比的增大,Rake接收机的误码率迅速下降
clear all;%接收部分code
user = 1;
Nc = 16; %扩频因子
ISI_Length = 1;%码元长度
SNR = [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_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( user, Tlen*Nc );
ray1 = zeros( user, 2*Tlen*Nc );
ray2 = zeros( user, 2*Tlen*Nc );
ray3 = zeros( user, 2*Tlen*Nc );
for i = 1:user
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( SNR )
EbNo = 10^( SNR(nEN)/10 ); %信噪比转换
pn = sqrt( 32/(2*EbNo) );
%接收信号
demp = power_unitary_factor1*ray1+...
power_unitary_factor2*ray2+...
power_unitary_factor3*ray3+...
( randn( 1,2*Tlen*Nc )+randn( 1,2*Tlen*Nc )*i )*pn;
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 = dt*wal16_d(1,:).';
%WALSH码延迟
wal16_delay1(1,2:32) = wal16_d(1,1:31);
%解扩第二径输出
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 = 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
plot(SNR,Bit_Error_Rate1,'r*-');hold on;
plot(SNR,Bit_Error_Rate2,'bo-');hold on;
plot(SNR,Bit_Error_Rate3,'g.-');
legend('最大比合并','选择式合并','等增益合并');
xlabel('信噪比');
ylabel('误比特率');
title('三种分集合并方式性能比较');
clear all;%发送部分代码
Numusers = 1;
Nc = 16; %扩频因子
ISI_Length = 1; %码元长度
EbN0db = 2; %信噪比
Tlen = 10; %数据长度
%我们假设出的每径功率因子
power_unitary_factor1 = sqrt( 6/9 );
power_unitary_factor2 = sqrt( 2/9 );
power_unitary_factor3 = sqrt( 1/9 );
s_initial = randsrc( 1, Tlen ); %产生数据源1或-1
figure(1)
stairs(s_initial,'Color','r');
ylim([-1.5,1.5]);
title('信号源序列','Color','b')
xlabel('数据长度')
ylabel('幅度')
%根据扩频因子等于16产生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 );
ray2( 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 ); %信噪比转换
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 )';
%将walsh码重复
wal16_d(1:2:31) = wal16(8,1:16);
wal16_d(2:2:32) = wal16(8,1:16);
%解扩后第一径输出
rdata1 = dt*wal16_d(1,:).';
%将walsh码延迟
wal16_delay1(1,2:32) = wal16_d(1,1:31);
%解扩后第二径输出
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 = 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)';
end
figure(2)
subplot(3,1,1)
stairs(r_Data1,'Color','r');
ylim([-1.5,1.5]);
title('最大比合并方式','Color','b')
xlabel('码元长度')
ylabel('幅度')
subplot(3,1,2)
stairs(r_Data2,'Color','r');
ylim([-1.5,1.5]);
title('等增益合并','Color','b')
xlabel('码元长度')
ylabel('幅度')
subplot(3,1,3)
stairs(r_Data3,'Color','r');
ylim([-1.5,1.5]);
title('选择式合并','Color','b')
xlabel('码元长度')
ylabel('幅度')
figure(3)
subplot(3,1,1)
stairs(s_initial,'Color','r');
ylim([-1.5,1.5]);
hold on
stairs(r_Data3,'Color','g');
ylim([-1.5,1.5]);
title('最大比值合并和初始序列比较','Color','b')
xlabel('码元长度')
ylabel('幅度')
ylim([-1.5,1.5]);
subplot(3,1,2)
stairs(s_initial,'Color','r');
ylim([-1.5,1.5]);
hold on
stairs(r_Data2,'Color','g');
ylim([-1.5,1.5]);
title('等增益合并和初始序列比较','Color','b')
xlabel('码元长度')
ylabel('幅度')
ylim([-1.5,1.5]);
subplot(3,1,3)
stairs(s_initial,'Color','r');
ylim([-1.5,1.5]);
hold on
stairs(r_Data1,'Color','g');
ylim([-1.5,1.5]);
title('选择式合并和初始序列比较','Color','b')
xlabel('码元长度')
ylabel('幅度')
ylim([-1.5,1.5]);
Walsh码是一种二进制线性码,它可以用来编码信息序列。Walsh码的名字来源于美国数学家J. L. Walsh,他在1923年提出了这种码制。
Walsh码有许多不同的变体,其中最常用的是Hadamard码。Hadamard码是Walsh码的一种特殊形式,它的码字长度为2的整数次幂。Hadamard码具有良好的正交性,因此在信道均衡、多径抵消和调制解调等应用中广泛使用。
Walsh码在通信领域中有许多应用。例如,它可以用来进行频谱分配,以便在无线电频谱中有效地传输信息。Walsh码也可以用来进行信号加扰,以增强信息的安全性。此外,Walsh码还可以用来进行无线局域网(WLAN)的调制解调,以及CDMA手机通信系统的多址传输。