function Sig_dist = Sig_dist_comp(Sig_Y) inner_product = Sig_Y * Sig_Y'; C = size(Sig_Y, 1); Sig_dist = max(diag(inner_product) * ones(1, C) + ones(C, 1) * diag(inner_product)' - 2 * inner_product, 0); Sig_dist = sqrt(Sig_dist); end

这段代码定义了一个名为 `Sig_dist_comp` 的函数，用于计算样本特征矩阵 `Sig_Y` 的样本间距离矩阵 `Sig_dist`。

首先，计算 `Sig_Y` 的内积矩阵 `inner_product`，即将 `Sig_Y` 乘以其转置。

然后，获取矩阵 `Sig_Y` 的行数 `C`。

接下来，通过以下公式计算样本间距离矩阵 `Sig_dist`：

Sig_dist = max(diag(inner_product) * ones(1, C) + ones(C, 1) * diag(inner_product)' - 2 * inner_product, 0);

其中，`diag(inner_product)` 返回 `inner_product` 的对角线元素，`ones(1, C)` 返回一个大小为 `(1, C)` 的全1矩阵，`ones(C, 1)` 返回一个大小为 `(C, 1)` 的全1矩阵。

最后，将 `Sig_dist` 中的每个元素开平方根，并将结果赋给 `Sig_dist`。

函数执行完毕后，会返回计算得到的样本间距离矩阵 `Sig_dist`。