Avaliação da operação dos relés digitais sob os efeitos aliasing e leakage

utilizando transformada de Fourier

Elian Salazar dos Reis*. Raiany Oliveira Freitas.**
Antonio Richard H. Carvalho*** Domingos Robson Silva Costa****

*Graduando em engenharia elétrica do Instituto Federal de Educação, Ciência e Tecnologia do Maranhão - Campus
Imperatriz, EEI170319,
BR (e-mail: elian.reis@acad.ifma.edu.br).
**Graduando em engenharia elétrica do Instituto Federal de Educação, Ciência e Tecnologia do Maranhão - Campus
Imperatriz, EEI170344,
BR (e-mail: raiany.oliveiras@acad.ifma.edu.br).
***Graduando em engenharia elétrica do Instituto Federal de Educação, Ciência e Tecnologia do Maranhão - Campus
Imperatriz, EEI170312,
BR (e-mail: antonio.richard@acad.ifma.edu.br)
.**** Professor orientador do curso de engenharia elétrica do Instituto Federal de Educação, Ciência e Tecnologia do
Maranhão - Campus Imperatriz,
BR (e-mail: domingos.costa@ifma.edu.br)

Abstract: The present work presents a study on the protection of electrical power systems whose
function is to cause the immediate removal of operation of elements affected by short circuits, as well as
elements that operate abnormally, interfering in the operation of the electrical system as a whole. Thus,
digital protection relays are an important component of these systems, and require voltage and current
signals present in the protected equipment. For this, they are fed through the associated PTs (potential
transformers) and CTs (current transformers), which are responsible for guaranteeing the fidelity of the
values ​to the relay so that it can process them and decide whether or not there is a need of acting. Thus,
digital relays are able to identify defects in these instruments, such as a lack of voltage in only one phase,
for example, make angular compensations in the cases of power transformers and treat other variations in
order to avoid wrong decisions at the time of their operation as protection element. Thus, this work seeks
to understand, through computer simulations carried out in matlab software, the basic theory of digital
protection relays, as well as to analyze and evaluate the impact that aliasing and leakage effects impose
on phasor measurements computed by the discrete transform of Fourier.
Resumo: O presente trabalho apresenta um estudo sobre a proteção dos sistemas elétricos de potência,
que possui por função provocar a remoção imediata de operação dos elementos acometidos por
curtos-circuitos, assim como dos elementos que operam de modo anormal, interferindo na operação do
sistema elétrico como um todo. Com isso, os relés digitais de proteção são um importante componente
desses sistemas, e necessitam dos sinais de tensão e corrente presentes no equipamento protegido. Para
isso, eles são alimentados através dos TPs (transformadores de potencial) e dos TCs (transformadores de
corrente) associados, que são os responsáveis por garantir ao relé a fidelidade dos valores para que este
possa processá-los e decidir se há ou não necessidade de atuação. Assim, os relés digitais são capazes de
identificar defeitos nestes instrumentos, como falta de tensão em apenas uma fase, por exemplo, fazer
compensações angulares nos casos dos transformadores de potência e tratar outras variações com intuito
de evitar decisões equivocadas na hora de sua operação como elemento de proteção. Desta forma, neste
trabalho procura-se compreender através de simulações computacionais realizadas no software matlab, a
teoria básica sobre os relés digitais de proteção, assim como, analisar e avaliar o impacto que os efeitos
aliasing e leakage impõe nas medições fasoriais computadas pela transformada discreta de Fourier.
Keywords: Protection; Systems; Reles; Digital; Fourier.
Palavras-chaves: Proteção; Sistemas; Réles; Digitais; Fourier.
 1. INTRODUÇÃO A proteção de sistemas elétricos de potência tem
