Как преобразовать итерационный цикл в векторизированную операцию в matlab?

Question

[mart kurtz]

Уважаемые коллеги, прошу ваших советов!

есть предельно простой цикл, написанный в matlab:

Fi(1) = 0
for k = 2:N
Fi(k) ~ F(k-1)......
end

проблема в том, что N = 10^7...10^8 и все считается очень долго.
Для ускорения пробую parfor, но возникла проблема: невозможно вычислить текущее значения без вычисления предшествующего результата Fi(k) ~ F(k-1). Поэтому возникает вопрос: возможно ли итерационный цикл переписать в векторном виде? И если да, то как?

Comments

[res2001]

Напишите целое выржение, по обрубку
Fi(k) ~ F(k-1)......
не возможно понять что вы хотите получить.
Тильда - это же логическое отрицание - унарный оператор, а у вас он употребляется с двумя аргументами.

[res2001]

Если бы тильда была бинарной операцией, то для векторизации можно было бы написать так:
Fi(2:N) ~ Fi(1:N-1)

[mart kurtz]

res2001!
~ это не отрицание, это я так неудачно обозначил, что вычислить последующее значение можно только зная предыдущее.

Fi = 0

for n = 1:NN

tblock = t(((n-1)*block + 1):1:n*block);
siblock = Si(((n-1)*block + 1):1:n*block);

SLFM = S *siblock.* exp(1j*((f0+n*dfN)*tblock + Fi));

if (n < NN)
nn = n*block + 1;
tblocklast = t(nn);
Fi = rem(Fi - dfN*tblocklast, 2*pi);
end

if n == 1
kniWave = SLFM;
else
kniWave = horzcat(kniWave, SLFM);
end
end

[res2001]

mart kurtz,

вычислить последующее значение можно только зная предыдущее

В такой формулировке, очевидно, что никакая векторизация не поможет.
Но судя по коду выше у вас не совсем такая ситуация.
Сделайте tblock и siblock двумерными (может быть лучше для экономии памяти сделать двумерными t и Si) и все остальные операции производите разом над двумерными массивами.

[mart kurtz]

res2001, признателен вам за совет.
Сам также подумывал в этом направлении.