序:
感觉三者在一定程度上有相似性,故将其整理至此文
瞬态声波方程与稳态声波方程
- 瞬态声波方程
概念:描述时域空间域(波动方程)
∇ 2 p = 1 c 0 2 ∂ 2 p ∂ t 2 \nabla^2 p=\frac{1}{c_0^2}\frac{\partial^2p}{\partial t^2} ∇2p=c021∂t2∂2p - 稳态声波方程
概念:描述频域空间域(亥姆霍兹方程)
∇ 2 P + k 2 P = 0 \nabla^2 P+k^2P=0 ∇2P+k2P=0 - 瞬态与稳态的关系
通过傅里叶变换联结(稳态声波方程假设声压单频简谐变化,声波方程化简为亥姆霍兹方程,而通过傅里叶变化可以将一个非周期的时域信号分解成不同频率的简谐信号的叠加)
p ( r ⃗ , t ) = 1 2 π ∫ − ∞ + ∞ P ( r ⃗ , ω ) e j ω t d ω = F − 1 [ P ( r ⃗ , ω ) ] p(\vec{r},t)=\frac{1}{2\pi}\int_{-\infty}^{+\infty} P(\vec{r},\omega)e^{j\omega t}{\rm d}\omega=\mathscr{F}^{-1}[P(\vec{r},\omega)] p(r,t)=2π1∫−∞+∞P(r,ω)ejωtdω=F−1[P(r,ω)]
角谱理论
概述:将一般的稳态声场表示为平面波的叠加
P
(
r
⃗
)
=
∬
A
(
k
⃗
)
e
j
k
⃗
⋅
r
⃗
d
k
⃗
P(\vec{r})=\iint A(\vec{k})e^{j\vec{k}\cdot \vec{r}} {\rm d} \vec{k}
P(r)=∬A(k)ejk⋅rdk
其中
k
⃗
\vec{k}
k为波矢,代表平面波传播的方向,且
∣
k
⃗
∣
=
ω
c
0
|\vec{k}|=\frac{\omega}{c_0}
∣k∣=c0ω,即所有频率相同但传播方向不同的平面波的叠加。(注:积分时只需要积两个变量即可,因为方向仅需两个变量即可表示,不需要径长)
P
(
r
⃗
)
=
∫
−
∞
+
∞
d
k
x
∫
−
∞
+
∞
d
k
y
A
(
k
⃗
)
e
j
(
k
x
x
+
k
y
y
+
k
z
z
)
P(\vec{r})=\int_{-\infty}^{+\infty}{\rm d} k_x \int_{-\infty}^{+\infty}{\rm d} k_y A(\vec{k})e^{j(k_xx+k_yy+k_zz)}
P(r)=∫−∞+∞dkx∫−∞+∞dkyA(k)ej(kxx+kyy+kzz)
分析:
物理意义与傅里叶变换类似,将空间域变换到波矢域,对传播方向不同但频率相同的平面波分别分析,然后再叠加。
模式理论
分析无限空间中声波的模式理论,设给定时刻的声压为
p
(
r
⃗
)
p(\vec{r})
p(r),稳态振动的复振幅为
P
(
r
⃗
)
P(\vec{r})
P(r),稳态波动方程的解即为各个模式,称其为振型(可理解为模式的振幅):
P
(
k
)
(
r
⃗
)
=
(
1
2
π
)
3
2
e
j
k
⃗
⋅
r
⃗
=
(
1
2
π
)
3
2
e
j
(
k
x
x
+
k
y
y
+
k
z
z
)
P^{(k)}(\vec{r})=(\frac{1}{2\pi})^{\frac{3}{2}}e^{j\vec{k}\cdot \vec{r}}=(\frac{1}{2\pi})^{\frac{3}{2}}e^{j(k_xx+k_yy+k_zz)}
P(k)(r)=(2π1)23ejk⋅r=(2π1)23ej(kxx+kyy+kzz)
其对应的共振频率是:
ω
(
k
)
=
c
0
∣
k
⃗
∣
\omega^{(k)}=c_0|\vec{k}|
ω(k)=c0∣k∣
1.任意的正实数都可以是共振频率,无限大空间的声场的谱是连续谱。(这里说的是波矢域吗?)
2.每个共振频率对应有无数个不同的波矢,他们的模相等,此时不同的振型代表不同传播方向的平面波。
故
k
⃗
\vec{k}
k模式即为:
P
(
k
)
(
r
⃗
)
e
−
j
ω
(
k
)
t
=
(
1
2
π
)
3
2
e
j
(
k
⃗
⋅
r
⃗
−
c
0
∣
k
⃗
∣
t
)
P^{(k)}(\vec{r})e^{-j\omega^{(k)}t}=(\frac{1}{2\pi})^{\frac{3}{2}}e^{j(\vec{k}\cdot \vec{r}-c_0|\vec{k}|t)}
P(k)(r)e−jω(k)t=(2π1)23ej(k⋅r−c0∣k∣t)
可以证明:不同振型满足正交归一性
三维傅里叶变换
在全空间满足一定条件的函数
p
(
r
⃗
)
p(\vec{r})
p(r)可以表示成三维傅里叶变换:
p
(
r
⃗
)
=
∭
a
(
k
⃗
)
(
1
2
π
)
3
2
e
j
k
⃗
⋅
r
⃗
d
k
⃗
=
∭
a
(
k
⃗
)
P
(
k
)
(
r
⃗
)
d
k
⃗
p(\vec{r})=\iiint a(\vec{k})(\frac{1}{2\pi})^{\frac{3}{2}}e^{j\vec{k}\cdot \vec{r}} {\rm d} \vec{k}=\iiint a(\vec{k})P^{(k)}(\vec{r}) {\rm d}\vec{k}
p(r)=∭a(k)(2π1)23ejk⋅rdk=∭a(k)P(k)(r)dk
物理意义:将空间域
r
⃗
\vec{r}
r和波矢域
k
⃗
\vec{k}
k建立联系,式中
a
(
k
⃗
)
a(\vec{k})
a(k)是各个振型的系数,不同的
a
(
k
⃗
)
a(\vec{k})
a(k)组成了一个线性的函数空间
记空间的三维傅里叶变换为
T
\mathcal{T}
T(实在找不到书上的那个符号。。。)
T
=
∭
d
k
⃗
(
1
2
π
)
3
2
e
j
k
⃗
⋅
r
⃗
=
∭
d
k
⃗
P
(
k
)
(
r
⃗
)
\mathcal{T}=\iiint {\rm d} \vec{k} (\frac{1}{2\pi})^{\frac{3}{2}}e^{j\vec{k}\cdot \vec{r}}=\iiint {\rm d}\vec{k} P^{(k)}(\vec{r})
T=∭dk(2π1)23ejk⋅r=∭dkP(k)(r)
T
−
1
=
∭
d
k
⃗
(
1
2
π
)
3
2
e
−
j
k
⃗
⋅
r
⃗
=
∭
d
k
⃗
[
P
(
k
)
(
r
⃗
)
]
∗
\mathcal{T}^{-1}=\iiint {\rm d} \vec{k} (\frac{1}{2\pi})^{\frac{3}{2}}e^{-j\vec{k}\cdot \vec{r}}=\iiint {\rm d}\vec{k} [P^{(k)}(\vec{r})]^{*}
T−1=∭dk(2π1)23e−jk⋅r=∭dk[P(k)(r)]∗
将瞬态声场
p
(
r
⃗
,
t
)
p(\vec{r},t)
p(r,t)和稳态声场
P
(
r
⃗
,
ω
)
P(\vec{r},\omega)
P(r,ω)用简正坐标表示:
p
(
r
⃗
,
t
)
=
T
[
a
(
k
⃗
,
t
)
]
p(\vec{r},t)=\mathcal{T}[a(\vec{k},t)]
p(r,t)=T[a(k,t)]
P
(
r
⃗
,
ω
)
=
T
[
A
(
k
⃗
,
ω
)
]
P(\vec{r},\omega)=\mathcal{T}[A(\vec{k},\omega)]
P(r,ω)=T[A(k,ω)]
p
(
r
⃗
,
t
)
=
F
−
1
[
T
[
A
(
k
⃗
,
ω
)
]
]
p(\vec{r},t)= \mathscr{F}^{-1}[\mathcal{T}[A(\vec{k},\omega)]]
p(r,t)=F−1[T[A(k,ω)]]
总结:物理坐标和简正坐标都能确定声场,瞬态和稳态的物理坐标
p
(
r
⃗
,
t
)
p(\vec{r},t)
p(r,t)和
P
(
r
⃗
,
ω
)
P(\vec{r},\omega)
P(r,ω)分别在时间空间域和频率空间域描述声场。瞬态和稳态的简正坐标
a
(
k
⃗
,
t
)
a(\vec{k},t)
a(k,t)和
A
(
k
⃗
,
ω
)
A(\vec{k},\omega)
A(k,ω)分别在时间波矢域和频率波矢域描述声场。
注意:
1.在做稳态声场的三维傅里叶变换时,由于简正坐标
A
(
k
⃗
,
ω
)
A(\vec{k},\omega)
A(k,ω)的波矢和频率满足关系
ω
=
c
0
∣
k
⃗
∣
\omega=c_0|\vec{k}|
ω=c0∣k∣,故三重积分退化成二重积分,此时的三维傅里叶变换就是角谱理论。
2.角谱理论的应用(此时声场形状不规则,难以通过坐标分离求得声场解析式):
- 将声场(待求)通过角谱理论分解成一系列传播方向不同的平面波,从空间域 r ⃗ \vec{r} r变换到波矢域 k ⃗ \vec{k} k
- 将边界条件也通过角谱理论分解,从空间域 r ⃗ \vec{r} r变换到波矢域 k ⃗ \vec{k} k
- 此时,每一束平面波(每一个模式 k ⃗ \vec{k} k)都有其在波矢域对应的边界条件,代入后即可求得每一束平面波的表达式
- 利用角谱理论叠加即可得到声场的表达式
Reference:
张海澜.理论声学[M].北京:高等教育出版社,2012.
徐佳新.二维阵列换能器声辐射力分布的计算分析[J].声学技术.2018(4).398-403