音频解码SBC_立体声音频编解码芯片

音频解码SBC_立体声音频编解码芯片SBC音频编解码算法浅析1.SBC算法简介SBC是subbandcode的缩写,也可称为子带编码在A2DP协议中,SBC算法是默认支持的蓝牙SBC算法是一种以中等比特率传递高质量音频数据的低计算复杂度的音频编码算法1.1算法基本框图SBC系统使用一个余弦调制的滤波器组,用来解析和同步。滤波器组可设定成4或8个子带子带信号的量化采用比特分配器和自适应脉冲编码器组调制可用的比特位数

大家好,又见面了,我是你们的朋友全栈君。如果您正在找激活码,请点击查看最新教程,关注关注公众号 “全栈程序员社区” 获取激活教程,可能之前旧版本教程已经失效.最新Idea2022.1教程亲测有效,一键激活。

Jetbrains全系列IDE使用 1年只要46元 售后保障 童叟无欺

SBC音频编解码算法浅析

1. SBC算法简介

  • SBC是sub band code的缩写,也可称为子带编码
  • 在A2DP协议中,SBC算法是默认支持的
  • 蓝牙SBC算法是一种以中等比特率传递高质量音频数据的低计算复杂度的音频编码算法

1.1 算法基本框图

  • SBC系统使用一个余弦调制的滤波器组,用来解析和同步。滤波器组可设定成4或8个子带
  • 子带信号的量化采用比特分配器和自适应脉冲编码器组调制
  • 可用的比特位数和编码器的块数都是可配置的
  • 编码后的数据由比特流打包器打包,以有线或无线方式传输
  • 解码是编码的逆过程
    这里写图片描述

1.1 算法所用参数

  • sampling_frequency:采样频率。常用的是16KHz,32KHz,44.1KHz和48KHz
  • channel mode:通道模式。可以是单声道,双声道,立体声和联合立体声
  • block length:块长度。可以是4,8,12,16
  • subbands:子带数量。可以是4或8
  • allocation method:分配方式。可以是SNR或loudness
  • bitpool:比特池。范围是2-250。此值越大,编码产生的数据越长

2 SBC算法实现

2.1 SBC编码算法实现

这里写图片描述

  • PCM格式的左右声道进入多相解析器,输出尺度因子和子带采样数据
  • 每个尺度因子分别对应一个子带
  • 量化后的子带采样数据需要进行打包,打包方式可以是分段或不分段
    这里写图片描述
  • 多相解析器的代码实现较为复杂,流程图如下[具体请看参考文献2的Appendix B]:
    这里写图片描述

2.2 SBC解码算法实现

这里写图片描述

  • 解码过程是编码过程的逆
  • 多相综合器的代码实现较为复杂,流程图如下[具体请看参考文献2的Appendix B]:
    这里写图片描述

3 SBC解码算法在某蓝牙主设备上的应用

3.1 帧格式

这里写图片描述

3.2 参数选择

  • sampling frequency:16
  • channel mode:单声道
  • number of subbands:8
  • number of channels:1
  • number of blocks:15
  • allocation method:SNR
  • bitpool:26

3.3 代码示例

