rx_radon_hurwitz.m

MATLAB Simulation of OFDM System

function decoupled_syms = rx_radon_hurwitz(rh_syms, channel_est)

decoupled_syms = zeros(size(rh_syms));

n_syms = size(rh_syms,2);
%移出可能多余的数据
if rem(n_syms,2) ~= 0
   rh_syms(:,n_syms) = [];
   n_syms = n_syms - 1;
end

%decoupled_syms = zeros(size(rh_syms));
%在发射端，有R1=H1*S1+H2*S2,   R2=-H1*conj(S2)+H2*conj(S1).
%故估计值为 S1_e=conj(H1)*R1+H2*conj(R2)=(abs(H1)^2+abs(H2)^2)*S1; (此处_e表示估计值）
%          S2_e=conj(H2)*R1-H1*conj(R2)=(abs(H1)^2+abs(H2)^2)*S2.
decoupled_syms(:,1:2:n_syms) = repmat(conj(channel_est(:,1)),1, n_syms/2).*rh_syms(:,1:2:n_syms) + ...
   repmat(channel_est(:,2),1,n_syms/2).*conj(rh_syms(:,2:2:n_syms));

decoupled_syms(:,2:2:n_syms) = repmat(conj(channel_est(:,2)),1, n_syms/2).*rh_syms(:,1:2:n_syms) - ...
   repmat(channel_est(:,1),1,n_syms/2).*conj(rh_syms(:,2:2:n_syms));