指数正弦扫频

最新推荐文章于 2022-08-30 17:32:33 发布

dbdxnuliba

最新推荐文章于 2022-08-30 17:32:33 发布

阅读量933

点赞数

分类专栏：数字信号处理文章标签：矩阵几何学线性代数

原文链接：https://blog.csdn.net/fzh2712/article/details/85289139?utm_source=app&app_version=4.21.0&code=app_1562916241&uLinkId=usr1mkqgl919blen

版权

数字信号处理专栏收录该内容

4 篇文章 1 订阅

订阅专栏

本文介绍了如何使用Python生成正弦扫频信号，包括sinesweep函数的数学定义和实现代码。通过示例展示了从生成信号到绘制波形以及保存为WAV文件的完整过程。同时提到了synchronized swept-sine方法及其相关论文，该方法在非线性系统识别中具有应用价值。

摘要由CSDN通过智能技术生成

本文主要参考：Sine Sweep(正弦扫频信号)_Recoding...-CSDN博客_正弦扫频信号

Sine Sweep(正弦扫频信号)_fzh2712的专栏-CSDN博客_扫频信号

读后记录如下：

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()

3.运行结果

syncsweptsine方法：

github:

https://github.com/SiggiGue/syncsweptsine

安装方法：pip install git+https://github.com/SiggiGue/syncsweptsine.git

Novak, A., Lotton, P., & Simon, L. (2015). Synchronized Swept-Sine: Theory, Application, and Implementation. Journal of the Audio Engineering Society, 63(10), 786–798.doi:10.17743/jaes.2015.0071

Novák, A., Simon, L., Kadlec, F., & Lotton, P. (2010). Nonlinear System Identification Using Exponential Swept-Sine Signal. IEEE Transactions on Instrumentation and Measurement, 59(8), 2220–2229.doi:10.1109/tim.2009.2031836