Я думаю, что вам нужно немного реорганизовать, чтобы ваша функция выводила 3 значения для ваших переменных phaseSystem
и phaseFailureTimes
как n на 3 матрицы.
Позвольте мне вызвать эту функцию calculate_phase_failure
.Тогда скрипт / функция, которая вызывает calculate_phase_failure
, может накапливать результаты.Наконец, отдельный цикл в конце может генерировать ваши графики.Если количество элементов отличается для каждой итерации ваших циклов, вам может потребоваться использовать массив ячеек для накопления ваших результатов.
Вот пример для простейшего случая, когда количество элементов согласовано между итерациямииз ваших петель.
for i=1:n
[phaseSystem(:,:,i), phaseFailureTime(:,:,i)] = calculate_phase_failure( <input variables> );
end
% now generate your plot
for i=1:n
for phase=1:size(phaseFailureTime,2)
figure(phase);
stem(1:size(phaseSystem,1), phaseFailureTime(:,phase,i))
end
end