转自:https://blog.csdn.net/wzz4420381/article/details/48676921

原作者:wzz4420381

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 代码示例

  1. void SBC_Decode(uint8_t * DataIn, FILE * fOutput)

  2. {

  3. #define SBC_SAMPLING_FREQ 16

  4. #define SBC_CHANNEL_MODE 0

  5. #define SBC_NUM_OF_SUBBANDS 8

  6. #define SBC_NUM_OF_CHANNELS 1

  7. #define SBC_NUM_OF_BLOCKS 15

  8. #define SBC_ALLOC_METHOD 0

  9. #define SBC_BITPOOL 26

  10. #define SBC_DECODED_BUFFER_SIZE (16*8)

  11. uint8_t blocks_per_packet = SBC_NUM_OF_BLOCKS;

  12. uint8_t num_bits = SBC_BITPOOL;

  13. const uint8_t * buf = (DataIn+1);//ignore CRC byte

  14. uint16_t len = SBC_GROUP_SIZE;

  15. uint16_t usDecodedBuffer[SBC_DECODED_BUFFER_SIZE];

  16. /* convenience */

  17. const uint8_t * end = buf + len;

  18. #define left (end - buf)

  19. uint16_t * outBufPtr = usDecodedBuffer;

  20. /* workspace */

  21. static INSAMPLE samples[16][8]; /* We blow the stack if this is not static. */

  22. ITER i, j, k;

  23. uint32_t scaleFactors[8]; //= {0x0f, 0x0c, 0x0b, 0x0b, 0x0a, 0x0a, 0x09, 0x09};

  24. int32_t bitneed[8];

  25. uint32_t bits[8];

  26. int32_t bitcount, slicecount, bitslice;

  27. uint8_t samplingRate, blocks, snr, numSubbands, bitpoolSz, bitpos = 0x80;

  28. int8_t max_bitneed = 0;

  29. #ifndef SPEED_OVER_ACCURACY

  30. int32_t levels[8];

  31. #endif

  32. #if (DEBUG_DECODING == 1)

  33. const uint8_t *start_buf = buf;

  34. pr_info("%s: blocks_per_packet = %d, num_bits = %d, buf = %p, len = %d\n",

  35. __func__, blocks_per_packet, num_bits, buf, len);

  36. for (i = 0; i < len; i++) {

  37. pr_info("buf[%d] = 0x%02x\n", i, buf[i]);

  38. }

  39. #endif

  40. /* look into the frame header */

  41. if (left < SBC_GROUP_SIZE) goto out;/* too short a frame header */

  42. /* use Bemote specific constants */

  43. samplingRate = 0; /* always 16000 Hz */

  44. blocks = blocks_per_packet;

  45. snr = 0;

  46. numSubbands = SBC_NUM_OF_SUBBANDS;

  47. bitpoolSz = num_bits;

  48. /* read scale factors */

  49. /* pr_info("sbc_decode: read scale factors, numSubbands = %d\n", numSubbands); */

  50. /**/

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

  52. if(bitpos == 0x80){

  53. scaleFactors[i] = (*buf) >> 4;

  54. bitpos = 0x08;

  55. }

  56. else{

  57. scaleFactors[i] = (*buf++) & 0x0F;

  58. bitpos = 0x80;

  59. }

  60. }

  61. /* calculate bitneed table and max_bitneed value (A2DP 12.6.3.1) */

  62. if(snr){

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

  64. bitneed[i] = scaleFactors[i];

  65. if(bitneed[i] > max_bitneed) max_bitneed = bitneed[i];

  66. }

  67. }

  68. else{

  69. const signed char* tbl;

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

  71. else tbl = (const signed char*)loudness_8[samplingRate];

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

  73. if(scaleFactors[i]){

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

  75. if(loudness > 0) loudness /= 2;

  76. bitneed[i] = loudness;

  77. }

  78. else bitneed[i] = -5;

  79. if(bitneed[i] > max_bitneed) max_bitneed = bitneed[i];

  80. }

  81. }

  82. /* fit bitslices into the bitpool */

  83. bitcount = 0;

  84. slicecount = 0;

  85. bitslice = max_bitneed + 1;

  86. /* pr_info("sbc_decode: fit bitslices into the bitpool, bitslice = %d\n", bitslice ); */

  87. do{

  88. bitslice--;

  89. bitcount += slicecount;

  90. slicecount = 0;

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

  92. if(bitneed[i] > bitslice + 1 && bitneed[i] < bitslice + 16) slicecount++;

  93. else if(bitneed[i] == bitslice + 1) slicecount += 2;

  94. }

  95. }while(bitcount + slicecount < bitpoolSz);

  96. /* distribute bits */

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

  98. if(bitneed[i] < bitslice + 2) bits[i] = 0;

  99. else{

  100. int8_t v = bitneed[i] - bitslice;

  101. if(v > 16) v = 16;

  102. bits[i] = v;

  103. }

  104. }

  105. /* allocate remaining bits */

  106. for(i = 0; i < numSubbands && bitcount < bitpoolSz; i++){

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

  108. bits[i]++;

  109. bitcount++;

  110. }

  111. else if(bitneed[i] == bitslice + 1 && bitpoolSz > bitcount + 1){

  112. bits[i] = 2;

  113. bitcount += 2;

  114. }

  115. }

  116. for(i = 0; i < numSubbands && bitcount < bitpoolSz; i++){

  117. if(bits[i] < 16){

  118. bits[i]++;

  119. bitcount++;

  120. }

  121. }

  122. /* reconstruct subband samples (A2DP 12.6.4) */

  123. #ifndef SPEED_OVER_ACCURACY

  124. for(i = 0; i < numSubbands; i++) levels[i] = (1 << bits[i]) - 1;

  125. #endif

  126. /* pr_info("sbc_decode: reconstruct subband samples, blocks = %d\n", blocks ); */

  127. for(j = 0; j < blocks; j++){

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

  129. if(bits[i]){

  130. uint32_t val = 0;

  131. k = bits[i];

  132. do{

  133. val <<= 1;

  134. #if (DEBUG_DECODING == 1)

  135. pr_info("%s: buf = %p, offset %d\n",

  136. __func__, buf, buf-start_buf);

  137. #endif

  138. if(*buf & bitpos) val++;

  139. if(!(bitpos >>= 1)){

  140. bitpos = 0x80;

  141. buf++;

  142. }

  143. }while(--k);

  144. val = (val << 1) | 1;

  145. val <<= scaleFactors[i];

  146. #ifdef SPEED_OVER_ACCURACY

  147. val = mulshift(val, bits[i]);

  148. #else

  149. val /= levels[i];

  150. #endif

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

  152. samples[j][i] = SAMPLE_CVT(val);

  153. }

  154. else samples[j][i] = SAMPLE_CVT(0);

  155. }

  156. }

  157. //sbc_decoder_reset();

  158. for(j = 0; j < blocks; j++){

  159. synth(outBufPtr, samples[j], numSubbands, gV);

  160. outBufPtr += numSubbands;

  161. }

  162. /* if we used a byte partially, skip the rest of it, it is "padding" */

  163. if(bitpos != 0x80) buf++;

  164. out:

  165. #if (DEBUG_DECODING == 1)

  166. if(left < 0)

  167. pr_err("SBC: buffer over-read by %d bytes.\n", -left);

  168. if(left > 0)

  169. pr_err("SBC: buffer under-read by %d bytes.\n", left);

  170. #endif

  171. fwrite(usDecodedBuffer, sizeof(uint16_t), 120, fOutput);

  172. fflush(fOutput);

  173. memset(usDecodedBuffer, 0, sizeof(usDecodedBuffer));

  174. }

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

