channel_estimation_rotate.m

来自「基于EMCA368国际标准的物理层超宽带协议的基本仿真程序,都是自己编写的」· M 代码 · 共 32 行

function [iv,qv]=channel_estimation_rotate(datai,dataq);


channel_est_r_fren=[0,1,1,-1,1,-1,-1,1,-1,1,-1,1,1,-1,-1,-1,1,1,-1,1,-1,1,1,1,-1,1,1,1,-1,-1,1,1,1,1,...
-1,-1,1,-1,1,1,1,-1,-1,1,1,-1,-1,1,-1,1,1,1,1,1,-1,1,1,-1,-1,-1,-1,-1,0,0,0,0,0,-1,...
-1,-1,-1,-1,1,1,-1,1,1,1,1,1,-1,1,-1,-1,1,1,-1,-1,1,1,1,-1,1,-1,-1,1,1,1,1,-1,-1,1,1,1,...
-1,1,1,1,-1,1,-1,1,1,-1,-1,-1,1,1,-1,1,-1,1,-1,-1,1,-1,1,1];

channel_est_i_fren=[0,1,1,-1,1,-1,-1,1,-1,1,-1,1,1,-1,-1,-1,1,1,-1,1,-1,1,1,1,-1,1,1,1,-1,-1,1,1,1,1,...
-1,-1,1,-1,1,1,1,-1,-1,1,1,-1,-1,1,-1,1,1,1,1,1,-1,1,1,-1,-1,-1,-1,-1,0,0,0,0,0,1,...
1,1,1,1,-1,-1,1,-1,-1,-1,-1,-1,1,-1,1,1,-1,-1,1,1,-1,-1,-1,1,-1,1,1,-1,-1,-1,-1,1,1,...
-1,-1,-1,1,-1,-1,-1,1,-1,1,-1,-1,1,1,1,-1,-1,1,-1,1,-1,1,1,-1,1,-1,-1]; 

tempi=zeros(1,128*6);
tempq=zeros(1,128*6);

for s=1:6;
    
    for k=2:62
        tempi((s-1)*128+k)=1./(channel_est_r_fren(k).*datai((s-1)*128+k)+channel_est_i_fren(k).*dataq((s-1)*128+k));
        tempq((s-1)*128+k)=1./(channel_est_r_fren(k).*dataq((s-1)*128+k)-channel_est_i_fren(k).*datai((s-1)*128+k));
    end
    for k=68:128
        tempi((s-1)*128+k)=1./(channel_est_r_fren(k).*datai((s-1)*128+k)+channel_est_i_fren(k).*dataq((s-1)*128+k));
        tempq((s-1)*128+k)=1./(channel_est_r_fren(k).*dataq((s-1)*128+k)-channel_est_i_fren(k).*datai((s-1)*128+k));
    end   
end
    
iv=tempi/sqrt(2);
qv=tempq/sqrt(2);