v_pll_dsc0.m

This file is used to develop Phase locked loop.

% **************************************************************************
% * CONFIGURACION DEL pll_dsc.mdl
% *
% * Articulo empleado: 
% *  "Phase-locked loop for grid-connected three-phase power conversion
% *   system" S.-K.Chung
% *************************************************************************
% Inicializacion
clear all;
% close all;

Ts=200e-6;
% Em=400*sqrt(2/3);
Em=400;
epsilon=0.707;
% wn=2*pi*1000;
wn=2*pi*100; % Prueba hecha el 18-09-2004


% Estas constantes se emplean para el PLL continuo y PLL discreto 1. Se
% trata de un controlador discreto por emulaci髇
Kp=2*epsilon*wn/Em;
tau=Kp*Em/wn^2;
Ki=Kp/tau;
Finit=50;
Phase_Init=0;


wn1=2*pi*100;
epsilon=0.707;
% Estas constantes se emplean para el PLL discreto 2. Las entradas al PLL
% est醤 sin normalizar respecto a la amplitud de la tensi髇 de entrada.
% Transformaci髇 invariante en potencia, por eso se divide por Em. Esto
% no es muy correcto porque la respuesta es sensible a la amplitud
kpJ=(2*(1-exp(-epsilon*wn1*Ts)*cos(wn1*Ts*sqrt(1-epsilon^2))))/(Ts*Em);
alfaJ=(1-exp(-2*epsilon*wn1*Ts))/(2*(1-exp(-epsilon*wn1*Ts)*cos(wn1*Ts*sqrt(1-epsilon^2))));
kiJ=(1-alfaJ)*kpJ/Ts;

% Estas constantes se emplean para el PLL discreto 3 normalizado con amplitud. 
% Las entradas al PLL est醤 normalizadas con la amplitud de la tensi髇 de entrada.
% De esta forma el PLL es independiente de hacer una transformaci髇
% invariante en potencia o invariante en amplitud
kp_a=(2*(1-exp(-epsilon*wn1*Ts)*cos(wn1*Ts*sqrt(1-epsilon^2))))/(Ts*1);
alfa_a=(1-exp(-2*epsilon*wn1*Ts))/(2*(1-exp(-epsilon*wn1*Ts)*cos(wn1*Ts*sqrt(1-epsilon^2))));
ki_a=(1-alfa_a)*kp_a/Ts;