C#获取声音信号并通过FFT得到声音频谱

最近研究了下用C#通过麦克风获取声音信号,再通过快速傅里叶变换,得到声音的频谱。

想应用到实际的led灯上,放歌曲,led灯展示频谱,现在用C#先模拟实现下,如下图。

上面的表展示的是声音的电压信号,下图是频谱信息

一、获取声音信号

通过麦克风获取声音用的是NAudio这个库,可以通过NuGet获取

具体使用方式是:

cs 复制代码
        private NAudio.Wave.WaveIn sourceStream = null;
        private NAudio.Wave.DirectSoundOut waveOut = null;

获取麦克风列表

cs 复制代码
List<NAudio.Wave.WaveInCapabilities> sources = new List<NAudio.Wave.WaveInCapabilities>();
for (int i = 0; i < NAudio.Wave.WaveIn.DeviceCount; i++)
{
    sources.Add(NAudio.Wave.WaveIn.GetCapabilities(i));
}
sourceListView.Items.Clear();
foreach (var source in sources)
{
    ListViewItem item = new ListViewItem(source.ProductName);
    item.SubItems.Add(new ListViewItem.ListViewSubItem(item, source.Channels.ToString()));
    sourceListView.Items.Add(item);
}

初始化麦克风

cs 复制代码
 if (sourceListView.SelectedItems.Count == 0) return;
 int deviceNumber = sourceListView.SelectedItems[0].Index;
 sourceStream = new NAudio.Wave.WaveIn();
 sourceStream.DeviceNumber = deviceNumber;

 int a=WaveIn.GetCapabilities(0).Channels;
 sourceStream.WaveFormat = new NAudio.Wave.WaveFormat(9600, 16, 1);
 sourceStream.DataAvailable += new EventHandler<NAudio.Wave.WaveInEventArgs>(sourceStream_DataAvailable);
 sourceStream.BufferMilliseconds = 20;
 sourceStream.StartRecording();

回调函数,lastBuffer就存储的是声音信号

cs 复制代码
private void sourceStream_DataAvailable(object sender, WaveInEventArgs e)
{    
    int samplesRecord = e.BytesRecorded / 2;
    if (lastBuffer is null||lastBuffer.Length!=samplesRecord)
        lastBuffer = new double[samplesRecord];
    for (int i = 0; i < samplesRecord; i++)
        lastBuffer[i] = BitConverter.ToInt16(e.Buffer, i * 2);
        
}

二、FFT变换

FFT变换原理之前文章我用C++写过

C++实现二维快速傅里叶变换(FFT)_c++ fft-CSDN博客

这里我用C#再实现下,复数可以直接using System.Numerics里面的Complex,而不用自己定义了。

FFT源码如下:

cs 复制代码
int ReverseBin(int a, int n)
{
    int ret = 0;
    for (int i = 0; i < n; i++)
    {
        if ((a & (1 << i)) != 0) ret |= (1 << (n - 1 - i));
    }
    return ret;
}

public void FFT(ref Complex[] a, int log2_N, int opt)
{
    int index;
    int length = 1 << log2_N;
    Complex[] tempA = new Complex[length];
    for (int i = 0; i < length; i++)
    {
        index = ReverseBin(i, log2_N);
        tempA[i] = a[index];
    }

    //生成WN表,减少重复计算
    Complex[] WN = new Complex[length / 2];
    for (int i = 0; i < length / 2; i++)
    {
        WN[i] = new Complex(Math.Cos(2 * Math.PI * i  / length), opt * -1 * Math.Sin(2 * Math.PI*i / length));
    }

    //蝶形运算
    int Index0, Index1;
    Complex temp;
    for (int steplenght = 2; steplenght <= length; steplenght *= 2)
    {
        for (int step = 0; step < length / steplenght; step++)
        {
            for (int i = 0; i < steplenght / 2; i++)
            {
                Index0 = steplenght * step + i;
                Index1 = steplenght * step + i + steplenght / 2;

                temp = tempA[Index1] * WN[length / steplenght * i];
                tempA[Index1] = tempA[Index0] - temp;
                tempA[Index0] = tempA[Index0] + temp;
            }
        }
    }

    for (int i = 0; i < length; i++)
    {
        if (opt == -1)
        {
            a[i] = tempA[i] / length;
        }
        else
        {
            a[i] = tempA[i];
        }
    }

}

三、总结

这样可以实时得到了一段声音信号,可以加个汉明窗处理下,再FFT下,就得到声音的频谱信息了。后续计划,可以通过这种方式,控制led灯,展示音乐的频谱跳动。

相关推荐
向宇it5 小时前
【unity小技巧】unity 什么是反射?反射的作用?反射的使用场景?反射的缺点?常用的反射操作?反射常见示例
开发语言·游戏·unity·c#·游戏引擎
九鼎科技-Leo5 小时前
什么是 WPF 中的依赖属性?有什么作用?
windows·c#·.net·wpf
Heaphaestus,RC6 小时前
【Unity3D】获取 GameObject 的完整层级结构
unity·c#
baivfhpwxf20236 小时前
C# 5000 转16进制 字节(激光器串口通讯生成指定格式命令)
开发语言·c#
直裾7 小时前
Scala全文单词统计
开发语言·c#·scala
ZwaterZ8 小时前
vue el-table表格点击某行触发事件&&操作栏点击和row-click冲突问题
前端·vue.js·elementui·c#·vue
ZwaterZ10 小时前
el-table-column自动生成序号&&在序号前插入图标
前端·javascript·c#·vue
SRC_BLUE_1713 小时前
SQLI LABS | Less-55 GET-Challenge-Union-14 Queries Allowed-Variation 2
oracle·c#·less
yngsqq13 小时前
037集——JoinEntities连接多段线polyline和圆弧arc(CAD—C#二次开发入门)
开发语言·c#·swift
Zԅ(¯ㅂ¯ԅ)14 小时前
C#桌面应用制作计算器进阶版01
开发语言·c#