Реализация LMS для QAM-4 Опции

[Wanderer.I]

Добрый день! Мне нужно реализовать модель приёмника с LMS эквалайзером при КАМ модуляции. Проблема состоит в следующем, мне не совсем понятно как быть в данном случае с квадратурной и синфазной составляющей, должно ли в данном случае присутствовать два эквалайзера один для квадратурной, а другой для синфазной составляющей ( т.е. для каждой составляющей импульсная характеристика будет линейной) или же данный алгоритм должен работать в комплексном виде( импульсная характеристика комплексная)? Дело в том, что все найденные мной описания оперирует с действительными числами.

[thermit]

Что-то вроде этого должно быть:

Код

clear all;
c4=[1+j -1+j 1-j -1-j];
HI = [ 0.00001   0.0004    0.0002   -0.0010    0.0010        0.0039      -0.0120 ...
      -0.0300    -0.1800  -0.1000   -0.0300    0.0020        0.0200       0.0500 ...
       0.0950    -0.0200  -0.3500   -0.3000   -0.1000        0.3970       0.3500 ...
      -0.0300    -0.1800  -0.0223   -0.0100    0.0000205     -0.00000039 -0.000000017...
       0.000001   0.00000008];
  
HR = [ 0.00001   0.0008    0.0002   -0.0010    0.0010        0.0039      -0.0120 ...
      -0.0300    -0.1900   0.1000   -0.0300    0.0020        0.0200       0.0500 ...
       0.0950    -0.0200  -0.3500   -0.3400   -0.1000        0.470       0.3500 ...
      -0.0300    -0.1800  -0.0223   -0.0100    0.000205      0.00000039  -0.000000117...
       0.000001   0.00000008];

EL=60;
NP=4;
tx=c4(randint(1,10000,NP)+1);

h=complex(HR,HI);


rx=filter(h,1,tx);

u=0.5/(EL*std(rx)^2);
y=lms([zeros(1,(length(h)/2)+EL/2) tx],rx,EL,u);

plot(y(5000:end),'.');

lms.m

Код

function cs = lms(r,t,N,mu)

mem=zeros(1,N);
coeff=zeros(N,1);
cs=[];
for i=1:length(t)
    mem=[t(i) mem(1:end-1)];
    ms=mem*coeff;
    e=r(i)-ms;
    cs=[cs ms];
    coeff=coeff+mem'*mu*e;
end;

return;