自适应增益AGC算法

最近有同事研究语音通话功能，发现嵌入式设备上录制音量过小。他跟我交流的时候说想要达到自适应调节音量的作用，然后我觉得这个点蛮有趣的，就百度搜索了大量PCM或wav的音量大小的资料，自己在matlab上做仿真验证自适应增益AGC算法。

测试语音文件一：

AGC算法处理后的测试语音文件一：

测试语音文件二：

AGC算法处理后的测试语音文件二：

然后我又分别计算响度了

format long;

[x,fs,nbits]=wavread('old.wav');
[N col]=size(x);
Px=zeros(size(x));
Py=zeros(size(x));
y=zeros(size(x));
g=zeros(size(x));
a=0.2;
mu=0.5;
Pref=0.4;
g(1,1)=1;
y(1,1)=x(1,1);
Py(1,1)=(1-a)*x(1,1)*x(1,1);
Px(1,1)=(1-a)*x(1,1)*x(1,1);
for n=2:N
Px(n,1)=a*Px(n-1,1)+(1-a)*x(n,1)*x(n,1);
g(n,1)=g(n-1,1)*(1+mu*Px(n,1)*(Pref-Py(n-1,1)));
y(n,1)=g(n,1)*x(n,1);
if y(n,1)>2
y(n,1)=2;
end
if y(n,1)<-2
y(n,1)=-2;
end
Py(n,1)=g(n,1)*g(n,1)*Px(n,1);
end

subplot(3,1,1);
plot(x);xlabel('number of sample');
title('Input signals'); 

wavwrite(y,fs,nbits,'result_old.wav');
subplot(3,1,2);
plot(y);xlabel('number of sample');

title('Output signals'); 
subplot(3,1,3);
plot(g);xlabel('number of sample');title('gain');legend('gain');
figure;plot(y,'r'); hold on;plot(x); 
title('Output signals & Input signals ');xlabel('number of sample');
legend('Output signals','Input signals');