ATIVIDADE1 Implementação da Proteção Diferencial Percentual no Software Matlab

Prévia do material em texto

1 
 
UNIVERSIDADE FEDERAL DE SANTA MARIA 
CENTRO DE TECNOLOGIA 
PROGRAMA DE PÓS-GRADUAÇÃO EM 
ENGENHARIA ELÉTRICA 
 
 
 
 
 
 
 
 
IMPLEMENTAÇÃO NO MATLAB DA PROTEÇÃO 
DIFERENCIAL PERCENTUAL APLICADA A 
GERADORES 
 
 
 
 
ATIVIDADE 1 
 
 
 
 
 
José Oizimas Junior 
 
 
 
 
 
Santa Maria, RS 
2016 
2 
 
IMPLEMENTAÇÃO NO MATLAB DA PROTEÇÃO 
DIFERENCIAL PERCENTUAL APLICADA A 
GERADORES 
 
 
 
 
 
José Oizimas Junior 
 
 
 
 
 
Atividade 1 apresentada a disciplina de Proteção de Geradores em Sistemas 
Elétricos de Potência da Pós-Graduação em Engenharia Elétrica, Área de 
Concentração em Processamento de Energia, da Universidade Federal de 
Santa Maria (UFSM, RS), como requisito parcial para aprovação curricular. 
 
 
 
 
 Prof°. Dr. Adriano Peres de Morais 
 
 
 
Santa Maria, RS, Brasil 
2016 
3 
 
LISTA DE FIGURAS 
 
Figura 1: Componente da Corrente de Falta ............................................................................................................... 7 
Figura 2: Gerador de Pólos Salientes ........................................................................................................................... 8 
Figura 3: Ajuste de Proteção Percentual .................................................................................................................... 13 
Figura 4: Coordenograma do Relé Percentual .......................................................................................................... 18 
Figura 5: Falta Externa a Zona de Proteção Diferencial .......................................................................................... 19 
Figura 6: Falta Monofásica-Terra Interna a zona de Proteção ............................................................................... 20 
Figura 7: Falta Bifásica ................................................................................................................................................. 23 
Figura 8: Falta Trifásica ................................................................................................................................................ 23 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
4 
 
 
LISTA DE TABELAS 
 
Tabela 1: Falhas em Usinas Hidrelétricas ........................................................................... 6 
Tabela 2: Dados do Gerador ............................................................................................. 15 
 
 
 
 
 
5 
 
SUMÁRIO 
 
1. INTRODUÇÃO ........................................................................................................................................ 6 
2. FUNDAMENTAÇÃO TEÓRICA ........................................................................................................... 7 
2.1 Geradores Síncronos de Pólos Salientes .................................................................................. 7 
2.2 Cálculo das Correntes de Falta do Gerador Síncrono ............................................................ 8 
2.2.1 Definição das Componentes de Quadratura e Eixo Direto do Gerador ................................. 9 
2.3 Proteção Diferencial ...................................................................................................................... 12 
2.4 Proteção Diferencial Percentual ................................................................................................. 13 
2.4.1 Ajuste da Corrente Mínima de Operação ................................................................................. 14 
2.4.2 Ajuste do Slope Mínimo de Operação ....................................................................................... 14 
3. IMPLEMENTAÇÃO DA PROTEÇÃO DIFERENCIAL PERCENTUAL NO MATLAB .............. 15 
3.1 Dados do Gerador em Estudo ..................................................................................................... 15 
3.2 Definição da Corrente de Pick-Up .............................................................................................. 15 
3.3 Definição do Slope Mínimo .......................................................................................................... 16 
4. CONCLUSÃO ....................................................................................................................................... 19 
REFERÊNCIAS BIBLIOGRÁFICAS ......................................................................................................... 21 
ANEXO A – Valor RMS do Sinal de Falta .............................................................................................. 22 
ANEXO B – Falta Bifásica e Trifásica Interna ..................................................................................... 23 
ANEXO C – Filtro Fourier de Ciclo Completo...................................................................................... 24 
ANEXO D – Filtro do Sinal ..................................................................................................................... 26 
ANEXO E – Simulação da Proteção Diferencial Percentual ............................................................ 28 
 
 
 
 
 
 
 
 
 
 
 
 
 
6 
 
 
1. INTRODUÇÃO 
 
 
 A grande parcela de energia elétrica consumida no Brasil é gerada através de 
fontes hidráulicas instaladas em todo território nacional, nas mais variadas potência e 
características operacionais. Além da responsabilidade de suprir valores acima de 70% da 
demanda de energia elétrica, as fontes geradoras hídricas têm a incumbência de garantir 
a estabilidade e qualidade do Sistema Elétrico de Potência (SEP). 
 Analisando o panorama nacional, as fontes geradoras hídricas possuem papel 
protagonista perante o SEP. Desta forma, é fundamental detectar e mitigar as principais 
falhas da planta geradora que causa a inoperância técnica. De acordo com a tabela 1, a 
principal causa responsável pela paralização da produção de energia elétrica das 
hidrelétricas está associada ao gerador. 
 
 
 
Tabela 1: Falhas em Usinas Hidrelétricas 
 
Fonte : Kinderman, 2008 
 
 Gerador síncrono é considerado a parte substancial da planta de uma usina 
hidráulica, é responsável pela conversão da energia mecânica em energia elétrica. Desta 
forma, o trabalho tem por objetivo estudar a proteção diferencial percentual aplicada a 
geradores utilizando o software Matlab® com o intuito de garantir a detecção de falhas 
internas e a não ocorrência de atuações intempestivas para perturbações externas. 
 
 
 
 
7 
 
 
2. FUNDAMENTAÇÃO TEÓRICA 
 
A elaboração da proteção diferencial percentual aplicada a geradores inicia-se com a 
modelagem matemática de seus parâmetros. As grandezas elétricas obtidas na 
modelagem são fundamentais à parametrização do sistema de proteção, além da 
veracidade destes valores influência acintosamente na qualidade do ajuste do relé. 
Analisando o momento da falta, a corrente elétrica disponibilizada pelo gerador terá 
comportamento inicial assimétrico formado pelo somatório do componente de Corrente 
Alternado (CA) e Corrente Continua (DC), como demostrado na figura 1. 
 
Figura 1: Componente da Corrente de Falta 
 
Fonte: Reimert, 2006 
 
 Examinando apenas o sinal AC da perturbação, verifica-se uma redução da 
amplitude, tendendo a um valor simétrico com o passar do tempo, formando um envoltório 
decrescente, subdividido no período subtransitório, transitório e permanente. 
 
2.1 Geradores Síncronos de Pólos SalientesO estudo é direcionado aos geradores síncronos de pólos salientes em virtude da sua 
utilização perante o SEP. Este modelo de gerador possui característica própria ligada a 
desconformidade na periferia do núcleo de ferro, que por sua vez resulta em variações 
cíclicas no entreferro, tornado visível à saliência dos pólos. 
8 
 
A variação cíclica do entreferro resulta em alternâncias na relutância nos eixos diretos 
“d” e quadratura “q”, originando a reatância de eixo direto e quadratura distintas, 
representado pela figura 2. 
 
 
Figura 2: Gerador de Pólos Salientes 
 
Fonte: Reimert, 2006 
 
 
2.2 Cálculo das Correntes de Falta do Gerador Síncrono 
O modelo de estudo considera a falta trifásica como o caso mais severo em que o 
gerador pode ser submetido. Nesta perturbação, a corrente tende a assumir valores 
de pico característicos as seguintes reatâncias do gerador: 
 Reatância subtransitória de eixo direto; 
 Reatância transitória de eixo direto; 
 Reatância síncrona de eixo direto; 
 Reatância subtransitória de eixo de quadratura; 
 Reatância transitória de eixo de quadratura; 
 Reatância síncrona de quadratura. 
 
 
9 
 
2.2.1 Definição das Componentes de Quadratura e Eixo Direto do Gerador 
Para perturbação trifásica, as correntes de falta de eixo direto são calculadas da 
seguinte maneira: 
a) Decomposição das componentes de quadratura e eixo direto antes da falta: 
Estes valores tendem a variar de acordo com a corrente passante na 
armadura momento antes da falta. 
 
 (
 ( )
 ( )
) 
 
 ( ) 
 ( ) 
 ( ) 
 ( ) 
 
 
Onde: 
 
 : ângulo de carga; 
 : corrente de armadura antes da falta; 
 reatância de eixo de quadratura; 
 : ângulo o fator de potência do gerador; 
 : tensão no terminal do gerador; 
 tensão de eixo de quadratura antes da falta; 
 corrente de quadratura antes da falta; 
 tensão de eixo direto antes da falta. 
 
 
b) Corrente síncrona de eixo direto: Corrente em regime permanente após 
decorrido o regime subtransitório e transitório. 
 
 
 
 
Onde: 
 Corrente síncrona de eixo direto; 
 Tensão de linha nominal; 
 Reatância de eixo direto; 
 Reatância de falta do sistema. 
 
10 
 
c) Corrente transitório de eixo direto: corrente que circula na armadura após o 
período subtransitório do curto-circuito, perdurando por um número maior de 
ciclos. 
 
 
 
 
 
 
 
Onde: 
 Corrente transitória de eixo direto; 
 Tensão transitória de eixo de quadratura antes da falta; 
 Reatância transitória de eixo direto. 
 
d) Corrente subtransitória de eixo direto: corrente que circula na armadura 
durante os primeiros ciclos da falta. 
 
 
 
 
 
 
 
Onde: 
 Corrente subtransitória de eixo direto; 
 Tensão subtransitória de eixo de quadratura antes da falta; 
 Reatância transitória de eixo direto. 
 
 
e) Corrente síncrona de eixo de quadratura: Corrente em regime permanente 
após decorrido o regime subtransitório e transitório. 
 
 
 
 
 
