为什么我的正弦波频率扫描不正确？

Question

我正在尝试创建一个简单的WAV文件，它会发出一个不断变化的音调。但是，写入文件的波形与get_sample返回的数据不对应。在

我希望音调以对数方式变化，从A10（28160赫兹）到A0（27.5赫兹）结束。每过一秒，音高应平稳下降一个倍频程。在

实际发生的事情很难解释。语调变了，但方式却出乎意料。而让我的问题更奇怪的是，降低采样率会使问题恶化。在每秒48000个样本的this output中，音高迅速下降，然后再次上升，只是为了再次缓慢下降。在每秒3000个样本的this output中，类似的效果也会发生，但它更加极端和混乱。我做错什么了？在

from math import pi, sin
from sys import byteorder
import wave

def get_sample(time):
    frequency = a10 / 2.0 ** time
    # print('{:.15f} {:.15f} {:.15f}'.format(time, frequency, sin(pi2 * frequency * time)))
    return sin(pi2 * frequency * time)

pi2 = 2 * pi
a10 = 28160.0

NUMBER_OF_CHANNELS = 1
SAMPLE_RATE = 48000  # samples per second
SAMPLE_WIDTH = 3  # bytes
DURATION = 10  # seconds

MAX_SAMPLE_VALUE = 2 ** (SAMPLE_WIDTH * 8 - 1)

samples = bytearray()

for i in range(SAMPLE_RATE * DURATION):
    time = i / SAMPLE_RATE
    sample = round(get_sample(time) * MAX_SAMPLE_VALUE)

    if sample == MAX_SAMPLE_VALUE:
        sample -= 1

    samples.extend(sample.to_bytes(SAMPLE_WIDTH, byteorder, signed=True))

with wave.open('output.wav', 'wb') as output:
    output.setnchannels(NUMBER_OF_CHANNELS)
    output.setsampwidth(SAMPLE_WIDTH)
    output.setframerate(SAMPLE_RATE)
    output.setnframes(NUMBER_OF_CHANNELS * SAMPLE_RATE * DURATION)
    output.setcomptype('NONE', 'not compressed')

    output.writeframes(samples)

Answer 1

28160Hz的频率对于48000的采样率来说太高了。在

当采样频率为3000hz时，最大频率小于1.5KHz

这与奈奎斯特采样率有关。简而言之，在给定的采样率下，您可以采样的最大频率是采样率的1/2。实际上，它不到采样率的1/2。在

请看：

https://en.wikipedia.org/wiki/Nyquist_frequency

0https://dsp.stackexchange.com/

给定48000Hz的采样率，您可以采样的最大频率是24000hz。这个最大频率是理想化的，它会少得多。在

要捕获28160Hz的频率，您需要大于56320Hz的采样率。比如64000Hz，或者更好 96000Hz采样率。在

编辑：顺便问一下，为什么频率函数会上升到时间的幂次？**

这会引起一些奇怪的锯齿效应

我认为应该是：

^{pr2}$

在 我懂了。。。你在做频率扫描。从而调整每个采样时间的频率。在