SBC音频编解码算法(转载)相关推荐

  1. 音频编解码基础知识(一)

    音频编解码常用的三种实现方案 1)采用专用的音频芯片对语音信号进行采集和处理,音频编解码算法集成在硬件内部,如 MP3 编解码芯片.语音合成分析芯片等.使用这种方案的优点就是处理速度块,设计周期短:缺 ...

  2. 蓝牙音频双剑客(二)--高质量音频分布协议(A2DP) SBC编解码算法

    零. 概述 主要介绍下蓝牙协议栈(bluetooth stack)传统蓝牙音频协议之高质量音频分布协议(A2DP) SBC编解码算法 Codec Specific Infomation Element ...

  3. 1.[总结]视音频编解码技术零基础学习方法【转载】【致敬雷神】

    一直想把视音频编解码技术做一个简单的总结,可是苦于时间不充裕,一直没能完成.今天有着很大的空闲,终于可以总结一个有关视音频技术的入门教程,可以方便更多的人学习从零开始学习视音频技术.需要注意的是,本文 ...

  4. 基于 FFMPEG 的音频编解码(二):音频解码

    音频解码 基于 FFMPEG 的音频编解码(一):Hello FFMPEG,安装与编译 基于 FFMPEG 的音频编解码(二):音频解码 基于 FFMPEG 的音频编解码(三):音频编码 在 Hell ...

  5. ffmpeg / avconv是通用的视频/音频编解码命令行工具

    转载自:http://blog.openlg.net/index.php/archives/767 ffmpeg / avconv是通用的视频/音频编解码命令行工具. 通用是既指他们可以处理各种各样的 ...

  6. 即时通讯音视频开发(六):如何开始音频编解码技术的学习

    前言 即时通讯应用中的实时音视频技术,几乎是IM开发中的最后一道高墙.原因在于:实时音视频技术 = 音视频处理技术 + 网络传输技术 的横向技术应用集合体,而公共互联网不是为了实时通信设计的. 系列文 ...

  7. 从编解码算法到全链路RTC架构,揭秘淘系直播技术演进之路

    从2016年直播元年至今,纯粹的直播已经逐渐失去竞争力,越来越多形式创新映入眼帘,而众多企业开始走向内容垂直化--秀场.游戏.电商.广电等内容特点深度结合.伴随2020年疫情爆发,电商为人们日常生活提 ...

  8. 视音频编解码技术及其实现

    核心提示:一.视音频编码国际标准化组织及其压缩标准介绍 国际上有两个负责视音频编码的标准化组织,一个是VCEG(VideocodeExpertGroup),是国际电信联合会下的视频编码专家组,一个是M ...

  9. 语音编码 c语言,语音编解码算法G.723.1在DSP - 嵌入式新闻 - 电子发烧友网

    1 引言 G.723.1是删组织于1996年推出的一种低码率的语音编码算法标准,也是目前该组织颁布的语音压缩标准中码率最低的一种标准.G.723.1主要用于对语音及其它多媒体声音信号的压缩,目前在一些 ...

  10. 基于 FFMPEG 的音频编解码(一):Hello FFMPEG,安装与编译

    Hello FFMPEG 基于 FFMPEG 的音频编解码(一):Hello FFMPEG,安装与编译 基于 FFMPEG 的音频编解码(二):音频解码 基于 FFMPEG 的音频编解码(三):音频编 ...

