% 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: IP Lazo Q_12.

% Parámetros constantes para el diseño del controlador.
t_e = 1e-2;                 % Tiempo establecimiento (Segundos).
v_d_12 = 100*sqrt(2);       % Tension pico AC VSC1 y VSC2 (Voltios).
                            % V_pico = V_rms * sqrt(2).

% Constante de tiempo de un sistema de primer orden:
% t_e = 5*tau ---> tau = t_e/5
tau = t_e/5;

% Constante de tiempo del sistema de control para la potencia reactiva:
% T_Q (s) = 1/(((3*v_d_12*Kp_Q-2)/(3*v_d_12*Ki_Q)) s + 1)
% tau = ((3*v_d_12*Kp_Q-2)/(3*v_d_12*Ki_Q))
% Ki_Q = ((3*v_d_12*Kp_Q-2)/(3*v_d_12*tau_p))
% Calculo de Kp_Q y Ki_Q.
Kp_Q = 1;
Ki_Q = ((3*v_d_12*Kp_Q-2)/(3*v_d_12*tau));

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

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

% Limites para las señales de control en amperios.
% VSC1.
lim_i_q_1_upper = 393.72177;
lim_i_q_1_lower = -576.70686;

% VSC2.
lim_i_q_2_upper = 304.57446;
lim_i_q_2_lower = -446.13252;

% Set-point de Q_1 y Q_2 en dos instantes diferentes.
% VSC1.
Q_1_1 = 0;
Q_1_2 = -2500;

% VSC2.
Q_2_1 = -3000;
Q_2_2 = 1000;

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

% Ruido en la medida de los sensores.
% Ruido sensor Q_1.
D_R_Q_1 = 30e-3;                 % Desviación tipica ruido Q_1.
V_R_Q_1 = D_R_Q_1^2;             % Varianza ruido Q_1.

% Ruido sensor Q_2.
D_R_Q_2 = 30e-3;                 % Desviación tipica ruido Q_2.
V_R_Q_2 = D_R_Q_2^2;             % Varianza ruido Q_2.

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

% Desviacion respecto al valor nominal.
I_Delta_v_d_1 = 0.15;               % Instante de la desviación.
P_Delta_v_d_1 = 0.03;               % Porcentaje de desviación.
Delta_v_d_1 = P_Delta_v_d_1*v_d_12; % Desviación.

% Desviacion respecto al valor nominal.
I_Delta_v_d_2 = 0.15;               % Instante de la desviación.
P_Delta_v_d_2 = 0.03;               % Porcentaje de desviación.
Delta_v_d_2 = P_Delta_v_d_2*v_d_12; % Desviación.
                                                            % Prueba 2
%--------------------------------------------------------------------------