dois principais objetivos: distribuir energia com o mínimo de
A energia move o mundo. Dentro desse contexto, o
interrupções e garantir segurança dos equipamentos e
sistema de energia elétrica composta por geração, transmissão
instalações. Segundo Goes (2013, p.03) “um sistema de
e distribuição é o principal responsável a suprir de forma
proteção deve garantir o processamento correto de um
segura e confiável a energia aos consumidores. No entanto, o
possível defeito ou falta (termo usado para designar um
sistema elétrico tem um alto grau de complexidade de
afastamento acidental das condições normais de operação)
operação, estando sujeito a vários tipos de perturbações que
bem como rapidez na resposta, iniciando ações corretivas
levam a falta de controle de operação do mesmo, resultando
para retornar o mais rápido possível para uma condição
na interrupção do fornecimento de energia afetando uma
normal de operação”.
grande parte dos consumidores.
Frazão (2019, p.11) afirma que “a função do sistema
Na sociedade atual, a interrupção do sistema de
de proteção é provocar a remoção imediata de operação dos
energia pode impactar negativamente tanto o sistema
elementos acometidos por curtos-circuitos, assim como dos
econômico como o social do país, pois está diretamente
elementos que operam de modo anormal, interferindo na
ligado às atividades econômicas e de lazer do dia a dia do
operação do sistema elétrico como um todo”.
consumidor. Portanto, a proteção dos sistemas elétricos vem
No trabalho de Oliveira (2018), para atuar de
com o objetivo de distribuir essa energia sem interrupções,
maneira rápida e precisa, o sistema de proteção deve atender
além de garantir a segurança.
as seguintes características:
De acordo Frazão (2019), em 1960 surgiram várias
- Sensibilidade: Entende-se por sensibilidade a capacidade de
pesquisas que visavam implementar funções de segurança
a proteção responder às anormalidades nas condições de
nos equipamentos, principalmente os computadores, foi
operação, e aos curtos circuitos para os quais foi projetada;
quando deu início ao uso de relés digitais na proteção dos
- Velocidade: é a rapidez em que o sistema de proteção atua
sistemas elétricos. Os relés digitais analisam os sinais
após a falta. Essa é uma das principais características do
analógicos de tensão e corrente, verificando se existe um
sistema, pois uma corrente de curto circuito circulando por
curto circuito no sistema.
um tempo longo pode prejudicar a integridade do isolamento
Esse artigo tem como finalidade avaliar a operação
dos equipamentos;
dos relés digitais sob os efeitos aliasing e leakage, utilizando
- Confiabilidade: define-se como a probabilidade de um
transformada de fourier. O seu propósito é fazer a análise de
componente, um equipamento ou um sistema satisfazer a
todos os fatores relevantes à proteção do sistema de energia,
função prevista, sob dadas circunstâncias;
além de expor os resultados da simulação feitas no programa
- Segurança: é a capacidade do sistema não operar para faltas
Matlab.
que ocorrem fora de sua zona de proteção;
- Seletividade: É a propriedade da proteção em reconhecer e
2. REFERENCIAL TEÓRICO selecionar entre aquelas condições para as quais uma
imediata operação é requerida, e aquelas para as quais
nenhuma operação ou um retardo de atuação é exigido.
2.1 Proteção dos sistemas elétricos de potência Desta forma, o sistema básico de proteção consiste
em um disjuntor, transformadores de corrente (TC),
O termo proteção de sistemas elétricos de transformadores de potencial (TP), relés auxiliares e o relés
transmissão se refere aos dispositivos que são responsáveis de proteção. Os conjuntos de proteção são compostos pelos
por garantir a continuidade do suporte de cargas, e evitar que relés de funções principais e os dispositivos auxiliares, que
os distúrbios que podem acontecer no sistema causem danos servem para sinalização, expansão das capacidades dos relés,
às instalações das subestações. Esses distúrbios podem ser circuitos de bloqueio e intertravamentos, entre outras funções
curtos-circuitos, sobrecargas ou sobretensões vindas de (OLIVEIRA, 2018).
rejeição de carga, ou perda de sincronismo em caso de
sistema de geração de energia (OLIVEIRA, 2018).
Os sistemas elétricos estão sujeitos a falhas de 2.2 Relés digitais de proteção
diversas naturezas. Nos casos das linhas de transmissão, estes
O uso de relés digitais na proteção dos sistemas
são os componentes mais susceptíveis de desarmes, por haver
elétricos se iniciou na década de 1960 com pesquisas que
a possibilidade de danos por vandalismo, descargas
estudavam a implementação das funções de segurança em
atmosféricas, acidentes com veículos, entre outros
computadores. Os problemas de curto-circuito, fluxo de carga
(OLIVEIRA, 2018).
e estabilidade estavam sendo convertidos em algoritmos
Nos demais equipamentos do SEP, como
computacionais, apesar dos computadores da época
transformadores e geradores, por exemplo, os defeitos podem
apresentarem baixa capacidade de processamento (FRAZÃO,
surgir por falhas na isolação elétrica, sobrecargas elevadas,
2019).
animais em contato com as partes energizadas, etc. Os
sistemas de proteção têm que atuar com segurança para todos “O relé digital foi pensado como ferramenta de auxílio desses
esses casos. computadores. Infelizmente, o início dessas novas
tecnologias foi um tanto quanto perturbado, pois não havia frequência máxima alcançada pelo sinal analógico para
incentivos econômicos para tal. No entanto, a busca pelo possibilitar o registro digital de todas as frequências
desenvolvimento de algoritmos para computadores e relés analisadas (LAZZETTA, 2009). Essa taxa de amostragem
pareceu atraente para vários pesquisadores, ou seja, através mínima desejada é conhecida como frequência de Nyquist.
da curiosidade acadêmica, esse fértil campo foi estimulado. A Quando a taxa de amostragem para determinado
evolução dos computadores e relés foi tão rápida que a sinal é menor que a frequência de Nyquist, acontece o efeito
robustez algorítmica exigida por esses equipamentos aliasing ou foldover, em que uma alta frequência é medida
finalmente encontrou correspondência na velocidade e erroneamente como sendo de frequência mais baixa. Ou seja,
capacidade de processamento, de modo que, atualmente, o um sinal cuja máxima frequência é 60 Hz, tem que ser
binômio computador-relé oferece a melhor solução amostrado com uma frequência de no mínimo 120 Hz. Caso a
econômica e técnica para os problemas de proteção – em frequência de amostragem não obedeça a essa relação, ocorre
muitos casos, a única solução viável” (PHADKE; THORP, o efeito aliasing, ou recobrimento de frequência.
apud FRAZÃO, 2019).
Segundo um estudo de Oliveira (2018) “os relés Fig.1- Efeito aliasing
digitais assemelham-se aos computadores pelo ponto de vista
do processamento das informações, visto que possuem
microprocessadores, memórias e a capacidade de executar
programas”. Também são compostos por componentes
eletrônicos e contam com os recursos de comunicação entre
dispositivos, podendo ser utilizado em rede ou isolados.
Para Goes (2013, p. 22) “esse tipo de relé é
gerenciado por um microprocessador e controlado por
software. Os princípios de funcionamento dos relés
convencionais são uma referência para seu
desenvolvimento”. A utilização de um processador para
tarefas de proteção possibilitou soluções para limitações
encontradas nas tecnologias analógicas.
Este relé é superior aos outros tipos pelo fato de
possuir uma melhor tecnologia de processamento sendo
capaz de processar mais de uma função de proteção ao
mesmo tempo, além de fazer o registro das faltas e dos Fonte: HIGUTI, 2012.
eventos, o que ajuda os profissionais na procura pelas causas
de determinadas ocorrências.
Com isso, para verificar se existe um curto circuito Quando ocorre este efeito, frequências originalmente
no sistema, o relé digital analisa sinais analógicos de tensão e altas no espectro do sinal aparecem em regiões de mais baixa
corrente da subestação. Essas tensões e correntes são dadas frequência, impossibilitando a recuperação do sinal original.
por meio de transformadores de potencial (TP) e “Caso a frequência de amostragem não obedeça o
transformadores de corrente (TC), respectivamente. Os sinais teorema de Nyquist, deve-se colocar um filtro passa-baixas
analógicos passam por um condicionamento para que suas antes da aquisição de dados, que limite a máxima frequência
magnitudes fiquem iguais ao nível de tensão da placa do sinal a metade da frequência de amostragem. Este filtro é
eletrônica do relé, em seguida, eles são convertidos em chamado de filtro anti-aliasing” (HIGUTI, 2012).
digitais pelo conversor analógico-digital (CAD) que faz o Desta forma, “o sinal analógico, antes da
processo de amostragem e retenção para aquisição de amostragem, é passado por um filtro anti-alising, com o
amostras do sinal analógico (FRAZÃO, 2019). objetivo de reduzir ou evitar o efeito aliasing” (Frazão, 2019,
p.30). Assim, quanto maior for a taxa de amostragem, mais
O CAD possui um erro de quantização também
medidas do sinal serão realizadas em um mesmo intervalo de
chamado de ruído de quantização. A equação do erro de
tempo, e assim, maior é a fidelidade do sinal digital a respeito
quantização pode ser escrita como:
do analógico.
𝑉 −𝑁
𝑄𝑒 = 𝑁−1 =2 𝑉 (1)
2𝑥2
Em que V é a tensão máxima de entrada em Volts e 2.4. Efeito Leakage
N é o número de bits.
Se um sinal estiver operando na frequência 𝑓0 e, por
2.3 Efeito aliasing
exemplo, N=36 amostras são retiradas do seu período T, a
Segundo o Teorema de Nyquist, é necessário uma frequência de amostragem será igual a 2160 Hz (valor
taxa de amostragem de no mínimo duas vezes o valor da inteiro), assim, para essa condição, os cálculos executados
pelo processador com base nessas amostras estarão livres de
erros. Contudo, se o sinal estiver operando na frequência, por Onde,
exemplo, de 60 + 0,02 Hz, a frequência de amostragem será
igual a 2160,72 Hz (valor racional), que não pode ser obtida 𝑇
2
−𝑗𝑛𝑤𝑜𝑡
1
fielmente por um CAD com clock de amostragem fixo. 𝑐𝑛 = ∫ 𝑓(𝑡)𝑒 (3)
𝑇
Portanto, o resíduo de 0,72 Hz irá impor erros nos cálculos −𝑇
2
executados pelo processador. Tal erro é denominado de efeito
leakage, que em essência provoca o que denominamos de Para transformar uma função periódica em
espalhamento espectral de frequência (FRAZÃO, 2019). aperiódica podemos fazer 𝑇 → ∞ o que significa que a
Este efeito pode ser visualizado no sinal a seguir, função não se repete. Também devemos observar que quando
que apresenta a mesma forma de onda mostrada T aumenta então a distância entre duas frequências
anteriormente, mas com um número não inteiro de ciclos, harmônicas (𝑛𝑤0𝑒 (𝑛 + 1)𝑤0) consecutivas se torna cada
bem como o respectivo espectro. Notar que a identificação da
vez menor. Assim, matematicamente temos o seguinte:
frequência ainda se mostrou correta (o que nem sempre
acontece, pois pode haver um deslocamento do pico), mas 1 1 ∆𝑤
com grande erro de amplitude. (S.M.Deckmann e ∆𝑤 = (𝑛 + 1)𝑤0 − 𝑛𝑤0 = 𝑤0 = 𝑇
⇒ 𝑇
= 2π
J.A.Pomilio, 2015) Da relação anterior, quando 𝑇 → ∞ temos o

