error_prop.m

基于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(cnstl_type,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;
   r=filter(c,1,a)+n;
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;
   ra=filter(ca,1,a)+na;
   rb=filter(cb,1,a)+nb;
end;


z=zeros(Lep,1);
res=zeros(Lep,1);
hard=zeros(Lep,1);
offset=max([Nh Nf Nd pdelay delta]);
n=offset:Lep;
ahat=zeros(Lep,1);
ahat(n)=a(n-delta+1);

if (spacing==1)
   [soft,dec,res,hard]=ep(spacing,a,r,ahat,offset+1,f,d,...
                     delta-1,pdelay,cnstl_type,Ms,csig);
elseif (spacing==1/2)
   [soft,dec,res,hard]=ep(spacing,a,ra,rb,ahat,offset+1,fa,fb,d,...
               delta-1,pdelay,cnstl_type,Ms,csig);
end;


if (real_source)
   soft=real(soft);
   dec=real(dec);
   res=real(res);
   res=real(res);
   hard=real(hard);
end;