Matlab - For para ode45

function ode45caballero12

%  Estoy tratando de resolver un circuito electrónico llamado elevador de 

%  tensión. Este dispositivo está descrito por las ecuaciones presentes en 

%  los links de arriba. Este posee dos estados uno es On y el otro es Off y

%  cada estado presenta ecuaciones diferenciales un poco diferentes, las 

%  que he llamado como función 'boost1' y 'boost2' . Por ejemplo este es 

%  boost2 (Off state):

% 

%  function F=boost2(t,x)

%  C= 1.16e-3; %Valor del condensador

%  R= 20e3; %Valor de la resistencia

%  L= 0.16e-3; %Valor de la Inductancia

%  Vs= 12; %tensión de entrada

%  A=[0,-(1/L) ; (1/C),-(1/(R*C))];

%  B=[(1/L) ; 0];

%  u=Vs;

%  F=A*x+B*u;

%  

%  Boost 1 es similar. Estas ecuaciones fueron acomodadas como una ecuación

%  de estados. Las matrices A y B son la matriz de estados y la de entrada

%  y F es la ecuación de estados. 

% 

%  Este circuito está gobernado por un control, el cual según su salida

%  0 ó 24 debo elegir desarrollar el estado On ó el estado Off, por eso 

%  tengo ese IF. 

% 

%  Sé que con ODE45 puedo resolver una ecuación diferencial, una sola vez.

%  Pero yo necesito que se calcule varias veces según se escoja el IF y

%  graficar eso en el tiempo. Es por eso que creo que necesito un FOR. 

%  No sé como tomar el resultado del ode45 para que en el siguiente paso

%  del FOR esa sea la nueva condición inicial del ode45.

global C R L  Vin

C= 1.16e-3; %Valor del condensador

R= 20e-3; %Valor de la resistencia

L= 0.16e-3; %Valor de la Inductancia

Vin= 12; %tensión de entrada

I= [0 1];

y0= [0 0];

fprintf('CONDICIONES:\n')

fprintf('\nEJEMPLO: C=1.16e-3   R= 20e-3  L= 0.16e-3   Vn= 12   I= [0 10]    y0= [0 0]\n\n')

% C=input('Ingrese el valor de Condensador C: ');

% R=input('Ingrese el valor de Resistencia R: ');

% L=input('Ingrese el valor de la Inductancia L: ');

% Vn=input('Ingrese el valor de tension de entrada Vn: ');

% I=input('Ingrese el intervalo I=');

% y0=input('Ingrese condicion inicial y0=');

[t,y]=ode45(@elevando_tension1,[I(1) I(2)],[y0(1) y0(2)]);

figure(gcf)

plot(t,y(:,1), t,y(:,2));

end

function dydt=elevando_tension1(t,y)

global L  Vin C R

dydt=zeros(2,1);

dydt(1)=-(1/L)*y(1)+ (1/L)*Vin;    %  diL/dt= -(1/L)*Vc +(1/L)*Vin 

dydt(2)=(1/C)*y(1)-(1/(R*C))*y(2);  %  dVc/dt=(1/C)*iL-(1/(R*C))*Vc

end

>> ode45caballero12

CONDICIONES:

EJEMPLO: C=1.16e-3   R= 20e-3  L= 0.16e-3   Vn= 12   I= [0 10]    y0= [0 0]