Onde: 
 Corrente síncrona de eixo de quadratura; 
 Reatância de eixo de quadratura. 
 
 
11 
 
f) Corrente transitório de eixo de quadratura: corrente que circula na armadura 
após o período subtransitório do curto, perdurando por um número maior de 
ciclos. 
 
 
 
 
 
 
 
 
Onde: 
 Corrente transitória de eixo de quadratura; 
 Tensão transitória de eixo direto antes da falta; 
 Reatância transitória de eixo de quadratura. 
 
g) Corrente subtransitória de eixo de quadratura: corrente que circula na 
armadura durante os primeiros ciclos. 
 
 
 
 
 
 
 
Onde: 
 Corrente subtransitória de eixo de quadratura; 
 Tensão subtransitória de eixo direto antes da falta; 
 Reatância transitória de eixo de quadratura; 
 
h) Definição dos períodos de quadratura e eixo direto no momento da falta: As 
reatâncias do gerador somado a impedância equivalente do sistema até o 
ponto da falta possibilitam a definição das constantes de tempo de cada 
etapa da perturbação. 
 
 
 (
 
 
) 
 (
 
 
) 
12 
 
 (
 
 
) 
 (
 
 
) 
 (
 
 
) 
 
 
 
Onde: 
 
 : Período transitório de eixo direto associado á falta; 
 : Período subtransitório de eixo direto associado á falta; 
 : Período transitório de eixo de quadratura associado á falta; 
 : Período subtransitório de eixo de quadratura associado á falta; 
 : Período transitório de eixo direto com circuito aberto; 
 : Período subtransitório de eixo direto com circuito aberto; 
 : Período transitório de eixo de quadratura com circuito aberto; 
 : Período subtransitório de eixo de quadratura com circuito aberto. 
 
 
i) Corrente RMS máxima de falta: O valor máximo da corrente de falta 
assimétrica referente ao gerador síncrono está associado ao componente 
AC e DC. 
 
 √ 
 
 ( ) √ 
 √ * ( ) 
 
 
Onde: 
 
 : Corrente máxima RMS; 
 ( ): Corrente alternada máxima RMS; 
 : Corrente máxima DC. 
 
2.3 Proteção Diferencial 
 
A proteção diferencial empregada nos geradores tem a finalidade de proteger contra 
curto-circuito produzido no enrolamento do estator. O sistema de proteção é similar ao 
13 
 
aplicado em transformadores (TR’s), porém com uma vantagem relacionada aos 
transformadores de corrente (TC’s), não necessitando utilizar equipamentos com relação 
de transformação (RTC) diferente, garantindo maior precisão comparada ao sistema 
aplicado aos TR’s. 
Este sistema de proteção é eficaz para correntes de curto-circuito elevadas, 
geralmente produzidas por C-C trifásico, C-C bifásico e C-C monofásico com a carcaça 
onde o sistema de aterramento do neutro é sólido. Já o sistema de aterramento do neutro 
com alta impedância tem característica de limitar as correntes de Curto- Circuito, fazendo 
com que algumas falhas monofásicas não sensibilizam a proteção diferencial. 
Outro fator importante relacionado a proteção diferencial está associado a geradores 
de grande porte, onde utilizam mais de uma espira por fase. Para está situação, a 
proteção não detecta falhas entre espiras da mesma fase, em função de não ocorrer 
divergência de valores entre os TC’s. 
 
 
2.4 Proteção Diferencial Percentual 
 
 A proteção diferencial percentual consiste no monitoramento de duas variáveis, 
denominadas corrente de operação (Iop) e corrente de restrição (IRT). O sistema de 
proteção irá atuar quando o valor do torque produzido pelo Iop ultrapassar o valor do 
torque gerado pela IRT, de acordo com o percentual de inclinação da reta de ajuste. A 
figura 3 demostra o ajuste de proteção percentual qualquer. 
 
 
Figura 3: Ajuste de Proteção Percentual 
 
Fonte: Reimert, 2006 
14 
 
2.4.1 Ajuste da Corrente Mínima de Operação 
 
