[求助][讨论][原创]声波的积分漂移

ALUCARD_1797

本人目前正在开发一种新的音频编码器，在滤波器上遇到了问题，希望懂信号处理方面的大侠帮小弟一把。
4 w6 U7 D) W  j! D8 D, x3 N, q信号的图释
* m1 A% t" P' h" mhttp://www.piccube.de/uploads/d650f8f2d2.gif
% a% X! H6 v/ s
http://www.piccube.de/uploads/cd0002d0e9.gif

信号频率范围（44100，10）HZ
% u: d# Q; B. f% O信号振幅（-32768，+32767）
$ S: i* j2 j, o5 E, [* B: K5 C  Y6 G2 t信号都是N个由未知振幅\频率\相位的正弦信号叠加而成的
' B5 n1 w$ A( l% v8 ?N是有限数量，也许是10个也许是50个也许是2个
特别的是，信号的振幅不是固定不变的，也就是说在t1时刻是x1而在t2时刻是x2
+ {/ y7 b) ^7 H% U: m图中红色的是原始声音信号，白色的是经过算法处理后的信号。原始信号是完全没问题的，声音信号是从CD里面抓取的。
- y; a) u$ u/ l0 O红色原始信号经过了积分后变成白色信号，原本在0的位置上下震荡的信号漂移到了n的位置。现在要把信号矫正过来。绿色的线是0，有没有一种算法可以求出n。
. B2 J9 ?, b, b2 i  q; P相关的算法：
& T* N/ Q/ r7 \0 \9 |! q- z滤波公式http://www.piccube.de/uploads/ef462128ea.gif
4 Q+ _, B6 N/ {1 }! l/ s+ @
double m_u_f=0;
5 f, c; S( f  k$ [* Q1 Z; lfor(x=0;x<size;x++)
( V, c# d' m# r9 X  k8 @{
9 Z' _. [' f# d( n! M* L m_u_f=pBuffer[x]+m_u_f;
0 a# s8 |6 ~* B) ] pBuffer[x]=m_u_f;
$ I6 ^7 e' C9 g. w# ]) x! v}
算法很简单，就是对信号进行累加，但是实际上信号是经过了一个积分的过程。
$ M/ U: M  Q8 ]公式里面的X是时间，N是频率，a是初始时间。m_u_f里面的数值就是在x时刻等式右边的值。信号漂移的原因是N*cos(a/N)（就是画红圈的）这个在滤波过程中的定值所引起的。实际上只要求出这个定值就可以解决漂移的问题。但是N和a(a还包含了相位角的问题)在滤波之前是未知的，所以根本无发通过具体的参数来求出N*cos(a/n)。
: u. E; M6 e3 P1 \1 b7 o目前我可以想到的方法就是，对滤波后的信号（白色的）通过某种算法求出振荡中心位置。
& l# |+ f5 e" Z# i有一点我要说明的，我的采样缓冲区长度不可能很长，只有1/10秒，就是4410个采样点。于是频率的范围是(10~22050)HZ。这是编码器的外在条件。
- q, Y, Q4 C  J9 |4 m  b若想联系本人请加QQ286726462，验证：声波的积分漂移

malei1984mm

<>小弟谢谢了</P>