简单的音频重采样

转载 2013年12月05日 19:19:05
Performance of audio resampling software:一些重采样的方法的测试比较

void _x_audio_out_resample_mono(short* input_samples, int in_samples,
    short* output_samples, int out_samples)
{
  int osample;
  /* 16+16 fixed point math */
  uint32_t isample = 0;
  uint32_t istep = ((in_samples-2) << 16)/(out_samples-2);

#ifdef VERBOSE
  printf ("Audio : resample %d samples to %d\n",
          in_samples, out_samples);
#endif

  for (osample = 0; osample < out_samples - 1; osample++) {
    int s1;
    int s2;
    int16_t os;
    uint32_t t = isample&0xffff;
   
    /* don't "optimize" the (isample >> 16)*2 to (isample >> 15) */
    s1 = input_samples[(isample >> 16)];
    s2 = input_samples[(isample >> 16)+1];
   
    os = (s1 * (0x10000-t)+ s2 * t) >> 16;
    output_samples[osample] = os;

    isample += istep;
  }
  output_samples[out_samples-1] = input_samples[in_samples-1];
}



使用方法:
如22050Hz-->8000Hz,参数为
_x_audio_out_resample_mono(inbuf, in_sample_num, outbuf, in_sample_num*8000/(22050*2));
inbuf为待采样音频,outbuf为输出缓存,保证尺寸不小于
in_sample_num*8000/(22050*2)


另一个函数,比上面的效果要好。

typedef short HWORD;
typedef int WORD;
typedef unsigned int UWORD;
typedef unsigned short UHWORD;
#define MAX_HWORD (32767)
#define MIN_HWORD (-32767)
#define Np 15
#define Pmask ((1<<Np)-1)

static inline HWORD WordToHword(WORD v, int scl)
{
    HWORD out;
    WORD llsb = (1<<(scl-1));
    v += llsb; /* round */
    v >>= scl;
    if (v>MAX_HWORD) {
#ifdef DEBUG
        if (pof == 0)
          fprintf(stderr, "*** resample: sound sample overflow\n");
        else if ((pof % 10000) == 0)
          fprintf(stderr, "*** resample: another ten thousand overflows\n");
        pof++;
#endif
        v = MAX_HWORD;
    } else if (< MIN_HWORD) {
#ifdef DEBUG
        if (nof == 0)
          fprintf(stderr, "*** resample: sound sample (-) overflow\n");
        else if ((nof % 1000) == 0)
          fprintf(stderr, "*** resample: another thousand (-) overflows\n");
        nof++;
#endif
        v = MIN_HWORD;
    } 
    out = (HWORD) v;
    return out;
}


static int 
  SrcLinear(HWORD X[], HWORD Y[], double factor, UWORD *Time, UHWORD Nx)
{
    HWORD iconst;
    HWORD *Xp, *Ystart;
    WORD v,x1,x2;
    
    double dt; /* Step through input signal */ 
    UWORD dtb; /* Fixed-point version of Dt */
    UWORD endTime; /* When Time reaches EndTime, return to user */
    
    dt = 1.0/factor; /* Output sampling period */
    dtb = dt*(1<<Np) + 0.5; /* Fixed-point representation */
    
    Ystart = Y;
    endTime = *Time + (1<<Np)*(WORD)Nx;
    while (*Time < endTime)
    {
        iconst = (*Time) & Pmask;
        Xp = &X[(*Time)>>Np]; /* Ptr to current input sample */
        x1 = *Xp++;
        x2 = *Xp;
        x1 *= ((1<<Np)-iconst);
        x2 *= iconst;
        v = x1 + x2;
        *Y++ = WordToHword(v,Np); /* Deposit output */
        *Time += dtb; /* Move to next sample by time increment */
    }
    return (- Ystart); /* Return number of output samples */
}


相关文章推荐

算法系列之二十四:离散傅立叶变换之音频播放与均衡器

导语 在算法系列的第二十二篇,我们介绍了离散傅立叶变换算法的实现,将时域的音频信号转换到频域进行分析,获取拨号音频的频率特征。这一篇我们将介绍一种频域均衡器的实现方法,所谓的频域均衡器,就是在频域...
  • orbit
  • orbit
  • 2015-05-04 20:06
  • 11598

音频重采样原理及技术实现

转载请注明出处! 音频重采样作为一个独立模块蕴含了数字信号处理理论的多方面内容,综合起来其物理原理及滤波器的实现优化可以作为一个独立的项目做较深入的研究,可谓是一门学问。推荐大家研读国外作者的理论...

音频重采样

很多时候在播放音频时,音频的采样率、通道数等并不能满足我们播放得需求,这就需要我们对音频进行重采样。 FFmpeg提供的方法libswresample.a库为我们很好的解决重采样的问题,下面就对如何使...

数字音频重采样

1. 什么是有限带宽插值(Bandlimited Interpolation)        在数字信号处理中,作用于离散时间信号的有限带宽插值是一个得到广泛应用的基本工具。一般而言,我们所面临的问题...

语音去噪LMS算法释义

如需转载请注明出处! LMS算法是音频去噪方法中的经典之作,在此基础上衍生出改进算法如NLMS、FxLMS算法等。 图1. LMS算法框架 在做算法框架搭建及传递函数估计仿真的时...

音频重采样的实现

在声纹识别中,为了满足对不同采样率的要求,常需要对语音进行重采样。重采样即将原始的采样频率变换为新的采样频率以适应不同采样率的要求。实现重采样的传统方法有三种:一是若原模拟信号x(t)可以再生,或是已...

算法系列之九:计算几何与图形学有关的几种常用算法(一)

我的专业是计算机辅助设计(CAD),算是一半机械一半软件,《计算机图形学》是必修课,也是我最喜欢的课程。热衷于用代码摆平一切的我几乎将这本教科书上的每种算法都实现了一遍,这种重复劳动虽然意义不大,但是...
  • orbit
  • orbit
  • 2011-12-18 23:13
  • 40979

《算法的乐趣》终于出版了

经过一年多的编写和准备,《算法的乐趣》终于出版了,4月中旬开始可以在京东、当当等网店购买,实体书店上架可能会慢一点。感谢图灵社区的各位老师的辛勤耕耘,也感谢《啊哈!算法》作者纪磊先生(啊哈磊),Lin...
  • orbit
  • orbit
  • 2015-04-22 23:23
  • 13430
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:深度学习:神经网络中的前向传播和反向传播算法推导
举报原因:
原因补充:

(最多只允许输入30个字)