Matlab运算性能不稳定如何解决？

Question

这里是我的核心代码，我提供loc文件的大小为46312，我使用了size（loc，1）取出loc文件的行数并赋值给n，希望能生成nn的矩阵dist，但是使用dist(1:n, 1:n)生成的dist大小是468*468，请问这种情况要如何解决呢？

Answer 1

引自免费微信小程序：皆我百晓生

从您提供的代码中，我观察到您在计算距离时使用了两个嵌套的循环，这是导致您得到的 dist 矩阵大小不稳定的原因。在您给出的代码中，dist 矩阵的大小是 n*n，但是在计算过程中，由于使用了两个嵌套的循环，dist 矩阵的大小会随着计算的进行而不断增加。

要解决这个问题，我建议您将计算距离的过程拆分成两个步骤。首先，计算每一对点之间的欧几里得距离，然后对每个距离进行平滑处理。这样，在计算每一对点之间的距离时，不会影响到其他点的计算。

以下是解决此问题的步骤：

1. 计算每一对点之间的欧几里得距离。

% 计算欧几里得距离
dist_euclidean = pdist2(loc, 'euclidean');

1. 对每个距离进行平滑处理。

% 平滑距离
dist_smoothed = dist_euclidean smoothtype('nearest');

其中，smoothtype 可以是 'nearest'（最邻近）、'bicubic'（三次多项式）或者 'spline'（三次样条），具体选择哪种平滑方法取决于您的需求。

这样，您就可以避免在计算过程中改变 dist 矩阵的大小。希望这能帮助您解决问题。