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]);