Приложение
Программа для сравнения материалов:
clear;
clc;
A = 10;
B = 10;
H = 5;
V = A*B*H;
S = 2*(A*B + B*H + A*H);
S0 = A*B;
w = [63 125 250 500 1000 2000 4000 8000];
L_p = [82 88 95 90 91 98 81 73];
B_p = [17 19 23 29 36 50 65 100];
alpha1 = [.15 .28 .50 .50 .70 .70 .70 .65]; %минеральная вата
alpha2 = [.10 .12 .35 .40 .82 .90 .92 .88]; %базальтовое волокно
a_n = [];
G1 = []; G2 = [];
B_pa1 = []; B_pa2 = [];
d_L1 = []; d_L2 = [];
L_r1 = []; L_r2 = [];
for i = 1:8
a_n(i) = B_p(i)/(B_p(i) + S);
G1(i) = S*a_n(i) - S0*(a_n(i) - alpha1(i));
G2(i) = S*a_n(i) - S0*(a_n(i) - alpha2(i));
B_pa1(i) = S*G1(i)/(S - G1(i));
B_pa2(i) = S*G2(i)/(S - G2(i));
d_L1(i) = 10*log10(B_pa1(i)/B_p(i));
d_L2(i) = 10*log10(B_pa2(i)/B_p(i));
L_r1(i) = L_p(i) - d_L1(i);
L_r2(i) = L_p(i) - d_L2(i);
end
figure(1);
semilogx(w, L_r1, 'ro', w, L_r2, 'bo', w, L_p, 'mx');
grid on;
xlabel('w, Гц');
ylabel('L_r(w), дБ');
title('Уровень шума после обработки');
legend({'Минеральная вата', 'Базальтовое волокно', 'До обработки'}, 'Location', 'best');
Программа для сравнения площадей облицовки:
clear;
clc;
A = 10;
B = 10;
H = 5;
V = A*B*H;
S = 2*(A*B + B*H + A*H);
S0 = A*B;
S01 = A*B + 2*B*H;
S02 = A*B + 2*B*H + 2*A*H;
w = [63 125 250 500 1000 2000 4000 8000];
L_p = [82 88 95 90 91 98 81 73];
B_p = [17 19 23 29 36 50 65 100];
alpha1 = [.15 .28 .50 .50 .70 .70 .70 .65]; %минеральная вата
a_n = [];
G1 = []; G2 = []; G3 = [];
B_pa1 = []; B_pa2 = []; B_pa3 = [];
d_L1 = []; d_L2 = []; d_L3 = [];
L_r1 = []; L_r2 = []; L_r3 = [];
for i = 1:8
a_n(i) = B_p(i)/(B_p(i) + S);
G1(i) = S*a_n(i) - S0*(a_n(i) - alpha1(i));
G2(i) = S*a_n(i) - S01*(a_n(i) - alpha1(i));
G3(i) = S*a_n(i) - S02*(a_n(i) - alpha1(i));
B_pa1(i) = S*G1(i)/(S - G1(i));
B_pa2(i) = S*G2(i)/(S - G2(i));
B_pa3(i) = S*G3(i)/(S - G3(i));
d_L1(i) = 10*log10(B_pa1(i)/B_p(i));
d_L2(i) = 10*log10(B_pa2(i)/B_p(i));
d_L3(i) = 10*log10(B_pa3(i)/B_p(i));
L_r1(i) = L_p(i) - d_L1(i);
L_r2(i) = L_p(i) - d_L2(i);
L_r3(i) = L_p(i) - d_L3(i);
end
figure(1);
semilogx(w, L_r1, 'ro', w, L_r2, 'bo', w, L_r3, 'mo', w, L_p, 'kx');
grid on;
xlabel('w, Гц');
ylabel('L_r(w), дБ');
title('Уровень шума после обработки');
legend({'Потолок', 'Потолок + 2 стены', 'Потолок + 4 стены','До обработки'}, 'Location', 'best');