%?Constants?used?in?the?MPEG-1?psycho-acoustic?model?#1
%?（for?Fs=44100?Hz）
%
%?Copyright?N.?Moreau?ENST?Paris?19/3/02

global?LTq_i?LTq_k?Table_z?Frontieres_i?Frontieres_k?Larg_f

%?“i“?to?“k“
i1?=?1:48;
i2?=?（49:72）-48;
i3?=?（73:108）-72;
i_to_k?=?[i1?（2*i2+48）?（4*i3+96）];

%?Absolute?auditory?threshold?as?a?function?of?“i“
LTq_i?=?[25.87?14.85?10.72??8.50??7.10??6.11??5.37??4.79?...
?????????????4.32??3.92??3.57??3.25??2.95??2.67??2.39??2.11?...
?????????????1.83??1.53??1.23??0.90??0.56??0.21?-0.17?-0.56?...
????????????-0.96?-1.38?-1.79?-2.21?-2.63?-3.03?-3.41?-3.77?...
????????????-4.09?-4.37?-4.60?-4.78?-4.91?-4.97?-4.98?-4.92?...
????????????-4.81?-4.65?-4.43?-4.17?-3.87?-3.54?-3.19?-2.82?...
????????????-2.06?-1.32?-0.64?-0.04??0.47??0.89??1.23??1.51?...
?????????????1.74??1.93??2.11??2.28??2.46??2.63??2.82??3.03?...
?????????????3.25??3.49??3.74??4.02??4.32??4.64??4.98??5.35?...
?????????????6.15??7.07??8.10??9.25?10.54?11.97?13.56?15.31?...
????????????17.23?19.34?21.64?24.15?26.88?29.84?33.05?36.52?...
????????????40.25?44.27?48.59?53.22?58.18?63.49?68.00?68.00?...
????????????68.00?68.00?68.00?68.00?68.00?68.00?68.00?68.00?...
????????????68.00?68.00?68.00?68.00];

%?Absolute?auditory?threshold?as?a?function?of?“k“
LTq_k?=?zeros（250?1）;
for?i?=?1:length（LTq_i）
????LTq_k（i_to_k（i））?=?LTq_i（i）;
end
for?k?=?1:250
????if?LTq_k（k）?==?0
????????LTq_k（k）?=?last_nonzero;
????else
????????last_nonzero?=?LTq_k（k）;
????end
end

%?axe_freq?=?（0:249）*Fe/512/1000;
%?figure（1）;?hold?off
%?plot（axe_freq?LTq_k）;?hold?on
%?axis（[0?20?-10?70]）

%?Bark?frequencies?as?a?fucntion?of?“i“
Table_z?=?[?.850?1.694?2.525?3.337?4.124?4.882?5.608?6.301?...
?????????????6.959??7.581??8.169??8.723??9.244??9.734?10.195?10.629?...
????????????11.037?11.421?11.783?12.125?12.448?12.753?13.042?13.317?...
????????????13.578?13.826?14.062?14.288?14.504?14.711?14.