Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
акустика / xie_bosun_spatial_sound_principles_and_applications.pdf
Скачиваний:
158
Добавлен:
04.05.2023
Размер:
28.62 Mб
Скачать

Spatial sound reproduction by wave field synthesis  463

of a secondary source in Section 10.1.4, the spatial window in Equation (10.2.26) causes an edge effect, but spatial windows with a smooth transition reduce the edge effect.

For a target incident plane wave, the pressure at an arbitrary receiver position is expressed in Equation (1.2.6). Substituting Equation (1.2.6) into Equation (10.2.27) yields

Epl3D r , f 2w r PA f j k n exp jk r

(10.2.28)

2w r PA f jkcos sn exp jk r ,

 

where the subscript “pl” denotes the target plane wave, and θsn = θS is the angle between the direction of the incident plane wave and the outward-normal direction at r′ of the boundary surface.

For a target point (monopole) source at rS, the pressure at an arbitrary receiver position is presented in Equation (1.2.3). Substituting Equation (1.2.3) into Equation (10.2.27) yields

3D

 

r , f

 

2w

 

r

 

p

f

 

1 jk|r rS |

r rS n

 

jk

 

S

p

2

 

 

E

 

 

 

 

Q

 

4

 

r rS

 

 

 

|r rS |

exp

 

r r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(10.2.29)

 

 

 

 

 

 

 

 

 

 

 

1 jk|r rS |

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2w r Qp f

4

 

r rS

 

2

cos sn exp jk | r rS | ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Where subscript “p” denotes the target point source, and θsn is the angle between the target source direction and the outward-normal direction at r′ of the boundary surface.

10.2.3  General theory of two-dimensional WFS

Secondary sources are infinite straight-line sources perpendicular to the horizontal plane. They are arranged uniformly and continuously to form a vertical cylindrical surface. LΣ is the intersecting curve line between the cylindrical surface and the horizontal plane z = 0. In the subsequent discussion, this case is simply termed “secondary straight-line sources arranged in a horizontal curve LΣ″. Thus, the target source is outside the horizontal region SΣ closed by LΣ and the receiver position is in SΣ.

For an ideal acoustical holography, the reconstructed pressure in the receiver region is calculated with the two-dimensional Kirchhoff–Helmholtz integral in Equation (10.2.16). In this case, the reconstructed pressure distribution is independent from the vertical (z) coordinate and thus identical in each plane parallel to the horizontal plane. Equation (10.2.16) also indicates that arrays of secondary monopole and dipole straight-line sources are generally needed for an ideal acoustical holography. If an array of secondary monopole straight-line sources only is used, the reconstructed pressure at a horizontal receiver position r is calculated as

P

 

 

 

 

 

free

r, r , f

 

 

 

r , f

 

 

r, f

 

 

 

G2D

 

E2D

 

dL ,

 

 

 

 

 

L

 

 

 

 

 

 

 

where

E2D r , f 2w r P r , f .

n

(10.2.30)

(10.2.31)

464  Spatial Sound

Equation (10.2.31) describes the driving signals of secondary monopole straight-line source. The spatial window similar to that in Equation (10.2.26) is chosen, but in this case, n′ is the unit vector at the inward-normal direction of LΣ. Equations (10.2.30) and (10.2.31) indicate that an array of secondary monopole straight-line sources can reconstruct the target sound field approximately in a horizontal receiver region, or more exactly, can reconstruct the vertical direction-independent sound field in a vertical cylindrical region inside the array of secondary sources, by appropriately choosing the driving signals.

For a target incident plane wave, the driving signals of secondary straight-line sources are identical to those expressed in Equation (10.2.28)

Epl2D r , f 2w r PA f j k n exp jk r

(10.2.32)

2w r PA f jkcos sn exp jk r .

 

For a target straight-line source, the radiated pressure in Equation (1.2.7) is described as

P r , rS , f Qli f

j

H0 k|r rS | .

(10.2.33)

 

4

 

 

Driving signals are derived by substituting Equation (10.2.33) into Equation (10.2.31) and using the derivative formula dH0(ξ)/dξ = −H1(ξ) of the zero-order Hankel function of the second kind as

Eli2D r ,f w r Qli f

jk

 

r rS n

H1 k|r rS |

2 |r rS |

 

(10.2.34)

w r Qli f jk2 cos snH1 k|r rS | ,

where H1(k |r′−rS|) is the first-order Hankel function of the second kind, and θsn is the angle between the target source direction and the outward-normal direction at r′ of the boundary surface.

According to the asymptotic formula of the Hankel function of the second kind, when k |r r′| >> 1, the following equation is obtained:

 

k|r r |

2

 

 

 

jk|r r | j

 

 

H0

 

 

 

exp

 

.

(10.2.35)

k |r

r |

4

 

 

 

 

 

 

As in Section 10.1.3, at the far-field, the magnitude of cylindrical wave created by a straight-line secondary source is inversely proportional to the square root of the distance between the secondary source and receiver position (−3 dB law rather than −6 dB law for a point source). In practice, WFS can be conveniently implemented using secondary point sources. Substituting Equation (9.2.35) into (9.2.5) and using Equation (10.1.2) yield

Gfree2D r, r , f

2 |r r | 1

 

exp jk|r r |

 

2 |r r |

Gfree3D r, r , f .(10.2.36)

jk

4

 

|r r |

jk

 

 

 

 

