首页->【超声成像实验室】

7 1

function [MV] = ui_MV(cha)
%chx=cha';
[row,col]=size(cha);
%sum=0;
%sum1=0;

w=zeros(row,col);
MV=zeros(row,1);

if col == 15
    L=7;
else
    L=8;
end
C=col-L+1;
a=ones(L,1);
for m=1:row
     x=cha(m,:);
     %采用平滑机制 防止出现 非奇异矩阵
     mv=0;
     %X=[x(1,);];
     for cnti=1:C
         cntj=cnti:cnti+L-1;
         y=x(1,cntj);
         X(cnti,:)=y;
     end
     
     R=X'*X/C;
     %R=inv(Rc);          %covariance matrix      
     %Ri=inv(R);         %Ri  inversion matrix
     num=R\a;
     den=(a'/R)*a;
     w=num./den;
     
     for cnti=1:C
         cntj=cnti:cnti+L-1;
         y=x(1,cntj);
         mv=mv+y*w;
     end  
        MV(m,:)=mv;
end
end

作者于2018-05-19 10:11:39修改!
楼主可见

  1. shwnyoo 1#

    最小方差波束合成算法例程