%該程序用來完成對輸入信號進行4QAM解調，屬于《鏈路級仿真軟件設計》程序二的4QAM解調模塊
%作者：趙亞利  
%編程日期：2005－3－7


function [llr]=deqam4(x,h,SNR_db)
%[llr]=deqam4(x)
%x是1*(Lp+8)的復向量，其中向量元素表示輸入到解調器中的信號
%h為信道估值向量，是1*(Lp+8)的復向量，其中向量元素表示對應采樣時刻的信道狀態。
%SNR_db為信道噪比。
%llr為x通過解調器后的輸出信號,是1*2(Lp+8) 的實向量，表示經過相關解調后的數據“似然比“信息。
%SNR_db
%SNR_db為信道信噪比，信號功率為1，噪聲功率為n_power
SNR_linr=exp(SNR_db*log(10)/10);
%噪聲方差
n_power=1/SNR_linr;
%n_power=1;
%len_input為輸入信號的長度
len_input=length(x);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
L=length(h)/4;
%均衡權值向量為q,由信道模塊給出的h來計算
%q=abs(h(4,:)).^2;
q=abs(h(1,:)).^2+abs(h(2,:)).^2+abs(h(3,:)).^2+abs(h(4,:)).^2;
%q=ones(1,375);
%%%%%%%%%下面是為單徑SCM服務的
%q=abs(h(1:L/4)).^2+abs(h(L/4+1:L/2)).^2+abs(h(L/2+1:3*L/4)).^2+abs(h(3*L/4+1:L)).^2;
%均衡輸出為r
%r=[q q q q];
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5
%計算低位的LLR信息
temp_1=real(x)>0;
temp_2=2*temp_1-1;
s0_L=(temp_2+j)*sqrt(2)/2;
s1_L=(temp_2-j)*sqrt(2)/2;

%循環判斷輸入信號x的實部是否大于0，從而確定相應的s0、s1值
%for I=1:len_input
%    if (real(x(I))>0)
%        s0_L(I)=sqrt(2)/2+j*sqrt(2)/2;%第一象限
%        s1_L(I)=sqrt(2)/2-j*sqrt(2)/2;%第四象限
%    else
%        s0_L(I)=-sqrt(2)/2+j*sqrt(2)/2;%第二象限
%        s1_L(I)=-sqrt(2)/2-j*sqrt(2)/2;%第三象限
%    end 
%end
%計算信號點與標準信號星座點之間的距離
         d0_L=abs(x-s0_L);
         d1_L=abs(x-s1_L);
%計算似然比信息llr
%llr_L=1.*(d1_L.^2-d0_L.^2)./(n_power.*(1./((abs(h)).^2)));
%llr_L=1.*(d1_L.^2-d0_L.^2)./(n_power.*(1./q));
llr_L=q.*(d1_L.^2-d0_L.^2);
%計算高位的LLR信息
temp_1=imag(x)>0;
temp_2=2*temp_1-1;
s0_H=(1+j*temp_2)*sqrt(2)/2;
s1_H=(-1+j*temp_2)*sqrt(2)/2;
%循環判斷輸入信號x的虛部是否大于0，從而確定相應的s0、s1值
%for I=1:len_input
%    if (imag(x(I))>0)
%        s0_H(I)=sqrt(2)/2+j*sqrt(2)/2;%第一象限
%        s1_H(I)=-sqrt(2)/2+j*sqrt(2)/2;%第二象限
%    else
%        s0_H(I)=sqrt(2)/2-j*sqrt(2)/2;%第三象限
%        s1_H(I)=-sqrt(2)/2-j*sqrt(2)/2;%第四象限
%    end   
%end
%%計算信號點與標準信號星座點之間的距離
       d0_H=abs(x-s0_H);
       d1_H=abs(x-s1_H);
%計算似然比信息llr
%llr_H=1.*(d1_H.^2-d0_H.^2)./(n_power.*(1./((abs(h)).^2)));
llr_H=(d1_H.^2-d0_H.^2).*q/n_power;
%將高低位排序，輸出解調結果的似然比序列 
llr_output=[];
for I=1:len_input
       llr_output=[llr_output,llr_L(I),llr_H(I)];
end
llr= llr_output;