Matlab---微穿孔板吸声系数MATLAB仿真代码.rar

f=0:10:1000; w=2*pi*f; c=340; d=0.8; t=2; D=0.06; p=0.78; u1=1.56; k=d.*sqrt(f./10); kr=sqrt(1+k.*k./32)+sqrt(2).*d./(8.*t); r1=32.*u1.*t.*kr./(p.*c.*d.*d); km=1+1./sqrt(9+k.*k/2)+0.85.*d./t; m=0.294e-3.*t.*km./p; a1=4.*r1./((1+r1).^2+(w.*m-cot(w.*D./c)).^2); plot(f,a1,'r','LineWidth',3); xlabel('频率 (Hz)'); ylabel('垂直入射吸声系数'); hold on