Fig.2- Amostragem de ciclos fracionários e seguinte:

respectivo espectro. 1
→
𝑑𝑤
𝑇 2π

Portanto, quando 𝑇 → ∞ a frequência deixa de ser uma

variável discreta e se torna uma variável contínua. Assim,
temos:
𝑛𝑤0 quando 𝑇 → ∞
Observando a relação (3) verificamos que 𝑐𝑛
Fonte: S.M.Deckmann e J.A.Pomilio, 2015. também varia inversamente com T. Assim, quando 𝑇 → ∞
então 𝑐𝑛 → 0 . Entretanto, o valor limite do produto 𝑐𝑛𝑇
Este efeito ocorre devido ao truncamento da série de
assume a seguinte forma:
sinais discretos no tempo quando estes são coletados em uma ∞
−𝑗𝑤𝑡
condição na qual a descontinuidade imposta pelo corte faz 𝑐𝑛𝑇 → ∫ 𝑓(𝑡)𝑒 𝑑𝑡 quando 𝑇 → ∞ (4)
surgir conteúdos de frequência inexistentes no sinal de −∞
resposta coletado. Dificilmente o sinal é de tal forma que um A integral mostrada na relação acima é chamada de
número inteiro de ciclos se encaixa no tempo de medição, transformada de Fourier de f(t) e pode ser representada da
gerando o espectro de frequências esperado. seguinte forma:
∞
Assim, o espectro de frequências é ampliado pelas −𝑗𝑤𝑡
"falhas" nas extremidades na grande maioria dos casos. O 𝐹(𝑤) = ፑ{𝑓(𝑡)} = ∫ 𝑓(𝑡)𝑒 𝑑𝑡 (5)
−∞
problema é agravado pelo fato de o tamanho do pulso
aleatório depender de quando a primeira medição ocorreu no
ciclo, assim a ampliação pode se alterar se a medição for Também podemos encontrar uma relação para a
repetida. Ao analisar uma sequência de dados finita através transformada inversa de Fourier analisando a forma limite da
de Fourier ou outro método de análise espectral, o relação (2) quando 𝑇 → ∞ . Assim, multiplicamos e
janelamento pode ser utilizado no intuito de se minimizar o dividimos o segundo membro por T e temos o seguinte:
efeito das formas de onda truncadas, reduzindo dessa forma a ∞ −𝑗𝑛𝑤0𝑡 1
perda espectral (PINTO apud ANDRADE;SOARES, 2000). 𝑓(𝑡) = ∑ [𝑐𝑛𝑇]𝑒
−∞
( )
𝑇

