归一化.m

用MATLAB实现了该算法

[F,fs,bits]=wavread('h_0b');
F=F(:,1);
F1=wavread('yy2nxx_pn_65s_1e-1_0db_0213.wav');
F1=F1(:,1);  %第一声道
N=length(F);
N1=length(F1);
if N>N1          %求min(N,N1)
    min=N1;
else
    min=N;
end
sum=0;                    %根据归一化公式编写代码
sum1=0;
for i=1:min
    sum=sum+abs( F(i) );
    sum1=sum1+abs( F1(i) );
end
fenmu=sum1/sum;
F2=F1/fenmu;%%%

F=F';

F_L=length(F);
t=(0:F_L-1)/fs;

F_L_1=length(F1);
t1=(0:F_L_1-1)/fs;

F_L_2=length(F2);
t2=(0:F_L_2-1)/fs;


subplot(3,1,1);            %画时域图
plot(t,F);xlabel('Time(s)');title('原始音频时域图');
subplot(3,1,2);plot(t1,F1);xlabel('Time(s)');title('录制音频时域图');
subplot(3,1,3);plot(t2,F2);xlabel('Time(s)');title('录制音频归一化处理后时域图');