Moor ANSI S3.4-2005响度计算程序

根据Moor ANSI S3.4-2005编制的响度计算程序。另增加了扬声器频响对响度影响的评估功能

function E_ft = E_EP_filter（x0）
Filter_n=（1.8:0.1:38.9）;
a=size（Filter_n）;
E_ft（::）=zeros（a（2）2）;
for i=1:a（2）;
% i=139; 
E_ft（i1）=Filter_n（i）;
    E_ft（i2）=ft_n_l（Filter_n（i）x0）;
end
% ===response for excitation of every filter
function filter_L = ft_n_l（nL_ERB）
f_n=（10^（n/21.366）-1）/0.004368;
ERB_n=24.673*（0.004368*f_n+1）;
a=size（L_ERB）;
for i=1:a（1）;
    x=L_ERB（i2）;    
    f=（10^（（L_ERB（i1））/21.366）-1）/0.004368;
    g0=abs（f-f_n）/f_n;
    W_a_n_t=（L_ERB（i2）-94）/20;
    W_a_n=10^（2*W_a_n_t）;
    if f       p（i）=4*f_n/ERB_n-（0.0115889*（4*f_n/ERB_n））*（x-51）; 
       if g0<1
          g（i）=g0;
       else
          g（i）=1;
       end 
    else
        p（i）=4*f_n/ERB_n; 
       if g0<4
          g（i）=g0;
       else
          g（i）=4;
       end    
    end   
    filter_L_t（i）=（1+p（i）*g（i））*exp（-p（i）*g（i））*W_a_n;
end
filter_L=sum（filter_L_t）/（4*10^-10）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件         920  2014-11-12 11:46  Moor ANSI S3.4-2005\E_EP_filter.m
     文件        1340  2014-11-11 15:52  Moor ANSI S3.4-2005\L_ERB.m
     文件        4452  2014-11-12 23:00  Moor ANSI S3.4-2005\func_drp_free.m
     文件        3079  2014-11-04 14:04  Moor ANSI S3.4-2005\func_mid_ear.m
     文件        1953  2014-11-12 22:36  Moor ANSI S3.4-2005\func_spec_l.m
     文件         199  2014-11-12 23:12  Moor ANSI S3.4-2005\main_draft.m
     文件        2007  2014-11-12 21:52  Moor ANSI S3.4-2005\spec_signal_3rd_oct.m
     目录           0  2014-11-13 10:17  Moor ANSI S3.4-2005\