fanxiangjifeng.m

脉冲反向积分发测试混响程序

a=load('RT_36_1.lvm');%加载文件
n=length(a);
y=a(1:n,2);
t=a(1:n,1);

%满场取值
yy=y;
tt=t;

y=y.^2;
for i=n-1:-1:1%反向积分
    y(i)=y(i)+y(i+1);
end
%plot(t,y);
y=10*log10(y);%求出声强级
y=rot90(y);
y=flipud(y);

%plot(t,y,'b--','LineWidth',1.5);%绘出声压级脉冲响应衰减曲线
grid on;
%Ylabel('dB');Xlabel('t /s');




fix=-1.2; %设置修正系数
yy1=yy.*exp(fix*tt);
yy2=yy1.^2;
for i=n-1:-1:1%反向积分
    yy2(i)=yy2(i)+yy2(i+1);
end
%plot(tt,yy2);
yy2=10*log10(yy2);%求出声强级
yy2=rot90(yy2);
yy2=flipud(yy2);



%Ylabel('dB');Xlabel('t /s');

%修正上移
d=y(1)-yy2(1);
yy2=yy2+d;
plot(t,yy2,'r','LineWidth',1.5);grid on;%绘出声压级脉冲响应衰减曲线
hold on;plot(t,y,'b--','LineWidth',1.5);%绘出声压级脉冲响应衰减曲线
h = legend('推算满场','实测空场');
Ylabel('dB');Xlabel('t /s');
title('脉冲反向积分法声压级衰减曲线')