2.5 Transformada de Fourier

Quando 𝑇 → ∞ então o somatório tende a uma
1 𝑑𝑤
Segundo Oppenhiem e Willsky (2010) pode-se encontrar integral 𝑐𝑛𝑇 → 𝐹(𝑤), 𝑛𝑤0 → 𝑤 𝑒 𝑇 ⇒ 2π . Assim, temos o
uma definição da transformada de Fourier como sendo um seguinte:
caso limite da série de Fourier. Lembrando que a série de
Fourier na forma exponencial assume a seguinte forma: ∞
𝑗𝑤𝑡 𝑑𝑤
∞ 𝑗𝑛𝑤𝑜𝑡 𝑓(𝑡) = ∫ 𝐹(𝑤)𝑒 2π
⇒
𝑓(𝑡) = ∑ 𝑐𝑛𝑒 (2) −∞
𝑛=−∞
 ∞ π
𝑗𝑤𝑡
𝑓(𝑡) =
1
2π
∫ 𝐹(𝑤)𝑒 𝑑𝑤 (6) 𝑥[𝑛] =
1
2π
∫𝑋 𝑒 ( 𝑗𝑤)𝑒𝑗𝑤𝑛𝑑𝑤
−∞ −π
(8)
Assim, a relação (5) transforma a expressão no
domínio do tempo f(t) em uma relação no domínio da As relações (7) e (8) foram deduzidas supondo que x[n] é
frequencia F(w). A relação (6) faz a operação inversa e aperiódica (sinal com duração finita).
transforma F(w) em f(t).

