function?[mmsemsef_optdelta]?=?MMSE_fse?（hhrecsnrLf）

%?=========================================================
%?function?[mmsemsef_optdelta]?=?MMSE_fse?（hhrecsnrLf）
%?---------------------------------------------------------
%?Find?the?MMSE?for?a?specific?estimated?channel.
%?The?two?channels?（hhrec）?must?be?of?the?same?length.
%?
%?hrec??:?the?reconstructed?channel
%?snr???:?the?SNR?at?the?equalizer?input
%?Lf????:?length?of?the?equalizer?to?use
%
%?mmse??:?minimum-MSE?achieved?over?all?delays
%?f_opt?:?the?equalizer?that?achieves?the?mmse
%?delta?:?the?optimal?delay
%
%?Author:?H.?Pozidis???September?23?1998
%?=========================================================

if?（length（h）-length（hrec）?~=?0）
??error（‘The?channels?must?be?of?the?same?length‘）;
end

if?（snr?>=?0）
??snr?=?10^（snr/10）;
??lambda?=?1/snr;
else
??lambda?=?0;
end
h?=?h（:）;
hrec?=?hrec（:）;
H?=?convmtx（hLf）;
C?=?convmtx（hrecLf）;
P?=?size（C1）;
He?=?H（1:2:P:）;
Ce?=?C（1:2:P:）;
P?=?size（Ce1）;

A?=?inv（Ce‘*Ce?+?lambda*eye（Lf））;
f_eq?=?A*Ce‘;
h_del?=?eye（P）;

mse_mat?=?（h_del-He*f_eq）‘*（h_del-He*f_eq）?+?f_eq‘*f_eq*lambda;
mse?=?diag（mse_mat）;
[mmsedelta]?=?min（mse）;
f_opt?=?f_eq（:delta）;

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件???????1195??2005-09-24?10:41??MMSE_fse.m

-----------?---------??----------?-----??----

?????????????????1195????????????????????1