masquage.m

来自「A partir d un son, permet de reperer les」· M 代码 · 共 46 行

%fonction renvoyant la puissance de masquage

function [Pm] = Masquage (f2,f1,P1) %f1 et f2 en hz

load TonaleTraited.mat;

f1B=bark(f1);

f2B=bark(f2);

if (((f1B-f2B)<-3) ||((f1B-f2B)>=8) )
    Pm=0;
    return
   
end
if (TonaleAtonale(f1))
    a=-1.525-0.275*f1B-4.5;
   
else
    a=-1.525-0.175*f1B-0.5;
end

f=f1B-f2B;

v=0;

if( (f>=-3) & (f<-1) )
    v=(17*(f+1)-(0.4*P1+6));
end

if( (f>=-1) & (f<0) )
    v=((0.4*P1+6)*f);
end

if( (f>=0) & (f<1) )
    v=(-17*f);
end

if( (f>=1) & (f<8) )
    v=(-(f-1)*(17-0.15*P1)-17);
end
    

Pm=P1+a+v;