Приложение

Листинг программы для расчетов.

clear; clc;

V = 480;

S0 = 32;

w = [63 125 250 500 1000 2000 4000 8000];

Lp = [76 81 75 71 68 66 60 54];

Ln = [71 61 54 49 45 42 40 38];

Bp = [15 18 24 28 34 48 51 96];

%Подходящие конструкции

R0_1 = [19 23 27 31 35 37 36 39]; %сталь толщиной 3 мм

R0_2 = [23 27 31 35 37 32 39 43]; %сталь толщиной 5 мм

R0_3 = [27 30 35 37 37 39 43 47]; %сталь толщиной 10 мм

R0_4 = [21 25 27 30 32 30 36 40]; %стекло силикатное толщиной 6 мм

Rt = [];

Lr_1 = []; Lr_2 = []; Lr_3 = []; Lr_4 = [];

count1 = 0; count2 = 0; count3 = 0; count4 = 0;

for i = 1:8

Rt(i) = Lp(i) - Ln(i) - 10*log10(Bp(i)) + 10*log10(S0);

if Rt(i) < R0_1(i)

count1 = count1 + 1;

end

if Rt(i) < R0_2(i)

count2 = count2 + 1;

end

if Rt(i) < R0_3(i)

count3 = count3 + 1;

end

if Rt(i) < R0_4(i)

count4 = count4 + 1;

end

Lr_1(i) = Lp(i) - 10*log10(Bp(i)) + 10*log10(S0) - R0_1(i);

Lr_2(i) = Lp(i) - 10*log10(Bp(i)) + 10*log10(S0) - R0_2(i);

Lr_3(i) = Lp(i) - 10*log10(Bp(i)) + 10*log10(S0) - R0_3(i);

Lr_4(i) = Lp(i) - 10*log10(Bp(i)) + 10*log10(S0) - R0_4(i);

end

figure(1);

semilogx(w, Lp, '-ro', w, Lr_1, '-bo', w, Lr_2, '-co', w, Lr_3, '-mo', w, Lr_4, '-yo', w, Ln, '-ko');

grid on;

xlabel('w, Гц');

ylabel('L_r(w), дБ');

title('Уровень шума до и после изоляции');

legend({'До изоляции', 'Сталь толщиной 3 мм', 'Сталь толщиной 5 мм', 'Сталь толщиной 10 мм', 'Стекло силикатное толщиной 6 мм', 'Предельный спектр уровня шума'}, 'Location', 'best');