xianxinnihe.m

建筑声学

clear all;
x=load('RT_21_4.lvm');
%设定数据提取时间区间进行拟合
BeginTime=1.77;
EndTime=3.7;
%读取文件
L=length(x);
t=x(1:L,1);
y=x(1:L,2);
%对原始数据在时间区间内进行提取，
for k=1:L
    if(t(k)>=BeginTime)
        break;
    end
end
for k1=k:L
    if(t(k1)>=EndTime)
        break;
    end
end

%获得提取后的数据
tt=t(k:k1);
yy=y(k:k1);
%对提取的数据进行线性拟合
p=polyfit(tt,yy,1);
yy1=polyval(p,tt);
%提取数据直接作高次曲线拟合
p=polyfit(tt,yy,15);
yy2=polyval(p,tt);
%对高次曲线在对数坐标里作线性拟合
yy2=log(yy2);
p=polyfit(tt,yy2,1);
yy2=polyval(p,tt);
yy2=exp(yy2);
%对提取后的数据在对数坐标里作线性拟合
%y1=log(y);
%p=polyfit(t,y1,1);
%y1=polyval(p,t);
%y1=exp(y1);
%绘出图形
plot(t,y,'b',tt,yy1,'r',tt,yy2,'y');grid on
%plot(t,y,'b',t,y4,'r',t,y2,'k',t,y3,'k',t,y1,'y');grid on
Ylabel('dB');Xlabel('t /s');
title('混响测量衰减曲线')
%axis([0.6 1.2,0 100]);