O processamento de sinal de processamento de sinal é a arte e a ciência de modificar dados de séries temporais adquiridos para fins de análise ou aprimoramento. Os exemplos incluem análise espectral (usando Fast Fourier ou outras transformações) e aprimorando dados adquiridos usando filtragem digital. O Igor é ideal para o processamento de sinal devido ao seu forte suporte para longos dados de séries temporais (ou quotwaveformquot). E porque seus muitos comandos internos de processamento de sinal podem ser facilmente usados em diálogos simples. Além disso, a linguagem de programação Igoracutes torna direto implementar qualquer tipo de algoritmo de processamento de sinal personalizado, muito auxiliado pelo poder de Igoracutes Fourier (e outras). Igor usa o algoritmo de Transformação de Fourier Rápida (FFT) para calcular uma Transformação de Fourier Discreta (DFT). A FFT pode ser usada para caracterizar a magnitude e a fase de um sinal, ou pode ser usada em combinação com outras operações para realizar cálculos mais envolvidos, como convolução ou correlação. A computação FFT presume que os dados de entrada se repetem repetidamente. Isso é importante quando os valores inicial e final de seus dados não são os mesmos: a descontinuidade causa aberrações no espectro computado pela FFT. QuotWindowingquot suaviza as extremidades dos dados para eliminar essas aberrações. QuotPower Spectraquot responde a pergunta que as frequências contêm o signalacutes powerquot. A resposta está na forma de uma distribuição de valores de potência como uma função de freqüência, onde quotpowerquot é considerado a média dos sinaisup2. No domínio da frequência, este é o quadrado da magnitude FFTacutes. Os espectros de potência podem ser calculados para todo o sinal ao mesmo tempo (um quotperiodogramquot) ou periodogramas de segmentos do sinal de tempo podem ser calculados em média para formar a densidade espectral do quotpower. A Hilbert Transform calcula um sinal de domínio do tempo que é 90 graus fora de fase com o sinal de entrada. As aplicações unidimensionais incluem a computação do envelope de um sinal modulado e a medição da taxa de decaimento de um sinusóide de decomposição exponencial freqüentemente encontrado em sistemas lineares e não lineares sem amortiguação. Quando você calcula o espectro de Fourier (ou Power Spectra) de um sinal, você descarta toda a informação de fase contida na transformada de Fourier. Você pode descobrir quais freqüências um sinal contém, mas você não sabe quando essas freqüências aparecem no sinal. Por exemplo, considere o sinal: a representação espectral de f (t) permanece essencialmente inalterada se trocamos as duas frequências f 1 e f 2. Claramente, o espectro de Fourier não é a melhor ferramenta de análise para sinais cujos espectros flutuam no tempo. Uma solução para este problema é o chamado quotShort-time Fourier Transformquot (ou quotSonogramquot) em que você pode calcular os espectros de Fourier usando uma janela temporal deslizante. Ao ajustar a largura da janela, você pode determinar a resolução do tempo dos espectros resultantes. Você pode usar a convolução para calcular a resposta de um sistema linear para um sinal de entrada. O sistema linear é definido por sua resposta de impulso. A convolução do sinal de entrada e a resposta ao impulso são a resposta do sinal de saída. A filtragem digital é realizada através da definição de uma resposta de impulso de sistema linear que, quando convolvida com o sinal, realiza o resultado desejado (filtro de passagem baixa ou passagem alta). O algoritmo de correlação é muito similar matematicamente à convolução, mas é usado para diferentes fins. É usado com mais freqüência para identificar o atraso no tempo em que dois sinais são subitáveis, ou são quase idênticos. O alisamento remove as variações de curto prazo ou quotnoisequot para revelar a importante forma subjacente dos dados. A forma mais simples de suavização é a quotmoving quot da média que simplesmente substitui cada valor de dados pela média de valores vizinhos. (Outros termos para esse tipo de alisamento são quotsliding averagequot, quotbox smoothingquot ou quotboxcar smoothingquot.) Igoracutes A operação suave executa o alisamento de caixas, o alinhamento quotbinomial (Gaussian) e o suavizado Savitzky-Golay (polinômico). Os diferentes algoritmos de suavização calculam médias ponderadas que multiplicam valores vizinhos por diferentes pesos ou quotcoeficientes para calcular o valor suavizado. Os filtros digitais são uma ferramenta natural quando os dados já estão digitalizados. As razões para a aplicação de filtragem digital para os dados incluem: Eliminação de componentes de sinal indesejados (quotnoisequot) Melhorando os componentes de sinal desejados Detectando a presença de certos sinais Simulação de sistemas lineares (computa o sinal de saída dado o sinal de entrada e o sistema possui uma função quottransfer). Filtros digitais geralmente Venha em dois sabores: Filtro de Impulso Finito (FIR) e Infinito de Impulso Resposta (IIR). Igor implementa o filtro digital FIR principalmente através da convolução do domínio do tempo usando os comandos Smooth ou SmoothCustom. (Apesar do nome de itacutes, o SmoothCustom convolve dados com coeficientes de filtro fornecidos pelo usuário para implementar qualquer tipo de filtro FIR, passagem baixa, passe alto, passagem de banda, etc.) O design dos coeficientes de filtro FIR usados com o SmoothCustom é mais Facilmente realizado usando o Igor Filter Design Laboratory (um produto separado que também requer Igor Pro). Os filtros digitais IIR são projetados e aplicados em dados usando IFDL. A detecção de nível é o processo de localização da coordenada X em que seus dados passam ou atingem um dado valor Y. Isso às vezes é chamado de interpolação quotinverse. Dito de outra forma, a detecção de nível responde a pergunta: quotgiven um nível Y, qual é o valor X correspondente Igor fornece dois tipos de respostas a essa pergunta. Uma resposta assume que seus dados Y são uma lista de valores Y únicos que aumentam ou diminuem monotonicamente. A outra resposta assume que seus dados Y variam de forma irregular, como seria com os dados adquiridos. Neste caso, pode haver múltiplos valores de X que atravessam o nível Y. Exemplos importantes disso são estatísticas de ponta e pulso. Uma questão relacionada, mas diferente, é dada uma função y f (x), encontre x onde y é zero (ou algum outro valor). Esta pergunta é respondida pela operação FindRoots. Documentação Este exemplo mostra como usar filtros médios móveis e reescrever para isolar o efeito de componentes periódicos da hora do dia em leituras horárias horárias, bem como remover o ruído indesejado da linha de um circuito aberto Medição de tensão. O exemplo também mostra como alisar os níveis de um sinal de relógio, preservando as bordas usando um filtro mediano. O exemplo também mostra como usar um filtro Hampel para remover grandes outliers. Motivation Smoothing é como descobrimos padrões importantes em nossos dados, deixando para fora coisas que não têm importância (ou seja, o ruído). Usamos filtragem para executar esse alisamento. O objetivo do alisamento é produzir mudanças lentas de valor para que seja mais fácil ver tendências em nossos dados. Às vezes, quando você examina dados de entrada, você deseja suavizar os dados para ver uma tendência no sinal. No nosso exemplo, temos um conjunto de leituras de temperatura em Celsius tomadas a cada hora no Aeroporto de Logan durante todo o mês de janeiro de 2017. Note que podemos visualizar visualmente o efeito que a hora do dia tem nas leituras de temperatura. Se você estiver interessado apenas na variação diária da temperatura ao longo do mês, as flutuações horárias só contribuem com o ruído, o que pode tornar as variações diárias difíceis de discernir. Para remover o efeito da hora do dia, gostaríamos agora de suavizar nossos dados usando um filtro de média móvel. Um filtro de média móvel Na sua forma mais simples, um filtro médio móvel de comprimento N leva a média de cada N amostras consecutivas da forma de onda. Para aplicar um filtro de média móvel a cada ponto de dados, construímos nossos coeficientes de nosso filtro de modo que cada ponto seja ponderado e contribua 124 para a média total. Isso nos dá a temperatura média em cada período de 24 horas. Retardamento do filtro Observe que a saída filtrada está atrasada em cerca de doze horas. Isto é devido ao fato de nosso filtro de média móvel ter um atraso. Qualquer filtro simétrico de comprimento N terá um atraso de (N-1) 2 amostras. Podemos explicar esse atraso manualmente. Extraindo diferenças médias Alternativamente, também podemos usar o filtro de média móvel para obter uma melhor estimativa de como a hora do dia afeta a temperatura geral. Para fazer isso, primeiro, subtrair os dados suavizados das medidas de temperatura horária. Em seguida, segmentar os dados diferenciados em dias e ter a média durante os 31 dias do mês. Extraindo o envelope máximo Às vezes, também gostaríamos de ter uma estimativa de variação suave de como os altos e baixos do nosso sinal de temperatura mudam diariamente. Para fazer isso, podemos usar a função envelope para conectar altas e baixas extremas detectadas em um subconjunto do período de 24 horas. Neste exemplo, garantimos que haja pelo menos 16 horas entre cada extremo alto e extremo baixo. Nós também podemos ter uma sensação de como os altos e baixos estão tendendo tomando a média entre os dois extremos. Filtros médios em movimento ponderados Outros tipos de filtros médios móveis não pesam cada amostra de forma igual. Outro filtro comum segue a expansão binomial de (12,12) n Este tipo de filtro se aproxima de uma curva normal para valores grandes de n. É útil para filtrar o ruído de alta freqüência para pequenos n. Para encontrar os coeficientes para o filtro binomial, convulsionar 12 12 com ele próprio e, então, combinar iterativamente a saída com 12 12 um número de vezes prescrito. Neste exemplo, use cinco iterações totais. Outro filtro um pouco semelhante ao filtro de expansão gaussiano é o filtro exponencial de média móvel. Este tipo de filtro de média móvel ponderada é fácil de construir e não requer um grande tamanho de janela. Você ajusta um filtro de média móvel ponderada exponencialmente por um parâmetro alfa entre zero e um. Um valor maior de alfa terá menos suavização. Amplie as leituras por um dia. Escolha o seu país
No comments:
Post a Comment