请为这段代码添加中文注释if ~existbSqrtvar bSqrt = 1;endif ~existfea_bvar isemptyfea_b aa = sumfea_afea_a2; ab = fea_afea_a; if issparseaa aa = fullaa; end D = bsxfunplusaaaa - 2a

% 判断变量bSqrt是否存在，如果不存在则赋值为1 if ~exist('bSqrt','var') bSqrt = 1; end

% 判断变量fea_b是否存在或为空 if (~exist('fea_b','var')) || isempty(fea_b) % 计算fea_a矩阵每行的平方和 aa = sum(fea_a.fea_a,2); % 计算fea_a矩阵的转置乘积 ab = fea_afea_a';

% 如果aa是一个稀疏矩阵，则将其转换为密集矩阵
if issparse(aa)
    aa = full(aa);
end

% 计算距离矩阵D
D = bsxfun(@plus,aa,aa') - 2*ab;
% 将小于0的值赋为0
D(D<0) = 0;
% 如果bSqrt为真，则对距离矩阵D进行开方
if bSqrt
    D = sqrt(D);
end
% 将距离矩阵D的上三角矩阵与下三角矩阵对称化
D = max(D,D');

else % 计算fea_a矩阵每行的平方和和fea_b矩阵每行的平方和 aa = sum(fea_a.*fea_a,2); bb = sum(fea_b.fea_b,2); % 计算fea_a矩阵和fea_b矩阵的转置乘积 ab = fea_afea_b';

% 如果aa是一个稀疏矩阵，则将其转换为密集矩阵
if issparse(aa)
    aa = full(aa);
    bb = full(bb);
end

% 计算距离矩阵D
D = bsxfun(@plus,aa,bb') - 2*ab;
% 将小于0的值赋为0
D(D<0) = 0;
% 如果bSqrt为真，则对距离矩阵D进行开方
if bSqrt
    D = sqrt(D);
end

end