声学波导管

声学波导管 ------------------------------------------------------------------------------------------------ 声学波导管 食不厌精 脍不厌细 1、恒定截面波导内的声传播 1.1、矩形波导管 1.2 、圆柱形波导管 设有一半径为a的圆柱形管，一端延伸到无限远。圆柱形管的声 波方程应以柱坐标系来描述。设管的径向坐标为r，极角为?，管轴用 z来示。直角坐标与柱坐标之间有如下关系 ?x?rcos???y?rsin? ?z?z? 而柱坐标系的拉普拉斯算符可表示为 1??1?2?2 (r)?2?2 (1-2-1) ??2r?r?rr???z2 于是三维声波动方程就可变换为: 1??p1?2p?2p1?2p(r)?2?? (1-2-2) r?r?rr??2?z2c2?t2 根据分离变量法，令解 p(r,?,z,t)?R(r)?(?)Z(z)ej?t, 将其代入(1-2-2)式可得如下三个常微分方 程 ?d2Z?kz2Z?0?2dz?d2?? ??m2??0 (1-2-3) 2d??22?dR?1dR?(k2?m)R?0r22?rdrdrr? 其中 k?2?2 —————————————————————————————————————— ------------------------------------------------------------------------------------------------ c2?kz2?kr2. (1-2-4) 由于圆柱管道向无限远处延伸，对于Z的方程可取行波解: Z(z)?Aze?jkzz; (1-2-5) 对于?的方程可取解为 ?(?)?A?cos(m???m), (1-2-6) 因为?(?)??(??2?)的关 系应该满足，所以式中m一定要为正整数。 对于R的方程我们作一适当变换，令krr?x,则方程就化为 d2R1dRm2 ?(1?2)R?0. (1-2-7) 2?dxxdxx 这是一个标准的m解贝塞尔方程，其一般解可表示为 R(krr)?ArJm(krr)?BrNm(krr), (1-2-8) 这里Jm(krr)与 Nm(krr)分别代表宗量为(krr)的m阶柱贝塞尔函数与柱诺伊曼函数。 按照柱诺伊曼函数在零点发散的性质，式中应取Br?0,于是(1-2-8) 式简化为 R(krr)?ArJm(krr)， (1-2-9) 由此求得 管中声压解为: pm?AmJm(krr)cos(m???m)ej(?t?kzz), (1-2-10) 由运动 方程j??Ur?? vrm??p可求得对应的径向速度 为: ?rj?pmjkdJ(kr)?Amr[mr]cos(m???m)ej(?t?kzz), (1-2-11) ?0??r?0?d(krr) —————————————————————————————————————— ------------------------------------------------------------------------------------------------ 设管壁为刚性，即在r?a处有vr?0,由此条件可得知如下关系: [ 按照贝塞尔函数的递推关系 dJm(krr)](r?a)?0, d(krr) ?dJm(x)1?dx?2[Jm?1(x)?Jm?1(x)] ?dJ0(x)???J1(x)dx? 可得到如下圆柱声波导的本征方程: Jm?1(kra)?Jm?1(kra) (m?0) J1(kra)?0 (m?0) 利用MATLAB可从这些方程解得一系列根植，部分根植列于下表 表 1.圆柱声波导本征值 在刚性壁条件下，kr应有一系列特定的数值，此特定值可用下标 m与n两个正整数表示，我们写成kr?kmn.在k?kmn时声压解可写成 如下形式 pmn?Amncos(m???m)Jm(kmnr)ej(?t?kzz), (1-2-12) 其中 2 kz?k2?kmn (1-2-13) 当k?kmn时，圆柱管中存在非传播形式的高次模式，这些高次 模式会随距离衰减，此时声压解可写成如下形式 pmn?Amncos(m???m)e??mnzJm(kmnr)ej?t, (1-2-14) 其中 222??kmn?k2 kz??jkmn?k ， mn —————————————————————————————————————— ------------------------------------------------------------------------------------------------ 当波导管的声源进行极轴对称振动时，即波导管中的声压与极角? 无关，因此我们可以取m?0,当k?kn时得到声压解为 pn?AnJ0(knr)ej(?t?kzz) (1-2-15) 其中 22 kz?k?kn (1-2-16) 同理，当k?kn时声压解可表示为 pn?Ane??nzJ0(knr)ej?t, (1-2-17) 其中 22kz??jkn?k2 ，?n?kn?k2 根据上表，我们可以与矩形管类似地得到圆柱形声波导管的截止 频率为 fc?f10?1.c0 (1-2-18) 2?a c0 (1-2-19) 2?a如果已知声源做极轴对称的 振动，则m?0,于是可以确定 fc?f01? 考 虑有限长圆环形声波导管的情况。设圆环的内径为a1，外径为a,长 度为l.假设入射声源做极轴对称的振动，则柱坐标系下的声波方程为 1??p?2p1?2p(r)?2?22 (1-2-20) r?r?r?zc?t 当圆环状波导管无限长时，根据分离变量法，令解 p?R(r)Z(z)ej?t， (1-2-21) 代入(1-2-20) 式可得到如下两个微分方 程 ?d2Z?kz2Z?0,?2?dz ?2 (1-2-22) dR1dR2???kR?0,r2?drrdr? —————————————————————————————————————— ------------------------------------------------------------------------------------------------ 如同半无限长圆柱形声波导管一样，我们可以求得: Z(z)?Aze?jkzz, R(krr)?ArJ0(krr)?BrN0(krr), 由于我们求解的是圆环形声波导管，我们无法根据柱诺伊曼函数 的零点发散性质使其系数Br?0,因此我们可以将其声压解表示为 ??? pn?[AnJ0(knr)?BnN0(knr)]ej(?t?kzz), (1-2-23) 其径向速度可表示为 v??rn 设内外管壁均为刚性，即在r?a,r?a1处有vr?0,由此条件可得到如 下关系 ?jkn[A?J1(knr)?B?N1(knr)]ej(?t?kzz), (1-2-24) nn?? Bn???An? 此时式(1-2-23)可简化为 J1(kna) N1(kna) ?j(?t?kzz)p??Ae?n(r), nn 其中 ?n(r)?J0(knr)?J1(kna)Y0(knr) (1-2-24) N1(kna) 在考虑圆环形声波导管的长度为有限长l时，由于管末端突变界 面的影响，此时管中将存在沿z轴负方向传播的反射波，根据以上对 入射波声压的求解，我们可以类似的得到反射声波声压 ?j[?t?kz(z?l)] p??Ae?n(r), (1-2-25) nn 由此，我们可以得到圆环形声波导管内的声压解为 ????jkzz?jkz(z?l)p(r,z)??[Ane?Ane]?n(r)ej?t??n?0 ? —————————————————————————————————————— ------------------------------------------------------------------------------------------------ (1-2-26) ??jkzzjkz(z?l)???P(r,z)??[Ane?Ane]?n(r)?n?0? 轴向声速解为 1????jkzz?jkz(z?l)v(r,z)?kz[Ane?Ane]?n(r)ej?t?????n?0 ? (1-2-27) ?1?jkzjk(z?l)???U(r,z)?kz[Anez?Anez]?n(r)????n?0? 根据径向声速在r?a1时为0，我们可以得到 J1(kna1)? 其中kz表示水平波数，并且有 J1(kna)Y1(kna1)?0 (1-2-28) N1(kna) kz?[k?k] 2122n 当k?kn时，当k?kn时，kz为实数，e?jkzz和ejkz(z?l)均不会沿着 传播方向衰减， 取kz为负虚数(此处不取正虚数的原因在矩形声波导管中已经 作了解释)，即 2 kz?-jkn?k2 (1-2-29) 此时，e?jkzz和ejkz(z?l)均会沿着传播方向衰减。因此在圆环形 声波导管中，声压是能够在管中正常传播模式和非正常传播模式的叠 加。 接下来，我们考虑圆环内径a1?0的情况，此时即为有限长的圆 柱形声波导管，结合半无限长圆柱形声波导管和有限长圆环形声波导 管的分析，我们可以求得管中的声压解为 ??jkzz-jkz(z-l) P(r,z)??[Ane?Ane]?n(r) (1-2-30) —————————————————————————————————————— ------------------------------------------------------------------------------------------------ n?0? 轴向声速解为 U(r,z)? 其中 ?n(r)?J0(knr) (1-2-32) 设波导管的管壁为刚性，即径向速度在r?a处为0，得其本征方程 J1(kna)?0 (1-2-33) 其轴向波数的分析同圆环形波导管。 计算单向传播的某高次波声场的声强， 其中 p??AnejkAnz?Am(r)ejwt,v?? I??1??n?0??jkzz-jkz(z-l)k[Ae?A]?n(r) (1-2-31) ?znne?1kAnAnejkAnz?Am(r)ejwt ??w111*2Ids?Re[pv]ds??Re?(k An)*A2J0(?n) ??SS2SS2?w 下面考虑一段是振动活塞，另一端为无限长度的柱形波导管模型，在已知波导开口处活塞速度分布的情况下，求解管内的声场。 设振速是轴对称分布的，由傅里叶-贝塞尔展开，把表面声压或振速函数 f(r)展成以特征函数Jm( ?mn a r)为特征函数的级数- ??0 f(r)??bmnJn( —————————————————————————————————————— ------------------------------------------------------------------------------------------------ ?mnr a ) ?mn由边界条件，应用柱贝塞尔函数的正交性求解bmn- m?n a ?rJ m ( ?mnr a )Jm( ?mnr a )dr 换元使得 ?mnr a ?x,?mn?x0 上式变为 a ?rJm( —————————————————————————————————————— ------------------------------------------------------------------------------------------------ ?mnra )Jm( x0 ?mnra )dr ?2(?2(?2(?2( 1 mna1 )2 22 [J(x)]dxm?0 mna1 )2 xJ(x0)? ( 22x0Jm(x0)? 202m 1 x0 mna1 )2 2'''2'[xJ(x)?xJ(x)?mJ(x)]J(x)dxmmmm?0 —————————————————————————————————————— ------------------------------------------------------------------------------------------------ x0 mna1 )2 ( xJ(x0)? 2 2m mna1 )2 2''''[xJ(x)?xJ(x)]J(x)dx?mmm?0 m2( x0 mna )2 ?J m (x)dJm(x) x0 2 mna )2 2( —————————————————————————————————————— ------------------------------------------------------------------------------------------------ mna )2 ?d[xJ(x)]? '2 m m22( mna )2 x02 Jm(x)|0 12m21' ?(a?)Jm(a)2?a2[Jm(a)]2 22(mn)2a 常有以下三类边界条件 ' 第一类，硬边界条件下，已知声压的分布函数，Jn(?mn)?0，此时，为数学 上的诺依曼问题，bmn?2 a2(1?m2 2?mna2)Jm(?mn)0?f(r)Jm(?mnra)rdr。 第二类，软边界条件下，已知振速的分布函数，Jm(?mn)?0，此时，为数学2?mnr上的狄利克雷问题，bmn?22f(r)J()rdr。 maJm?1(?mn)?a0 —————————————————————————————————————— ------------------------------------------------------------------------------------------------ 第三类，相当于弹性支撑下的边界条件，a?mn a'HJm(?mn)?Jm(?mn)?0。此时 bmn?12a4 (a?2?22)2?mn?mnH2a2m2aJ(?mn)?f(r)Jm(02m?mnra)rdr。 1.4、波导管的T型网络类比 2、缓变截面波导内的声传播 2.1、缓变截面管内的声场解析表达式 2.2、指数型、悬链型、锥形缓变截面号筒 3、突变截面波导内的声传播 3.1、平面波假设下的常见突变截面管模型 考虑波导管是轴对称分布下的刚性壁管道，沿着轴向不同位置的垂直横截面是相似的，横截面积之时轴向的一维函数s?s(x)，突变截面管s(x)为分段连续函数，以在间断点处的截面分界，分别分析和联系，对传播特性进行研究。 3.1.1、模型一 图一 tI?It4s2?2(s?) 122Ii(s12?1)s1 2(s12?1)2 隔声量nI?采用收缩式截面管，nI?1，突变的截面相当于介质特4 性阻抗的突变，起到反射作用。 3.1.2、模型二 图二 —————————————————————————————————————— ------------------------------------------------------------------------------------------------ 由边界上的声压连续和体积速度连续 ?pa?pb?pc1?pd1?Ua?Ub?Uc1?Uc2???pc2?pd2?pe?Uc2?Ud2?Ue?? pc2?pc1e?jkl ?jkl?p?pec2?d2 ? pe2S1S2?2pa2S1S2coskl?j(S12?S2)sinkl pe24tI?2?pa4cos2(kl)?(S21?S12)2sinkl 1(S21?S12)2sinkl2nI??cos(kl)?tI4 (S21?S12)2 nImax? 4 nImax?1,无论扩张式亦或收缩式，其相当于中间插入过度介质。 nImin?1,此时 f?ncn?,l?(n?N?) 2l2 这种结构可以实现对声波的全投射却无法实现全阻隔。 在 中间管道中的声场 p?p3?p4? ?2S1(S2?S1)p1cos[k(x?l)]ejwt(驻波项)222S1S2coskl?j(S1?S2)sinkl2 12Sp1ej[wt?k(x?l)](行波项)222S1S2cokls?j(S1?S2)sinkl 可见在突变截面下，驻波项和行波项是永恒存在的，而当S2?S1 时，没有了反射，驻波项消失，投射系数为1。 考虑这种波导隔声量的Q值，使得nI?nImax 2 S1S22 ?)S2S1 ?co2skl —————————————————————————————————————— ------------------------------------------------------------------------------------------------ 2[(1?2)2?4]S2S1 ( 2 c2(S12?S2)carcc2)，以第一处峰值f?计算 f?022?l2(S1?S2)4l 2 cc2(S12?S2)f0 带宽?f??arcc2)，Q??2 2l?l2(S1?S2)?f 1 2 2(S12?S2) 2?arcc2)2 ?2(S1?S2) 4 在 S1 ?1的前提下，提高该比例，会使得带宽加大，同时隔声量nI减 小，改S2 善消声带宽和提高消声量是一对矛盾。 模型二,平面波假设 NI/dB 200400600800 —————————————————————————————————————— ------------------------------------------------------------------------------------------------ 1000f/hz 12001400160018002000 图三 3.1.3、模型三 图四 ?p1?p2?p3?p4?p?p?p?125??UA?UB?Uc ??v|x?d?0 p2SA?5?p1SA?SC?jSBtankd 22p54SA(SA?SC)2(SBtankd)2tI?2?,nI??22p1(SA?SC)2?(SBtankd)24SA4SA 当kd?? 2?n?(n?N)时，tImin?0,nImax??? 24SAtImax?,2(SA?SC) 和模型二相反，SB越大，共振空气柱储能越多，消声量提高，消声带宽减小，显然其对消声量和带宽的影响与模型二相反，考虑将二者综合来讨论消声情况。 model4,planewave hypothesis 0NI/dB2004006008001000 f/hz12001400160018002000 3.1.4、模型四 图六 ?pa1?pa2?pb1?pb2?p?p?p?pa2c1c2?a1 —————————————————————————————————————— ------------------------------------------------------------------------------------------------ ??pc1?pc2?pd1??U?U?U?U?U?Ua2b1b2c1c2?a1 ?Uc1?Uc2?Ud1???v|(x1?d)?0 令Sa?Sd pd1(2SaSccoskl?SaSbsinkltankd)?j[Ssinkl?Ssinkl?SbSccoskltankd]?pa12Sa Sc tI?|pd12|pa12c2a 模型三，平面波假设 图七 二者都是在插管长度为四分之一波长的奇数倍对应的频率处达到最大的消声量，显然通过有限个插管B的组合是不能够完成对模型三所有的谷点进行补偿的。但是可以通过应用有限个插管，来对低于某个频率的传播损失曲线进行一个均衡。 3.2、考虑高次模式下的突变截面管内的声传播 3.2.1、 高次模式简正波声压系数的确定 为确定高次波的系数，需要联立声压连续和体积速度连续的方程和方程组，根据1中有关高次模式下声场的解，可以看到当声波频率小于某号简正波的截止频率时，该号简正波不再有行波特性，而呈现随距离的衰减特性，而且可以看到，2?k2越大，该号简正波对声场的影响简正波号序号数越大，衰减系数?mn?kmn 越小。故实际操作中可以取有限号简正波来近似对声场进行描述，显然，取得总的号数越多，与实际声场越逼近。通过柱贝塞尔函数的—————————————————————————————————————— ------------------------------------------------------------------------------------------------ 正交归一性，得到系数矩阵，用高斯消元法求解系数矩阵，通过编程 进行实现。 针对柱贝塞尔函数的正交归一性，做出以下的归纳和证明，设 B0(x)是任意可能的柱贝塞尔函数，如J0(x)或Y0(x)(又记为N0(x))， 那么有 ?B0(?x)B0(?x)xdx? 证明如下: r[?B0(?x)B1(?x)??B0(?x)B1(?x)](???) ?2??2 ?B(?x)B(?x)xdx00 ??B0(?x)d(1 ? xB(?x)B1(?x)xB(?x)?0??1dB0(?x)?? xB(?x)B1(?x)?xB1(?x)B1(?x)?0??dx?? xB(?x)B1(?x)?xB1(?x)?0??dB0(?x)??2 ?xB0(?x)B1(?x)xB1(?x)) ??xB1(?x)B0(?x)?2 ???2B0(?x)B0(?x)xdx2?? ? ?2xB(?x)B1(?x)?xB1(?x)B0(?x)(1?2)?B0(?x)B0(?x)xdx?0????2 ? 原式?r[?B0(?x)B1(?x)??B0(?x)B1(?x)](???)22??? 当???时，对上式取极限， r[?B0(?x)B1(?x)??B0(?x)B1(?x)] ????2??2lim —————————————————————————————————————— ------------------------------------------------------------------------------------------------ ?r2B0(?x)B0(?x)??r2B1(?x)B1(?x) ?lim ???2? r22?(B0(?x)?B12(?x))2 B0(x)和B1(x)可以是不同类型的柱贝塞尔函数，如- r2 ?[J0(?x)Y0(?x)?J1(?x)Y1(?x)]2?J0(?x)Y0(?x)xdx 根据此正交性公式，并利用边界条件，可以列出系数方程组。 3.2.2、考虑高次模式下常见突变截面管模型的消声量随频率的变 化 3.2.2.1、考虑模型二，以最常见的收缩式突变截面管模型为例。 各部分的声压和振动速度的表达式为: ?PA(r,z1)?e?jkz1??Ane ?PB(r,z1)??(Bne?jkB,nz1jkA,nz1?A,njkB,n(z1?Lb)??Bne)?B,n PC(r,z2)??Cne?jkC,nz2?C,n ?P由j??U??得:?z 1?jkA,nz1UA(r,z1)?[ke?jkz1??kA,nAne?A,n] ?? 1UB(r,z1)????kB,n(Bn?e?jkB,nz1??BnejkB,n(z1?Lb))?B,n 1?jkC,nZ2UC(r,z2)?kCe?C,n?C,nn?w 利用截面上的声压和体积速度连续- (Bs??Be(Bs?e ??jkB,sLbs )??B,s?B,s?B??r?B,sdr??An???A,n?B,s?B a1 ?jkB,sLb ?Bs?)??B,s?B,s?B??Cn??C,n?B,s?C —————————————————————————————————————— ------------------------------------------------------------------------------------------------ ??jkB,nLbn a ?k?k ? sA,s B,n (B?Be ?n )??B,n?A,s?B?k?r?A,sdr?Ak ??A,s?A,s?A ?(BB,nne ?jkB,nLb ?Bn?)??B,n?C,s?B?kC,sCs??C,s?C,s?C 其中 ?A,n?J0( ?B,n?J0( ?nr a ),其中J1(?n)?0),其中J1(?n)?0 ),其中J1(?n)?0 ?nr a1 ?C,n?J0( —————————————————————————————————————— ------------------------------------------------------------------------------------------------ ?nr a1 利用柱贝塞尔函数的正交性，上式可展成: a1?n?3 a?a?J(?)J()1n0s1??12212??jkB,sLb?)?a1J0(?s)?a1s?0??An??(Bs?Bse2222 22a1?n?a?s ??a1?n3 a?a?J(?)J()?1n0s1122??jkB,sLb??Bs)?a1J0(?s)??Cn?(Bse2 2a12?n?a?s2? ? a1?s3?a?a?J(?)J()1s0n1 ?k(B??B?e?jkB,sLb)??k1a2s?0?A?k1a2J2(?) nsA,s0s2??B,nna12?s2?a2?n22 ? a1?s?3 a?a?J(?)J()1s0n1?jkB,sLb?122??k(Be?B)??kCaJ0(?s)nC,ss??B,nn2222 a??a?21sn? 共有4×5=20个位置数，每个方程可列出5个等式，有20个等 式，设20×20的矩阵A，每个方程占5行，第一个方程1-5行，第 二个方程6-10行...对应每个方程的每一行代表s为1-5，每个系数占 5列，A占据1-5列，B?占据6-10列...对应的每个系数每一列代表不 同n.X为20×1矩阵，从上到下代表 —————————————————————————————————————— ------------------------------------------------------------------------------------------------ A0?,A1?,...B0?,B1?,...B0?,B1?,...,C0?,...C4?。AX=b.A是稀疏矩阵，用LU分解的方法求解方程组，得到更精确的解。 mode2 with higher modes being considered NL/dbf/hz 可见，在频率小于半径最大的圆柱管(即A管和C管)的(0,1)号简正波的截止频率(f01??0c?2997Hz)之前，主要是波导内除了平面波外主要是非2?a 传播的高次模式在起作用，而且随着频率接近f01，峰值向左偏移量越大，峰值的高度也有所增加。在频率大于f01时，可传播的行波模式的高次波开始起作用，平面波下的吸声量周期规律不复存在，出现的吸声峰的“破碎”。 可以得到以下结论，在可传播高次波出现之前，根据峰值处fn? 可以计算重复的峰的个数约为为(2n?1)c?f01， 4ll2?01l(d?2a)，约为2.440，由此可知当d?d l?0.41时，这种重复的峰会不再出现。如下图所示。 d mode2 with higher modes being considered,Lb/(2*a)=0.07 NL/dbf/hz 考虑方矩形波导，计算模型二传播损失随声波频率的关系。 管壁仍为刚性,依次列写的声压和速度的表达式与上式一致，此处 m?xn?y))aa m?xn?y?Bmn?)) bb m?xn?y?cmn?))aa?Amn? 由声压和体积速度连续- —————————————————————————————————————— ------------------------------------------------------------------------------------------------ ??(Bpq?Bpqe ?(Bpqe?jkBpqLb?jkBpqLb)??Bpq?Bpq?B???Bpqdxdy??Amn??Amn?Bpq?BS??Bpq)??Bpq?Bpq?B??Cmn??Cmn?Bpq?B ???jkBmnLb??Bmne)??Bmn?Apq?B?k??Apqdxdy?ApqkApq??Apq?Apq?A?kBmn(Bmn S ??jkBmnLb?k(Be?B?Bmnmnmn)??Bmn?Cpq?B?kCpqCpq??Cpq?Cpq?C bm?bn??24p?q4abmn(?1))sin()2??b(2p??sin(2p?))(2q??sin(2q?))??jkBpqLb2?)??4bp,q?0??Amn??(Bpq?Bpqe22222222224pq??(bm?ap)(bn?aq)??bm?bn?24p?q4abmn(?1)))2???jkLb(2p??sin(2p?))(2q??sin(2q?))?Bpqb?Bpq)???Cmn?(Bpqe24pq??2(b2m2?a2p2)(b2n2?a2q2)??bp?bq?24m?n?4abpq(?1)sin()sin()2?k(B??B?e?jkBmnLb)??4ka2p,q?0?A?ka(2p??sin(2p?))(2q??sin(2q?)) pqApq??Bmnmnmn?2(b2p2?a2m2)(b2q2?a2n2)4pq?2 ?bp?bq?24m?n?4abpq(?1)sin()sin()?a2(2p??sin(2p?))(2q??sin(2q?))??jkBmnLb??Bmn)??kCpqCpq222222222??kBmn(Bmne?(bp?am)(bq?an)4pq?2? 3.2.2.2、在考虑高次模式下重新计算模型四- 列出各部分的声压和振速- ?PA(r,z1)?e?jkz1??Ane ?PB(r,z1)??(Bne?jkB,nz1 —————————————————————————————————————— ------------------------------------------------------------------------------------------------ ?jkC,nz1jkA,nz1?A,njkB,n(z1?Lb)??Bne??Cne)?B,n)?C,n?PC(r,z1)??(CnejkC,n(z1?LC) PD(r,z2)??Dne 由j??U???jkD,nz2?D,n ?P得:?z 1?jkA,nz1UA(r,z1)?[ke?jkz1??kA,nAne?A,n] ?? 1UB(r,z1)? UC(r,z1)? UD(r,z2)???1?k(B?B,nne?jkB,nz1??BnejkB,n(z1?Lb))?B,n)?C,n??1?k(C?C,nne?jkC,nz1??CnejkC,n(z1?Lc) ???kD,nDn?e?jkC,nz2?D,n 根据声压和体积速度连续条件- ? 根据B管的刚性边界条件，Bne?jkB,nLb????Bn?0,得 Bn?BnejkB,nLb。 (Cs??Ce (B?Be (Cse a?s??jkC,sLbs?)??C,s?C,s?C??r?C,sdr??An??A,n?C,s?C0?)??B,s?B,s?B??r?B,sdr??An??A,n?B,s?Ba1aa2??jkB,sLbs??jkC,sLc?Cs?)??C,s?C,s?C??Dn??D,n?C,s?C?jkB,nLb??)??B,n?A,s?B??kC,n(Cn?Cne?jkC,nLb )??C,n?A,—————————————————————————————————————— ------------------------------------------------------------------------------------------------ s?C?kB,n(Bn??Bn?e ?k(C?C,nne?k?r?A,sdr?As?kA,s??A,s?A,s?A?jkC,nLc??Cn)??C,n?D,s?C?kD,sDs??D,s?D,s?D ?A,n?J0( ?B,n?J0( ?C,n?J0(?nra),其中J1(?n)?0)?J1(?n)?r?aJ(?)?aY0(n),其中 J1(n1)?1nY1(n1)?0Y1(?n)aaY1(?n)a?nra?nr a2 ),其中J1(?n)?0 ),其中J1(?n)?0?D,n?J0(?nr a a2?n?3a?a?J(?)J()2n0s1??12212??jkC,sLc?)?a2J0(?s)?a2s?0??An??(Cs?Cse222222a??a?s2n??a1?sa1?s222a(J(?)Y(?)?J(?)Y(?))?a(J(?)Y()?J()Y1(?s))2?1s0s0s1s11s00jkL?jkL?Bs?(eB,sb?eB,sb)22Y1(?S)??a?a?a??aa1?nJ1(1n)(J1(?s)Y0(1s)?J0(1s)Y1(?s))??1(a2?a2)s?0?A??1n22?2(?n??s)Y1(?s)?a2?n?3a?a?J(?)J()2n0s1???jkC,sLC122??Cs)?a2J0(?s)??Dn?(Cse22222a2?n?a?s? a1?sa1?na1?n?aa?J()[J(?)Y()?J()Y1(?n)]1s11n00?jkB,nLb?jkB,nLb??e)???kB,nBn(e22(???)Y(?)sn1n??a2?s3aa?J(?)J()?s0n12???jkC,nLb)???kC,n(Cn?Cne2222 a??a?n?2s?11??k?a2s?0?As?kA,s?a2J02(?S)22??a2?s3aa?J(?)J()?2s0n1??—————————————————————————————————————— ------------------------------------------------------------------------------------------------ jkC,nLC??kD?1a2J2(?)?Cn)??kC,n(CneD,ss0S2222?a??a?22sn? model 4 NL/db f/hz NL/db f/hzmodel 4 NL/db f/hz 3.2.2.3、模型四的异构体 将模型四的B管进行改进，从内嵌型和后置型两种异构体模型上 对吸声效果进行比较。 内嵌型 各部分的声压和振速: PA(r,z1)?e?jkz1? PB(r,z1)? PC PD PE1?n?n ?A?nejkA,nz1?A,njkB,n(z1?Lb)?(Be(r,z)??(Ce(r,z)??(De(r,z)??Ee1 ?n3n?jkB,nz1?jkC,nz1?jkD,nz2??Bne??Cne)?B,n)?C,n)?C,njkC,n(z1?Lb)?? DnejkD,n(z2?Ld)?jkE,nz3 ?E,n 由j??U?? UA(r,z1)? —————————————————————————————————————— ------------------------------------------------------------------------------------------------ UB(r,z1)? UC(r,z1)? UD(r,z2)? UE(r,z3)???1?P得:?z1jkz?[ke?jkz1??kA,nAneA,n1?A,n] B,n?(Bne?jkB,nz1??BnejkB,n(z1?Lb)k???k???11)?B,n)?C,n )?D,nC,n?(Cne ?jkC,nz1??CnejkC,n(z1?Lb)??1?e?kD,n(Dn?jkD,nz2??Dne?jkD,n(z2?Ld)k? E,nEne?jkE,nz3?E,n 根据刚性边界条件的得到系数B的关系。 利用声压和体积速度连续 (C?Ce(Bs??Be ?s??jkC,sLbs ? )??C,s?C,s?C??r?C,sdr??An??A,n?C,s?C ?)??B,s?B,s?B??r?B,sdr??An??A,n?B,s?B a1 ?)??D,s?D,s?D??(Cne ?jkC,nLb ? ?Cn)??C,n?D,s?D a a1 ??jkB,sLb s —————————————————————————————————————— ------------------------------------------------------------------------------------------------ (Ds??Ds?e(Ds?e ?jkD,sLd ?jkD,sLd ? ?Ds?)??D,s?D,s?D??En??D,s?E,n?D ?jkB,nLb ?? )??B,n?A,s?B??kC,n(Cn?Cne ?jkC,nLb ?k B,n ?? (Bn?Bne )??C,n?A,s?C a ?k?r?A,sdr?As?kA,s??A,s?A,s?A kC,s(Cs?e ?jkC,sLb ?? ?Cs?)??C,s?C,s?C??kD,n(Dn?Dne ?jkD,nLd )??D,n?C,s?D —————————————————————————————————————— ------------------------------------------------------------------------------------------------ ?k D,n ? (Dne ?jkD,nLd ? ?Dn)??D,n?E,s?D?kE,sEs???E,s?E,s?E ?A,n?J0(?B,n?J0(?C,n?J0( ?nr a ),其中J1(?n)?0)? J1(?n)?r?aJ(?)?a Y0(n),其中J1(n1)?1nY1(n1)?0Y1(?n)aaY1(?n)a ?nr a ?nr a1 ),其中J1(?n)?0),其中J1(?n)?0 ?D,n?J0(?E,n?J0( ?nr a2 ?nr —————————————————————————————————————— ------------------------------------------------------------------------------------------------ a ),其中J1(?n)?0 11?)?a12J02(?s)?a12s?0??An?22) a?a13?nJ0(?s)J1( 2 a12?n?a2?s2 (Cs??Cs?e ?jkC,sLc a1?n )Bs?(e jkB,sLb ?e ?jkB,sLb a2(J1(?s)Y0(?s)?J0(?s)Y1(?s))2?a12(J1(?s)Y0( 2Y12(?S) aa1?nJ1( a1?sa? )?J0(1s)Y1(?s))21? ?(a2?a12)s?0??An2(Ds??Ds?e ?jkD,sLd a1?na?a? )(J1(?s)Y0(1s)?J0(1s)Y1(?s))22 —————————————————————————————————————— ------------------------------------------------------------------------------------------------ (?n??s)Y1(?s) 3 a1?a2?nJ0(?s)J1( 22 a2?n?a12?s2 122??jkC,sLb? )?a2J0(?s)??(Cne?Cn)2 a2?n )a1 (Ds?e ?jkD,sLd 122? ?Ds?)?a2J0(?s)??En 2 3 a?a2?nJ0(?s)J1( 22 a2?n?a2?s2 a2?n ) ?kB,nBn?(e?k C,n —————————————————————————————————————— ------------------------------------------------------------------------------------------------ jkB,nLb aaa aa1?sJ1(1?s)[J1(?n)Y0(1?n)?Y1(?n)J0(1?n)]?jkL?eB,nb)?22 (?s??n)Y1(?n) ?jkC,nLb ??(Cn?Cne ) aa3?sJ0(?n)J1(1 2 a12?s2?a2?n a1?s ) 3a1a2?sJ0(?n)J1( 11 ?k?a2s?0?As?kA,s?a2J02(?S) 22kC,s(Cs?e ?jkC,sLc 1???jkD,nLd ?Cs?)?a12J02(?S)??kD,n(Dn?Dne) 2 ??Dn) 3 —————————————————————————————————————— ------------------------------------------------------------------------------------------------ aa2?sJ0(?n)J1( a2?s )a1 22a2?s?a12?n2 ?kD,n(Dn?e ?jkD,nLd a2 ?s) 122?kE?aJ0(?S)E,ss22 a2?s?a2?n22 mode4 with inlet NL/db frequency/hz 内嵌型(inlet)和外置型(outlet)波导管相比，内置型的Lb作用了两次，即起到了共振腔的作用，又起到了突变截面管的作用，所以，Lb引起的共振峰和外置型的相比更加尖锐。 后置型 ? PA(r,z1)?e?jkz1??Ane?PB(r,z2)??(Bne?PC(r,z1)??(Cne ?jkB,nz2 jkA,nz1 ?A,n jkB,n(z2?Lb)jkC,n(z1?LC) —————————————————————————————————————— ------------------------------------------------------------------------------------------------ ? ?Bne??Cne )?B,n)?C,n ?jkC,nz1 PD(r,z2)??Dne ?jkD,nz2 ?D,n ?P 由j??U??得: ?z1?jkA,nz1 UA(r,z1)?[ke?jkz1??kA,nAne?A,n] ?? 1 UB(r,z2)?UC(r,z1)?UD(r,z2)? ?? 1 ?kB,n(Bn?e ? k(C?C,nne?kD?D,nne ?jkB,nz2 ? ?Bne —————————————————————————————————————— ------------------------------------------------------------------------------------------------ jkB,n(z2?Lb) )?B,n ?jkC,nz1 ?? 1 ??Cne jkC,n(z1?Lc) )?C,n ?jkC,nz2 ?? ?D,n jkB,nLb ??? ?Bn?0,Bn?Bne ?jkB,nLb ? 根据B管的刚性边界条件，Bne a2 。 (Cs??Ce(Bs?e ??jkC,sLc s —————————————————————————————————————— ------------------------------------------------------------------------------------------------ ? )??C,s?C,s?C??r?C,sdr??An??A,n?C,s?C jkB,sLb ? )??B,s?B,s?B??Dn??D,n?B,s?B ?jkB,sLb?jkC,sLc ?Bs?e (Cs?e ?Cs?)??C,s?C,s?C??Dn??D,n?C,s?C ?Be ?n jkB,nLb ?k?k B,n (Be ??jkB,nLbn )??B,n?D,s?B??kC,n(Ce a2 ??jkC,nLcn ?C)??C,n?D,s?C ?n ?Ds?kD,s??D,s?D,s?D —————————————————————————————————————— ------------------------------------------------------------------------------------------------ C,n (C?Ce ?n??jkC,nLcn )??C,n?A,s?C?k?r?A,sdr?kA,sAs???A,s?A,s?A ?A,n?J0( ?B,n?J0( ?C,n?J0(?nra),其中J1(?n)?0)?J1(?n)?r?aJ(?)?aY0(n),其中 J1(n1)?1nY1(n1)?0Y1(?n)aaY1(?n)a?nra?nr a2 ),其中J1(?n)?0 ),其中J1(?n)?0?D,n?J0(?nr a a2?n?3a?a?J(?)J()2n0s1??12212??jkC,sLc?)?a2J0(?s)?a2s?0??An??(Cs?Cse222222a??a?s2n??a?a?a2(J1(?s)Y0(?s)?J0(?s)Y1(?s))2?a12(J1(?s)Y0(1s)?J0(1s)Y1(?s))2?jkL?jkL?Bs?(eB,sb?eB,sb)22Y1(?S)??a?a?a??aa1?nJ1(1n)(J1(?s)Y0(1s)?J0(1s)Y1(?s))??D??n22?(?n??s)Y1(?s)?a2?n?3a?a?J(?)J()2n0s1???jkC,sLC122??Cs)?a2J0(?s)??Dn?(Cse222a2?n?a2?s2? a1?sa1?na1?n?aa?J()[J(?)Y()?J()Y1(?n)]1s11n00?jkB,nLb?jkB,nLb??e)???kB,nBn(e22(?s??n)Y1(?n)??a2?saa3?J(?)J()?n12s0??jkC,nLb??Cn)???kC,n(Cne22a2?s?a2?n2 ??1??Ds?kD,s?a2J02(?S)2??a2?s3aa?J(?)J()?2s0n112122???jkC,nLC?)?ka|?kA?aJ0(?S)??kC,n(Cn?Cnes?0A,ss2222?a??a?222sn? —————————————————————————————————————— ------------------------------------------------------------------------------------------------ model 4 with outlet in the back NL/dbf/hz compara beteen model 4 with outlet in forehead and back0100020003000400050006000 与前置性内插管的消声规律相比较，该模型相当于把前者颠倒，输入作输出，输出作输入。二者的在低频消声峰值基本相同，但图像有很大差别。这是由于此波导具有非互易性的特点造成的。 波导内的声信号 相速度、群速度 频散现象 有关声学波导的声学计量和实验方法 四传声器法测量隔声模型的传输损失 3.3.1 四传声器法的原理 四传声器法的原理图如图13所示，信号发生器发出信号，经扬声器变为声波进入声源管后产生平面波。隔声构件置于声源管和接收管之间，构件的前方为入射波部分，构件的后方为透射波部分。图中pA,pB分别代表声源管中向右和向左传播的平面波，pC,pD分别代表接收管中向右和向左传播的平面波。用t1,t2 表示测试构件前后两个面的声压透射系数,用r1,r2 表示测试构件前后两个面的声压反射系数。X1,X2,X3,X4分别表示四个传声器的位置坐标。 图表 1 四传声器法原理图 由于样品表面的反射 ,入射部分形成驻波场 ,采用驻波分离方—————————————————————————————————————— ------------------------------------------------------------------------------------------------ 法 ,即用两个传声器把入射波与反射波分开. 在接收管部分也采用双 传声器把正向透射波与末端的反射波分开 ,便组成了四传声器测试 方法。此时如果忽略掉接收管中pD透过构件对pB的影响，忽略掉 pD在构件后表面发生反射时对pC的影响，即理想状态为pD为零(吸 声材料使pC不存在回波)。根据管内平面波传播公式, 4个传声器位置处的声压可分别表示为 ?p1?Ae?jk1X?Bejk1X??jkXjkX?p2?Ae2?Be2 ? (3-3-1) ?jk3Xjk3Xp?Ce?De?3 ?jk4X??Dejk4X?p4?Ce 进而可以求得 ?1A??2j???C?1 ?2j?p1ejkX2?p2ejkX1sin[k(X1?X2)] p3ejkX4?p4ejkX3 sin[k(X3?X4)] 根据声压透射系数的计算公式可以得到 pCsin[k(X1?X2)]p3ejk(X4?X3)?p4jk(X2?X3)|tp|?||???epAsin[k(X4?X3)]p1 ?p2ejk(X1?X2) 进而可以得到传输损失为 TL??20lg|tp| 在上述的分析中，我们认为吸声材料的反射系数为0，而实际上 吸声材料是存在声反射的。也就是说，接收管中向右传播的波不仅包 括入射波中透过材料的那部分声能 , 还包括末端反射再经材料反射 —————————————————————————————————————— ------------------------------------------------------------------------------------------------ 的二次反射波 , 将隔声管中向右传播的波作为透射波会引起计算偏 差. 此时我们可以采用改进的四传声器法来消除尾端影响吸声材料 效果不好的影响。 在考虑构件前后表面的发射和透射现象后，隔声管中pC可表示 为pA经过样品正向透射与pD经过样品反向反射的和.pB可表示为pA经过样品正向反射与 pD经过样品反向透射的和。由此，我们可以写出管内平面波关系矩阵 ?C??tr??A? ????12??? (3-3-2) ?B??r1t2??D? 其中t1为构件正想透射系数，是我们欲求的物理量。 通过(3-3-1)我们可以求得四个平面波的声压值为 ?1p1ejkX2?p2ejkX11A?， B??2jsin[k(X1?X2)]2j??jkX3jkX4?p4e1?C?1p3e， D??2jsin[k(X3?X4)]2j?p1e?jkX1?p2e?jkX2sin[k(X1?X2)] p3e?jkX3?p4e?jkX4 sin[k(X3?X4)] 这样通过测量4个传声器的声压值 ,求出4个平面波声压值,即可 带入式 (3-3-2)进行计算。在式3-3-2)中,被求矩阵有4 个未知数，而 进行1次测量只能确定2个方程 ,故需测量2次得到4个方程才能求解.2次测量分别用下标a和 b表示 ,如Ca表示第1次测得的隔声管 中右向传播的声压值。则t1为 —————————————————————————————————————— ------------------------------------------------------------------------------------------------ t1?CaDb?CbDa (3-2-3) AaDb?AbDa 如果方程(3-2-2)中分母在不同的数值变得比减数的绝对值更小， 那么这个方程的解是不稳定的。显然，两次测试的边界条件变化明显 有利于计算。一种有效的方法是改变隔声管末端的吸声系数，即一次 测试打开末端使隔声管中的声波近似于行波，二次测试为吸声系数较 高的吸声末端。这样就可以消除吸声材料性能不足带来的误差。 考虑到传声器的性能，不同传声器的幅值灵敏度和相位灵敏度可 能都会存在差异，并且传声器的相位灵敏度并不便于测量。因此，我 们寻找了一种利用数学处理消除不同传声器之间相位灵敏度差异的 方法。这种方法只需要传声器3对传声器1的幅值灵敏度比，但是需 要在传声器按照1，2，3，4号顺序排放测量完毕后，将1，2号传声 器和3，4号传声器交换位置做再次测量。此时我们可以写出第一次 测量时管道中的声压解 ?p1?(Iue?jkx1?Ruejkx1)ej?t ??jk(x1??x)?Ruejk(x1??x))ej?tJ2ej?2?p2?(Iue ? (3-2-4) ?jkx3jkx3j?3j?t?Rde)eJ3e?p3?(Ide ?jk(x3??x)??Rdejk(x3??x))ej?tJ4ej?4?p4?(Ide 其中I,R分别表示正向波和反向波的复数幅值。 将传声器1，2和传声器3，4位置交换后，我们得到第二次测量 时管道中的声压解 ?p1??(Iue?jk(x1??x)?Ruejk(x1??x))ej(?t??) ???(Iue?jkx1?Ruejkx1)ej(?t??)J2ej?2?p2 ? —————————————————————————————————————— ------------------------------------------------------------------------------------------------ (3-2-5) ?jk(x3??x)jk(x3??x)j?3j(?t??)??(Ide?Rde)eJ3e?p3 ???(Ide?jkx3?Rdejkx3)ej(?t??)J4ej?4?p4 根据(3-2-4)和(3-2-5)式，我们可以用两次测量得到的声压 ?,p2?,p3?,p4?来表示Iu,Ru,Id,Rd以及他们之间的关系 p1,p2,p3,p4,p1 ?RuH12?e?jk?x ?2jkx1H21??e?jk?xIH12?eH21?u ?jk?x ?R (3-2-6) ?HH?e4334?2jkx3?d??ejk?x?Id?HH?e4334? ?,H12,H43?,H34分别表示其中H21?p?p2p3p1,,,4。 ?p1p4?p3p2 将(3-2-6)式代入(3-2-4)式，可分别得到p1,p3为 ??H12?e?jk?xH21?jkx1(1?)?p1?Iuejk?xHH?e?2112 ? (3-2-7) ?jk?x?H34?eH43??jkx3(1?)J3jk?x?p3?Ide?HH?e4334? 此时，利用p3便可以得到 p1 H34?e?jk?xH12?e?jk?xH43H21Idjk(x3?x1)?e(1?)/[H31J3(1?)] jk?xjk?xIuH12?eH34?eH21H43 其中H31表示p1。根据传输损失的计算公式可以得到 p3 ?H34?e?jk?x?H12?e?jk?xH43H21Iu?20log[H31J3(1?)/(1?)] TL?20log(3-2-8) jk?xjk?xId?H34?e?H12?eH43H21 至此，我们一共提出了三种利用四传声器测量传输损失的方法， 在具体的实验过程中，我们会根据设备的实际情况来选取数据处理方 —————————————————————————————————————— ------------------------------------------------------------------------------------------------ 法。 3.3.2 实验模型传输损失测量 在实际的模型传输损失测量实验时，我们各加工了图1和图3中的模型。模型的具体参数为 a?7cm,a1?1.3cm,a2?2.1cm,Lc?7.7cm,Lb?3.3cm。在极轴对称声源情况下，模型的截止频率为2997Hz。按照之前的分析，此时的TL峰值恰好可以消除其谷点的存在,如图11。模型材料为有机玻璃，管A,D的厚度为1cm,管B的厚度是8mm。由于模型的加工存在一定的误差，实际得到的模型的a1?1.26cm,a2?2.12cm,管B的厚度为8.6mm，长度为3.2cm，管C长度为7.8cm。 这对理论计算结果的影响，如图 11. 我们首先对吸声材料的声压反射系数利用传递函数的方法进行了测试。实验原理图如下: 吸声材料 图表 2 rp测试原理图 设传声器1，2的位置坐标分别为x1,x2，当入射波的频率低于管道的截止频率时，我们认为管道内传播的是平面波。那么传声器1，2接收到的声压可以分 ?ikxikx??p1?pie1?pre1 别表示为 ? ?ik2xik2x??pre?p2?pie 其中pi,pr分别表示入射波和反射波在x?0处的幅值，并且存在pr?r?pi，r是末端吸声材料的声压反射系数。由此，我们可以得到 p2e?jkx2?re?jkx2H12???jkx1 p1e?re?jkx1 —————————————————————————————————————— ------------------------------------------------------------------------------------------------ 整理后可以求得 e?jk(x2?x1)?H122jkx1e r? (3-3-2-1) jk(x2?x1)H12?e 进而可以得到吸声棉的吸声系数 ??1-|r|2 如果考虑到传声器之间的灵敏度差异(包括相位和幅值)，我们可以利用交换传声器位置，做两次测量的方法加以纠正。此时，可以将传声器两次接收到的声压表示为 ?p1?pie?ikx1?preikx1 ??ikxikxj??p2?(pie2?pre2)J2e2 ??ikx2ikx2??pie?pre?p1?p??(pe?ikx1?peikx1)Jej?2ir2?2 ?,p2?分别表示传声器1，2前后两次测得的声压。此时我们取其中p1,p2,p1 ?H?H12H21?p2p1，用H替换(3-3-2-1)式中的H12，这样就可以消除传??p1p2 H1222j?2?J2e可以求得传声器2对传声器1的 H21声器之间的灵敏度差异。另外利用 在完成对吸声材料吸声效果的测定和传声器灵敏度的校准后，由图17、18可以看出，传声器灵敏度随频率的变化情况较大。我们选取了300Hz-2950Hz的频率段进行测量，并方法一、三对实验数据进行了处理。 图 19 隔声模型传输损失测量原理图 —————————————————————————————————————— ------------------------------------------------------------------------------------------------ 实验数据采集的过程中，我们选用了SKC的传声器和TopView2000采集系统，采样频率为100K，采样长度为50*1024点。我们将传声器编号为1、2、3、4，按照1、2、3、4的顺序将四个传声器从左向右摆放，并将采集到的数据记为:数据1234_1;然后再按照2、1、4、3的顺序摆放传声器，将 在此基础上，我们进行了第二组的测量，即将传声器1、2、3、4依序从左到右摆放，记测量数据为:数据1234_2;然后交换传声器顺序，记录数据为: 通过上图可以看出，实验数据在相对较低的频率段与理论计算结果符合较好，但在高频段存在一定的偏差。这主要是因为管道在极轴对称声源下的截止频率为2997Hz，在非对称声源下的截止频率为c0 ，即1439Hz。当入射声波2?a 频率大于管道的截止频率时，管道中便会存在能够正常传播的高次模式，这与四传声器法原理中管道内传播的是平面波的假设不符。在实际的实验中，如果扬声 器的摆放存在偏差，便会出现上述情况，使实验数据存在误差。另外，利用方法1处理得到的结果相对方法3处理得到的结果较好，这主要是因为方法3对噪声敏感。在噪声环境下，第一测量得到的声压应表示如下 ?p1?(Iue?jkx1?Ruejkx1?noise1)ej?t ??jk(x1??x)?Ruejk(x1??x)?noise1)ej?tJ2ej?2?p2?(Iue ??jkx3jkx3j?3j? tp?(Ie?Re?noise1)eJedd3?3 —————————————————————————————————————— ------------------------------------------------------------------------------------------------ ?jk(x3??x)??Rdejk(x3??x)?noise1)ej?tJ4ej?4?p4?(Ide 交换顺序后测量得到的声压与前一次测量得到的声压中噪声不 同 ?p1??(Iue?jk(x1??x)?Ruejk(x1??x)?noise2)ej(?t??) ???(Iue?jkx1?Ruejkx1?noise2)ej(?t??)J2ej?2?p2 ??jk(x3??x)jk(x3??x)j?3j(?t??)?p?(Ie?Re?noise2)eJedd3?3 ???(Ide?jkx3?Rdejkx3?noise2)ej(?t??)J4ej?4?p4 在这种情况下，求解得到的声压反射系数和声压透射系数均存在 较大误差。 ——————————————————————————————————————