- 博客(11)
- 资源 (7)
- 收藏
- 关注
RSS订阅原创 matlab自带的VAD检测程序
在matlab的DSP工具箱中，自带了VAD程序，本文主要记录如何使用这个程序。clear allclose allclc%% 设置音频源，读入一个speaker.wav文件，并对读入的文件进行分帧，每帧80个采样点audioSource = dsp.AudioFileReader('SamplesPerFrame',80,... ...
2018-07-31 14:50:39
6753
6
原创 正交性原理与维纳霍夫（正则）方程
有期望信号d（n），纯净信号x（n），以及噪声信号g（n）;有滤波器h（m）,以及滤波器输出信号y（n），滤波器输出纯净信号x的估计值y，并且使y与x误差的期望为0. （相关性）根据最小均方差准则设计滤波器,滤波器设计问题转为函数求驻点的问题，在驻点处能够取得极值，如果能证明该式有连续二阶偏导，且得到的hesse矩阵正定...
2018-07-31 10:33:49
12886
1
原创 快速查找某个范围内的所有素数
观察以下素数表，不能被小于素数n的素数整除的最小数为下一个最小素数。换句话说所谓素数可以看做是坐标中不断使用更小素数的倍数填充后剩余的最小数。例：1】2是最小素数，使用2的倍数不断填充坐标后，未被填充的最小数就是3，于是3就是大于2的下一个素数2】然后使用3的倍数不断填充坐标，剩余的未被填充的最小数就是5（4已经被2的倍数填充），可以发现5是大于3的下一个素数3】依次类推可以顺序的...
2018-07-30 14:35:33
7573
1
原创 卷积与圆周卷积
关于卷积的两种手算方法见链接matlab中fft函数是对序列做DFT，第二个参数是采样点的数量，当缺省时，默认为输入信号序列的长度；当大于信号序列的长度时候，对输入序列后面补零计算，频域信号的分辨率因此而变细，但是这时候的分辨率是伪分辨率。设定信号序列s1，length（s1） = 9s1 = [1,7,8,9,5,4,6,3,2]调用fft函数，对s1进行傅里叶变换plo...
2018-07-25 10:57:30
8751
原创 枪手博弈
博弈论中的枪手模型说，三个枪手中最弱的生存下来的可能性最大，这其实挺反直觉的。模型中设置三者同时开枪、且开枪只有一轮，我觉得这是不合理的。我觉得现实中应该是，一个裁判喊开始，但是枪手的出枪速度有块有慢，而且决斗必须决定出谁是最后那一个站着的人。我觉得最弱的那个人，一定是生存几率最低的那个人，能力越强生存下去的几率就越大。 有了这个想法，我们就需要来验证一下。找三个人来站着对着开枪...
2018-07-24 18:07:17
3298
1
原创 webRTC AEC 添加舒适噪声算法
webRTC的舒适噪声实际上就是，一个随机噪声经过同一帧消回声输出NLP参数的加权处理产生的噪声。如果是32K信号，则对产生的高频段随机噪声统一加权，加权值是同一帧消回声输出NLP参数的平均值。static void ComfortNoise(AecCore* aec, float efw[2][PART_LEN1], ...
2018-07-17 15:48:12
2922
1
原创 webRTC AEC 非线性滤波处理算法
先把代码贴上，有空时候回来注释static void NonLinearProcessing(AecCore* aec, short* output, short* outputH) { float efw[2][PART_LEN1], dfw[2][PART_LEN1], xfw[2][PART_LEN1]; complex_t comfortNoiseHband[PART_LEN...
2018-07-17 15:06:27
4607
4
原创 webRTC AEC 发散处理与抑制算法
所谓发散问题就是在自适应滤波的过程中出现了权值不收敛的问题，一般是由于参考信号与近场数据时间点没有对齐引起的。当滤波器的权值不收敛的时候，自适应滤波器就不能正常的跟踪参考信号产生正确的的合成回声，以至于无法完成回声消除的功能，甚至自适应滤波器会产生一些乱七八糟的信号，淹没近端信号或产生音爆。综上，回声消除过程中要对回声消除的残差进行监控，如果残差信号的能量超过了近场信号的数据，就说明出现了发散...
2018-07-17 14:08:25
2183
4
原创 webRTC AEC 远场滤波器
最后几句代码将远场经过滤波之后加到了传入的yf（合成回声信号）中，但是在调用本滤波器之前已经将yf全部设置为0了，因此输出的合成回声信号就是滤波器的输出。static void FilterFarSSE2(AecCore* aec, float yf[2][PART_LEN1]) { int i; // 设置分块数 const int num_partitions = aec-&...
2018-07-17 14:04:37
850
原创 webRTC AEC 归一化误差信号
因为误差信号的大小变化差异较大，不便于设置门限阈值，因此将误差信号e对远端信号x进行归一化，便于对与阈值比较。static void ScaleErrorSignalSSE2(AecCore* aec, float ef[2][PART_LEN1]) { // extern __m128 _mm_set_ps1(float _W); //返回一个__m128的寄存器，Sets the f...
2018-07-17 14:01:38
1044
原创 webRTC AEC 滤波器的自适应
权值更新就是使用SSE指令集加速完成NLMS算法的频域计算，具体实现细节见代码注释。static void FilterAdaptationSSE2(AecCore* aec, float* fft, float ef[2][PART_LEN1]) { in...
2018-07-17 13:56:44
1491
脉冲时不变法设计数字滤波器的matlab实现.m
2019-11-23
Optimal Time-Domain Noise Reduction Filters--A Theoretical Study.pdf
2018-10-17
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人