24小时热门版块排行榜

返回列表

thisisit

铜虫 (小有名气)

应助: 0 (幼儿园)
金币: 280.7
帖子: 124
在线: 11.8小时
虫号: 996869

[交流] 【讨论】线性预测编码中Levinson-Durbin算法在音频信号去相关应用中的问题

很多论文中指出对音频激励信号进行预白处理后得到近似白噪声，将与处理后的音频信号输入NLMS自适应滤波器，得到的残差相比未预白的原始信号所得残差要小，而且NLMS自适应滤波器的收敛速度也会加快，而且常用的方法有使用线性预测去除原始音频信号的自相关性，通过解格型预测器可以得到预测器的系数，常用的方法为Levinson-Durbin算法。

于是我查阅了Levinson-Durbin算法原理及其实现的资料，并将该算法用c语言实现，对于某一帧音频信号，我求得p阶预测器的系数为a0,a1,...,ap，其中a0=1，由此可知，如果该组系数对的话所得的预测语音和原始语音将非常相似，而且其差值e=s(n)+a1s(n-1)+a2s(n-2)+...+aps(n-p)，将非常的小，而且相关性很小，近似为高斯白噪声。而且根据该组系数和所得残差序列将能够获得原始信号，s(n)=e-(a1s(n-1)+a2s(n-2)+...+aps(n-p))。

可是仿真的结果却并非如此，残差序列比较大，而且通过系数和残差序列返回的原始信号与原始信号有较大差异，尤其是音频信号值较大时。

是不是Levinson-Durbin算法不适合音频信号的去相关，可是有的资料却介绍使用该算法，而且在返回时的误差可能是因为截断和类型转换时产生的，可是按理说该值将很小，可是在音频数据较大时，该值却非常大，有时从某一正值编程某一负值，这样的误差太离谱了。

有研究该邻域的虫友用到过Levinson-Durbin算法吗，它有些什么样的特性，而误差有时如何产生的，请虫友帮忙回答，感谢万分。

[ Last edited by mze04532 on 2011-2-22 at 08:25 ]

回复此楼

» 猜你喜欢

售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急已经有4人回复
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急已经有4人回复
b口会评已经有6人回复
信息学部会评时间已经有3人回复
祈祷青基必中已经有9人回复
在职考研已经有7人回复
MSER送审了还被拒稿已经有5人回复
关于如何从代码看上不上会已经有14人回复
硝基还原已经有5人回复
化学口会评？应用基础研究和基础研究，希望能分开已经有3人回复

» 本主题相关价值贴推荐，对您同样有帮助:

云南省应用基础研究面上项目与结果了吗！已经有15人回复
收集的天大、哈工大的MATLAB讲义，还有几个应用实例，分享一下已经有239人回复
求助，两组数据的相关性分析已经有19人回复
帮忙分析下雷达测向中辐射源信号的相关性对测向性能的影响具有研究价值吗？已经有4人回复
DVD音频提取工具：DVD Audio Extractor绿色版【转载】已经有9人回复
非线性最小二乘拟合（初值问题）lsqcurvefit 函数的设定已经有14人回复
关于单取代苯的HMBC相关信号的疑问已经有10人回复
蛋白质信号肽预测的相关软件已经有8人回复
【重金求助】同伦延拓算法是如何具体实现的已经有6人回复
怎么才能做出三个9的线性方程啊？？！！愁死我了已经有11人回复
微波会影响红外信号吗已经有4人回复
南大英语分数线预测贴！已经有33人回复
【求助】2010-2015年中国芳纶纤维行业投资分析及前景预测报告已经有17人回复

» 抢金币啦！回帖就可以得到:

查看全部散金贴

1楼 2011-02-21 19:11:20

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

thisisit

铜虫 (小有名气)

应助: 0 (幼儿园)
金币: 280.7
帖子: 124
在线: 11.8小时
虫号: 996869

误差是由于在LD算法的求自相关系数是产生的，由于其大幅度音频信号下该帧的自相关系数值非常的大，32bit的存储空间远远不够，所以产生了溢出，以至于LD算法产生的去相关滤波器系数不准确，而同时在返回时也产生了同样的问题，将其改为64bit数据类型后，该误差消失，同时却在返回时产生了+-1的误差，该误差是由于float型强制转换为short int型产生的，做以下处理可以解决：

out = (short int)buffer;
if(buffer-(float)out>0.5) out +=1;
else if(buffer-(float)out<-0.5) out -=1;

同时另一解决方案为：
#define WORD2INT(x) ((x) < -32767.5f ? -32768 : ((x) > 32766.5f ? 32767 : floor(.5+(x))))

out = WORD2INT(buffer);

欢迎做音频信号处理的虫友们一起讨论。。