1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29
| fid=fopen('045639 T01 02 1028710.txt','r');
data = strrep(fread(fid, [1 inf], '*char'), ',', '.');
fclose(fid);
T = textscan(data,'%f%f%f%f%f%f%f', 'CollectOutput', 1);
t=0:1:length(T{1}(5:end,1))-1;
t=t';
y=[];
Tmax=[];
x=[];
ENERGIE=[];
RMS=[];
i=[];
l=length(x);
j=1;
for i=5:1:7
y(i)=max(max(T{1}(5:end,i))).*0.1;
Tmax(i)=max(max(T{1}(5:end,i)));
x(i)=max(y(i).*ones(length(T{1}(5:end,i)),1),T{1}(5:end,i))-y(i);
ENERGIE(i)=trapz(x(i),[-t*1.195]);
while x(j,i)<=0
j=j+1;
end
while x(l,i)<=0
l=l-1;
end
RMS(i)=sqrt((1000/((l-j)*1.195))*norm(T{1}(j:l,i-4),2));
end
La=20*log10(RMS/10^-6);
Lv=20*log10(RMS/10^-9); |
Partager