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

function [llr]=debpsk(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);

%首先判斷是奇數位還是偶數位
for i=1:len_input
    if (mod(i,2)==1)
        s1=1;
        s2=-1;
        %計算信號點與標準信號星座點之間的距離
        d1=abs(x(i)-s1);
        d2=abs(x(i)-s2);
    else 
        s1=j;
        s2=-j;
         %計算信號點與標準信號星座點之間的距離
        d1=abs(x(i)-s1);
        d2=abs(x(i)-s2);
    end
%計算該位的似然比信息,越0越正，越1越負
llr(i)=1*(d2^2-d1^2)/(n_power*(1/((abs(h(i)))^2)));    
end