Question

我已经为我的程序调整了一些现有的代码，但是我遇到了一个我不知道原因的错误。我有N个观察数据，我的目标是将数据分解为增加的较小子样本，并对每个子样本进行计算。为了确定子样本大小将如何变化，程序找到N的除数并将其存储到数组OptN中。

dmin = 2;
% Find  OptN such that it has the largest number of 
% divisors among all natural numbers in the interval [0.99*N,N] 
N = length(x); 
N0 = floor(0.99*N);
dv = zeros(N-N0+1,1);
for i = N0:N,
    dv(i-N0+1) = length(divisors(i,dmin));
end
OptN = N0 + find(max(dv)==dv) - 1;
% Use the first OptN values of x for further analysis
x = x(1:OptN);
% Find the divisors >= dmin for OptN
d = divisors(OptN,dmin);  

function d = divisors(n,n0)
% Find all divisors of the natural number N greater or equal to N0
i = n0:floor(n/2);
d = find((n./i)==floor(n./i))' + n0 - 1; % Problem line

在功能除数中出现问题。我有'使用错误./矩阵尺寸必须同意。'但是，这适用于长度为60的输入数据，但是当我尝试长度为1058的数据时，它会给出上述错误。

Answer 1

我认为对于大型数据集，find(max(dv)==dv)可能会返回多个数字。所以OptN将成为一个向量，而不是标量。

然后i的长度（BTW不是MATLAB中变量的好名称，它也是一个复数i）将是不可预测的，可能与n不同，导致下一个语句中的维度错误

您可以尝试使用find(max(dv)==dv,1)来获得第一场比赛。或者添加一个循环。

矩阵维必须在matlab中同意错误？

1 个答案: