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 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51
| % Calcul du périmètre
% d= 24 et w=60
function [P] = perimetre (d, w)
P = 2*(d+w)
% Calcul de la superficie
% h=24 et w=60
function [S] = superficie (h, w)
S = w*h
% calcul de TL [i]
n = input (' donner le nombre des itérations n');
switch n
case 1
disp ('f=125 Hz et alpha= 0.11')
TL[1]= 1.5 * (P/S) * alpha * l
break
case 2
disp ('f=250 Hz et alpha = 0.19')
TL[2]= 1.5 * (P/S) * alpha * l
break
case 3
disp ('f=500 Hz et alpha= 0.41')
TL[3]= 1.5 * (P/S) * alpha * l
break
case 4
disp ('f=1000 Hz et alpha= 0.54')
TL[4]= 1.5 * (P/S) * alpha * l
break
case 5
disp ('f=2000 Hz et alpha= 0.6')
TL[5]= 1.5 * (P/S) * alpha * l
break
case 6
disp ('f=4000 Hz et alpha= 0.75')
TL[6]= 1.5 * (P/S) * alpha * l
end
% Traçage de la courbe TL = f (f)
% TL (dB) :transmission Loss
% f (Hz) : Fréquence en bonde d'octave
figure;
plot (f, TL)
title ( 'Transmission Loss TL (dB) en fonction de la fréquence f (Hz)');
xlabel ('fréquence (Hz)');
ylabel ('TL (dB)');
% Echelle logarithmique sur x
x = [125 250 500 1000 2000 4000];
semilogx (x, exp(x)); |
Partager