% Trabajo de Fin de Grado.
% Grado en Ingeniería Electrónica Industrial.
% Título: Diseño y control de un convertidor Back-to-Back: Análisis teórico
% y validación mediante simulación.
% Autor: José Antonio Ramírez Jurado.
% Tutor: Dr. Manuel Jesús López Sánchez.
% Cádiz, junio 2024.

%--------------------------------------------------------------------------

% Subtítulo: PI Lazo V_DC.

% Parámetros constantes para el diseño del controlador.
C_DC = 6.8e-3;              % Capacitancia condensador (Faradios).
Xi = 1;                     % Factor amortiguamiento.
t_e = 1e-2;                 % Tiempo establecimiento (Segundos).
omega_n = 5.86/t_e;         % Frecuencia natural (Radianes/Segundos).
v_d_1 = 100*sqrt(2);        % Tension pico AC VSC1 (Voltios).
                            % V_pico = V_rms * sqrt(2).

% Condiciones de operación.
V_DC0 = 500;                % Tension en enlace DC (Voltios).
i_d0_1 = 14.14213562;       % Corriente AC VSC1 (Amperios).
I_DC0_2 = -6;                % Corriente DC VSC2 (Amperios).


% Ecuación característica estándar de un sistema de segundo orden:
% s^2 + 2*Xi*omega_n s + omega_n^2 = 0
q_11 = 2*Xi*omega_n;
q_12 = omega_n^2;

% Ecuación característica del sistema en lazo cerrado:
% s^2 + (3*v_d_1)/(2*V_DC0*C_DC) * Kp s + (3*v_d_1)/(2*V_DC0*C_DC) * Ki = 0
q_21 = (3*v_d_1)/(2*V_DC0*C_DC);
q_22 = (3*v_d_1)/(2*V_DC0*C_DC);

% Calculo de Kp_DC y Ki_DC.
Kp_DC = q_11/q_21;
Ki_DC = q_12/q_22;

% Verificacion del controlador.
disp('Accion proporcional Kp:');
disp(Kp_DC);

% Verificacion del controlador.
disp('Accion integral Ki:');
disp(Ki_DC);

% Limites para la señal de control en amperios.
lim_i_d_1_upper = 75.07117;
lim_i_d_1_lower = -75.07117;

                                                            % Prueba 1
%--------------------------------------------------------------------------

% Desviación de I_DC_2 respecto al valor nominal.
Delta_I_DC_2 = -2*I_DC0_2;

                                                            % Prueba 2
%--------------------------------------------------------------------------

% Ruido en la medida de los sensores.
% Ruido sensor V_DC.
D_R_V = 30e-3;                 % Desviación tipica ruido voltaje.
V_R_V = D_R_V^2;               % Varianza ruido voltaje.

                                                            % Prueba 3
%--------------------------------------------------------------------------


% Para introducir filtro al ruido.
A = -1/(1e-4);
B = 1/(1e-4);
C = 1;
D = 0;

Gf1_ss = ss(A,B,C,D);