void SBC_Decode(uint8_t * DataIn, FILE * fOutput)
{
#define SBC_SAMPLING_FREQ 16
#define SBC_CHANNEL_MODE 0
#define SBC_NUM_OF_SUBBANDS 8
#define SBC_NUM_OF_CHANNELS 1
#define SBC_NUM_OF_BLOCKS 15
#define SBC_ALLOC_METHOD 0
#define SBC_BITPOOL 26 
#define SBC_DECODED_BUFFER_SIZE (16*8)

    uint8_t blocks_per_packet = SBC_NUM_OF_BLOCKS;
    uint8_t num_bits = SBC_BITPOOL;
    const uint8_t * buf = (DataIn+1);//ignore CRC byte
    uint16_t len = SBC_GROUP_SIZE;
    uint16_t usDecodedBuffer[SBC_DECODED_BUFFER_SIZE];

    /* convenience */
    const uint8_t * end = buf + len;
#define left (end - buf)
    uint16_t * outBufPtr = usDecodedBuffer;

    /* workspace */
    static INSAMPLE samples[16][8]; /* We blow the stack if this is not static. */             
    ITER i, j, k;
    uint32_t scaleFactors[8]; //= {0x0f, 0x0c, 0x0b, 0x0b, 0x0a, 0x0a, 0x09, 0x09};
    int32_t bitneed[8];
    uint32_t bits[8];
    int32_t bitcount, slicecount, bitslice;
    uint8_t samplingRate, blocks, snr, numSubbands, bitpoolSz, bitpos = 0x80;
    int8_t max_bitneed = 0; 
#ifndef SPEED_OVER_ACCURACY
    int32_t levels[8];
#endif

#if (DEBUG_DECODING == 1)
    const uint8_t *start_buf = buf;
    pr_info("%s: blocks_per_packet = %d, num_bits = %d, buf = %p, len = %d\n",
        __func__, blocks_per_packet, num_bits, buf, len);
    for (i = 0; i < len; i++) {
        pr_info("buf[%d] = 0x%02x\n", i, buf[i]);
    }
#endif

    /* look into the frame header */
    if (left < SBC_GROUP_SIZE) goto out;/* too short a frame header */

    /* use Bemote specific constants */
    samplingRate = 0; /* always 16000 Hz */
    blocks = blocks_per_packet;
    snr = 0;
    numSubbands = SBC_NUM_OF_SUBBANDS;
    bitpoolSz = num_bits;               

    /* read scale factors */
    /* pr_info("sbc_decode: read scale factors, numSubbands = %d\n", numSubbands); */
    /**/
    for(i = 0; i < numSubbands; i++){

        if(bitpos == 0x80){

            scaleFactors[i] = (*buf) >> 4;
            bitpos = 0x08;
        }
        else{

            scaleFactors[i] = (*buf++) & 0x0F;
            bitpos = 0x80;
        }
    }

    /* calculate bitneed table and max_bitneed value (A2DP 12.6.3.1) */
    if(snr){

        for(i = 0; i < numSubbands; i++){

            bitneed[i] = scaleFactors[i];
            if(bitneed[i] > max_bitneed) max_bitneed = bitneed[i];
        }
    }
    else{

        const signed char* tbl;

        if(numSubbands == 4) tbl = (const signed char*)loudness_4[samplingRate];
        else tbl = (const signed char*)loudness_8[samplingRate];

        for(i = 0; i < numSubbands; i++){

            if(scaleFactors[i]){

                int loudness = scaleFactors[i] - tbl[i];

                if(loudness > 0) loudness /= 2;
                bitneed[i] = loudness;
            }
            else bitneed[i] = -5;
            if(bitneed[i] > max_bitneed) max_bitneed = bitneed[i];
        }
    }               

    /* fit bitslices into the bitpool */
    bitcount = 0;
    slicecount = 0;
    bitslice = max_bitneed + 1;
    /* pr_info("sbc_decode: fit bitslices into the bitpool, bitslice = %d\n", bitslice ); */
    do{
        bitslice--;
        bitcount += slicecount;
        slicecount = 0;
        for(i = 0; i < numSubbands; i++){

            if(bitneed[i] > bitslice + 1 && bitneed[i] < bitslice + 16) slicecount++;
            else if(bitneed[i] == bitslice + 1) slicecount += 2;
        }

    }while(bitcount + slicecount < bitpoolSz);              

    /* distribute bits */
    for(i = 0; i < numSubbands; i++){

        if(bitneed[i] < bitslice + 2) bits[i] = 0;
        else{

            int8_t v = bitneed[i] - bitslice;
            if(v > 16) v = 16;
            bits[i] = v;
        }
    }       

    /* allocate remaining bits */
    for(i = 0; i < numSubbands && bitcount < bitpoolSz; i++){

        if(bits[i] >= 2 && bits[i] < 16){

            bits[i]++;
            bitcount++;
        }
        else if(bitneed[i] == bitslice + 1 && bitpoolSz > bitcount + 1){

            bits[i] = 2;
            bitcount += 2;
        }
    }
    for(i = 0; i < numSubbands && bitcount < bitpoolSz; i++){

        if(bits[i] < 16){

            bits[i]++;
            bitcount++;
        }
    }               

    /* reconstruct subband samples (A2DP 12.6.4) */
#ifndef SPEED_OVER_ACCURACY
        for(i = 0; i < numSubbands; i++) levels[i] = (1 << bits[i]) - 1;
#endif

    /* pr_info("sbc_decode: reconstruct subband samples, blocks = %d\n", blocks ); */
    for(j = 0; j < blocks; j++){

        for(i = 0; i < numSubbands; i++){

            if(bits[i]){

                uint32_t val = 0;
                k = bits[i];
                do{

                    val <<= 1;
#if (DEBUG_DECODING == 1)
                    pr_info("%s: buf = %p, offset %d\n",
                        __func__, buf, buf-start_buf);
#endif
                    if(*buf & bitpos) val++;
                    if(!(bitpos >>= 1)){
                        bitpos = 0x80;
                        buf++;
                    }
                }while(--k);

                val = (val << 1) | 1;
                val <<= scaleFactors[i];

                #ifdef SPEED_OVER_ACCURACY
                    val = mulshift(val, bits[i]);
                #else
                    val /= levels[i];
                #endif

                val -= (1 << scaleFactors[i]);

                samples[j][i] = SAMPLE_CVT(val);
            }
            else samples[j][i] = SAMPLE_CVT(0);
        }
    }       

    //sbc_decoder_reset();

    for(j = 0; j < blocks; j++){
        synth(outBufPtr, samples[j], numSubbands, gV);
        outBufPtr += numSubbands;
    }

    /* if we used a byte partially, skip the rest of it, it is "padding" */
    if(bitpos != 0x80) buf++;

    out:
#if (DEBUG_DECODING == 1)
        if(left < 0)
            pr_err("SBC: buffer over-read by %d bytes.\n", -left);
        if(left > 0)
            pr_err("SBC: buffer under-read by %d bytes.\n", left);
#endif 

    fwrite(usDecodedBuffer, sizeof(uint16_t), 120, fOutput);
    fflush(fOutput);
    memset(usDecodedBuffer, 0, sizeof(usDecodedBuffer));

}

