rou_mask.m

B3g_phase2_C语言_Matlab程序及说明

function [rou,h]=Rou_Mask(h_half,Ns)

Lh=length(h_half);
rou=abs(diag(h_half(:,1:Ns)));
h=zeros(Ns,2*Lh-Ns);
h(:,Lh-Ns+1:Lh)=h_half(:,1:Ns).'-diag(rou);

N=Lh/Ns-1;

for l=1:N
    hl=h_half(:,l*Ns+(1:Ns));
    h(:,Lh-(l+1)*Ns+(1:Ns))=hl.';
    h(:,Lh+(l-1)*Ns+(1:Ns))=conj(hl);
end