allign3.m

来自「基于卷积信号的MIMO系统盲信号估计」· M 代码 · 共 43 行

function [PERM]=allign3(estC2,estC0,n)

%%% solve for the permutation ambiguioty between two PARAFAC iteration

[m,n]=size(estC0);

estC2_norm=zeros(m,n);

estC0_norm=zeros(m,n);

PERM=zeros(n,n);

for jj=1:n
    estC2_norm(:,jj)=estC2(:,jj)/std(estC2(:,jj));
    estC0_norm(:,jj)=estC0(:,jj)/std(estC0(:,jj));
end

for ii=1:n
    for jj=1:n
        PERM(ii,jj)=squeeze(estC2_norm(:,jj)).'*squeeze(estC0_norm(:,ii));
    end
end

PERM1=abs(PERM);

tmp=zeros(n,n);
            
for ii=1:n
                
[V1,I1]=max(abs(PERM1));
            
[V2,I2]=max(V1);
                
tmp(I1(I2),I2)=1;
               
PERM1(I1(I2),:)=zeros(1,n);
                
PERM1(:,I2)=zeros(n,1);
                
end
            
PERM=tmp;