|
|
почему результат фильтрации массива x фильтром h рассчитанным по d совершенно не совпадает с d, хотя вроде бы должен быть похож на d значительно более чем x. x это бывшее d гауссово сглаженное, промодулированное gmsk c BT=0.3, добавлен аддитивный шум и затем продемодулированное обычным частотным детектором.
d-известная тренинг последовательность.
x-то что пришло из канала.
Пытаюсь рассчитать h канала по d и x.
Затем притменяю рассчитанную h к x и мечтается увидеть нечто похожее на d, но получается результат, который интерпретировать не могу.
Где залет?
d=[-1 1 1 -1 -1 1 1 1 1 -1 -1 -1 -1 1 -1 1 1];
x=[-0.3257 0.3746 -0.7216 -0.6853 0.2987 -0.3676 0.7288 0.6018 -0.6440 -0.9535 -0.9883 -0.6290 0.3427 -0.3447 0.6477 1.0000 0.9612];
M=length(d);
N=3;
y=zeros(1,M);
eg=zeros(1,M);
h=zeros(1,N);
h(round(N/2)+0)=1;
Px=sum(x.^2)/(M+1);
delta=1/(1*Px*N);
for n=N:M
x1=x(n:-1:n-N+1);
y1=h*x1';
e=d(n)-y1;
eg(n)=e;
h=h+delta*e*x1;
y(n)=y1;
end;
figure, plot(fftfilt(h,x),'-*'), hold on, plot(d,'-r*') , hold off, grid on
|
|
E-mail: 
info@telesys.ru
