Sine Sweep(正弦扫频信号)

最新推荐文章于 2024-10-31 16:04:14 发布

fzh2712

最新推荐文章于 2024-10-31 16:04:14 发布

阅读量1.9w

点赞数 7

分类专栏：音频处理

音频处理专栏收录该内容

7 篇文章

订阅专栏

本文介绍Sinesweep信号的基本函数表达式及其在Python中的实现方法，包括使用numpy和matplotlib库绘制信号波形，以及如何将生成的信号保存为WAV文件。通过具体代码示例展示了从信号生成到音频文件输出的全过程。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

本文主要参考：https://blog.youkuaiyun.com/okobul1/article/details/46415287

读后记录如下：

1. Sine sweep波的函数

sweep(t) = sin(2*PI*f0*(k^t-1)/log(k))

2. python的代码

from math import * 
import array
import matplotlib.pyplot as plt
import wave

def sinesweep(f0, f1, sweeptime, samplingrate, peak): 
    k = exp(log(float(f1)/f0)/sweeptime)
    data_len = sweeptime * samplingrate
    data = array.array('i',[0]*data_len)
    dt = 1.0/samplingrate
    t = 0.0 
    p = 2*pi*f0/log(k)
    for i in range(data_len):
        data[data_len-i-1] = int(peak*sin( p*(pow(k,t)-1) ))
        t += dt 
    return data

def sindata(f0, time, samplingrate, peak): 
    data = [int(peak*sin(2*pi*f0*x/samplingrate)) for x in xrange(int(time*samplingrate))] 
    return array.array('h', data)

if 1:
    SAMPLING_RATE = 44100 
    SWEEP_TIME = 2
    F0 = 20 
    F1 = 20000 
    PEAK = 0x2000 
    data = sinesweep(F0, F1, SWEEP_TIME, SAMPLING_RATE, PEAK) 
    plt.title("sweep data")
    plt.plot(data)
    plt.show()
if 1 :
    f = wave.open("sweep.wav","wb")
    f.setnchannels(1) # mono wave 
    f.setsampwidth(2) # 16bit 
    f.setframerate(SAMPLING_RATE) #sampling rate 
    f.writeframes(data)
    f.close()