3. RESULTADOS E DISCUSSÃO
Fig.3- Gráfico da transformada contínua de Fourier do sinal
de 5 Hz.
A primeira simulação é feita sem os efeitos, cujo
código se encontra no anexo A. Sendo o número de
amostras de N = 32, a frequência fundamental fo=60 Hz, a
frequência de amostragem 1920 Hz, para mais detalhes
observar o código em anexo.

Fig.4 - Coeficientes de Fourier

Fonte: Cesar da Costa, 2017

Na imagem acima, é mostrado um sinal no tempo e

abaixo a transformada de Fourier deste sinal. Assim, é visto o
sinal no domínio do tempo, então é feito a amplitude do sinal
conforme varia o tempo, e onde foi feita a transformada tem o
módulo da amplitude desse sinal variando conforme a
frequência.
Fonte: Autoral
2.6 Transformada discreta de Fourier
Fig.5 - Fasores harmônico
A transformada discreta de Fourier é usada quando
pretendemos analisar sinais não periódicos e de tempo
discreto.
A transformada de Fourier (DTFT) do sinal 𝑥 𝑛 assume
a seguinte forma:

∞
𝑗Ω −𝑗𝑤𝑛
( )
𝑋 𝑒 = ∑ 𝑥[𝑛]𝑒 (7)
𝑛=−∞

Por outro lado, a transformada discreta de Fourier

inversa assume a seguinte forma:

Fonte: Autoral
 efeito aliasing, a qual pode observar também alterações em
As Figuras 6 e 7 são com acréscimo de desvio de relação às figuras 4 e 5 que não tem efeitos. Portanto pode-se
frequência de 0.92 Hz para a avaliação do efeito de leakage, concluir a parte da simulação que os efeitos aliasing e leakage
utilizando ainda o código no Anexo A com essa alteração faria com que os relés de proteção atuasse de maneira
específica. imprecisa com o esperado, seja atuando quando não deveria
. ou ausência de tomada de decisão para proteger o sistema.
Fig.6 - Coeficiente de Fourier com o efeito leakage
Fig.8 - Coeficientes de Fourier com efeito aliasing

Fonte: Autoral
Fonte: Autoral

Fig.7 - Fasores harmônicos com efeito leakage Fig.9 - Fasores harmônico com efeito aliasing

Fonte: Autoral Fonte: Autoral

A parte do tempo 1s os valores começa a diferenciar
do coeficiente de fourier sem efeito leakage de forma que o
4. CONCLUSÕES
relés irá ler com a distorção do sinal analógico para o sinal
digital, observando as figuras 4 e 6 que são do coeficiente de O desenvolvimento do presente estudo, possibilitou
fourier, já nas figuras 5 e 7 que são dos harmônicos a uma análise e melhor compreensão através do uso de
distorção é mais significativas tanto para os coeficientes de simulações feitas sobre os relés digitais de proteção e como
fourier quanto para o fasores harmônico a parte do tempo de ele funcionam. Além disso, propiciou observar e avaliar o
14s, percebe a presença da frequência de espúrios, que é uma impacto causado com o uso dos efeitos aliasing e leakage nas
frequência inexistente no sinal análogico, mas que aparece simulações, e a utilização da transformada discreta de
por esse efeito neste caso. fourier nas medições fasorial computadas.
Por último, fazendo alteração na taxa de amostragem
para 30 Hz no código em Anexo A, gera-se os gráficos do
 REFERÊNCIAS

