butterfly.asv

VC6.0开发环境下开发的ATA接口控制程序

function [p,w,dam,dem]=butterfly(num1,num2,num3,num4,pathm,wm,datam,detecm,inp,seqn,N)
detvec=[0,1;2,3;0,1;2,3;-2,-1;0,1;-2,-1;0,1;-1,0;1,2;-1,0;1,2;-3,-2;-1,0;-3,-2;-1,0];
datvec1=zeros(16,1);
datvec2=ones(16,1);
datvec=[datvec1 datvec2];

a1=wm(num1,seqn)-2*detvec(num1,1)*inp(seqn);
a2=wm(num2,seqn)-2*detvec(num2,1)*inp(seqn);
if a1<=a2;
    wm(num3,seqn+1)=a1;
    pathm(num3,seqn)=num1;
    detecm(num3,seqn)=detvec(num1,1);
    datam(num3,seqn)=datvec(num1,1);
else 
      wm(num3,seqn+1)=a2;
    pathm(num3,seqn)=num2;
    detecm(num3,seqn)=detvec(num2,1);
    datam(num3,seqn)=datvec(num2,1);
end
a3=wm(num1,seqn)-2*detvec(num1,2)*inp(seqn);
a4=wm(num2,seqn)-2*detvec(num2,2)*inp(seqn);
if a3<=a4;
    wm(num4,seqn+1)=a3;
    pathm(num4,seqn)=num1;
    detecm(num4,seqn)=detvec(num1,2);
    datam(num4,seqn)=datvec(num1,2);
else 
      wm(num4,seqn+1)=a4;
    pathm(num4,seqn)=num2;
    detecm(num4,seqn)=detvec(num2,2);
    datam(num4,seqn)=datvec(num2,2);
end
p=pathm;
w=wm;
dam=datam;
dem=detecm;