clc; close all
clear
P = bodeoptions; 
P.FreqUnits = 'Hz'; 
P.grid='on'; 

tfinal=0.006;
paso=1e-6;
t=0:paso:tfinal-paso;

s=tf('s');

% Parametros del convertidor 
Pn=6000;                % VALOR DE POTENCIA NOMINAL EN VA
Vpico=311;              % VALOR DE PICO MAXIMO DE LA TENSION DE SALIDA EN VOLTS
Vrms=220;               % VALOR EFICAZ DE LA TENSION DE SALIDA EN VOLTS
Vcc=350;                % VALOR DE LA FUENTE DE TENSION CC DE ENTRADA EN VOLTS 
In=Pn/Vrms;             % VALOR NOMINAL DE LA CORRIENTE DE CARGA EN AMPERES
Inpico=In*sqrt(2);      % VALOR DE PICO DE LA CORRIENTE DE CARGA NOMINAL
Rc=Vrms/In;             % Resistencia de carga para la potencia nominal en Ohms

% Parametros de la Planta. Filtro LC y Carga Resistiva Pura
L=1500e-6;              % Indutancia en uH
C=22e-6;                % Capacitancia en uF
rL=0.15;                % Resistencia serie del inductor en Ohm
wr=1/sqrt(L*C);         % Frecuencia angular de resonancia del filtro LC en rad/seg
fr=wr/2/pi;             % Frecuencia de resonancia del filtro LC en Hz

% Función de transferencia Ga=IL/(Vab - Vo) 
Ga=1/(s*L + rL);
% Función de transferencia Gb=Vo/Ic 
Gb=1/(s*C);
% Función de transferencia de la carga conocida Ioc/Vo=1/Zc, Zc=Rc
Gc=1/Rc;

%% Funcion de transferencia de la planta de tension
%% Relacion entre Vo y la tension de entrada al inversor Vab con W = 0

% Gv=minreal((Ga*Gb)/(1 + (Ga*Gb) + (Gb*Gc)));
Gv=(1/(L*C))/(s^2 + (1/Rc/C + rL/L)*s + ((Rc + rL)/(Rc*L*C)));

zpk(Gv)
[numGv,denGv] = tfdata(Gv,'v');
figure; step(Vcc*Gv,t)
legend('Respuesta al Escalon Entrada-Salida')

%% Funcion de transferencia de perturbación Gw con Vab = 0
Gw = -(s*L + rL)/(s^2 + (1/Rc/C + rL/L)*s + ((Rc + rL)/(Rc*L*C)));

zpk(Gw)
[numGw,denGw] = tfdata(Gw,'v');
figure; step(8.38*Gw,t)
legend('Respuesta al Escalon Perturbacion-Salida')

figure; bode(Gv,P);
title('Respuesta en Frecuencia Entrada-Salida')
figure; bode(Gw,P);
title('Respuesta en Frecuencia Perturbacion-Salida')