Você realmente precisa de uma taxa de amostragem de pouco mais de 2 kHz para amostrar ondas senoidais de 1 kHz corretamente. É
não
f N ≤ f S / 2
fN<fS/2
fN≤fS/2
PS Se você levou seu sinal para um espaço complexo, onde um sinusóide é da forma
onde t é tempo, A é amplitude, f é frequência e θ é deslocamento de fase,
f N
v(t)=Aej(2πft−θ)=A(cos(2πft−θ)+jsin(2πft−θ))
é o ponto em que a frequência "se dobra", ou seja, você não pode distinguir
fde
-f. Após a amostragem, outros aumentos na frequência parecerão subtraídos da frequência de amostragem, no caso de um sinusóide puro.
fN=fS/2
Não Sinusoides
No caso de uma onda quadrada a 1 kHz com um ciclo de trabalho menor ou igual a 10%, amostrado a 10 kHz, você está entendendo mal a entrada.
Primeiro, você precisa decompor sua forma de onda em uma série de Fourier para descobrir quais são as amplitudes das harmônicas dos componentes. Você provavelmente ficará surpreso ao saber que os harmônicos desse sinal ultrapassam os 5 kHz! (A regra geral do terceiro harmônico é 1/3 tão forte quanto o fundamental, e o quinto é 1/5 do fundamental, somente se aplica a ondas quadradas de 50% do ciclo de serviço .)
A regra básica para um sinal de comunicação é que sua largura de banda complexa é igual à inversa do tempo do menor pulso; portanto, neste caso, você está procurando uma largura de banda mínima de 10 kHz (-5 kHz a 5 kHz) para um ciclo de trabalho de 10% com o fundamental em 1 kHz (ou seja, 10 kbps).
Então, o que vai arruiná-lo é que esses fortes harmônicos de ordem superior se dobram e interferem (construtiva ou destrutivamente) com os harmônicos da banda; portanto, é perfeitamente esperado que você não obtenha uma boa amostra porque há muita informação fora do Nyquist banda.