Question

使用 Run＆amp;我的算法上的时间我发现在向整数添加标准偏差时有点慢。首先，我创建了大整数矩阵：

NumeroCestelli = 5;

lover_bound = 0;
upper_bound = 250;
steps = 10 ;

Alpha = 0.123

livello = [lover_bound:steps:upper_bound];
L = length(livello);

[PianoSperimentale] = combinator(L,NumeroCestelli,'c','r');
for i=1:L
    PianoSperimentale(PianoSperimentale==i)=livello(i);
end

然后我添加标准偏差（sigma = alpha * mu）和错误（称重）：

%Standard Deviation
NumeroEsperimenti = size(PianoSperimentale,1);

PesoCestelli = randn(NumeroEsperimenti,NumeroCestelli)*Alfa;
PesoCestelli = PesoCestelli.*PianoSperimentale + PianoSperimentale;
random = randn(NumeroEsperimenti,NumeroCestelli);
PesoCestelli(PesoCestelli<0) = random(PesoCestelli<0).*(Alfa.*PianoSperimentale(PesoCestelli<0) + PianoSperimentale(PesoCestelli<0));

%Error
IncertezzaCella = 0.5*10^(-6);
Incertezza = randn(NumeroEsperimenti,NumeroCestelli)*IncertezzaCella;
PesoIncertezza = PesoCestelli.*Incertezza+PesoCestelli;
PesoIncertezza = (PesoIncertezza<0).*(-PesoIncertezza)+PesoIncertezza;

有更快的方法吗？

Answer 1

您多次计算PesoCestelli<0。你可以只计算一次并保存值。您还可以创建一整套随机数，但只使用其中PesoCestelli<0的子集。您可以通过仅创建所需的随机数来加快速度。

目前尚不清楚Alfa是什么，但如果它是标量，而不是

Alfa.*PianoSperimentale(PesoCestelli<0) + PianoSperimentale(PesoCestelli<0)

可能会更快

(1+Alfa).*PianoSperimentale(PesoCestelli<0)

Answer 2

我没有足够的信息来测试它，但我敢打赌，消除你所做的所有重复计算将导致加速。我试图删除其中一些：

PesoCestelli = randn(NumeroEsperimenti,NumeroCestelli)*Alfa;
PesoCestelli = (1+PesoCestelli).*PianoSperimentale;
random = randn(NumeroEsperimenti,NumeroCestelli);
idx = PesoCestelli<0;
PesoCestelli(idx) = random(idx).*(1+Alfa).*PianoSperimentale(idx);

%Error
IncertezzaCella = 0.5*10^(-6);
Incertezza = randn(NumeroEsperimenti,NumeroCestelli)*IncertezzaCella;
PesoIncertezza = abs((1+PesoCestelli).*Incertezza);

请注意，我将最后两行缩减为一行。

生成随机数...更快的方式？

2 个答案: