Relatorio Filtros Digitais - Analise e Implementacao com Scilab

Processamento Digital de Sinais Prof Eng Viviana R Zurro MSc M11 Filtros digitais Roteiro 1 Atividade Prática Filtros Digitais OBJETIVO Verificar o sinal de saída a partir de um sinal de entrada determinado MATERIAL UTILIZADO Ambiente matemático Scilab Apostilas disponíveis no AVA na Aula Ambiente Matemático Scilab Listas de Exercícios disponíveis em todas as aulas do AVA ORIENTAÇÕES Para realizar esta atividade leia atentamente e estude todo o material disponível no AVA principalmente as apostilas e listas de exercícios Atenção Coloque no relatório todo o desenvolvimento matemático prévio ao desenvolvimento do algoritmo Se não houver desenvolvimento matemático como indicado no vídeo do experimento será descontada nota Inclua imagens de todos os procedimentos solicitados Nas imagens não se esqueça de colocar nomes nos eixos xlabel e ylabel Será descontada nota Para facilitar o desenvolvimento da atividade use o aplicativo SciNotes que permite gravar sua atividade como um programa página 7 da Apostila 1 Introdução ao Scilab Coloque o algoritmo completo no relatório com o detalhe de cada uma das linhas como o exemplo indicado Será descontada nota Trabalhos iguais serão considerados plágio e a nota será zero para todos os alunos que entregarem o mesmo trabalho ATIVIDADE Um determinado sistema digital tem a seguinte estrutura Figura 1 Estrutura do filtro Processamento Digital de Sinais Prof Eng Viviana R Zurro MSc M11 Filtros digitais Roteiro 2 Para chegar na resposta ao impulso será necessário calcular a função do sistema 𝐻𝑧 a partir da estrutura do filtro Esta atividade deve ser desenvolvida considerando um RU de 7 números Se seu RU tiver menos de 7 números deverá preencher com zeros os últimos números Exemplo RU 12345 RU1 RU2 RU3 RU4 RU5 RU6 RU7 1 2 3 4 5 0 0 Se seu RU tiver mais de 7 números deverá desconsiderar os últimos números Exemplo RU 123456789 RU1 RU2 RU3 RU4 RU5 RU6 RU7 1 2 3 4 5 6 7 Sendo 𝑎 𝑅𝑈1 𝑏 𝑅𝑈210 se 𝑅𝑈2 0 𝑏 02 𝑐 𝑅𝑈3 se 𝑅𝑈3 0 𝑐 3 𝑑 𝑅𝑈4 se 𝑅𝑈4 0 𝑑 4 Vetor de entrada sinal de entrada 𝑥𝑛 𝑒𝑅𝑈3𝑛 10 2𝑛 10 𝑐𝑜𝑠 04𝑘𝑛 𝜋 4 𝑅𝑈4 𝑛 4 𝜋 número 𝜋 sintaxe no Scinotes pi 𝑘 𝑅𝑈2 9 se 𝑅𝑈2 0 adotar 𝑘 022 Usando o ambiente matemático Scilab calcular o sinal de saída 𝑦𝑛 definido por 𝑦𝑛 𝑥𝑛 ℎ𝑛 1 Onde ℎ𝑛 é a resposta ao impulso domínio do tempo correspondente à 𝐻𝑧 domínio da frequência E o sinal de saída 𝑦𝑛 é resultante da convolução entre o sinal de entrada 𝑥𝑛 e a resposta ao impulso do sistema ℎ𝑛 PROCEDIMENTO É conveniente usar o aplicativo SciNotes para escrever os comandos As funções impulso unitário e degrau unitário explicadas na Apostila 1 Introdução ao Scilab serão fundamentais para esta atividade Elas deverão ser definidas no início da série de comandos Processamento Digital de Sinais Prof Eng Viviana R Zurro MSc M11 Filtros digitais Roteiro 3 1 6 pontos A partir da estrutura do filtro mostrada na Figura 1 calcular a resposta ao impulso ℎ𝑛 do sistema Todos os cálculos deverão ser apresentados detalhadamente no relatório estes cálculos serão realizados no caderno As listas de exercícios mencionadas no MATERIAL UTILIZADO têm vários problemas e resoluções similares a 05p Estrutura do filtro b 15p Procedimento c 3p Cálculo dos parâmetros d 1p Resposta ao impulso ℎ𝑛 correta 2 Scilab a Gerar um vetor 𝑛 entre 20 e 20 para 𝑥𝑛 e ℎ𝑛 b Gerar um vetor 𝑛1 entre 40 e 40 para 𝑦𝑛 c 3 pontos Algoritmo matemático detalhado de preferência no SciNotes i 15p Vetor de entrada 𝑥𝑛 correto ii 15p Resposta ao impulso ℎ𝑛 correta 3 1 ponto Usando os comandos subplot e plot2d3 Apostila 2 plotar 𝑥𝑛 ℎ𝑛 e 𝑦𝑛 no mesmo gráfico a Gráficos corretos b Nomes nos eixos dos gráficos será descontada nota se os gráficos não tiverem nome nos eixos Exemplo RU 1234567 RU1 RU2 RU3 RU4 RU5 RU6 RU7 1 2 3 4 5 6 7 𝑥𝑛 𝑒03𝑛 2𝑛 10 𝑐𝑜𝑠 0088𝑛 𝜋 4 4 𝑛 4 Processamento Digital de Sinais Prof Eng Viviana R Zurro MSc M11 Filtros digitais Roteiro 4 Gráficos Processamento Digital de sinais Prof Eng Viviana R Zurro MSc Atividade Prática Filtros Digitais Célio 1660242 ATIVIDADE Um determinado sistema digital tem a seguinte estrutura Figura 1 Estrutura do filtro Calcular a função do sistema 𝐻𝑧 a partir da estrutura do filtro RU1 RU2 RU3 RU4 RU5 RU6 RU7 1 6 6 0 2 4 2 x n 1660242 Sendo o aRU 1 o bRU 210 se for zero adotar 02 o cRU 3 se for zero adotar 3 o dRU 4 se for zero adotar 4 Vetor de entrada sinal de entrada o x n e RU 3n102 n10cos04kn π 4RU 4n4 TrabalhoDigitaldocx 1 Processamento Digital de sinais Prof Eng Viviana R Zurro MSc kRU 2 9 se for zeroadotar k022 Função de transferência do filtro Desenvolvimento matemático A partir da estrutura do filtro mostrada na Figura 1 calcular a resposta ao impulso ℎ𝑛 do sistema Todos os cálculos deverão ser apresentados detalhadamente no relatório estes cálculos serão realizados no caderno As listas de exercícios mencionadas no MATERIAL UTILIZADO têm vários problemas e resoluções similares a Estrutura do filtro Inicialmente é importante verificar o caminho percorrido pelos sinais e os locais onde estes se somam TrabalhoDigitaldocx 2 Processamento Digital de sinais Prof Eng Viviana R Zurro MSc b Procedimento Observando o diagrama de blocos podese observar o sinal W z W zX zaX zb z 1W z z 1 W zX za z 1X zb z 2W z W zb z 2W zX za z 1 X z W z1b z 2X z1a z 1 W zX z 1a z 1 1b z 2 Observando o sinal R z temse R z cX z dX z bR z z 1 R zcX z d z 1 X zb z 1R z R zb z 1R zcX zd z 1 X z R z1bz 1X zcd z 1 R zX z cd z 1 1b z 1 TrabalhoDigitaldocx 3 Processamento Digital de sinais Prof Eng Viviana R Zurro MSc A saída é dada pela soma dos dois sinais Y z W zR z Y z X z 1a z 1 1b z 2X z cd z 1 1b z 1 Y z X z 1a z 1 1b z 2 cd z 1 1b z 1 Y z X z 1a z 1 1b z 2 cd z 1 1b z 1 H z 1a z 1 1b z 2 cd z 1 1b z 1 c Calculo dos parâmetros H z z 2az z 2b czd zb H z z za z 2b c z zb d zb H z z za zb zb c z zb d zb Por frações parciais temse za zb zb A zb B zb A zbB zbza z AB b Ab Bza AB1 b Ab Ba AB1 ABab TrabalhoDigitaldocx 4 Processamento Digital de sinais Prof Eng Viviana R Zurro MSc 2 A1 a b A05 a 2b A05ab 2b B1A B105 ab 2b B05 ab 2b H z z 05ab 2b zb 05 ab 2b zb c z zb d zb H z 0506 12 z z06 05 06 12 z z06 6 z z06 4 1 z06 H z 01455 z z06 11455 z z06 6 z z06 4 1 z06 d Resposta ao impulso h n correta h n01455 06 n11455 06 n606 n4 06 n1 h n01455 07746 n11455 07746 n606 n406 n1 Algoritmo Limpase as variáveis comandos e janelas gráficas para garantir que o programa será implementado sem sujeiras de outros códigos e inicializase a variável para manipular os gráficos TrabalhoDigitaldocx 5 Processamento Digital de sinais Prof Eng Viviana R Zurro MSc clear Limpa as variáveis clc Limpa a janela de comandos clf Limpa as janelas gráficas fgcf Manipulador da figura As funções de impulso e degrau ajudam a implementar a entrada saída e função hn do sistema As funções podem ser vistas a seguir function yimpx Função impulso yzeros 1lengthx yfindx01 endfunction function ydegx Função degrau yzeros1lengthx yfindx01 endfunction Para inicializar as variáveis é importante escrever os valores de matrícula com os valores de matrícula inicializados podese calcular as variáveis a b c e d do exercício RU11 Inicia os valores de matricula RU26 RU36 RU40 RU52 RU64 RU72 aRU1 Inicia variáveis do exercício bRU210 if b0 then b02 end cRU3 if c0 then c3 end dRU4 if d0 then d4 end O sistema necessita de dois vetores temporais um para os sinais de entrada e hn e outro para a saída pois este é obtido a partir da convolução que retorna a quantidade de pontos de xn somado a quantidade de pontos de hn subtraído de um Além das variáveis temporais serão utilizadas as contantes e π e k dado no exercício TrabalhoDigitaldocx 6 Processamento Digital de sinais Prof Eng Viviana R Zurro MSc n20120 n para xn e hn n140140 n para yn pipi pi ee e kRU29 variável k if k0 then k022 end A função xn dada no exercício é válida no intervalo RU 4n4 para garantir isto escrevese a janela a seguir janeladegnRU4degn4 Janela em que xn é válido A função x n então é escrita e multiplicada pela janela para garantir o intervalo onde ela é válida xeRU310n2n10cos04knpi4 xn dado no exercício xxjanela xn no intervalo mencionado A função h n é escrita utilizando o degrau e os valores de variável encontrados matematicamente O y n então é calculado utilizando a função conv que retorna a convolução de dois sinais udegn un é a função degrau for i2020 Calcula hn hi2101455sqrt06iui21 11455sqrt06iui21 606iui21406i1ui21 end yconvxh Calcula yn Os gráficos do sinal de entrada x n da função h n e da saída y n são apresentados em uma mesma figura utilizando a função subplot e o plot2d3 A grossura da linha do gráfico é alterada para 2 e são escritos o título a legenda do eixo das ordenadas e do eixo das abscissas subplot311 subplot 1 plot2d3nx style2 plot discreto com cor azul fchildrenchildren1childrenthickness2linha grossura 2 titlexn xlabeln ylabelamplitude subplot312 subplot 2 TrabalhoDigitaldocx 7 Processamento Digital de sinais Prof Eng Viviana R Zurro MSc plot2d3nh style3 plot discreto com cor verde fchildrenchildren1childrenthickness2linha grossura 2 titlehn xlabeln ylabelamplitude subplot313 subplot 3 plot2d3n1y style5 plot discreto com cor vermelho fchildrenchildren1childrenthickness2linha grossura 2 titleyn xlabeln ylabelamplitude TrabalhoDigitaldocx 8 Processamento Digital de sinais Prof Eng Viviana R Zurro MSc Gráficos de entrada função de transferência e saída TrabalhoDigitaldocx 9 Processamento Digital de sinais Prof Eng Viviana R Zurro MSc Código completo clear Limpa as variáveis clc Limpa a janela de comandos clf Limpa as janelas gráficas fgcf Manipulador da figura function yimpx Função impulso yzeros 1lengthx yfindx01 endfunction function ydegx Função degrau yzeros1lengthx yfindx01 endfunction RU11 Inicia os valores de matricula RU26 RU36 RU40 RU52 RU64 RU72 aRU1 Inicia variáveis do exercício bRU210 if b0 then b02 end cRU3 if c0 then c3 end dRU4 if d0 then d4 end n20120 n para xn e hn n140140 n para yn pipi pi ee e kRU29 variável k if k0 then k022 end janeladegnRU4degn4 Janela em que xn é válido xn dado no exercício xeRU310n2n10cos04knpi4 xxjanela xn no intervalo mencionado udegn un é a função degrau for i2020 Calcula hn TrabalhoDigitaldocx 10 Processamento Digital de sinais Prof Eng Viviana R Zurro MSc hi2101455sqrt06iui21 11455sqrt06iui21 606iui21406i1ui21 end yconvxh Calcula yn subplot311 subplot 1 plot2d3nx style2 plot discreto com cor azul fchildrenchildren1childrenthickness2linha grossura 2 titlexn xlabeln ylabelamplitude subplot312 subplot 2 plot2d3nh style3 plot discreto com cor verde fchildrenchildren1childrenthickness2linha grossura 2 titlehn xlabeln ylabelamplitude subplot313 subplot 3 plot2d3n1y style5 plot discreto com cor vermelho fchildrenchildren1childrenthickness2linha grossura 2 titleyn xlabeln ylabelamplitude TrabalhoDigitaldocx 11 Processamento Digital de sinais Prof Eng Viviana R Zurro MSc TrabalhoDigital 1 Atividade Prática Filtros Digitais Célio 1660242 ATIVIDADE Um determinado sistema digital tem a seguinte estrutura Figura 1 Estrutura do filtro Calcular a função do sistema 𝐻𝑧 a partir da estrutura do filtro RU1 RU2 RU3 RU4 RU5 RU6 RU7 1 6 6 0 2 4 2 𝑥𝑛 1 6 6 0 2 4 2 Sendo o 𝑎 𝑅𝑈1 o 𝑏 𝑅𝑈210 se for zero adotar 02 o 𝑐 𝑅𝑈3 se for zero adotar 3 o 𝑑 𝑅𝑈4 se for zero adotar 4 Vetor de entrada sinal de entrada o 𝑥𝑛 𝑒𝑅𝑈3𝑛10 2𝑛10 cos 04𝑘𝑛 𝜋 4 𝑅𝑈4 𝑛 4 𝑘 𝑅𝑈2 9 𝑠𝑒 𝑓𝑜𝑟 𝑧𝑒𝑟𝑜 𝑎𝑑𝑜𝑡𝑎𝑟 𝑘 022 Processamento Digital de sinais Prof Eng Viviana R Zurro MSc TrabalhoDigital 2 Função de transferência do filtro Desenvolvimento matemático A partir da estrutura do filtro mostrada na Figura 1 calcular a resposta ao impulso ℎ𝑛 do sistema Todos os cálculos deverão ser apresentados detalhadamente no relatório estes cálculos serão realizados no caderno As listas de exercícios mencionadas no MATERIAL UTILIZADO têm vários problemas e resoluções similares a Estrutura do filtro Inicialmente é importante verificar o caminho percorrido pelos sinais e os locais onde estes se somam b Procedimento Processamento Digital de sinais Prof Eng Viviana R Zurro MSc TrabalhoDigital 3 Observando o diagrama de blocos podese observar o sinal 𝑊𝑧 𝑊𝑧 𝑋𝑧 𝑎𝑋𝑧 𝑏𝑧1𝑊𝑧𝑧1 𝑊𝑧 𝑋𝑧 𝑎𝑧1𝑋𝑧 𝑏𝑧2𝑊𝑧 𝑊𝑧 𝑏𝑧2𝑊𝑧 𝑋𝑧 𝑎𝑧1𝑋𝑧 𝑊𝑧1 𝑏𝑧2 𝑋𝑧1 𝑎𝑧1 𝑊𝑧 𝑋𝑧 1 𝑎𝑧1 1 𝑏𝑧2 Observando o sinal 𝑅𝑧 temse 𝑅𝑧 𝑐𝑋𝑧 𝑑𝑋𝑧 𝑏𝑅𝑧𝑧1 𝑅𝑧 𝑐𝑋𝑧 𝑑𝑧1𝑋𝑧 𝑏𝑧1𝑅𝑧 𝑅𝑧 𝑏𝑧1𝑅𝑧 𝑐𝑋𝑧 𝑑𝑧1𝑋𝑧 𝑅𝑧1 𝑏𝑧1 𝑋𝑧𝑐 𝑑𝑧1 𝑅𝑧 𝑋𝑧 𝑐 𝑑𝑧1 1 𝑏𝑧1 A saída é dada pela soma dos dois sinais 𝑌𝑧 𝑊𝑧 𝑅𝑧 𝑌𝑧 𝑋𝑧 1 𝑎𝑧1 1 𝑏𝑧2 𝑋𝑧 𝑐 𝑑𝑧1 1 𝑏𝑧1 𝑌𝑧 𝑋𝑧 1 𝑎𝑧1 1 𝑏𝑧2 𝑐 𝑑𝑧1 1 𝑏𝑧1 𝑌𝑧 𝑋𝑧 1 𝑎𝑧1 1 𝑏𝑧2 𝑐 𝑑𝑧1 1 𝑏𝑧1 𝐻𝑧 1 𝑎𝑧1 1 𝑏𝑧2 𝑐 𝑑𝑧1 1 𝑏𝑧1 c Calculo dos parâmetros Processamento Digital de sinais Prof Eng Viviana R Zurro MSc TrabalhoDigital 4 𝐻𝑧 𝑧2 𝑎𝑧 𝑧2 𝑏 𝑐𝑧 𝑑 𝑧 𝑏 𝐻𝑧 𝑧 𝑧 𝑎 𝑧2 𝑏 𝑐 𝑧 𝑧 𝑏 𝑑 𝑧 𝑏 𝐻𝑧 𝑧 𝑧 𝑎 𝑧 𝑏𝑧 𝑏 𝑐 𝑧 𝑧 𝑏 𝑑 𝑧 𝑏 Por frações parciais temse 𝑧 𝑎 𝑧 𝑏𝑧 𝑏 𝐴 𝑧 𝑏 𝐵 𝑧 𝑏 𝐴𝑧 𝑏 𝐵𝑧 𝑏 𝑧 𝑎 𝑧𝐴 𝐵 𝑏𝐴 𝑏𝐵 𝑧 𝑎 𝐴 𝐵 1 𝑏𝐴 𝑏𝐵 𝑎 𝐴 𝐵 1 𝐴 𝐵 𝑎𝑏 2𝐴 1 𝑎 𝑏 𝐴 05 𝑎 2𝑏 𝐴 05 𝑎𝑏 2𝑏 𝐵 1 𝐴 𝐵 1 05 𝑎𝑏 2𝑏 𝐵 05 𝑎𝑏 2𝑏 Processamento Digital de sinais Prof Eng Viviana R Zurro MSc TrabalhoDigital 5 𝐻𝑧 𝑧 05 𝑎𝑏 2𝑏 𝑧 𝑏 05 𝑎𝑏 2𝑏 𝑧 𝑏 𝑐 𝑧 𝑧 𝑏 𝑑 𝑧 𝑏 𝐻𝑧 05 06 12 𝑧 𝑧 06 05 06 12 𝑧 𝑧 06 6 𝑧 𝑧 06 4 1 𝑧 06 𝐻𝑧 01455 𝑧 𝑧 06 11455 𝑧 𝑧 06 6 𝑧 𝑧 06 4 1 𝑧 06 d Resposta ao impulso ℎ𝑛 correta 𝒉𝒏 𝟎 𝟏𝟒𝟓𝟓 𝟎 𝟔 𝒏 𝟏 𝟏𝟒𝟓𝟓 𝟎 𝟔 𝒏 𝟔 𝟎 𝟔𝒏 𝟒 𝟎 𝟔𝒏𝟏 𝒉𝒏 𝟎 𝟏𝟒𝟓𝟓 𝟎 𝟕𝟕𝟒𝟔𝒏 𝟏 𝟏𝟒𝟓𝟓 𝟎 𝟕𝟕𝟒𝟔𝒏 𝟔 𝟎 𝟔𝒏 𝟒 𝟎 𝟔𝒏𝟏 Algoritmo Limpase as variáveis comandos e janelas gráficas para garantir que o programa será implementado sem sujeiras de outros códigos e inicializase a variável para manipular os gráficos clear Limpa as variáveis clc Limpa a janela de comandos clf Limpa as janelas gráficas fgcf Manipulador da figura As funções de impulso e degrau ajudam a implementar a entrada saída e função ℎ𝑛 do sistema As funções podem ser vistas a seguir function yimpx Função impulso yzeros 1lengthx yfindx01 endfunction function ydegx Função degrau yzeros1lengthx yfindx01 endfunction Processamento Digital de sinais Prof Eng Viviana R Zurro MSc TrabalhoDigital 6 Para inicializar as variáveis é importante escrever os valores de matrícula com os valores de matrícula inicializados podese calcular as variáveis a b c e d do exercício RU11 Inicia os valores de matricula RU26 RU36 RU40 RU52 RU64 RU72 aRU1 Inicia variáveis do exercício bRU210 if b0 then b02 end cRU3 if c0 then c3 end dRU4 if d0 then d4 end O sistema necessita de dois vetores temporais um para os sinais de entrada e ℎ𝑛 e outro para a saída pois este é obtido a partir da convolução que retorna a quantidade de pontos de 𝑥𝑛 somado a quantidade de pontos de ℎ𝑛 subtraído de um Além das variáveis temporais serão utilizadas as contantes 𝑒 𝜋 e k dado no exercício n20120 n para xn e hn n140140 n para yn pipi pi ee e kRU29 variável k if k0 then k022 end A função 𝑥𝑛 dada no exercício é válida no intervalo 𝑅𝑈4 𝑛 4 para garantir isto escrevese a janela a seguir janeladegnRU4degn4 Janela em que xn é válido A função 𝑥𝑛 então é escrita e multiplicada pela janela para garantir o intervalo onde ela é válida xeRU310n2n10cos04knpi4 xn dado no exercício xxjanela xn no intervalo mencionado Processamento Digital de sinais Prof Eng Viviana R Zurro MSc TrabalhoDigital 7 A função ℎ𝑛 é escrita utilizando o degrau e os valores de variável encontrados matematicamente O 𝑦𝑛 então é calculado utilizando a função conv que retorna a convolução de dois sinais udegn un é a função degrau for i2020 Calcula hn hi2101455sqrt06iui21 11455sqrt06iui21 606iui21406i1ui21 end yconvxh Calcula yn Os gráficos do sinal de entrada 𝑥𝑛 da função ℎ𝑛 e da saída 𝑦𝑛 são apresentados em uma mesma figura utilizando a função subplot e o plot2d3 A grossura da linha do gráfico é alterada para 2 e são escritos o título a legenda do eixo das ordenadas e do eixo das abscissas subplot311 subplot 1 plot2d3nx style2 plot discreto com cor azul fchildrenchildren1childrenthickness2linha grossura 2 titlexn xlabeln ylabelamplitude subplot312 subplot 2 plot2d3nh style3 plot discreto com cor verde fchildrenchildren1childrenthickness2linha grossura 2 titlehn xlabeln ylabelamplitude subplot313 subplot 3 plot2d3n1y style5 plot discreto com cor vermelho fchildrenchildren1childrenthickness2linha grossura 2 titleyn xlabeln ylabelamplitude Processamento Digital de sinais Prof Eng Viviana R Zurro MSc TrabalhoDigital 8 Gráficos de entrada função de transferência e saída Processamento Digital de sinais Prof Eng Viviana R Zurro MSc TrabalhoDigital 9 Código completo clear Limpa as variáveis clc Limpa a janela de comandos clf Limpa as janelas gráficas fgcf Manipulador da figura function yimpx Função impulso yzeros 1lengthx yfindx01 endfunction function ydegx Função degrau yzeros1lengthx yfindx01 endfunction RU11 Inicia os valores de matricula RU26 RU36 RU40 RU52 RU64 RU72 aRU1 Inicia variáveis do exercício bRU210 if b0 then b02 end cRU3 if c0 then c3 end dRU4 if d0 then d4 end n20120 n para xn e hn n140140 n para yn pipi pi ee e kRU29 variável k if k0 then k022 end janeladegnRU4degn4 Janela em que xn é válido xn dado no exercício xeRU310n2n10cos04knpi4 xxjanela xn no intervalo mencionado udegn un é a função degrau for i2020 Calcula hn Processamento Digital de sinais Prof Eng Viviana R Zurro MSc TrabalhoDigital 10 hi2101455sqrt06iui21 11455sqrt06iui21 606iui21406i1ui21 end yconvxh Calcula yn subplot311 subplot 1 plot2d3nx style2 plot discreto com cor azul fchildrenchildren1childrenthickness2linha grossura 2 titlexn xlabeln ylabelamplitude subplot312 subplot 2 plot2d3nh style3 plot discreto com cor verde fchildrenchildren1childrenthickness2linha grossura 2 titlehn xlabeln ylabelamplitude subplot313 subplot 3 plot2d3n1y style5 plot discreto com cor vermelho fchildrenchildren1childrenthickness2linha grossura 2 titleyn xlabeln ylabelamplitude