Quiero encontrar la tensión de salida de un circuito como este:
yesteesmicódigodeMatlabpararesolverlo:
functionfind_xt=0:0.001:10;x0=0;R1=90000;R2=1000;C=0.001;[t,x]=ode45(@rhs,t,x0);plot(t,x);functiondxdt=rhs(t,x)dxdt=-(C/R1+C/R2)*x+C/R1*abs(sin(100*pi*t));endendperomedalarespuestaincorrecta,porquelarespuestacorrectaescasiunvoltajedeCC.
¿Puede alguien ayudarme?
Parece que tu ecuación diferencial no está configurada correctamente. La ecuación básica para empezar es
\ $ I = C \ dfrac {\ mathrm {d} v} {\ mathrm {d} t} \ $
y para las resistencias
\ $ I = \ dfrac {v} {R} \ $
Eso significa que deberías tener algo como
\ $ \ dfrac {\ mathrm {d} v} {\ mathrm {d} t} = \ dfrac {-v} {R_1 C} + \ dfrac {V_s (t) - v} {R_2 C} \ $
donde he enterrado la sinusoide rectificada en \ $ V_s (t) \ $.
Observe que C aparece en los denominadores, no en los numeradores.
Además, la fuente solo carga el condensador a través de la resistencia de 90 kOhm, por lo que con sus valores actuales la carga tiene una constante de tiempo de 90 s. Querrá ejecutar la simulación durante 5 o 10 veces esa constante de tiempo para ver cómo la salida se establece alrededor del valor de estado estable.