大家好,欢迎来到IT知识分享网。
概述
对于接受阵的某个阵元来说,其空间位置固定,接收到的是声场时域信号。
在理想波导中,我们假设声源是单频的,推导出了声场的表达式,是由可以传播的各号简正波叠加而成:
p ( z , r ; ω ) = ∑ n = 1 N 2 π k r n r φ n ( k z n z 0 ) φ n ( k z n z ) e − j k r n r p(z,r;\omega)=\sum_{n=1}^{N}\sqrt{\frac{2\pi}{k_{rn}r}}\varphi_{n}(k_{zn}z_0)\varphi_{n}(k_{zn}z)e^{-jk_{rn}r} p(z,r;ω)=n=1∑Nkrnr2πφn(kznz0)φn(kznz)e−jkrnr
其中 φ n \varphi_{n} φn是垂直方向上的模态函数,影响该垂直位置的柱面波的振幅。各号模态函数的分布如下图所示,n号简正波的模态函数有n-1个零点:
单频信号
我们取某一阵元的深度 z h z_h zh,考虑第一号简正波和第二号简正波的区别,首先写出它们的表达式(依旧考虑单频声源):
p 1 ( r ; t ) = B 1 φ 1 ( k z 1 z 0 ) φ 1 ( k z 1 z h ) e j ( ω t − k r 1 r ) p_1(r;t)=B_1\varphi_{1}(k_{z1}z_0)\varphi_{1}(k_{z1}z_h)e^{j(\omega t-k_{r1}r)} p1(r;t)=B1φ1(kz1z0)φ1(kz1zh)ej(ωt−kr1r) p 2 ( r ; t ) = B 2 φ 2 ( k z 2 z 0 ) φ 2 ( k z 2 z h ) e j ( ω t − k r 2 r ) p_2(r;t)=B_2\varphi_{2}(k_{z2}z_0)\varphi_{2}(k_{z2}z_h)e^{j(\omega t-k_{r2}r)} p2(r;t)=B2φ2(kz2z0)φ2(kz2zh)ej(ωt−kr2r)
可以发现,如果固定水平位置 r r r,即在一个阵元上接收到的时域信号的频率是相同的,不同号简正波的差别在振幅和初相位;而固定时间 t t t,即在某一时刻,不同号简正波在空间上的重复不一样(也就是水平波数),当然振幅和初相位也不同。
这里就可以引入相速度的概念,即同一相位的传播速度。对于固定深度的某一号简正波,其相位可以表示为 ω t − k r r \omega t-k_{r}r ωt−krr,故其相速度为:
c p n = w k r = c 0 1 − ( k z n k 0 ) 2 c_{pn}=\frac{w}{k_r}=\frac{c_0}{\sqrt{1-(\frac{k_{zn}}{k_0})^2}} cpn=krw=1−(k0kzn)2c0
对应的也就是该深度处该号简正波的传播速度。所以如果是单频声源,那么不同号简正波可以通过不同的传播速度区分开来。
宽带信号
如果声源的频谱是存在一定的带宽的,根据傅里叶变换对,信号可以写成多个单频的简谐波叠加的形式,对于简谐波依旧可以用简正波理论去进行研究。
p ( t ) = 1 2 π ∫ ω 0 ω 0 + △ ω p ( ω ) e j ω t d ω p(t)=\frac{1}{2\pi}\int_{\omega_0}^{\omega_0+\triangle \omega}p(\omega)e^{j\omega t}{\rm d \omega} p(t)=2π1∫ω0ω0+△ωp(ω)ejωtdω
所以声场是由不同号简正波叠加形成,而每一号简正波又是通过不同频率的简谐柱面波叠加而成,这里选取某一号简正波来进行研究,我们知道它是由不同频率的简谐柱面波叠加得到的,所以可以进行傅里叶分解。
p n ( t , r ) = 1 2 π ∫ ω 0 ω 0 + △ ω S ( ω ) p n ( ω ) e j ω t d ω p_n(t,r)=\frac{1}{2\pi}\int_{\omega_0}^{\omega_0+\triangle \omega}S(\omega)p_n(\omega)e^{j\omega t}{\rm d \omega} pn(t,r)=2π1∫ω0ω0+△ωS(ω)pn(ω)ejωtdω p n ( ω ) = 2 π k r n r φ n ( k z n z 0 ) φ n ( k z n z h ) e − j k r n r p_n(\omega)=\sqrt{\frac{2\pi}{k_{rn}r}}\varphi_{n}(k_{zn}z_0)\varphi_{n}(k_{zn}z_h)e^{-jk_{rn}r} pn(ω)=krnr2πφn(kznz0)φn(kznzh)e−jkrnr
因为是同一号简正波,所以 k z n k_{zn} kzn相同, k r n k_{rn} krn随频率而变,而 S ( ω ) S(\omega) S(ω)则是频谱密度,仅与频率有关,与简正波号数无关。不同频率的简谐波叠加会有拍现象,形成波包:
对叠加信号进行分解,将 ω ( k r ) \omega(k_r) ω(kr)在 k r 0 k_{r0} kr0处泰勒展开 ω ( k r ) = ω 0 + d ω d k r 0 ( k r − k r 0 ) \omega(k_r)=\omega_0+\frac{d \omega}{d k_{r0}}(k_r-k_{r0}) ω(kr)=ω0+dkr0dω(kr−kr0):
f ( t , r ) = ∫ c ( ω ) e j ( ω t − k r r ) d ω f(t,r)=\int c(\omega)e^{j(\omega t-k_r r)}{\rm d \omega} f(t,r)=∫c(ω)ej(ωt−krr)dω f ( t , r ) = ∫ c ( ω ) e j ( ω 0 + d ω d k r ( k r − k r 0 ) ) t e − j ( k r 0 + k r − k r 0 ) r d ω f(t,r)=\int c(\omega)e^{j(\omega_0+\frac{d \omega}{d k_r}(k_r-k_{r0})) t} e^{-j(k_{r0}+k_r-k_{r0} )r}{\rm d \omega} f(t,r)=∫c(ω)ej(ω0+dkrdω(kr−kr0))te−j(kr0+kr−kr0)rdω f ( t , r ) = c ( ω 0 ) e j ( ω 0 t − k r 0 r ) ∫ e j ( k r 0 − k r ) ( d ω d k r t − r ) d ω f(t,r)=c(\omega_0) e^{j(\omega_0 t-k_{r0}r)}\int e^{j(k_{r0}-k_r )(\frac{d \omega}{d k_r}t-r)}{\rm d \omega} f(t,r)=c(ω0)ej(ω0t−kr0r)∫ej(kr0−kr)(dkrdωt−r)dω d ω d k r ∣ k r = k r 0 = c o n s t \frac{d \omega}{d k_r}|_{k_r=k_{r0}}=const dkrdω∣kr=kr0=const
个人理解,上式是信号的频带在较窄的情况进行线性展开的近似,前一项是调制,后一项是波包。是对这一号简正波的近似。
由此引出群速度的定义:
c g n = d w d k r = c 0 1 − ( k z n k 0 ) 2 c_{gn}=\frac{dw}{dk_r}=c_0 \sqrt{1-(\frac{k_{zn}}{k_0})^2} cgn=dkrdw=c01−(k0kzn)2
物理意义是该号简正波能量的传播速度,也可以理解为振幅的传播速度or波包的传播速度。相速度大于群速度,如果相位传过去了,但振幅没有传过去,能量还是没有传过去。
模间频散与模内频散
模间频散:不同号简正波的群速度不同,号数低的跑得快
模内频散:同一号简正波,频率高的跑得快(频率高的群速度大)
个人理解:
1.在推导群速度时用到了小区间内的泰勒展开,此时的角频率范围应该在 [ ω 0 − △ ω , ω 0 + △ ω ] [\omega_0-\triangle \omega,\omega_0+\triangle \omega] [ω0−△ω,ω0+△ω],其中 △ ω ≪ ω 0 \triangle \omega \ll \omega_0 △ω≪ω0,而200-300Hz的频率区间显然过大,故将其分成很多个小区间,在小区间中进行泰勒展开,形成波包+调制的一个信号,整体的时域信号就是这些频率小区间上信号的叠加,所以会有模间频散的现象,即高频率的信号跑得快。
2.而又因为号数高的简正波在同一频率区间中群速度变化率大,所以模间频散的现象更加明显。
PS:做了一下仿真,频率区间范围越大,模间频散越厉害,那么上述理解应该对的。
免责声明:本站所有文章内容,图片,视频等均是来源于用户投稿和互联网及文摘转载整编而成,不代表本站观点,不承担相关法律责任。其著作权各归其原作者或其出版社所有。如发现本站有涉嫌抄袭侵权/违法违规的内容,侵犯到您的权益,请在线联系站长,一经查实,本站将立刻删除。 本文来自网络,若有侵权,请联系删除,如若转载,请注明出处:https://yundeesoft.com/27386.html