Spatial sound reproduction by wave field synthesis  465

The reconstructed pressure at r is approximated by substituting Equation (10.2.36) into Equation (10.2.30) as follows:

P

 

r, f

 

 

 

2 |r r |

E2D

 

r , f

 

G3D

r, r , f

 

dL .

(10.2.37)

 

 

 

jk

 

 

free

 

 

 

 

 

 

 

L

 

 

 

 

 

 

 

 

 

 

In Equation (10.2.37), when the target source and receiver position are restricted in the horizontal plane, a sound field can be reconstructed by a uniform and continuous array of secondary monopole point sources arranged on LΣ instead of a straight-line source. In this case, driving signals should be equalized; that is,

P

 

r, f

 

 

 

G3D

r, r , f

 

E2.5D

 

r , f

 

dL .

(10.2.38)

 

 

 

free

 

 

 

 

L

The equalized driving signals for secondary monopole point sources are related to those for secondary monopole straight-line sources by the following equation:

E2.5D r , f

2 |r r |

E2D r , f .

(10.2.39)

 

jk

 

 

Equation (10.2.38) specifies a 2.5-dimensional WFS or reproduction. Equation (10.2.39) describes 2.5-dimensional driving signals, which are obtained through frequencyand dis- tance-dependent equalization of two-dimensional driving signals. Equalization depends on the receiver position. In 2.5-dimensional reproduction, the reconstructed sound field deviates from the target sound field when either the source or receiver position deviates from the horizontal plane.

The driving signals in Equation (10.2.39) depend on the receiver position. As indicated in Section 10.1.3, the magnitude of driving signals is equalized at a given reference position rref in practice. Then, 2.5-dimensional driving signals in Equation (10.2.39) become

E2.5D r , f

2 |rref r

 

| E2D r , f .

(10.2.40)

 

jk

 

 

 

 

 

For a target plane wave, substituting Equation (10.2.32) into Equation (10.2.40) yields

2.5D

 

 

 

jk

2 |rref r

 

 

 

Epl

r , f 2w r

 

| PA f cos sn exp jk r

 

 

 

 

 

jk

 

 

 

 

 

(10.2.41)

 

 

4 w r

 

 

|rref r

 

|cos snP r , f ,

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

where P(r′, f) = PA(f) exp (−jk r′) is the pressure of the target plane wave at the position of secondary sources. The physical significance of Equation (10.2.41) is similar to that of Equations (10.1.20) and (10.1.21) except for the spatial window w(r′) of secondary sources arranged in LΣ.

466  Spatial Sound

Spors et al. (2008) suggested deriving 2.5-dimensional driving signals of secondary sources from the pressure of a three-dimensional target point source rather than a target straight-line source for a target source at the horizontal position rS because of the uneven frequencyspectral characteristics of the radiation of a target straight source in Equation (10.2.33). Substituting E2D(r, f) in Equation (10.2.40) with Ep3D r , f in Equation (10.2.29) yields

E2.5D

 

r , f

 

2w

 

r

 

Q

f

 

 

2 |rref

r |1 jk|r rS | r rS n

exp

jk

 

r

r

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

jk

 

4

r rS

2

 

|r rS |

 

 

 

 

 

 

 

S

(10.2.42)

For k |r′- rS| >>1, Equation (10.2.42) becomes

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2.5D

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

exp jk

 

r

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

rS n

 

 

 

 

S

 

 

Ep

 

r , f

2w r

 

Qp f

jk

2 |rref r

 

|

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

|r rS |

 

 

 

4 |r rS |

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2w r

jk

2 |rref r |cos snP r , rS, f

 

 

 

 

 

 

 

 

 

(10.2.43)

 

 

 

 

 

 

 

w r 4

 

 

jk

 

|rref r |cos snP r , rS, f ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

where θsn is the angle of the target source with respect to the outward-normal direction at rof the secondary source array. For a horizontal linear array of secondary monopole point sources in Section 10.1.3, w(r) = 1. Equation (10.2.43) is equivalent to Equation (10.1.21) when | rrS | >> 1. Therefore, a horizontal linear array of secondary monopole point sources is a special case of the discussion in this section.

If driving signals are directly derived from the pressure of a target straight-line source, the two-dimensional driving signals in Equation (10.2.34) should be equalized or multiplied with the factor jk so that the radiation of the target straight-line source has a frequency-spectral characteristic identical to that of a point source. Then, the equalized two-dimensional driving signals are converted into 2.5-dimensional driving signals by using Equation (10.2.40):

2.5D

 

 

 

 

 

 

jk r rS n

 

 

 

Eli

r , f w r

Qli f

2 |rref r

 

|

 

 

 

H1 k|r

 

rS | . (10.2.44)

 

2

 

|r rS |

 

 

 

 

 

 

 

 

 

 

 

 

When k |r-rS| >> 1, the asymptotic formula of the first-order Hankel function of the second kind yields

 

k|r rS|

2

 

 

jk|r rS | j

 

j

 

 

H1

 

 

exp

 

 

.

(10.2.45)

k |r rS |

2

4

 

 

 

 

 

 

 

In the case of Qli(f) = Qp(f), Equation (10.2.44) becomes

Eli2.5 r , f 2 |r rS |Ep2.5D

 

 

 

 

 

 

w r Qli f

jk

 

2 |rref r |

r rS n exp jk|r rS|

.

(10.2.46)

2

 

|r rS |

 

|r rS |