vv.m

自己做的一个室内无线通信仿真

function [c1]=vv(y)


s0=[0 0];
s1=[1 0];
s2=[0 1];
s3=[1 1];
%代表了状态机的四种状态，
s01=[0 0 0];
s02=[1 1 0];
s11=[1 1 1];
s12=[0 0 1];
s21=[1 0 0];
s22=[0 1 0];
s31=[0 1 1];
s32=[1 0 1];

%代表了格图中每一个部分路径的输出码子和，Si,j， i表示第几个状态，j表示上或者下面的路径。
Ps1=0;
Ps2=0;
Ps3=0;
Ps4=0;
%幸存路径的量度
ds01=0;
ds02=0;
ds11=0;
ds12=0;
ds21=0;
ds22=0;
ds31=0;
ds32=0;

%部分路径的量度
   
c1=0;
c2=0;
c3=0;
c4=0;
%时时的码子估计
us1=[];
us2=[];
us3=[];
us4=[];
%保存上一个时间的幸存路径度量
%第一次路径选择
temp=y([1 2]);
ds01=distance(temp,s01([1:2]));
ds11=distance(temp,s11([1:2]));
us1=[us1 s01(3)];
us2=[us2 s11(3)];
c1=us1;
c2=us2;
Ps1=ds01;
Ps2=ds11;
Ps3=0;
Ps4=0;


Ps1_=ds01; %保存上次的路径
Ps2_=ds11; %
Ps3_=0;
Ps4_=0;


%第二次路径选择
temp=y([3 4]);
ds01=distance(temp,s01([1:2]));
ds11=distance(temp,s11([1:2]));
ds21=distance(temp,s21([1:2]));
ds31=distance(temp,s31([1:2])); %分支度量

c1=[us1  s01(3)];
c2=[us1  s11(3)];
c3=[us2  s21(3)];
c4=[us2  s31(3)]; %预测的输入

Ps1=Ps1_+ds01;
Ps2=Ps1_+ds11;
Ps3=Ps2_+ds21;
Ps4=Ps2_+ds31; %路径长度

us1=c1;
us2=c2;
us3=c3;
us4=c4; %保存

Ps1_=Ps1;
Ps2_=Ps2;
Ps3_=Ps3;
Ps4_=Ps4; %保存

%进入全状态的情况
for ii=3:length(y)/2-2  
    temp=y([2*ii-1,2*ii]);
ds01=distance(temp,s01([1:2]));
ds11=distance(temp,s11([1:2]));
ds21=distance(temp,s21([1:2]));
ds31=distance(temp,s31([1:2]));
ds02=distance(temp,s02([1:2]));
ds12=distance(temp,s12([1:2]));
ds22=distance(temp,s22([1:2]));
ds32=distance(temp,s32([1:2])); % 每个分支都要度量
   
%进入选择比较的单元

%对于第一节点的处理
   if (ds01+Ps1_)<=(ds02+Ps3_)
       Ps1=ds01+Ps1_;
       c1=[us1  s01(3)];       
   else
        c1=[us3  s02(3)];
        Ps1=ds02+Ps3;
    end
    
    %对于第二个节点的处理
    if ds11+Ps1_<=ds12+Ps3_; %还是依据格图的情况，将本次的分支度量与上次的路径相加
         Ps2=ds11+Ps1_; 
        c2=[us1  s11(3)]; %依据
    else
         Ps2=ds12+Ps3_;
        c2=[us3  s12(3)];
    end
    
    %对于第三个节点的处理
    if ds21+Ps2_<=ds22+Ps4_;        
        c3=[us2  s21(3)];
         Ps3=ds21+Ps2_;
    else
         Ps3=ds22+Ps4_;
        c3=[us4  s22(3)];
    end
    
    %对于第四个节点的处理
    if ds31+Ps2_<=ds32+Ps4_;
        c4=[us2  s31(3)];
         Ps4=ds31+Ps2_;
    else
        Ps4=ds32+Ps4_;
        c4=[us4  s32(3)];
    end
    
    
us1=c1;
us2=c2;
us3=c3;
us4=c4;
  
Ps1_=Ps1;
Ps2_=Ps2;
Ps3_=Ps3;
Ps4_=Ps4; %保存选择出来的路径
    
end 
    
% 进入收尾的比较
if (ds01+Ps1_)<=(ds02+Ps3_)
       Ps1=ds01+Ps1_;
       c1=[us1  s01(3)];
       
   else
        c1=[us3  s02(3)];
        Ps1=ds02+Ps3;
 end

if ds21+Ps2_<=ds22+Ps4_
        
        c3=[us2  s21(3)];
         Ps3=ds21+Ps2_;
    else
         Ps3=ds22+Ps4_;
        c3=[us4  s22(3)];
    end
    
us1=c1;
us3=c3;
  
Ps1_=Ps1;
Ps3_=Ps3;

if (ds01+Ps1_)<=(ds02+Ps3_)
       Ps1=ds01+Ps1_;
       c1=[us1  s01(3)];
       
   else
        c1=[us3  s02(3)];
        Ps1=ds02+Ps3;
 end