Moving Average Fft Matlab

Um modo simples (ad hoc) é apenas tomar uma média ponderada (ajustável por alfa) em cada ponto com seus vizinhos: ou alguma variação do mesmo. Sim, para ser mais sofisticado, Fourier pode transformar seus dados primeiro, depois cortar as altas freqüências. Algo como: Isso corta as 20 freqüências mais altas. Tenha cuidado para cortá-los simetricamente, caso contrário, a transformada inversa não é mais real. Você precisa escolher cuidadosamente a freqüência de corte para o nível correto de suavização. Este é um tipo de filtragem muito simples (caixa de filtragem no domínio da frequência), para que você possa tentar suavemente atenuar as freqüências de alta ordem se a distorção for inaceitável. Respondeu 4 de outubro 09 às 9:16 FFT não é uma má idéia, mas provavelmente é exagerado aqui. As médias em execução ou em movimento dão resultados geralmente ruins e devem ser evitadas para qualquer coisa, além da lição de casa tardia (e ruído branco). Eu uso filtragem Savitzky-Golay (em Matlab sgolayfilt (.)). Isso lhe dará os melhores resultados para o que você procura - algum alisamento local, mantendo a forma da curva. Resposta de freqüência do filtro médio de corrida A resposta de freqüência de um sistema LTI é o DTFT da resposta de impulso, a resposta de impulso De uma média móvel em L é uma vez que o filtro de média móvel é FIR, a resposta de freqüência reduz-se à soma finita. Podemos usar a identidade muito útil para escrever a resposta de freqüência como onde nós deixamos ae menos jomega. N 0 e M L menos 1. Podemos estar interessados ​​na magnitude desta função, a fim de determinar quais freqüências obtêm o filtro desatualizado e atenuados. Abaixo está um gráfico da magnitude desta função para L 4 (vermelho), 8 (verde) e 16 (azul). O eixo horizontal varia de zero a pi radianes por amostra. Observe que em todos os três casos, a resposta de freqüência possui uma característica de passagem baixa. Um componente constante (zero freqüência) na entrada passa pelo filtro não atenuado. Certas freqüências mais altas, como pi 2, são completamente eliminadas pelo filtro. No entanto, se a intenção era projetar um filtro de passagem baixa, então não fizemos muito bem. Algumas das frequências mais altas são atenuadas apenas por um fator de cerca de 1 10 (para a média móvel de 16 pontos) ou 1 3 (para a média móvel de quatro pontos). Podemos fazer muito melhor do que isso. A trama acima foi criada pelo seguinte código Matlab: omega 0: pi 400: pi H4 (1 4) (1-exp (-omeome4)). (1-exp (-maomega)) H8 (1 8) (1-exp (-iomega8)). (1-exp (-maomega)) H16 (1 16) (1-exp (-omeome16)). (1-exp (-maome)) trama (omega, abs (H4) abs (H8) abs (H16)) eixo (0, pi, 0, 1) Copyright cópia 2000- - Universidade da Califórnia, Berkeley