sort_d_v.m

来自「这个第四代移动通信中」· M 代码 · 共 42 行

function [D_R_sort, V_R_sort] = sort_D_V(A, D_R, V_R);         
% this function is to reshape the matrix D_opt,V_opt,D_R,V_R,accordig to A,S

size_A = length(A);
size_D = length(D_R);

% to use angle of a Complex number
angle_A = diag(angle(A))/pi;       
angle_D_R = diag(angle(D_R(1:16,1:16)))/pi; 

i = 1;
while (i<= size_A)
    index = 1;
    while (index <= size_D)
        if abs(abs(angle_D_R(index)) - 1) < 0.0625
            angle_D_R(index) = 1;       %to identify the over search angle_D_R
        end;
%         if abs(abs(angle_D_R(index)) - 1) < 0.0625  % 1/MG
%             %D_R_new(3,3) = A(3,3);
%             D_R_new(3,3) = D_R(index,index);
% 		    V_R_new(:,3) = V_R(:,index);  
%             angle_D_R(index) = 100;     %to identify the over search angle_D_R
%             index = index + 1;
%             break;
%         end;
        if abs(angle_D_R(index) - angle_A(i)) < 0.0625  % 1/MG
		    D_R_new(i,i) = D_R(index,index);
		    V_R_new(:,i) = V_R(:,index);  
            angle_D_R(index) = 100;     %to identify the over search angle_D_R
            %index = index + 1;
            break;
        end;
        index = index + 1;
    end;
    i = i + 1;
end;

D_R_sort = D_R_new;
V_R_sort = V_R_new;