clear
clc
%%% 插入损失的计算(离散频率100Hz到1500Hz)
fs=10000;
sensitivity=[42.44 40.11 35.67 35.94]*10^(-3); %% mv/pa
%%% 3米处声压插入损失
gain=[10 10 10 50 50 50 10 10 50 50 10 10 10 10 10]; %% 3米处声压测点的放大倍数
p0=2.0e-5;
k=1;
for i=100:100:1500
filename=strcat('TL_OUT_',num2str(i),'_time.lvm');
data=textread(filename,'%f','headerlines',23);
AA=data(3:8:length(data))/gain(k)/sensitivity(3);
t_pressure(:,k)=AA;
k=k+1;
end
f_pressure=abs(fft(t_pressure,fs)*2/fs);
f_pressure_dB=20*log10(f_pressure/p0);
for i=1:15
SPL_OUT(i,1)=max(f_pressure_dB(i*100-5:i*100+5,i));
end
clear t_pressure
k=1;
fs=10000;
for i=100:100:1500
filename=strcat('TL_IN_',num2str(i),'_time.lvm');
data=textread(filename,'%f','headerlines',23);
AA=data(3:8:length(data))/50/sensitivity(3);
t_pressure(:,k)=AA;
k=k+1;
end
f_pressure=abs(fft(t_pressure,fs)*2/fs);
f_pressure_dB=20*log10(f_pressure/p0);
for i=1:15
SPL_IN(i,1)=max(f_pressure_dB(i*100-5:i*100+5,i));
end
IL = SPL_OUT-SPL_IN;
figure(1)
stem(100:100:1500,IL);
figure(2)
plot(f_pressure_dB(2:1000,2))
% f_p_2=f_pressure.^2;
%
% f_pressure_mean=sqrt(mean(f_p_2,2));
%
% f_pressure_mean2=mean(f_pressure,2);
%
% SPL_Z=20*log10(f_pressure_mean(1:3000,1)/p0);
%
% N=3000;
% f=1:N; %frequency
% for i=1 : N
% A(i) = 10*log10(3.5*(10^16)*(f(i)^8)/(((20.6^2 + f(i)^2)^2)*(107.7^2 + f(i)^2)*(737^2+f(i)^2)*((12194.2^2 + f(i)^2)^2))) ; %A weighting
% end
%
% SPLA_data=SPL_Z+A';
%
% % SPLA_total=10*log10(sum(10.^(SPLA_data(1:N/2)./10)))
%
% SPLA_total=10*log10(sum(10.^(SPLA_data(20:5:N/2)./10)))
%
%
% %%% 计算辐射声功率级
%
% SPowerLA_=SPLA_total+10*log10(s/s0)
%
%
% figure(3)
% plot(21:1500,SPLA_data(21:1500),'r','linewidth',2)
% xlabel('Frequency (Hz)','fontsize',14)
% ylabel('Sound Power Level (dBA)','fontsize',14)
% title('Sound Power Level of A Phase Reactor','fontsize',14)
% set(gca,'fontsize',14)