COSTA, César. EXERCICIOS TRANSFORMADA DE

FOURIER COM MATLAB. 2017. Disponível ANEXO A
em:http://professorcesarcosta.com.br/upload/imagens_upload
/Apostila%202%20-%20MATLAB.pdf. Acesso em 01 out.
2021. N = 32; % Nº de amostras
em um Ciclo. Período Fundamental
DECKMANN, S.M. POMILIO, J.A. Avaliação da
fo = 60; % freq.
Qualidade da Energia Elétrica. 2015. Disponível em:
fundamental Hz (W = 2pi/N; tb é a menor frequência )
https://www.dsce.fee.unicamp.br/~antenor/pdffiles/qualidade/
b4.pdf. Acesso em 03 out. 2021. Df = 0; % define desvio de
frequência em Hz %
FRAZÃO, R. J. A. Proteção do sistema elétrico de encaixa sobre efeito leakage e fenômeno Aliasing
potência. 2019. Editora e Distribuidora Educacional S.A.
tsamp = linspace(0,1/fo,N+1); % Tempo
Disponível em:
de amostragem em Segundos a partir da frequência de
http://cm-kls-content.s3.amazonaws.com/201901/INTERATI
amostragem fa x[n]=x[n+N]
VAS_2_0/PROTECAO_DO_SISTEMA_ELETRICO_DE_P
OTENCIA/U1/LIVRO_UNICO.pdf. Acesso em 01 out. DC = 0;
2021.
Harm = 18;
HIGUTI, R. T. Amostragem e Reconstruçãoo de Sinais. Amp_harm = 2;
2012. Disponível
Amp_fund = 10;
em:https://www.feis.unesp.br/Home/departamentos/engenhari
aeletrica/ele1095_3_amostragem.pdf. Acesso em 06 nov, ysamp = DC + Amp_fund*cos(2*pi*(fo+Df).*tsamp + pi/7)
2021. + Amp_harm*cos(Harm*2*pi*(fo+Df).*tsamp+ pi/8); %
vetor de amostras
GOES, A. R. G. Modernização da proteção de sistemas
counter = 0; % Incrematador
elétricos de potência. 2013. Disponível em:
do valor da T. de F.
http://repositorio.poli.ufrj.br/monografias/monopoli10006073
.pdf. Acesso em 05 nov. 2021.
for k=0:N/2
LAZZETTA, F. Análise Dinâmica dos Sistemas elétricos de
potência. Escola de Comunicações e Artes da Universidade for n=0:N-1
de São Paulo. 2009. Disponível em: https://www.eca.usp.br/
counter=counter+ysamp(n+1)*(cos(2*pi*n*k/N) -
Acesso em 06 nov. 2021
1i*sin(2*pi*n*k/N));
OLIVEIRA, L. H. E. B. Ensaios de relés digitais de end
proteção utilizando o protocolo IEC 61850. 2018.
phasors(k+1,:)=counter;
Disponível em:
http://esc.com.br/artigos/ensaios-de-reles-digitais-de-protecao counter=0;
-utilizando-o-protocolo-iec-61850. Acesso em 01 out. 2021.

OPPENHIEM, A. V.; WILLSKY, A. S. Sinais e Sistemas. 2 end

ed. São Paulo: Pearson Education do Brasil LTDA. 2010.
DC_offset=(1/N)*phasors(1,:)
PINTO, T. H. L. Processamento de sinais aplicado à Fund_harmonic_phasors = (sqrt(2)/N)*phasors((2:end),:)
análise comparativa do timbre de instrumentos
Fourier_coefficients = conj(2*(1/N)*phasors((2:end),:))
produzidos com madeiras nacionais e tradicionais. Anais
do SIMMEC 2016 ABMEC. 2016. Disponível
em:http://acervo.ufvjm.edu.br/jspui/bitstream/1/1553/1/simm
stem(Fourier_coefficients)
ec_processamento_sinais.pdf. Acesso em 03 out. 2021.
figure
stem(Fund_harmonic_phasors)