Amostragem de sinais e teorema de Nyquist

Defina uma frequência de sinal e uma taxa de amostragem para ver se o sinal é capturado fielmente ou sofre aliasing.

A carregar a simulação interativa...

Duas frequências diferentes podem deixar amostras idênticas

A amostragem só funciona quando a taxa de amostragem captura a oscilação mais rápida presente. Se fs < 2 fsinal, frequências contínuas diferentes produzem amostras idênticas, e a reconstrução se fixa em uma frequência de alias falsa mais baixa.

Pelo menos dois pontos por onda

Para fixar uma onda é preciso captá-la subindo e descendo: pelo menos duas amostras em cada ciclo. Com menos, esses mesmos pontos poderiam pertencer a muitas ondas diferentes e a cópia digital adivinha errado. O áudio de CD usa 44,100 amostras por segundo porque a audição humana chega a cerca de 20 kHz, e 2 × 20 kHz deixa uma margem confortável.

Por que as rodas giram ao contrário no cinema

O efeito roda de carroça é o Nyquist no mundo visível. Uma câmera de cinema amostra o movimento a 24 quadros por segundo; se uma roda com raios gira um pouco menos que o vão entre raios entre dois quadros, cada tomada fica um pouco atrás da anterior e o olho vê a roda girando devagar para trás. A mesma matemática do aliasing de áudio: a taxa de quadros é a taxa de amostragem, e uma rotação rápida se dobra numa falsa rotação lenta.