O ajuste da corrente mínima de operação deve considerar o erro máximo dos 
transformadores de corrente, juntamente com o erro máximo do relé para corrente 
nominal do sistema. Desta maneira evitará atuações indesejadas da proteção quando o 
sistema estiver operando normalmente. 
 A definição da corrente mínima de operação ocorre através da seguinte equação:( (
 
 
) – ( (
 
 
)] 
 
Onde: 
 : corrente de pick-up do relé; 
 : corrente nominal do secundário do transformador de corrente; 
 : percentual de erro do transformador de corrente; 
 : percentual de erro do relé. 
 
 
2.4.2 Ajuste do Slope Mínimo de Operação 
 
O ajuste mínimo do slope deve prever o maior nível de corrente de curto-circuito fora 
da zona de proteção diferencial, deve-se levar em conta o somatório dos erros máximo 
dos TC’s e relé, a fim de garantir a não operação do sistema de proteção para falta 
externas. 
A condição mínima de operação deve atender a seguinte condição: 
 
 (
 
 
) 
 
| | | |
 
 
 
Onde: 
 : percentual mínima de inclinação da reta de operação; 
 : fator de segurança (geralmente 2); 
 : corrente de operação; 
 : corrente de restrição; 
 
 
 
 
 
 
 
 
15 
 
3. IMPLEMENTAÇÃO DA PROTEÇÃO DIFERENCIAL PERCENTUAL NO MATLAB 
 
O objetivo da atividade é implementar a proteção diferencial percentual no software 
Matlab® para um determinado gerador e testar algumas possibilidades de faltas em 
que o sistema pode ser submetido, a fim de validar o ajuste proposto. 
 
3.1 Dados do Gerador em estudo 
 
Tabela 2: Dados do Gerador 
Type of Rotor Salient Pole 
Frequency (Hz) 60 
Number of poles 48 
Rs (pu) 0.002636 
Xls (pu) 0.015 
Xd (pu) 0.985 
Xq (pu) 0.609 
X'd (pu) 0.340 
X'q (pu) 0.60 
X''d (pu) 0.249 
X''q (pu) 0.272 
Td0' 7.348 
Tq0' N/A 
Td0'' 0.075 
Tq0'' 0.14 
Sn 370 MVA 
FP 0.85 
Vn 13.8 KV 
Erro TC ±10% 
Erro Relé ±2% 
Fonte: Autoria Própria 
 
 No desenvolvimento dos cálculos foram utilizados valores de base iguais aos 
valores nominais do gerador. 
 
3.2 Definição da Corrente de Pick-Up 
 
 A definição da corrente de operação mínima será em função da capacidade 
nominal do gerador influenciada pelo erro máximos do TC’s e relé. 
 
 Corrente nominal do gerador 
 
 
 
√ 
 
 
 
16 
 
 Relação de transformação: 
 
 
 
 
 
 
 
 Corrente de operação nominal refletida no secundário do TC 
 
 
 
 
 
 Corrente de operação mínima 
 
 ( (
 
 
) – ( (
 
 
)] 
 
 A faixa de ajuste da corrente de pick-up do relé é de (0,5 – 5,0) A em degraus de 
0.5 A. Portanto o valor mínimo de operação é: 
 
 
 
 Para este caso valor da corrente mínima de operação assume valor constante de 
1A na faixa de (0 – 4,84) A no eixo representado pela corrente de restrição. 
 
 
3.3 Definição do Slope Mínimo 
 A característica dos transformadores de corrente é aumentar o erro de forma 
gradativa com o incremento da corrente em patamares que ultrapassa a capacidade 
nominal dos TC. Desta forma é necessário o acréscimo na corrente de operação mínima, 
na faixa superior a capacidade nominal do TC a fim de evitar atuações na proteção 
diferencial quando ocorrer falta externa a zona de atuação. 
 
 Corrente máxima de curto-circuito nos terminais do gerador 
 
 (
 ( )
 ( )
) 
 
 ( ) 
 ( ) 
 ( ) 
 ( ) 
 
 
17 
 
 
 
 
 
 
 
 
 
 ( ) √ 
 √ * ( ) 
 
 √ 
 
 IRMS é o maior valor de corrente de curto circuito que eventualmente poderá ocorrer 
fora da zona protegida pela proteção diferencial percentual. 
 
 % de Slope: A faixa de ajuste do slope é de (10 – 100)% em degraus de 5%. 
 
 
| | | |
 
 
 (
 
 
) 
 
 
 
 Função que representa o Slope 
 
{
 
 
 
 ( ) 
 ( ) {
 
 
 
Para esta situação, onde se desconhece a impedância do sistema, não haverá ajustes 
intermediários de Slope. A precisão de ajustes intermediários é dependente dos níveis de 
curto-circuito em pontos distantes do gerador, onde os níveis de C-C são relativamente 
inferior, comparados as corretes de curto-circuito nos terminais do gerador. A figura 4 
representa o ajuste para este gerador desconsiderando a impedância do sistema. 
18 
 
Figura 4: Coordenograma do Relé Percentual 
 
Fonte: Autoria Própria 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
19 
 
4. CONCLUSÃO 
 O trabalho descreveu o processo de implementação do ajuste da proteção 
diferencial percentual aplicado a geradores utilizando o software Matlab®. Através das 
simulações foi comprovado que o sistema, quando ajustado adequadamente é 
extremamente eficaz para faltas de grande magnitude (faltas bifásicas e trifásicas). 
 O estudo teve por objetivo testar sinais de corrente proveniente de faltas dentro da 
zona de atuação e externa a área de atuação do sistema de proteção. Os sinais de falta 
externa de maiores magnitudes propicia a validação da proteção relacionada ao erro 
máximo dos equipamentos, para estas perturbações o sistema de proteção não deve 
atuar, tal situação é comprovada pela figura 5, onde o sistema de proteção foi submetido 
a perturbação externa e comprovou a não sensibilização do relé. 
 
 
Figura 5: Falta Externa a Zona de Proteção Diferencial 
 
Fonte: Autoria Própria 
 
 Referente as faltas internas, um sistema ideal deveria detectar todos os tipos de 
falhas que podem vir a ocorrer na zona de proteção. Em função dos erros dos 
equipamentos isto não ocorre, dado o fato de que o sistema deve ser ajustado para o 
máximo erro em operação normal, este ajuste na grande maioria das vezes é superior aos 
20 
 
valores de corrente de curto circuito monofásico-terra, desta forma, a magnitude destas 
perturbações não sensibiliza a proteção diferencial percentual. 
 A figura 6 representa a situação onde ocorreu a falha interna de baixa magnitude, 
tendo por característica a não atuação do sistema de proteção, evidenciando a 
inoperância do sistema para correntes de curto-circuito de baixo valor. 
 
Figura 6: Falta Monofásica-Terra Interna a zona de Proteção 
 
Fonte: Autoria Própria 
 
 De base nos argumentos descrito, o sistema de proteção se mostra muito eficaz 
para perturbações que produzem elevado nível de corrente de curto-circuito, em 
contraponto, é um sistema inábil para faltas de baixa intensidade, necessitando sistema 
de proteção auxiliares para detecção destes eventos. 
 
 
 
 
21 
 
REFERÊNCIAS BIBLIOGRÁFICAS 
 
KINDERMANN, G. PROTEÇÃO DE SISTEMAS ELÉTRICOS DE POTÊNCIA. 
Florianópolis: Edição do autor, 2008. 
REIMERT, D. PROTECTIVE RELAYING FOR POWER GENERATION SYSTEMS. Boca 
Raton: CRC Press, 2006. 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
22 
 
ANEXO A – Valor RMS do Sinal de Falta 
 
 
 
 
 
 
23 
 
ANEXO B – Falta Bifásica e Trifásica Interna 
Figura 7: Falta Bifásica 
 
 
 
Figura 8: Falta Trifásica 
 
24 
 
ANEXO C – Filtro Fourier de Ciclo Completo 
 
function [Yc,Ys,soma1,mod,ang,angfinal]=fourier1(yq,N,ts,yfref) 
 
%% FILTRO FOURIER DE CICLO COMPLETO %% 
 
%% Novo yq (com zeros na janela inicial) 
N=96; 
jan=zeros((N-1),1); %% zeros para completar a janelayq=[jan;yq]; 
ref=[jan;yfref]; 
 
%% Coeficientes do filtro 
j=1; 
for k=1:N 
 Cfc(1,j)=(sqrt(2)/N)*cos(((k-1)*2*pi)/N); 
 Cfs(1,j)=-(sqrt(2)/N)*sin(((k-1)*2*pi)/N); 
 j=j+1; 
end 
 
%% Cálculos 
 
for k=1:length(yq)-N+1 
 janela(:,1)=yq(k:k-1+N); 
 janela2(:,1)=ref(k:k-1+N); 
 Yc(k,1)= Cfc*janela(1:N); 
 Ys(k,1)= Cfs*janela(1:N); 
 Yc1(k,1)= Cfc*janela2(1:N); 
 Ys1(k,1)= Cfs*janela2(1:N); 
end 
 
soma1=Yc+Ys; 
mod=abs(Yc+i*Ys); 
ang=(180/pi)*angle(Yc+i*Ys); 
 
modref=abs(Yc1+i*Ys1); 
angref=(180/pi)*angle(Yc1+i*Ys1); 
 
angfinal=ang-angref; 
 
for k=1:length(angfinal) 
 if angfinal(k)>=180 
 angfinal(k)=angfinal(k)-360; 
 end 
 if angfinal(k)<=-180 
 angfinal(k)=angfinal(k)+360; 
 end 
 if abs(angfinal(k))==angfinal(k) & angfinal(k)<=1e-4 
 angfinal(k)=0; 
 end 
 if abs(angfinal(k))~=angfinal(k) & angfinal(k)>=-1e-4 
 angfinal(k)=0; 
 end 
end 
a2 =1+1; 
 
% Gráfico do sinal de entrada no filtro + saída (módulo e ângulo) 
 
figure(1) 
25 
 
plot(ts,yq(N:length(yq))) 
hold on 
plot(ts,mod,'r') 
 
h = legend('Sinal quantizado','Sinal após DFT',1); 
title('Módulo'); 
grid on 
 
 
figure(4) 
tempo=[1:size(mod,1)]; 
plot(mod) 
 
hold on 
plot(ts,yq,'r') 
h = legend('Saída','Entrada',1); 
title('Módulo - Sinal de entrada e saída do filtro Fourier'); 
grid on 
 
figure(5) 
plot(ts,ang) 
hold on 
plot(ts,yq*20,'r') 
h = legend('Saída','Entrada*20',1); 
title('Fase - Sinal de entrada e saída do filtro Fourier'); 
grid on 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
26 
 
 
ANEXO D – Filtro do Sinal 
function [Val_Fas] = filtros (Input_Signals,fn,fm) 
 
% if flag==0 
 %% Dados %% 
 
 fc = 480; %frequencia de corte do filtro Butterworth 
% Q=sqrt(2)/2; %fator de qualidade do filtro Butterworth 
 %fn = 50; %frequencia fundamental 
 
 N = fm/fn; %número de amostras por ciclo ... En DS se le dió 1200Hz 
1200/60 
 bits = 8; %número de bits para a quantização 
 N1 = 300; %número de amostras para plotar a FFT 
 ord = 2; %ordem do filtro Butterworth 
 m = 49; %número da harmônica que se deseja obter com o FCDFT 
 t = Input_Signals(:,1); 
 step = t(5)-t(4); 
 ts = step; 
 
 %% Filtros 
 k = 0; 
 for k = 2:7 % 12 variáveis (Ia,b,c; Va,b,c) lado m e lado n 
 y = Input_Signals(:,k); 
 [yref]=referencia(y,t,fn); 
 [yfilt,coefb,coefa]=prefiltragem1(y,fc,ord,t); 
% [yfref,coefb1,coefa1]=prefiltragem1(yref,fc,ord,fs,t); 
 [Ycfor,Ysfor,soma,mod1,angulo,angfinal]=fourier1(yfilt,N,t,yref); % 
filtro fourier; 
 
 if k==2 
 Ian_mod_rms = mod1; 
 Ian_ang_rms1 = angfinal; 
 Ian_real = Ian_mod_rms.*cos(Ian_ang_rms1.*pi/180); 
 Ian_imag = Ian_mod_rms.*sin(Ian_ang_rms1.*pi/180); 
 %plot (Ian_mod_rms); 
 end 
 if k==3 
 Iat_mod_rms = mod1; 
 Iat_ang_rms1 = angfinal; 
 Iat_real = Iat_mod_rms.*cos(Iat_ang_rms1.*pi/180); 
 Iat_imag = Iat_mod_rms.*sin(Iat_ang_rms1.*pi/180); 
 %plot (IAn_mod_rms); 
 end 
 if k==4 
 Ibn_mod_rms = mod1; 
 Ibn_ang_rms1 = angfinal; 
 Ibn_real = Ibn_mod_rms.*cos(Ibn_ang_rms1.*pi/180); 
 Ibn_imag = Ibn_mod_rms.*sin(Ibn_ang_rms1.*pi/180); 
 
 end 
 if k==5 
 Ibt_mod_rms = mod1; 
 Ibt_ang_rms1 = angfinal; 
 Ibt_real = Ibt_mod_rms.*cos(Ibt_ang_rms1.*pi/180); 
 Ibt_imag = Ibt_mod_rms.*sin(Ibt_ang_rms1.*pi/180); 
 
27 
 
 end 
 if k==6 
 Icn_mod_rms = mod1; 
 Icn_ang_rms1 = angfinal; 
 Icn_real = Icn_mod_rms.*cos(Icn_ang_rms1.*pi/180); 
 Icn_imag = Icn_mod_rms.*sin(Icn_ang_rms1.*pi/180); 
 end 
 if k==7 
 Ict_mod_rms = mod1; 
 Ict_ang_rms1 = angfinal; 
 Ict_real = Ict_mod_rms.*cos(Ict_ang_rms1.*pi/180); 
 Ict_imag = Ict_mod_rms.*sin(Ict_ang_rms1.*pi/180); 
 
 end 
 
 end 
 
 Ianf = complex(Ian_real,Ian_imag); % forma complexa de IAm 
 Ibnf = complex(Ibn_real,Ibn_imag); % forma complexa de IBm 
 Icnf = complex(Icn_real,Icn_imag); % forma complexa de ICm 
 
 Iatf = complex(Iat_real,Iat_imag); % forma complexa de IAn 
 Ibtf = complex(Ibt_real,Ibt_imag); % forma complexa de IBn 
 Ictf = complex(Ict_real,Ict_imag); % forma complexa de ICn 
 
 
 
 Val_Fas=[Ianf Iatf Ibnf Ibtf Icnf Ictf]; 
end 
 
28 
 
ANEXO E – Simulação da Proteção Diferencial Percentual 
 
clc 
clear all 
% Implementação Relé Diferencial 
%% ///////////////////////////Dados de Entrada////////////////////////////// 
fn=60; 
fm=96*fn; 
Sn=370e+6; %Potência aparente nominal do gerador 
Vl=13.8e+3; % Tensão de linha do gerador 
n=48; % Número de pólos 
f=60; % Frequência do Gerador em HZ 
%Rs_pu=0.002636; % Resistência do sistema em pu 
%Xls_pu=0.015; % Reatância do sistema em pu 
Xd_pu=0.985; % Reatância de eixo direto pu 
Xq_pu=0.609; % Reatância de quadratura em pu 
X1d_pu=0.34; % Reatância transitória de eixo direto em pu 
X1q_pu=0.60; % Reatância transitória de quadratura em pu 
X2d_pu=0.249; % Reatância subtransitória de eixo direto em pu 
X2q_pu=0.272; % Reatância subtransitória de quadratura em pu 
T1do=7.348; % Período transitório de eixo direto em segundos 
T2do=0.075; % Período subtransitório de eixo direto em segundos 
T2qo=0.14; % Período subtransitório de quadratura em segundos 
fp=0.85; % Fator de potência do gerador 
Rs_pu=0;%input ('Entre com a resistência do sistema:'); 
Xls_pu=0;%input ('Entre com a reatância do sistema:'); 
Iln_pu=1; 
Vl_pu=1; 
Iln=Sn/(Vl*sqrt(3)); 
RTC_n=Iln/5; 
%input ([ 'O valor do RTC de carga maxíma real nominal é:', 
num2str(RTC_n),'A']); 
RTC=3200;%input(' Entre com a relação de correte do transformador comercial:'); 
I_base_P=Iln; % Corrente de base no Primário do TC 
I_base_S=I_base_P/RTC;% Corrente de base no secundário do TC 
 
%% Ajuste da corrente mínima de operação 
 
TC_erro_max =10;%input ('Entre com o valor do erro maxímo do TC:'); %Erro maxímo 
indidual do Tc 
RE_erro_max=2;% input ('Entre com o valor do erro maxímo do Relé:'); %Erro 
maxímo indidual do Relé 
I_tc1=Iln/RTC; 
I_tc2=Iln/RTC; 
 
Iop=[[I_tc1*(1+(TC_erro_max/100))]-[I_tc2*(1-
(TC_erro_max/100))]*(1+(RE_erro_max/100))]; 
 
 % Faixa de ajuste do pickup do relé Iop=(0,5 - 5,0)no secundário em degraus 
de 0.5 A. 
 
if ((Iop>0.01) && (Iop<=0.75)) % Ver valores maiores que 5 
 a=0.5; 
 %input ([ 'O valor da corrente minima de operação é:', num2str(a),'A']); 
elseif ((Iop>0.75) && (Iop<=1.25)) 
 a=1; 
 % input ([ 'O valor da corrente minima de operação é:', 
num2str(a),'A']); 
elseif ((Iop>1.25) && (Iop<=1.75)) 
29 
 
 a=1.5; 
 %input ([ 'O valor da corrente minima de operação é:', num2str(a),'A']); 
elseif((Iop>1.75) && (Iop<=2.25)) 
 a=2; 
 input ([ 'O valor da corrente minima de operação é:', num2str(a),'A']); 
elseif ((Iop>2.25) && (Iop<=2.75)) 
 a=2.5; 
 input ([ 'O valor da corrente minima de operação é:', num2str(a),'A']); 
elseif ((Iop>2.5) && (Iop<=2.75)) 
 a=2.5; 
 input ([ 'O valor da corrente minima de operação é:', num2str(a),'A']); 
elseif ((Iop>2.5) && (Iop<=3.25)) 
 a=3; 
 input ([ 'O valor da corrente minima de operação é:',num2str(a),'A']); 
elseif ((Iop>3.25) && (Iop<=3.75)) 
 a=3.5; 
 input ([ 'O valor da corrente minima de operação é:', num2str(a),'A']); 
elseif ((Iop>3.75) && (Iop<=4.25)) 
 a=4; 
 input ([ 'O valor da corrente minima de operação é:', num2str(a),'A']); 
elseif ((Iop>4.25) && (Iop<=4.75)) 
 a=4.5; 
 input ([ 'O valor da corrente minima de operação é:', num2str(a),'A']); 
elseif ((Iop>4.75) && (Iop<=5)) 
 a=5; 
 input ([ 'O valor da corrente minima de operação é:', num2str(a),'A']); 
elseif (Iop>5); 
end 
 
% Faixa nominal da corrente de restrição 
Irt_o= 0:0.01:I_tc1; 
Iop=a; 
%% Ajuste do Primeiro Slope 
 
ang_fp=acos(fp); %angulo do fator de potência; 
a=Iln_pu*Xq_pu*fp; 
b=Vl_pu + (Iln_pu*Xq_pu*sin(acos(fp))); 
ang_carga=atan(a/b); 
Eqo_pu= Vl_pu*cos(atan(ang_carga)); % Tensão de quadratura no instante que 
antecede a falta e p.u. 
Iqo_pu= Iln_pu*cos(ang_fp + ang_carga); 
Edo_pu= Vl_pu*sin(ang_carga); 
Ido_pu= Iln_pu*sin(ang_fp + ang_carga); 
 
% Tensões de Quadratura 
E2qo_pu= Eqo_pu + Ido_pu*X2d_pu; % Tensão subtransitória em p.u. 
E1qo_pu= Eqo_pu + Ido_pu*X1d_pu; % Tensão transitória em p.u. 
El_pu = Eqo_pu + Ido_pu*Xd_pu; % Tensão em regime de sincronismo em 
p.u. 
 
% Tensão de Eixo Direto 
E2do_pu= Edo_pu - Iqo_pu*X2d_pu; % Tensão subtransitória em p.u. 
E1do_pu= Edo_pu - Iqo_pu*X1d_pu; % Tensão transitória em p.u. 
 
% Constante de tempo considerando a reatância de falta 
T2df= T2do *[(X2d_pu + Xls_pu)/(X1d_pu + Xls_pu)]; 
T1df= T1do *[(X1d_pu + Xls_pu)/(Xd_pu + Xls_pu)]; 
T2qf= T2qo *[(X2q_pu + Xls_pu)/(X1q_pu + Xls_pu)]; 
%T1qf= T1qo *[(X1q_pu + Xls_pu)/(Xq_pu + Xls_pu)]; 
 
30 
 
% Correntes de eixo direto 
Id_pu= El_pu/(Xd_pu + Xls_pu); % Sincrona 
I1d= E1qo_pu/(X1d_pu + Xls_pu); % Transitória 
I2d_pu= E2qo_pu/(X2d_pu + Xls_pu); % Subtransitória 
 
% Correntes de Quadratura 
%Iq= Ed/(Xq+Xf); % Síncrona 
I1q_pu= E1do_pu/(X1q_pu + Xls_pu); % Transitória 
I2q_pu= E2do_pu/(X2q_pu + Xls_pu); %Subtransitória 
 
% Corrente RMS Assimétrica 
Iac_rms1_pu= sqrt(I2d_pu^2 + I2q_pu^2); 
Idc_pu= sqrt(2)*Iac_rms1_pu; 
I_rms1_pu= sqrt(Iac_rms1_pu^2 + Idc_pu^2); 
I_rms1=I_rms1_pu*I_base_S; 
%input(['O valor da corrente RMS é:', num2str(I_rms1), 'A']); 
 
TC_erro_slope1 =10;%input ('Entre com o valor de erro do TC (metade do erro 
maximo:'); %Erro 50% do valor maxímo individual do TC 
 
Iop_slope1=[[I_rms1*(1+(TC_erro_slope1/100))]-[I_rms1*(1-
(TC_erro_slope1/100))]*(1+(RE_erro_max/100))]; 
%input(['O valor da corrente de operção slope 1 é:', num2str(Iop_slope1), 
'A']); 
 
Irt1_max= ((I_rms1*2)+ Iop_slope1)/2; 
 
%Definição da equação do slope1////////Considera a impêdancia do sistema// 
 
a=(Iop-Iop_slope1)/(I_tc1-I_rms1); 
b= Iop -I_tc1*[(Iop - Iop_slope1)/(I_tc1-I_rms1)]; 
 
Irt_1= I_tc1:0.01:I_rms1; 
y= a*Irt_1 + b; 
K=2; 
min_slope_per=K*100*(Iop_slope1/I_rms1); 
%input(['O valor do mínimo slope 1 é:', num2str(min_slope_per), '%']); 
 
 
%% Ajuste do slope 2 
Eqo_pu= Vl_pu*cos(ang_carga); % Tensão de quadratura no instante que antecede 
a falta e p.u. 
Iqo_pu= Iln_pu*cos(ang_fp + ang_carga); 
Edo_pu= Vl_pu*sin(ang_carga); 
Ido_pu= Iln_pu*sin(ang_fp + ang_carga); 
 
% Tensões de Quadratura 
E2qo_pu= Eqo_pu + Ido_pu*X2d_pu; % Tensão subtransitória em p.u. 
E1qo_pu= Eqo_pu + Ido_pu*X1d_pu; % Tensão transitória em p.u. 
El_pu = Eqo_pu + Ido_pu*Xd_pu; % Tensão em regime de sincronismo em 
p.u. 
 
% Tensão de Eixo Direto 
E2do_pu= Edo_pu - Iqo_pu*X2d_pu; % Tensão subtransitória em p.u. 
E1do_pu= Edo_pu - Iqo_pu*X1d_pu; % Tensão transitória em p.u. 
 
% Constante de tempo considerando a reatância de falta 
T2df_2= T2do *(X2d_pu / X1d_pu); 
T1df_2= T1do *(X1d_pu / Xd_pu); 
31 
 
T2qf_2= T2qo *(X2q_pu / X1q_pu); 
%T1qf_2= T1qo *(X1q_pu/Xq_pu); 
 
% Correntes de eixo direto 
Id_2_pu= El_pu/Xd_pu; % Sincrona 
I1d_2_pu= E1qo_pu/X1d_pu; % Transitória 
I2d_2_pu= E2qo_pu/X2d_pu; % Subtransitória 
 
% Correntes de Quadratura 
%Iq= Ed/(Xq+Xf); % Síncrona 
I1q_2_pu= E1do_pu/X1q_pu; % Transitória 
I2q_2_pu= E2do_pu/X2q_pu; %Subtransitória 
 
% Corrente RMS Assimétrica 
Iac_rms2_pu= sqrt(I2d_2_pu^2 + I2q_2_pu^2); 
Idc_2_pu= sqrt(2)*Iac_rms2_pu; 
I_rms2_pu= sqrt(Iac_rms2_pu^2 + Idc_2_pu^2); 
I_rms2=I_rms2_pu*I_base_S; 
%input(['O valor da corrente RMS de curto-circuito maxímo é:', num2str(I_rms2), 
'A']); 
 
TC_erro_slope2 =10;%input ('Entre com o valor de erro do TC:'); %Erro 100% do 
valor maxímo individual do TC 
 
Iop_slope2=[[I_rms2*(1+(TC_erro_slope2/100))]-[I_rms2*(1-
(TC_erro_slope2/100))]*(1+(RE_erro_max/100))]; 
%input(['O valor da corrente de operção slope 2 para curto-circuito maxímo é:', 
num2str(Iop_slope2), 'A']); 
 
%Irt2_max= ((I_rms2*2)+ Iop_slope2)/2; 
 
%% Definição da equação do slope2 
c= (Iop_slope1-Iop_slope2)/(I_rms1-I_rms2); 
d= (Iop_slope1)-(I_rms1*c); 
%(Iop_slope1-Iop_slope2)/(I_rms1-I_rms2)) 
Irt_2=I_rms1:0.01:I_rms2; 
y2= c*Irt_2 + d; 
K=2; 
min_slope2_per=K*100*(Iop_slope2/I_rms2); 
 
%% Pontos de falta 
 disp('Codigos das faltas no gerador'); 
 disp('O codigo da falta F1_P15_75_1G é:1'); 
 disp('O codigo da falta F1_P15_99_1G é:2'); 
 disp('O codigo da falta F1_P170_5_1G é:3'); 
 disp('O codigo da falta F2_P170_E_BT é:4'); 
 disp('O codigo da falta F2_P330_1_70_1G é:5'); 
 disp('O codigo da falta F3_P15_1_1_1G é:6'); 
 disp('O codigo da falta F3_P15_E_BT é:7'); 
 disp('O codigo da falta F3_P170_1_95_1G é:8'); 
 disp('O codigo da falta F4_P15_90_70_80_1G é:9'); 
 disp('O codigo da falta F4_P15_E_AT é:9'); 
 disp('O codigo da falta F4_P170_1_5_3_1G é:10'); 
 disp('O codigo da falta F4_P170_95_1_1_1G é:11'); 
 disp('O codigo da falta F4_P330_90_70_80_1G é:12'); 
 disp('O codigo da falta F4_P15_90_70_80_1G é:13'); 
 disp('O codigo da falta F4_P15_E_AT é:14'); 
 disp('O codigo da falta F4_P170_1_5_3_1G é:15'); 
 disp('O codigo da falta F4_P170_95_1_1_1G é:16'); 
 disp('O codigo da falta F4_P330_90_70_80_1G é:17'); 
32 
 
 p_falta=input('Entre com o codigo da falta:'); 
 
 
% Aplicar o filtro de fourier e o usar o valor RMS para testar os sinais no 
% ajuste da proteção 
TC_erro=5; 
switch (p_falta); 
 %% Falta fase-terra 
case 1 
 load F1_P15_75_1G; 
 Ian=awgn([Gen_B_IF_Ian].*(1+(TC_erro/100)),45); 
 Iat=awgn([Gen_B_IF_Iat].*(1-(TC_erro/100)),45); 
 Ibn=awgn([Gen_B_IF_Ibn].*(1+(TC_erro/100)),45); 
 Ibt=awgn([Gen_B_IF_Ibt].*(1-(TC_erro/100)),45); 
 Icn=awgn([Gen_B_IF_Icn].*(1+(TC_erro/100)),45); 
 Ict=awgn([Gen_B_IF_Ict].*(1-(TC_erro/100)),45); 
[Inputs]=[x Ian/RTC Iat/RTC Ibn/RTC Ibt/RTC Icn/RTC Ict/RTC]; 
[Val_fas]=filtros(Inputs,fn,fm); 
IopA = abs (Val_fas(:,1)+Val_fas(:,2)); 
 IopB = abs (Val_fas(:,3)+Val_fas(:,4)); 
 IopC = abs (Val_fas(:,5)+Val_fas(:,6)); 
 IrtA = (abs (Val_fas(:,1))+abs(Val_fas(:,2)))/2; 
 IrtB = (abs (Val_fas(:,3))+abs(Val_fas(:,4)))/2; 
 IrtC = (abs (Val_fas(:,5))+abs(Val_fas(:,6)))/2; 
plot(IrtA,IopA,'k*','LineWidth',2) 
hold on 
plot(IrtB,IopB,'m+','LineWidth',2) 
hold on 
plot(IrtC,IopC,'yo','LineWidth',2) 
legend ('Trip fase A','Trip fase B','Trip fase C') 
plot(Irt_o,Iop,'b','LineWidth',2); 
ylim([-1,10]); 
hold on 
plot(Irt_1,y,'k','LineWidth',2); 
ylabel('Corrente de Operação'); 
title('Coordenograma do Relé Diferencial Percentual'); 
xlabel('Corrente de Restrição'); 
hold on 
plot(Irt_2,y2,'b','LineWidth',2); 
hold off 
%% Falta fase-terra 
 
case 2 
 load F1_P15_99_1G; 
 Ian=awgn([Gen_B_IF_Ian].*(1+(TC_erro/100)),45); 
 Iat=awgn([Gen_B_IF_Iat].*(1-(TC_erro/100)),45); 
 Ibn=awgn([Gen_B_IF_Ibn].*(1+(TC_erro/100)),45); 
 Ibt=awgn([Gen_B_IF_Ibt].*(1-(TC_erro/100)),45);Icn=awgn([Gen_B_IF_Icn].*(1+(TC_erro/100)),45); 
 Ict=awgn([Gen_B_IF_Ict].*(1-(TC_erro/100)),45); 
[Inputs]=[x Ian/RTC Iat/RTC Ibn/RTC Ibt/RTC Icn/RTC Ict/RTC]; 
[Val_fas]=filtros(Inputs,fn,fm); 
IopA = abs (Val_fas(:,1)+Val_fas(:,2)); 
 IopB = abs (Val_fas(:,3)+Val_fas(:,4)); 
 IopC = abs (Val_fas(:,5)+Val_fas(:,6)); 
 IrtA = (abs (Val_fas(:,1))+abs(Val_fas(:,2)))/2; 
 IrtB = (abs (Val_fas(:,3))+abs(Val_fas(:,4)))/2; 
 IrtC = (abs (Val_fas(:,5))+abs(Val_fas(:,6)))/2; 
plot(IrtA,IopA,'k*','LineWidth',2) 
hold on 
plot(IrtB,IopB,'m+','LineWidth',2) 
33 
 
hold on 
plot(IrtC,IopC,'yo','LineWidth',2) 
legend ('Trip fase A','Trip fase B','Trip fase C') 
plot(Irt_o,Iop,'b','LineWidth',2); 
ylim([-1,10]); 
hold on 
plot(Irt_1,y,'k','LineWidth',2); 
ylabel('Corrente de Operação'); 
title('Coordenograma do Relé Diferencial Percentual'); 
xlabel('Corrente de Restrição'); 
hold on 
plot(Irt_2,y2,'b','LineWidth',2); 
hold off 
%% Falta fase-terra 
 
case 3 
 load F1_P170_5_1G; 
 Ian=awgn([Gen_B_IF_Ian].*(1+(TC_erro/100)),45); 
 Iat=awgn([Gen_B_IF_Iat].*(1-(TC_erro/100)),45); 
 Ibn=awgn([Gen_B_IF_Ibn].*(1+(TC_erro/100)),45); 
 Ibt=awgn([Gen_B_IF_Ibt].*(1-(TC_erro/100)),45); 
 Icn=awgn([Gen_B_IF_Icn].*(1+(TC_erro/100)),45); 
 Ict=awgn([Gen_B_IF_Ict].*(1-(TC_erro/100)),45); 
[Inputs]=[x Ian/RTC Iat/RTC Ibn/RTC Ibt/RTC Icn/RTC Ict/RTC]; 
[Val_fas]=filtros(Inputs,fn,fm); 
IopA = abs (Val_fas(:,1)+Val_fas(:,2)); 
 IopB = abs (Val_fas(:,3)+Val_fas(:,4)); 
 IopC = abs (Val_fas(:,5)+Val_fas(:,6)); 
 IrtA = (abs (Val_fas(:,1))+abs(Val_fas(:,2)))/2; 
 IrtB = (abs (Val_fas(:,3))+abs(Val_fas(:,4)))/2; 
 IrtC = (abs (Val_fas(:,5))+abs(Val_fas(:,6)))/2; 
plot(IrtA,IopA,'k*','LineWidth',2) 
hold on 
plot(IrtB,IopB,'m+','LineWidth',2) 
hold on 
plot(IrtC,IopC,'yo','LineWidth',2) 
legend ('Trip fase A','Trip fase B','Trip fase C') 
plot(Irt_o,Iop,'b','LineWidth',2); 
ylim([-1,10]); 
hold on 
plot(Irt_1,y,'k','LineWidth',2); 
ylabel('Corrente de Operação'); 
title('Coordenograma do Relé Diferencial Percentual'); 
xlabel('Corrente de Restrição'); 
hold on 
plot(Irt_2,y2,'b','LineWidth',2); 
hold off 
%% Falta fase-terra 
 
case 4 
 load F1_P170_E_BT; 
 Ian=awgn([Gen_B_IF_Ian].*(1+(TC_erro/100)),45); 
 Iat=awgn([Gen_B_IF_Iat].*(1-(TC_erro/100)),45); 
 Ibn=awgn([Gen_B_IF_Ibn].*(1+(TC_erro/100)),45); 
 Ibt=awgn([Gen_B_IF_Ibt].*(1-(TC_erro/100)),45); 
 Icn=awgn([Gen_B_IF_Icn].*(1+(TC_erro/100)),45); 
 Ict=awgn([Gen_B_IF_Ict].*(1-(TC_erro/100)),45); 
[Inputs]=[x Ian/RTC Iat/RTC Ibn/RTC Ibt/RTC Icn/RTC Ict/RTC]; 
[Val_fas]=filtros(Inputs,fn,fm); 
IopA = abs (Val_fas(:,1)+Val_fas(:,2)); 
 IopB = abs (Val_fas(:,3)+Val_fas(:,4)); 
34 
 
 IopC = abs (Val_fas(:,5)+Val_fas(:,6)); 
 IrtA = (abs (Val_fas(:,1))+abs(Val_fas(:,2)))/2; 
 IrtB = (abs (Val_fas(:,3))+abs(Val_fas(:,4)))/2; 
 IrtC = (abs (Val_fas(:,5))+abs(Val_fas(:,6)))/2; 
plot(IrtA,IopA,'k*','LineWidth',2) 
hold on 
plot(IrtB,IopB,'m+','LineWidth',2) 
hold on 
plot(IrtC,IopC,'yo','LineWidth',2) 
legend ('Trip fase A','Trip fase B','Trip fase C') 
plot(Irt_o,Iop,'b','LineWidth',2); 
ylim([-1,10]); 
hold on 
plot(Irt_1,y,'k','LineWidth',2); 
ylabel('Corrente de Operação'); 
title('Coordenograma do Relé Diferencial Percentual'); 
xlabel('Corrente de Restrição'); 
hold on 
plot(Irt_2,y2,'b','LineWidth',2); 
hold off 
%% Falta fase-fase 
 
case 5 
 load F2_P15_1_30_1G; 
 Ian=awgn([Gen_B_IF_Ian].*(1+(TC_erro/100)),45); 
 Iat=awgn([Gen_B_IF_Iat].*(1-(TC_erro/100)),45); 
 Ibn=awgn([Gen_B_IF_Ibn].*(1+(TC_erro/100)),45); 
 Ibt=awgn([Gen_B_IF_Ibt].*(1-(TC_erro/100)),45); 
 Icn=awgn([Gen_B_IF_Icn].*(1+(TC_erro/100)),45); 
 Ict=awgn([Gen_B_IF_Ict].*(1-(TC_erro/100)),45); 
[Inputs]=[x Ian/RTC Iat/RTC Ibn/RTC Ibt/RTC Icn/RTC Ict/RTC]; 
[Val_fas]=filtros(Inputs,fn,fm); 
IopA = abs (Val_fas(:,1)+Val_fas(:,2)); 
 IopB = abs (Val_fas(:,3)+Val_fas(:,4)); 
 IopC = abs (Val_fas(:,5)+Val_fas(:,6)); 
 IrtA = (abs (Val_fas(:,1))+abs(Val_fas(:,2)))/2; 
 IrtB = (abs (Val_fas(:,3))+abs(Val_fas(:,4)))/2; 
 IrtC = (abs (Val_fas(:,5))+abs(Val_fas(:,6)))/2; 
plot(IrtA,IopA,'k*','LineWidth',2) 
hold on 
plot(IrtB,IopB,'m+','LineWidth',2) 
hold on 
plot(IrtC,IopC,'yo','LineWidth',2) 
legend ('Trip fase A','Trip fase B','Trip fase C') 
plot(Irt_o,Iop,'b','LineWidth',2); 
ylim([-1,10]); 
hold on 
plot(Irt_1,y,'k','LineWidth',2); 
ylabel('Corrente de Operação'); 
title('Coordenograma do Relé Diferencial Percentual'); 
xlabel('Corrente de Restrição'); 
hold on 
plot(Irt_2,y2,'b','LineWidth',2); 
hold off 
%% Falta fase-fase 
 
 case 6 
 load F2_P15_45_95_1G; 
 Ian=awgn([Gen_B_IF_Ian].*(1+(TC_erro/100)),45); 
 Iat=awgn([Gen_B_IF_Iat].*(1-(TC_erro/100)),45); 
 Ibn=awgn([Gen_B_IF_Ibn].*(1+(TC_erro/100)),45); 
35 
 
 Ibt=awgn([Gen_B_IF_Ibt].*(1-(TC_erro/100)),45); 
 Icn=awgn([Gen_B_IF_Icn].*(1+(TC_erro/100)),45); 
 Ict=awgn([Gen_B_IF_Ict].*(1-(TC_erro/100)),45); 
[Inputs]=[x Ian/RTC Iat/RTC Ibn/RTC Ibt/RTC Icn/RTC Ict/RTC]; 
[Val_fas]=filtros(Inputs,fn,fm); 
IopA = abs (Val_fas(:,1)+Val_fas(:,2)); 
 IopB = abs (Val_fas(:,3)+Val_fas(:,4)); 
 IopC = abs (Val_fas(:,5)+Val_fas(:,6)); 
 IrtA = (abs (Val_fas(:,1))+abs(Val_fas(:,2)))/2; 
 IrtB = (abs (Val_fas(:,3))+abs(Val_fas(:,4)))/2; 
 IrtC = (abs (Val_fas(:,5))+abs(Val_fas(:,6)))/2; 
plot(IrtA,IopA,'k*','LineWidth',2) 
hold on 
plot(IrtB,IopB,'m+','LineWidth',2) 
hold on 
plot(IrtC,IopC,'yo','LineWidth',2) 
legend ('Trip fase A','Trip fase B','Trip fase C') 
plot(Irt_o,Iop,'b','LineWidth',2); 
ylim([-1,10]); 
hold on 
plot(Irt_1,y,'k','LineWidth',2); 
ylabel('Corrente de Operação'); 
title('Coordenograma do Relé Diferencial Percentual'); 
xlabel('Corrente de Restrição'); 
hold on 
plot(Irt_2,y2,'b','LineWidth',2); 
hold off 
%% Falta fase-fase 
 
 case 7 
 load F2_P15_E_AT; 
 Ian=awgn([Gen_B_IF_Ian],45); 
 Iat=awgn([Gen_B_IF_Iat],45); 
 Ibn=awgn([Gen_B_IF_Ibn],45); 
 Ibt=awgn([Gen_B_IF_Ibt],45); 
 Icn=awgn([Gen_B_IF_Icn],45); 
 Ict=awgn([Gen_B_IF_Ict],45); 
[Inputs]=[x Ian/RTC Iat/RTC Ibn/RTC Ibt/RTC Icn/RTC Ict/RTC]; 
[Val_fas]=filtros(Inputs,fn,fm); 
IopA = abs (Val_fas(:,1).*(1+(TC_erro_slope1/100))+Val_fas(:,2).*(1-
(TC_erro_slope1/100))); 
IopB = abs (Val_fas(:,3).*(1+(TC_erro_slope1/100))+Val_fas(:,4).*(1-
(TC_erro_slope1/100))); 
IopC = abs (Val_fas(:,5).*(1+(TC_erro_slope1/100))+Val_fas(:,6).*(1-
(TC_erro_slope1/100))); 
IrtA = (abs (Val_fas(:,1).*(1+(TC_erro_slope1/100))+abs(Val_fas(:,2).*(1-
(TC_erro_slope1/100)))))/2; 
IrtB = (abs (Val_fas(:,3).*(1+(TC_erro_slope1/100))+abs(Val_fas(:,4).*(1-
(TC_erro_slope1/100)))))/2; 
IrtC = (abs (Val_fas(:,5).*(1+(TC_erro_slope1/100))+abs(Val_fas(:,6).*(1-
(TC_erro_slope1/100)))))/2; 
plot(IrtA,IopA,'k*','LineWidth',2) 
hold on 
plot(IrtB,IopB,'m+','LineWidth',2) 
hold on 
plot(IrtC,IopC,'yo','LineWidth',2) 
legend ('Trip fase A','Trip fase B','Trip fase C') 
plot(Irt_o,Iop,'b','LineWidth',2); 
ylim([-1,10]); 
hold on 
plot(Irt_1,y,'k','LineWidth',2); 
36 
 
ylabel('Corrente de Operação'); 
title('Coordenograma do Relé Diferencial Percentual'); 
xlabel('Corrente de Restrição'); 
hold on 
plot(Irt_2,y2,'b','LineWidth',2); 
hold off 
%% Falta fase-fase 
 
 case 8 
 load F2_P170_E_BT 
 Ian=awgn([Gen_B_IF_Ian].*(1+(TC_erro/100)),45); 
 Iat=awgn([Gen_B_IF_Iat].*(1-(TC_erro/100)),45); 
 Ibn=awgn([Gen_B_IF_Ibn].*(1+(TC_erro/100)),45); 
 Ibt=awgn([Gen_B_IF_Ibt].*(1-(TC_erro/100)),45); 
 Icn=awgn([Gen_B_IF_Icn].*(1+(TC_erro/100)),45); 
 Ict=awgn([Gen_B_IF_Ict].*(1-(TC_erro/100)),45); 
[Inputs]=[x Ian/RTC Iat/RTC Ibn/RTC Ibt/RTC Icn/RTC Ict/RTC]; 
[Val_fas]=filtros(Inputs,fn,fm); 
IopA = abs (Val_fas(:,1)+Val_fas(:,2)); 
 IopB = abs (Val_fas(:,3)+Val_fas(:,4)); 
 IopC = abs (Val_fas(:,5)+Val_fas(:,6)); 
 IrtA= (abs (Val_fas(:,1))+abs(Val_fas(:,2)))/2; 
 IrtB = (abs (Val_fas(:,3))+abs(Val_fas(:,4)))/2; 
 IrtC = (abs (Val_fas(:,5))+abs(Val_fas(:,6)))/2; 
plot(IrtA,IopA,'k*','LineWidth',2) 
hold on 
plot(IrtB,IopB,'m+','LineWidth',2) 
hold on 
plot(IrtC,IopC,'yo','LineWidth',2) 
legend ('Trip fase A','Trip fase B','Trip fase C') 
plot(Irt_o,Iop,'b','LineWidth',2); 
ylim([-1,10]); 
hold on 
plot(Irt_1,y,'k','LineWidth',2); 
ylabel('Corrente de Operação'); 
title('Coordenograma do Relé Diferencial Percentual'); 
xlabel('Corrente de Restrição'); 
hold on 
plot(Irt_2,y2,'b','LineWidth',2); 
hold off 
%% Falta fase-fase 
 
 case 9 
 load F2_P330_1_70_1G; 
 Ian=awgn([Gen_B_IF_Ian].*(1+(TC_erro/100)),45); 
 Iat=awgn([Gen_B_IF_Iat].*(1-(TC_erro/100)),45); 
 Ibn=awgn([Gen_B_IF_Ibn].*(1+(TC_erro/100)),45); 
 Ibt=awgn([Gen_B_IF_Ibt].*(1-(TC_erro/100)),45); 
 Icn=awgn([Gen_B_IF_Icn].*(1+(TC_erro/100)),45); 
 Ict=awgn([Gen_B_IF_Ict].*(1-(TC_erro/100)),45); 
[Inputs]=[x Ian/RTC Iat/RTC Ibn/RTC Ibt/RTC Icn/RTC Ict/RTC]; 
[Val_fas]=filtros(Inputs,fn,fm); 
IopA = abs (Val_fas(:,1)+Val_fas(:,2)); 
 IopB = abs (Val_fas(:,3)+Val_fas(:,4)); 
 IopC = abs (Val_fas(:,5)+Val_fas(:,6)); 
 IrtA = (abs (Val_fas(:,1))+abs(Val_fas(:,2)))/2; 
 IrtB = (abs (Val_fas(:,3))+abs(Val_fas(:,4)))/2; 
 IrtC = (abs (Val_fas(:,5))+abs(Val_fas(:,6)))/2; 
plot(IrtA,IopA,'k*','LineWidth',2) 
hold on 
plot(IrtB,IopB,'m+','LineWidth',2) 
37 
 
hold on 
plot(IrtC,IopC,'yo','LineWidth',2) 
legend ('Trip fase A','Trip fase B','Trip fase C') 
plot(Irt_o,Iop,'b','LineWidth',2); 
ylim([-1,10]); 
hold on 
plot(Irt_1,y,'k','LineWidth',2); 
ylabel('Corrente de Operação'); 
title('Coodenograma do Relé Diferencial Percentual'); 
xlabel('Corrente de Restrição'); 
hold on 
plot(Irt_2,y2,'b','LineWidth',2); 
hold off 
%% Falta trifásica 
 
 case 10 
 load F3_P15_1_1_1G; 
 Ian=awgn([Gen_B_IF_Ian].*(1+(TC_erro/100)),45); 
 Iat=awgn([Gen_B_IF_Iat].*(1-(TC_erro/100)),45); 
 Ibn=awgn([Gen_B_IF_Ibn].*(1+(TC_erro/100)),45); 
 Ibt=awgn([Gen_B_IF_Ibt].*(1-(TC_erro/100)),45); 
 Icn=awgn([Gen_B_IF_Icn].*(1+(TC_erro/100)),45); 
 Ict=awgn([Gen_B_IF_Ict].*(1-(TC_erro/100)),45); 
[Inputs]=[x Ian/RTC Iat/RTC Ibn/RTC Ibt/RTC Icn/RTC Ict/RTC]; 
[Val_fas]=filtros(Inputs,fn,fm); 
IopA = abs (Val_fas(:,1)+Val_fas(:,2)); 
 IopB = abs (Val_fas(:,3)+Val_fas(:,4)); 
 IopC = abs (Val_fas(:,5)+Val_fas(:,6)); 
 IrtA = (abs (Val_fas(:,1))+abs(Val_fas(:,2)))/2; 
 IrtB = (abs (Val_fas(:,3))+abs(Val_fas(:,4)))/2; 
 IrtC = (abs (Val_fas(:,5))+abs(Val_fas(:,6)))/2; 
plot(IrtA,IopA,'k*','LineWidth',2) 
hold on 
plot(IrtB,IopB,'m+','LineWidth',2) 
hold on 
plot(IrtC,IopC,'yo','LineWidth',2) 
legend ('Trip fase A','Trip fase B','Trip fase C') 
plot(Irt_o,Iop,'b','LineWidth',2); 
ylim([-1,10]); 
hold on 
plot(Irt_1,y,'k','LineWidth',2); 
ylabel('Corrente de Operação'); 
title('Coordenograma do Relé Diferencial Percentual'); 
xlabel('Corrente de Restrição'); 
hold on 
plot(Irt_2,y2,'b','LineWidth',2); 
hold off 
%% Falta trifásica 
 
 case 11 
 load F3_P15_E_BT; 
 Ian=awgn([Gen_B_IF_Ian].*(1+(TC_erro/100)),45); 
 Iat=awgn([Gen_B_IF_Iat].*(1-(TC_erro/100)),45); 
 Ibn=awgn([Gen_B_IF_Ibn].*(1+(TC_erro/100)),45); 
 Ibt=awgn([Gen_B_IF_Ibt].*(1-(TC_erro/100)),45); 
 Icn=awgn([Gen_B_IF_Icn].*(1+(TC_erro/100)),45); 
 Ict=awgn([Gen_B_IF_Ict].*(1-(TC_erro/100)),45); 
[Inputs]=[x Ian/RTC Iat/RTC Ibn/RTC Ibt/RTC Icn/RTC Ict/RTC]; 
[Val_fas]=filtros(Inputs,fn,fm); 
IopA = abs (Val_fas(:,1)+Val_fas(:,2)); 
 IopB = abs (Val_fas(:,3)+Val_fas(:,4)); 
38 
 
 IopC = abs (Val_fas(:,5)+Val_fas(:,6)); 
 IrtA = (abs (Val_fas(:,1))+abs(Val_fas(:,2)))/2; 
 IrtB = (abs (Val_fas(:,3))+abs(Val_fas(:,4)))/2; 
 IrtC = (abs (Val_fas(:,5))+abs(Val_fas(:,6)))/2; 
plot(IrtA,IopA,'k*','LineWidth',2) 
hold on 
plot(IrtB,IopB,'m+','LineWidth',2) 
hold on 
plot(IrtC,IopC,'yo','LineWidth',2) 
legend ('Trip fase A','Trip fase B','Trip fase C') 
plot(Irt_o,Iop,'b','LineWidth',2); 
ylim([-1,10]); 
hold on 
plot(Irt_1,y,'k','LineWidth',2); 
ylabel('Corrente de Operação'); 
title('Coordenograma do Relé Diferencial Percentual'); 
xlabel('Corrente de Restrição'); 
hold on 
plot(Irt_2,y2,'b','LineWidth',2); 
hold off 
%% Falta trifásica 
 
 case 12 
 load F3_P170_1_95_1G; 
 Ian=awgn([Gen_B_IF_Ian].*(1+(TC_erro/100)),45); 
 Iat=awgn([Gen_B_IF_Iat].*(1-(TC_erro/100)),45); 
 Ibn=awgn([Gen_B_IF_Ibn].*(1+(TC_erro/100)),45); 
 Ibt=awgn([Gen_B_IF_Ibt].*(1-(TC_erro/100)),45); 
 Icn=awgn([Gen_B_IF_Icn].*(1+(TC_erro/100)),45); 
 Ict=awgn([Gen_B_IF_Ict].*(1-(TC_erro/100)),45); 
[Inputs]=[x Ian/RTC Iat/RTC Ibn/RTC Ibt/RTC Icn/RTC Ict/RTC]; 
[Val_fas]=filtros(Inputs,fn,fm); 
IopA = abs (Val_fas(:,1)+Val_fas(:,2)); 
 IopB = abs (Val_fas(:,3)+Val_fas(:,4)); 
 IopC = abs (Val_fas(:,5)+Val_fas(:,6)); 
 IrtA = (abs (Val_fas(:,1))+abs(Val_fas(:,2)))/2; 
 IrtB = (abs (Val_fas(:,3))+abs(Val_fas(:,4)))/2; 
 IrtC = (abs (Val_fas(:,5))+abs(Val_fas(:,6)))/2; 
plot(IrtA,IopA,'k*','LineWidth',2) 
hold on 
plot(IrtB,IopB,'m+','LineWidth',2) 
hold on 
plot(IrtC,IopC,'yo','LineWidth',2) 
legend ('Trip fase A','Trip fase B','Trip fase C') 
plot(Irt_o,Iop,'b','LineWidth',2); 
ylim([-1,10]); 
hold on 
plot(Irt_1,y,'k','LineWidth',2); 
ylabel('Corrente de Operação'); 
title('Coordenograma do Relé Diferencial Percentual'); 
xlabel('Corrente de Restrição'); 
hold on 
plot(Irt_2,y2,'b','LineWidth',2); 
hold off 
%% Falta trifásica 
 
 case 13 
 load F4_P15_90_70_80_1G; 
 Ian=awgn([Gen_B_IF_Ian].*(1+(TC_erro/100)),45); 
 Iat=awgn([Gen_B_IF_Iat].*(1-(TC_erro/100)),45); 
 Ibn=awgn([Gen_B_IF_Ibn].*(1+(TC_erro/100)),45); 
39 
 
 Ibt=awgn([Gen_B_IF_Ibt].*(1-(TC_erro/100)),45); 
 Icn=awgn([Gen_B_IF_Icn].*(1+(TC_erro/100)),45); 
 Ict=awgn([Gen_B_IF_Ict].*(1-(TC_erro/100)),45); 
[Inputs]=[x Ian/RTC Iat/RTC Ibn/RTC Ibt/RTC Icn/RTC Ict/RTC]; 
[Val_fas]=filtros(Inputs,fn,fm); 
IopA = abs (Val_fas(:,1)+Val_fas(:,2)); 
 IopB = abs (Val_fas(:,3)+Val_fas(:,4)); 
 IopC = abs (Val_fas(:,5)+Val_fas(:,6)); 
 IrtA = (abs (Val_fas(:,1))+abs(Val_fas(:,2)))/2; 
 IrtB = (abs (Val_fas(:,3))+abs(Val_fas(:,4)))/2; 
 IrtC = (abs (Val_fas(:,5))+abs(Val_fas(:,6)))/2; 
plot(IrtA,IopA,'k*','LineWidth',2) 
hold on 
plot(IrtB,IopB,'m+','LineWidth',2) 
hold on 
plot(IrtC,IopC,'yo','LineWidth',2) 
legend ('Trip fase A','Trip fase B','Trip fase C') 
plot(Irt_o,Iop,'b','LineWidth',2); 
ylim([-1,10]); 
hold on 
plot(Irt_1,y,'k','LineWidth',2); 
ylabel('Corrente de Operação'); 
title('Coordenograma do Relé Diferencial Percentual'); 
xlabel('Corrente de Restrição'); 
hold on 
plot(Irt_2,y2,'b','LineWidth',2); 
hold off 
%% Falta trifásica 
 
 case 14 
 load F4_P15_E_AT; 
 Ian=awgn([Gen_B_IF_Ian].*(1+(TC_erro/100)),45); 
 Iat=awgn([Gen_B_IF_Iat].*(1-(TC_erro/100)),45); 
 Ibn=awgn([Gen_B_IF_Ibn].*(1+(TC_erro/100)),45); 
 Ibt=awgn([Gen_B_IF_Ibt].*(1-(TC_erro/100)),45); 
 Icn=awgn([Gen_B_IF_Icn].*(1+(TC_erro/100)),45); 
 Ict=awgn([Gen_B_IF_Ict].*(1-(TC_erro/100)),45); 
[Inputs]=[x Ian/RTC Iat/RTC Ibn/RTC Ibt/RTC Icn/RTC Ict/RTC]; 
[Val_fas]=filtros(Inputs,fn,fm); 
IopA = abs (Val_fas(:,1)+Val_fas(:,2)); 
 IopB = abs (Val_fas(:,3)+Val_fas(:,4)); 
 IopC = abs (Val_fas(:,5)+Val_fas(:,6)); 
 IrtA = (abs (Val_fas(:,1))+abs(Val_fas(:,2)))/2; 
 IrtB = (abs (Val_fas(:,3))+abs(Val_fas(:,4)))/2; 
 IrtC = (abs (Val_fas(:,5))+abs(Val_fas(:,6)))/2; 
plot(IrtA,IopA,'k*','LineWidth',2) 
hold on 
plot(IrtB,IopB,'m+','LineWidth',2) 
hold on 
plot(IrtC,IopC,'yo','LineWidth',2) 
legend ('Trip fase A','Trip fase B','Trip fase C') 
plot(Irt_o,Iop,'b','LineWidth',2); 
ylim([-1,10]); 
hold on 
plot(Irt_1,y,'k','LineWidth',2); 
ylabel('Corrente de Operação'); 
title('Coordenograma do Relé Diferencial Percentual'); 
xlabel('Corrente de Restrição');hold on 
plot(Irt_2,y2,'b','LineWidth',2); 
hold off 
40 
 
%% Falta trifásica 
 
case 15 
 load F4_P170_1_5_3_1G; 
 Ian=awgn([Gen_B_IF_Ian].*(1+(TC_erro/100)),45); 
 Iat=awgn([Gen_B_IF_Iat].*(1-(TC_erro/100)),45); 
 Ibn=awgn([Gen_B_IF_Ibn].*(1+(TC_erro/100)),45); 
 Ibt=awgn([Gen_B_IF_Ibt].*(1-(TC_erro/100)),45); 
 Icn=awgn([Gen_B_IF_Icn].*(1+(TC_erro/100)),45); 
 Ict=awgn([Gen_B_IF_Ict].*(1-(TC_erro/100)),45); 
[Inputs]=[x Ian/RTC Iat/RTC Ibn/RTC Ibt/RTC Icn/RTC Ict/RTC]; 
[Val_fas]=filtros(Inputs,fn,fm); 
IopA = abs (Val_fas(:,1)+Val_fas(:,2)); 
 IopB = abs (Val_fas(:,3)+Val_fas(:,4)); 
 IopC = abs (Val_fas(:,5)+Val_fas(:,6)); 
 IrtA = (abs (Val_fas(:,1))+abs(Val_fas(:,2)))/2; 
 IrtB = (abs (Val_fas(:,3))+abs(Val_fas(:,4)))/2; 
 IrtC = (abs (Val_fas(:,5))+abs(Val_fas(:,6)))/2; 
plot(IrtA,IopA,'k*','LineWidth',2) 
hold on 
plot(IrtB,IopB,'m+','LineWidth',2) 
hold on 
plot(IrtC,IopC,'yo','LineWidth',2) 
legend ('Trip fase A','Trip fase B','Trip fase C') 
plot(Irt_o,Iop,'b','LineWidth',2); 
ylim([-1,10]); 
hold on 
plot(Irt_1,y,'k','LineWidth',2); 
ylabel('Corrente de Operação'); 
title('Coordenograma do Relé Diferencial Percentual'); 
xlabel('Corrente de Restrição'); 
hold on 
plot(Irt_2,y2,'b','LineWidth',2); 
hold off 
%% Falta trifásica 
 
case 16 
 load F4_P170_95_1_1_1G; 
 Ian=awgn([Gen_B_IF_Ian].*(1+(TC_erro/100)),45); 
 Iat=awgn([Gen_B_IF_Iat].*(1-(TC_erro/100)),45); 
 Ibn=awgn([Gen_B_IF_Ibn].*(1+(TC_erro/100)),45); 
 Ibt=awgn([Gen_B_IF_Ibt].*(1-(TC_erro/100)),45); 
 Icn=awgn([Gen_B_IF_Icn].*(1+(TC_erro/100)),45); 
 Ict=awgn([Gen_B_IF_Ict].*(1-(TC_erro/100)),45); 
[Inputs]=[x Ian/RTC Iat/RTC Ibn/RTC Ibt/RTC Icn/RTC Ict/RTC]; 
[Val_fas]=filtros(Inputs,fn,fm); 
IopA = abs (Val_fas(:,1)+Val_fas(:,2)); 
 IopB = abs (Val_fas(:,3)+Val_fas(:,4)); 
 IopC = abs (Val_fas(:,5)+Val_fas(:,6)); 
 IrtA = (abs (Val_fas(:,1))+abs(Val_fas(:,2)))/2; 
 IrtB = (abs (Val_fas(:,3))+abs(Val_fas(:,4)))/2; 
 IrtC = (abs (Val_fas(:,5))+abs(Val_fas(:,6)))/2; 
plot(IrtA,IopA,'k*','LineWidth',2) 
hold on 
plot(IrtB,IopB,'m+','LineWidth',2) 
hold on 
plot(IrtC,IopC,'yo','LineWidth',2) 
legend ('Trip fase A','Trip fase B','Trip fase C') 
plot(Irt_o,Iop,'b','LineWidth',2); 
ylim([-1,10]); 
hold on 
41 
 
plot(Irt_1,y,'k','LineWidth',2); 
ylabel('Corrente de Operação'); 
title('Coordenograma do Relé Diferencial Percentual'); 
xlabel('Corrente de Restrição'); 
hold on 
plot(Irt_2,y2,'b','LineWidth',2); 
hold off 
%% Falta trifásica 
 
case 17 
 load F4_P330_90_70_80_1G; 
 Ian=awgn([Gen_B_IF_Ian].*(1+(TC_erro/100)),45); 
 Iat=awgn([Gen_B_IF_Iat].*(1-(TC_erro/100)),45); 
 Ibn=awgn([Gen_B_IF_Ibn].*(1+(TC_erro/100)),45); 
 Ibt=awgn([Gen_B_IF_Ibt].*(1-(TC_erro/100)),45); 
 Icn=awgn([Gen_B_IF_Icn].*(1+(TC_erro/100)),45); 
 Ict=awgn([Gen_B_IF_Ict].*(1-(TC_erro/100)),45); 
[Inputs]=[x Ian/RTC Iat/RTC Ibn/RTC Ibt/RTC Icn/RTC Ict/RTC]; 
[Val_fas]=filtros(Inputs,fn,fm); 
IopA = abs (Val_fas(:,1)+Val_fas(:,2)); 
 IopB = abs (Val_fas(:,3)+Val_fas(:,4)); 
 IopC = abs (Val_fas(:,5)+Val_fas(:,6)); 
 IrtA = (abs (Val_fas(:,1))+abs(Val_fas(:,2)))/2; 
 IrtB = (abs (Val_fas(:,3))+abs(Val_fas(:,4)))/2; 
 IrtC = (abs (Val_fas(:,5))+abs(Val_fas(:,6)))/2; 
plot(IrtA,IopA,'k*','LineWidth',2) 
hold on 
plot(IrtB,IopB,'m+','LineWidth',2) 
hold on 
plot(IrtC,IopC,'yo','LineWidth',2) 
legend ('Trip fase A','Trip fase B','Trip fase C') 
plot(Irt_o,Iop,'b','LineWidth',2); 
ylim([-1,10]); 
hold on 
plot(Irt_1,y,'k','LineWidth',2); 
ylabel('Corrente de Operação'); 
title('Coordenograma do Relé Diferencial Percentual'); 
xlabel('Corrente de Restrição'); 
hold on 
plot(Irt_2,y2,'b','LineWidth',2); 
hold off 
end