最新文章

  1. 一篇文章告诉你标准化和归一化的区别?
  2. 批处理start命令学习
  3. python 二维数组 长度_Python创建二维数组的正确姿势
  4. 洛谷 题解 P4955 【[USACO14JAN]Cross Country Skiing 越野滑雪】
  5. android的热点分配ip_Android 得到连接热点的ip的方法
  6. java面试题之hashCode与equals()方法
  7. 【服务器】【个人网盘】宝塔安装NextCloud
  8. python123 测验三_作业要求 20190919-3 效能分析
  9. java程序员中英文简历_Java程序员英文简历范文
  10. word排版教程-----------技巧
  11. 南琼考试系统APP服务器,南琼考试系统
  12. 康佳电视手机遥控器android版,康佳电视遥控器
  13. MongoDB——聚合管道之$limit$skip$sort操作
  14. IOI 2022 简要题解
  15. http://jingyan.baidu.com/article/a3761b2b66fe141577f9aa51.html
  16. 会议OA(会议排座送审)
  17. 2023最新SSM计算机毕业设计选题大全(附源码+LW)之java在线电影院售票系统5u8st
  18. linux 系统安装aria2以及配置web端
  19. 64位系统,定义int* a[2][3],占几个字节?
  20. Java Web --HTML(尚硅谷2022版Javaweb)

热门文章

  1. 保研之路——北邮网研院交换中心夏令营
  2. 爬取雪球网的股票信息评论
  3. java 加权平均_SPSS用加权平均值计算百分位数
  4. qq远程控制无法连接服务器,win10系统QQ远程协助,无法控制对方电脑,怎么办?...
  5. 网站跳出率过高的原因有哪些?怎么解决?
  6. Spring Boot整合Redis---学习09
  7. 为什么使用用友盟的摇一摇截屏功能,截出来的地图图片是黑屏的,渴求大师解答,感激不尽
  8. 出口托收结算方式DA、DP、CAD详解
  9. word中批量添加印章、签名图片
  10. Python 爬虫--下载音乐