Problemas resolvidos na íntegra

  1. A regra dos 6 dB por bit e o 1,76 isolado para um conversor de 8 bits 5 passos

    Deduza a regra dos 6 dB por bit de raiz para um conversor de 8 bits — de onde vem o 1,76 isolado? — e depois diga o que ela não promete.

    1. O arredondamento para o nível mais próximo deixa um erro não superior a meio passo e, num sinal complexo, esse erro distribui-se uniformemente pelo intervalo. A sua variância é o integral de e² sobre uma distribuição uniforme.

      Δ=VFS2N,e∼U ⁣(−Δ2,Δ2),σe2=1Δ∫−Δ/2Δ/2 ⁣e2 de=Δ212\Delta=\frac{V_{FS}}{2^{N}},\qquad e\sim\mathcal{U}\!\left(-\tfrac{\Delta}{2},\tfrac{\Delta}{2}\right),\qquad \sigma_e^{2}=\frac{1}{\Delta}\int_{-\Delta/2}^{\Delta/2}\!e^{2}\,de=\frac{\Delta^{2}}{12}
    2. O sinal de referência é uma senoide que preenche toda a gama de entrada, pelo que a sua potência é metade do quadrado da sua amplitude.

      Psig=(VFS/2)22=VFS28P_{\text{sig}}=\frac{(V_{FS}/2)^{2}}{2}=\frac{V_{FS}^{2}}{8}
    3. Divida os dois. A gama de fundo de escala cancela-se, razão pela qual o resultado é uma função pura do número de bits e não da tensão.

      SNR=VFS2/8Δ2/12=128(VFSΔ)2=32⋅22N\mathrm{SNR}=\frac{V_{FS}^{2}/8}{\Delta^{2}/12}=\frac{12}{8}\left(\frac{V_{FS}}{\Delta}\right)^{2}=\frac{3}{2}\cdot2^{2N}
    4. Aplique logaritmos. O 6,02 é 20 log₁₀2 — um bit duplica o número de níveis — e o 1,76 é o 3/2 restante da linha anterior.

      SNRdB=10log⁡10 ⁣32+20Nlog⁡102=1.76+6.02 N\mathrm{SNR_{dB}}=10\log_{10}\!\tfrac{3}{2}+20N\log_{10}2=1.76+6.02\,N
    5. Calcule para 8 bits, depois para 16. Cada bit adicional vale um pouco mais de 6 dB, pelo que os oito bits adicionais do áudio de CD compram 48 dB.

      N=8: 6.02(8)+1.76=49.92 dB;N=16: 98.08 dBN=8:\ 6.02(8)+1.76=49.92\ \mathrm{dB};\qquad N=16:\ 98.08\ \mathrm{dB}

    Resposta

    49,92 dB, e o 1,76 é 10 log₁₀(3/2). É a razão entre o valor médio quadrático de uma senoide e a variância do erro uniforme do quantizador, e é a única parte da fórmula que não é apenas uma contagem de bits. Note-se o que nunca aparece: a frequência de amostragem. Um conversor de 24 bits a funcionar a 8 kHz tem 146 dB de resolução em amplitude e, ainda assim, não consegue representar um tom de 5 kHz — os bits compram precisão vertical, a frequência de amostragem compra largura de banda, e nenhuma quantidade de um substitui o outro. É por isso que os dois seletores deslizantes nesta página são eixos genuinamente independentes, em vez de duas maneiras de dizer "qualidade".

    SNRdB=6.02 N+1.76=49.92 dB  (N=8)\mathrm{SNR_{dB}}=6.02\,N+1.76=49.92\ \mathrm{dB}\ \ (N=8)
  2. Um tom de 30 kHz num CD, e onde ele reaparece 7 passos

    Um microfone capta uma componente de 30 kHz, bem acima do que qualquer pessoa consegue ouvir, e esta chega a um conversor de 44,1 kHz sem passar por qualquer filtro. O que fica na gravação e de que outros tons passa a ser indistinguível?

    Abrir este problema na ferramenta

    1. Comece pelo limite superior. Um conversor a operar a 44 100 amostras por segundo consegue representar com fidelidade qualquer frequência abaixo de metade dessa taxa. O nosso tom está acima do limite, pelo que, independentemente do aspeto final das amostras, elas não podem significar aquilo que aparentam significar.

      fNyquist=fs2=44 1002=22 050 Hz30 000>22 050f_{Nyquist}=\frac{f_s}{2}=\frac{44\,100}{2}=22\,050\ \mathrm{Hz}\qquad 30\,000>22\,050
    2. Registe o que o conversor realmente armazena. Ele nunca vê a onda. Vê o valor da onda uma vez a cada 1/44 100 de segundo, e nada mais.

      x[n]=sin⁡ ⁣(2πf nfs)=sin⁡ ⁣(2π⋅30 00044 100 n)x[n]=\sin\!\left(2\pi f\,\frac{n}{f_s}\right)=\sin\!\left(2\pi\cdot\frac{30\,000}{44\,100}\,n\right)
    3. Adicione agora um número inteiro de ciclos entre uma amostra e a seguinte. Cada amostra ganha 2πk adicionais de fase, e uma sinusoide permanece inalterada por voltas completas. Os números armazenados saem idênticos, dígito por dígito.

      sin⁡ ⁣(2π f+kfsfs n)=sin⁡ ⁣(2πffsn+2πkn)=sin⁡ ⁣(2πffsn)k∈Z\sin\!\left(2\pi\,\frac{f+k f_s}{f_s}\,n\right)=\sin\!\left(2\pi\frac{f}{f_s}n+2\pi k n\right)=\sin\!\left(2\pi\frac{f}{f_s}n\right)\qquad k\in\mathbb{Z}
    4. Portanto, as amostras não identificam uma frequência. Identificam uma família, e a reprodução tem de escolher um elemento dessa família. Escolhe aquele que estiver mais próximo de zero, que é a distância de 30 000 Hz ao múltiplo mais próximo da taxa de amostragem.

      k=round⁡ ⁣(30 00044 100)=1,falias=∣30 000−1⋅44 100∣=14 100 Hzk=\operatorname{round}\!\left(\frac{30\,000}{44\,100}\right)=1,\qquad f_{alias}=\left|30\,000-1\cdot44\,100\right|=14\,100\ \mathrm{Hz}
    5. Esse é o valor que a ferramenta apresenta. 14 100 Hz é um silvo fino e agudo, audível para a maioria dos ouvintes jovens e uma das primeiras coisas a perder-se com a idade, e não estava na sala: uma componente inaudível foi transformada numa componente audível pelo simples ato de a gravar.

      OSR=fs2fsignal=44 1002⋅30 000=0.735<1\mathrm{OSR}=\frac{f_s}{2f_{signal}}=\frac{44\,100}{2\cdot30\,000}=0.735<1
    6. Que outros tons coincidem na mesma posição? Todos os membros da família. E existem infinitos. Aqueles que se encontram uma taxa de amostragem inteira acima reaparecem dígito por dígito; os espelhados apresentam todas as amostras de sinal trocado, o que equivale ao mesmo tom com a fase invertida. Em ambos os casos, a gravação preserva apenas um número por amostra e não guarda memória do tom que o produziu.

      ∣k fs±14 100∣={14 100, 30 000, 58 200, 74 100, 102 300, … }\left|k\,f_s\pm14\,100\right|=\{14\,100,\ 30\,000,\ 58\,200,\ 74\,100,\ 102\,300,\ \dots\}
    7. É por isso que o filtro tem de ser analógico e tem de ficar antes do conversor, o que fixa o quão abrupto deve ser: deixar passar 20 kHz, cortar aos 22 050 Hz, um intervalo de cerca de um sétimo de oitava. Construir isto em eletrónica analógica é difícil, razão prática pela qual existem conversores com sobreamostragem. Amostrando mais depressa, a banda de guarda alarga-se, o filtro analógico pode ser suave e a parte mais exigente da tarefa transfere-se para a aritmética, onde o custo é reduzido.

      log⁡222 05020 000=0.141 octave\log_{2}\frac{22\,050}{20\,000}=0.141\ \text{octave}

    Resposta

    14 100 Hz, com intensidade máxima. Um tom de 30 kHz amostrado a 44,1 kHz é indistinguível de um tom de 14 100 Hz, assim como o são 58 200 Hz, 74 100 Hz e qualquer outro membro de |k·fs ± 14 100|. O aliasing destrói a informação no momento da amostragem, pelo que nenhum processamento posterior consegue voltar a separar os tons rebatidos. O único local para o impedir é antes do conversor.

    falias=∣f−fsround⁡ ⁣(ffs)∣=14 100 Hzf_{alias}=\left|f-f_s\operatorname{round}\!\left(\tfrac{f}{f_s}\right)\right|=14\,100\ \mathrm{Hz}
Referências (2)