%Los que no estan muy familiarizados con MATLAB puede agarrar este "codigo"
%y tomarlo como base para hacer un RC, RL, LC, RLC paralelo, trata de calcular las fases en
%funcion de la frecuencia (que no hago) y lo que se les ocurra.

%Consideren que cuando hagan otro circuito tengan que definir algunas
%variables. Por ejemplo, ¿Cual puede ser un buen parametro para reemplazar
% w_0 en un RC?

%Otra opcion es pasar al dominio temporal (estacionario)

%Disclaimer: Algunas cosas no fueron hechas con mucha rigurosidad

%% RLC Serie

clear i; %borro en caso de que hayan usado la i como variable o parametro en otro codigo, asi es la unidad imaginaria

%Parametros del circuito. Todo en MKS 
R = 1000;
L = 0.25;
C = 2.0000e-09; %Matlab reconoce esta notacion automaticamente
V_0 = 1; 

%Frecuencia de resonancia y vector de frecuencias
w_0 = 1/sqrt(L*C); 
w = linspace(10,100*w_0,500000); %Crea un vector entre 10 y 100*frecuencia de resonancia(para el diagrama de Bode), con 500000 puntos intermedios (tal vez un poco exagerado)

%Impedancias Complejas
Z_L = i*w*L;    %reconoce automaticamente "i" como la unidad imaginaria. OJO si la reescriben (por ejemplo, en un for)
Z_C = -i./(w*C);   % ./ le indica a Matlab que hay que hacer "i" dividido cada elemento del vector w

%Corriente del circuito y caida de tension en cada componente
I = V_0./(R + Z_L + Z_C); %I: corriente como numero complejo; diferencia mayusculas d
Tens_R = I*R;
Tens_C = I.*Z_C;  % hay que poner ".*" porque es la multiplicacion de dos vectores elemento a elemento. NO ES PRODUCTO ESCALAR
Tens_L = I.*Z_L;

%Grafico el modulo de las tensiones(abs) en funcion de la frecuencia
plot(w,abs(Tens_R)) %
hold on;    %con este comando no reemplaza el grafico en la figura que ya hicimos
plot(w,abs(Tens_C))
plot(w,abs(Tens_L))
xlabel('2*Pi*Frecuencia(Hz)')
ylabel('Tension(V)')
legend('Tension Resistencia','Tension Capacitancia','Tension Inductancia')
grid on; %pone la grilla en el grafico
xlim([10, 2*w_0]) %le dice que muestre el grafico entre esos limites para el eje x. Sino ponen esto(prueben) no se ve la campana salvo que este en log(x)

%Diagramas de Bode - No visto en clase al momento de hacer este script
figure;   %abre otra figura, sino graficaria en la que abrimos y "holdeamos" antes
loglog(w,abs(Tens_R)/V_0)
hold on;
loglog(w,abs(Tens_C)/V_0)
loglog(w,abs(Tens_L)/V_0)
xlabel('2*Pi*Frecuencia(Hz)')
ylabel('Ganancia(B)')
grid on;

% ¿Puede la caida de tension en algun componente del RLC serie ser mayor
% que la tension de la fuente? En caso afirmativo, ¿en cuales? ¿Trae algun
% inconveniente esto?¿Se puede usar como amplificador?¿Hay algun
% inconveniente?

%¿El circuito RLC es un filtro pasabajos, pasaaltos, pasabanda?¿Todo junto?
