Storage of variables being generated in a loop

2 views (last 30 days)
Hy guys!
My code uses random variables to generate times of up and down of components and feed a state map. and im using the sum of the times to calculate indexs.
I need to store theses indexes to make an average of them(10 is good enough)
Here's my code:
U=rand(n,1);
for i=1:n
%Primeiro Tup
T7=-1*log(U)/((lambdaT+lambdaP));
end
%0 = Up
%1= down transitorio
%2= down permanente
% Primeiro tempo, todos os componentes ON
Tclock7=zeros(1,n);
%Calcular o menor T para ver o componente que ira transitar de estagio
comp=find(T7==min(T7));
%Mudando estado do componente envolve sortear se caso esteja transitando
%de ON para OFF sera a falha permanente ou transitoria
if Tclock7(comp)==0
if rand>0.111
Tclock7(comp)=1;
mi=mit(comp);
else
Tclock7(comp)=2;
mi=mip(comp);
end
else
Tclock7(comp)=0
end
% Calculando Tdown do primeiro Componente off
T7(comp)=T7(comp)-1*log(U2)/((mi*trecho(comp))/horas);
% Encontrando próximo componente a transitar
comp=find(T7==min(T7));
%Mudando estado do componente
if Tclock7(comp)==0
if rand>0.111
Tclock7(comp)=1;
mi=mit(comp);
contT7=1;
contP7=0
else
Tclock7(comp)=2;
mi=mip(comp);
contP7=1;
contT7=0;
end
else
Tclock7(comp)=0
end
Tclock7=Tclock7;
T7=T7;
% Calculando novo Tup ou Tdown
% Essa rotina tem que ser automatizada até que todos os componentes cheguem
% a 8760 para fazer a avaliação dos criterios de disponibilidade de 1 ano
% Iteração para todos os componentes(17) chegarem a 1 ano de planejamento
%Variaveis de contagem
iter7=1;
contON7=0;
w=10;
%Loop para analisar todos os componentes
while length(S7) < 10
for all=1:n
%Loop para chegar até analise de 1 ano
while T7(all)<horas
iter7=iter7+1;
%Encontrar qual componente transita
comp=find(T7==min(T7));
%Caso o componente esteja DOWN
if Tclock7(comp)==1 | Tclock7(comp)==2
%Transicoes para ON
contON6=contON6+1;
%Calculo do TUP
T7(comp)=T7(comp) -1*log(rand)/((lambdaT+lambdaP));
%Faz a transiçao no mapa de estados
Tclock7(comp)=0;
% Caso o componente esteja UP
else Tclock7(comp)=0;
%Sorteio para ver se a falha sera transitoria ou permanente
if rand>0.111
% Puxo o mi transitorio
mi=mit(comp);
T7(comp)=T7(comp)-tmorto;
%Faz a transicao no mapa de estados
Tclock7(comp)=1;
% Contador de numero de falhas transitorias por ano
contT7=contT7+1;
%Puxo o mi permanente
else mi=mip(comp)
%Faz a transicao
Tclock7(comp)=2;
%Variavel aleatoria para calcular TDOWN
U3=rand;
%Contador de numero de falhas permanentes por ano
contP7=contP7+1;
end
%Calculo do TDOWN
T7(comp)=T7(comp) -1*log(U3)/((mi)*trecho(comp)/horas);
TDT7=(-1*log(U3)/((mi)*trecho(comp)/horas));
TDT7=TDT7+(-1*log(U3))/((mi)*trecho(comp)/horas);
% Indice de Frequencia
% SAIFI
consumidores=sum(consumidores);
saifi7=contP7/consumidores;
%SAIDI - indice de duraçao das interrupcoes sustentadas
saidi7=TDT7/consumidores;
% Indice de Confiabildidade
%ASAI
for m=1:10
asai7=(horas*n-TDT7)/(horas*n);
S7=[asai7(m)]
  14 Comments
Rik
Rik on 4 Dec 2019
Since you refuse to provide working code and don't provide any details about what you mean with "Doesnt work" it is very hard to help you. I can't run your code, so I can't diagnose your problem for you. The choice is yours:
  1. provide a fully working example (attach data in a mat file if you need to)
  2. provide much more details about what you mean by "Doesnt work"
  3. only provide incomplete descriptions of failure and/or code that cannot be run independently
If you go for option 3 I will not be able to help you. Please let me know your choice, because if it option 3 I know I am wasting my time.
Daniel Ribeiro Gonçalves
Thank you a lot for your time Rik! I got it by other way. The way you were sugesting would'nt work cuz my while loop would stop every time it reached the sample(horas) I seted.
So I just raised my time of iteration in my main while loop, and divided my indexes vector to the new sample I wanted to set(ask for input), and also created a for only to build my indexes..

Sign in to comment.

Answers (1)

Rik
Rik on 25 Nov 2019
The general solution to a problem like this is to put some code around it. The most general solution is to use cell arrays, but if you already know the dimension and type it is a better idea to use an array.
%example of working code:
a=10;b=100;c=30;
some_output=sum(rand(a,b)*c,'all');
some_output2=a*b/c;
array_out_1=cell(10,1);
array_out_2=cell(10,1);
for n=1:numel(array_out)
%first run the working code
a=10;b=100;c=30;
some_output=sum(rand(a,b)*c,'all');
some_output2=a*b/c;
%then store your variables of interest to the cell arrays
array_out_1{n}=some_output;
array_out_2{n}=some_output2;
end

Categories

Find more on Loops and Conditional Statements in Help Center and File Exchange

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!