%該程序用來完成對輸入信號進行4QAM解調;
function [llr]=deqpsk(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為信道信噪比，信號功率為1，噪聲功率為n_power
SNR_linr=10^(SNR_db/10);
%噪聲功率
n_power=1/SNR_linr;

%len_input為輸入信號的長度
len_input=length(x);



%計算低位的LLR信息
temp_1=real(x)>0 ;%x的實部對應低位
temp_2=2*temp_1-1;
s0_L=(temp_2+j)*sqrt(2)/2;
s1_L=(temp_2-j)*sqrt(2)/2;%可以描述是一、三象限還是二、四，如果是一三就繼續根據下面的llr來區分具體是哪一個。

%循環判斷輸入信號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信息
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_output=[];
for I=1:len_input
       llr_output=[llr_output,llr_L(I),llr_H(I)];
end
llr= llr_output;