error_prop.bak

基于Matlab

% this simulates transmission for Lep symbols
% to look at error propagation events
function [] = error_prop(Lep);
dfec_global; 
ui_private_global;

a=make_src(constellation,Lep);
source=a;
sig=10^(-max(SNR)/10);

if (spacing==1),
   n=make_noise(sig,real_noise,Lep);
   w=filter(f,1,n);
   noise=n;
   cnoise=w;
elseif (spacing==1/2)
   na=make_noise(sig,real_noise,Lep);
   nb=make_noise(sig,real_noise,Lep);
   w=filter(fa,1,na)+filter(fb,1,nb);
   noise=[na
          nb];
   cnoise=w;
end;


z=zeros(1,Lep);
res=zeros(1,Lep);
q=zeros(1,Lep);
ahat=zeros(1,Lep);
offset=Nh+Nd+pdelay+delta;
n=offset-Nh:Lep;
ahat(n)=a(n-delta+1);
DELTA=delta+pdelay-1;
Ndtrunc=min(Nh-DELTA,Nd);

for n=offset:Lep,
      A=a(n:-1:n-Nh+1);
      Ahat=ahat(n-pdelay:-1:n-pdelay+1-Ndtrunc);
      z(n)=A.'*h-Ahat.'*d+w(n);
      ahat(n)=quantize(constellation,z(n));
      q(n)=a(n-delta+1)-ahat(n);
      res(n)=z(n)-a(n-delta+1);
end;
dec=ahat;