compensador proporcional

'Controlador Porporcional'

num=1.2*[0.25 1];

den=[0.5 1.5 1.44 0.89];

Ts=0.1;

G=tf(num,den) %Gerar Funcaao de transferencia no plado S

BoG=c2d(G,Ts) % Gera a F.T G(s) para BoG(s) com o período de amostragem definido ?Ts?

FTMAr=BoG %Funcao de transferencia de malha aberta

zpk(BoG) %BoG na forma de polos e zeros

'Controlador Proporcional'

%Controlador Proporcional

rlocus(FTMAr) %Grafico do lugar das raizes

zeta=(-log(10/100) ) /( sqrt ( pi^2+(log(10/100)^2)) )

axis ('equal')

hold on

zgrid(zeta,0)

pause

[K,polos_mf]=rlocfind(FTMAr)

[theta,rho]=cart2pol(real(polos_mf(1)),imag(polos_mf(1))); % Realiza a transformação da coordenadas

%cartesianas em polares, e resulta num vetor de amplitude rho com ângulo de inclinação de theta (radianos)

'No Grafico selecione o ponto no circulo azul mais proximo da cobrinha'

theta

rho

deg = (rho *180) / pi

%Obter a função de transferência em malha fechada,

%com H(s) = 1 e com ganho proporcional encontrado no gráfico,e realimentação unitária:

FTMF = feedback(K*FTMAr,1)

dcgain(FTMF) % Encontra o ganho DC (estático) para FTMF

G_grau=10/ans % Determina o ganho para fixar a saída em regime permanente

FTMF = (G_grau*FTMF) %FTMF com ganho dc em regime permanente

figure(2)

step(FTMF)

disp('a função de transferencia da planta G(s)')

Gs=tf(num,den)

T=0.1;

disp('A função de transferencia discreta da planta G(s)')

Gz=c2d(Gs,T,'zoh')

disp('A função de transferencia de retorno')

Gr=feedback(Gz,1)

[numz,denz]=tfdata(Gr,'v');

n = length(denz)-1;

'Controlador Porporcional'

num=1.2*[0.25 1];

den=[0.5 1.5 1.44 0.89];

Ts=0.1;

G=tf(num,den) %Gerar Funcaao de transferencia no plado S

BoG=c2d(G,Ts) % Gera a F.T G(s) para BoG(s) com o período de amostragem definido ?Ts?

FTMAr=BoG %Funcao de transferencia de malha aberta

zpk(BoG) %BoG na forma de polos e zeros

'Controlador Proporcional'