【投稿】微型器件仿真中微声管阻抗模型的建立： 考虑媒质粘滞运动方程，管末端修正与声管本身的体积顺性
【一位朋友的投稿。他是科班硕士，一线大公司研发经验。理论功底和仿真水平都很好。 欢迎各位来稿。(●'◡'●)】
【公式较多，直接粘贴会打乱格式，所以转换成图片。把“微管严格解与近似解绘图MATLAB代码”**在最后了，供参考。】
声学短管是设计动铁单元时经常会遇到的一种结构。无论是单体本身的导声管还是ITE或者BTE的模拟声管，我们都需要更为精确的结构模型以获得精确的模拟结果。对于声学管模型在ER 076B; ER122A; ER167中都有论述。本篇报告主要针对短声管（即声管长度小于波长的十分之一）的建模。建模时主要考虑三个方面：粘滞媒质的运动方程，管末端修正以及声管本身的微小体积顺性修正。
附录：
微管严格解与近似解绘图MATLAB代码
%%%%%Impedancefor micro-tube(strict solution&approximate solution)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%by Qing Wang
clear;
f=(10:10:20000);%frequency from 10Hz to 20kHz
k=(1-j)*(1.2*3.14*f/1.8*100000).^0.5;
a=besselj(0,k.*0.00055245);%0 bessel function
b=besselj(1,k.*0.00055245);%1 bessel function
t0=f./f;% 1
t1=(2*b)./(k.*a*0.00055245);
z=(-j*2*1.2*f.*0.0012/(0.00055245^2)).*((t0-t1).^-1);%acoustic impedance
z0=abs(real(z));%acoustic resistance
z1=abs(imag(z));%acoustic reactance
%%%%%drawing%%%%%
figure;h=plot(f,z0);gridon;
xlabel('Frequency');ylabel('Acoustic Resistance');title('Acoustic Resistance');
set(h,'LineSmoothing','on')
figure;h=plot(f,z1);gridon;
xlabel('Frequency');ylabel('Acoustic Reactance');title('Acoustic Reactance');
set(h,'LineSmoothing','on')
%%%%%dataexporting%%%%%
z0=z0';z1=z1';f=f';