4 总结

  • 在章节3中给出的例子中,压缩前的数据有120*2=240Byte,压缩后的数据有54Byte,压缩比接近4.4:1,压缩比是可调的,其与编码参数有关
  • 压缩后的帧数据长度可以由编码参数计算出来。此例中,
    帧长度=4+(4 * 子带数量 * 通道数量)/8+(块数量 * 通道数量 * bitpool)/8
    =4+(4*8*1)/8+(15*1*26)/8=8+48.75=57(加上CRC正好58字节)
  • 如果在传输过程中丢失了个别帧,解压后回放时会出现一小段音频的整体丢失,对整体的解压没有影响

参考文献

1. LOW-POWER IMPLEMENTATION OF THE BLUETOOTH SUBBAND AUDIO CODEC.pdf
2. A2DP_SPEC_V12.pdf
3. 音频数据SBC算法编码方法及蓝牙立体声子系统.pdf

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 举报,一经查实,本站将立刻删除。

发布者:全栈程序员-用户IM,转载请注明出处:https://javaforall.cn/194905.html原文链接:https://javaforall.cn

【正版授权,激活自己账号】: Jetbrains全家桶Ide使用,1年售后保障,每天仅需1毛

【官方授权 正版激活】: 官方授权 正版激活 支持Jetbrains家族下所有IDE 使用个人JB账号...

(0)


相关推荐

  • 图论简介[通俗易懂]

    图论简介[通俗易懂]这里介绍图论(GraphTheory),图论是计算机科学中非常重要的一部分内容,甚至可以单独划分成为一个领域。很多人第一次接触到图论这个词,就觉得图论是研究和图画相关的内容。不过当大家真的去学习图

  • vuetify富文本编辑器_vue富文本编辑器的使用

    vuetify富文本编辑器_vue富文本编辑器的使用最近再弄一个后台管理系统,挑选了不少的编辑器,最终选择了tinymce,UI精美,功能模块多,可按需加载配置来一张tinymce官网的完整功能的图(没梯子可能访问速度有点慢…)下面开始工作:插件安装tinymce官方提供了一个vue的组件tinymce-vue如果有注册或购买过服务的话,直接通过组件配置api-key直接使用,像我这种懒的注册或者购买的直接下载tinymce,自力更生…

  • IntelliJ IDEA常用设置和好用插件,不定时更新 2021-08-12更新

    IntelliJ IDEA常用设置和好用插件,不定时更新 2021-08-12更新一:设置方面1:改黑色背景主题2:改字体设置到手第一件事就是改字体,改idea显示字体,和代码显示字体首先改软件显示字体及大小,最喜欢使用的是CourierNew字体然后改代码编辑区字体及大小和字体间距3注释模板设置首先是新创建文件时的文档注释/***[description]*@Author:某某某*@Date:${DATE}*@versionv1….

  • pyhcarm激活码(最新序列号破解)[通俗易懂]

    pyhcarm激活码(最新序列号破解),https://javaforall.cn/100143.html。详细ieda激活码不妨到全栈程序员必看教程网一起来了解一下吧!

  • java连接ldap服务器_从Java应用程序连接LDAP服务器「建议收藏」

    java连接ldap服务器_从Java应用程序连接LDAP服务器「建议收藏」我正在构建一个基于GXT(J2EE)的应用程序.现在的问题是我必须将应用程序连接到LDAP服务器.您能告诉我如何从我们的Java应用程序连接LDAP服务器以及我将使用哪些库或API?解决方法:要连接到LDAP,请查看以下包/类:javax.naming.directory.*javax.naming.ladp.*com.sun.jndi.ldap.LdapCtxFactorycom.sun.jnd…

  • Android 基于TCP的 Socket 编程实现(结合 okio)

    Android 基于TCP的 Socket 编程实现(结合 okio)两个进程如果要进行通讯最基本的一个前提就是能够唯一的标识一个进程,在本地进程通讯中我们可以使用PID来唯一标识一个进程,但PID只在本地是唯一的,网络中两个进程PID冲突几率很大,这时我们就需要通过其他手段来唯一标识网络中的进程了,我们知道IP层的ip地址可以唯一标示主机,而TCP层协议和端口号结合就可以唯一标示主机的一个进程了。

发表回复

您的电子邮箱地址不会被公开。

关注全栈程序员社区公众号