eqccom.m

这是一个用于语音信号处理的工具箱

function   [L,R,C,LR,RR]=eqccom(af,rlen,nsec)

ro=0.00114;    % air density
cv=35300;      % sound velocity
mu=0.000186;   % shear viscosity
ro05=ro/2;
rocc=ro*cv*cv;
f_r=0.1531;
f_rr=57.99;
f_lr=0.000546;


for i=1:nsec
    L(i)=ro05*rlen/af(i);
    C(i)=af(i)*rlen/rocc;
    R(i)=f_r*rlen/af(i)/sqrt(af(i));    
end

RR=f_rr/af(nsec);
LR=f_lr/sqrt(af(nsec));