Sound field separation technique using the principle of double layer patch acoustic radiation modes

Abstract. In order to solve the problems of near-field acoustic holography in applications such as external interference and aperture effects, a sound field separation technique using the principle of double layer patch acoustic radiation modes is proposed in this paper. The radiated acoustic pressures over two planar surfaces at certain distances from the sources are calculated first. Then, the effects resulting from the backscattering interference in non-free sound fields can be eliminated by a double-layer sound field separation technique. Next, data interpolation and extrapolation are performed on the separated data to increase the sound source's pressures on the holographic plane equivalently for holographic images with higher spatial resolution. Simulation and experimental results demonstrate that good agreements can be obtained with few measuring points.


Introduction
Near-field acoustic holography (NAH) [1] as a new technique in source identification has been used to reconstruct and visualize the complete radiated sound field.Specifically, ideal NAH requires an anechoic room, and the sound sources should be located on the same side of the measurement plane.However, most testing sites in practical applications cannot meet these requirements, and thus sound field separation techniques are introduced for eliminating the effects of interference sources [2].Additionally, when NAH technique is applied, the holographic aperture should be greater than the actual size of the sound source in order to suppress the aperture effect.As complicated structures over large areas are analyzed, the number of measuring points should be increased, which significantly increases the workload and time requirements.
For solving the aforementioned problems, many advances have been made, such as the sound field separation method in noisy environments.The spatial Fourier transform method, which was used to separate incident and scattered fields has been proposed by Hallman and Bolton [3].Statistically optimized NAH for sound field reconstruction in non-free-field environments with double layer pressure measurements [4], double layer velocity transducer array [5], and an array of pressure-velocity probes [6] can also be used.Further, the boundary element method (BEM) based on double-layer and single-plane measurements [7,8] and the equivalent source method (ESM) using double-layer pressure measurements and particle velocity measurements [9][10][11] have been proposed for the same purpose.Recently, taking into account the scattering effect of disturbing sources, sound field separation methods based on single-plane and double-layer pressure measurements by using spherical wave superposition method were proposed [12,13].Based on the plane wave approximation and sound pressure superposition principles, a sound-field separation method for double arrays was also developed [14].
Aiming at solving the problems induced by aperture effects, data interpolation and extrapolation methods based on Patch NAH have been proposed.A data interpolation method based on the reconstruction of band-limited signals and data extrapolated by constructing the weighted normalization in the frequency domain was used by Xu et al. [15,16].A local NAH technique based on the equivalent source method and a distributed source-boundary point method were developed by Bi et al. [17,18].Moreover, a method based on wave-superposition-based method was proposed for enhancing the spatial resolution of NAH images [19,20] and a NAH resolution enhancement method by using the interpolation of orthogonal spherical waves were studied [21,22].Methods for reducing the measuring points of holographic planes based on support vector regression (SVR) were also proposed by Mao, Du, et al. [23,24].
However, these methods only focus on a single issue and lack comprehensive considerations to address multiple issues.When two issues exist simultaneously, we should combine the sound field separation technique with the Patch method based on data interpolation and extrapolation, and perform the optimization selection and control on the related parameters uniformly.Therefore, the disturbance of the interference source can be accurately eliminated by using fewer measurement data points, and the spatial resolution of the holographic plane can be further enhanced.Meanwhile, the workload in measurements can be significantly reduced.The combination of these two techniques will be of great theoretical value and significance when promoting the engineering application of NAH.
In the 1990s, Elliott, Borgiotti and Cunefare proposed an acoustic radiation mode analysis method that can serve as another approach for solving the acoustic radiation problems in low-to-medium frequency vibration structures [25][26][27].Recently, Nie et al. [28] introduced a sound-field reconstruction formula for complex structures into the field of NAH based on acoustic radiation modes, and many favourable results were acquired in sound reconstruction and spatial resolution enhancement.In this paper, using sound radiation modes analysis and the reconstruction formula for sound fields, we consider the relationship between the relevant parameters and the proposed problems in a comprehensive account and propose a method to accurately eliminate the effects of interference sources, while simultaneously enhancing the spatial resolution of the holographic plane when fewer measuring points are used.Using the proposed method, the sound field distribution on the holographic plane with high spatial resolution can be obtained, which can provide accurate holographic data for further reconstruction of the sound field and have favourable engineering applicability.

Double-layer sound field separation technique
A vibration structure is placed in a homogeneous fluid whose sound velocity and density are and , respectively.The angular frequency of vibrations is .By defining a structural surface , the whole space can be divided into inner space and outer space, respectively, denoted as and .The variable denotes the outer normal direction on the surface of the structure.In the outer space of the vibration structure, the pressure generated by the vibration follows the Helmholtz equation and satisfies the Neumann boundary conditions on the boundary surface.Moreover, the pressure satisfies the Sommerfeld radiation condition at infinite distance.From the Helmholtz integral equations for the inner space and outer space, the following relation is obtained: where ′ and represent the position vectors on the measurement and vibration structure surface, respectively, and denotes the free Green's function.Based on acoustic radiation mode theory, the pressure matrix on the measurement plane can be obtained by: where = ( ′ ) ( ′ ) ⋯ ( ′ ) ; the elements in the matrix are denoted as = − ( | ′ ), denotes the element area after equal area partitioning on the structure surface, denotes the column vector of the expansion coefficients of different acoustic radiation modes, and = denotes the mode matrix of the sound field distribution.If the pressure is measured in two parallel planes, as displayed in Fig. 1, it can be expressed as in Eq. ( 2) in each plane.Considering that the pressure is a scalar quantity, we can obtain the following expression: where the pressure responses on the holographic plane generated by the target source and the interference source are denoted as and , respectively and the pressure responses on the holographic plane generated by the target source and the interference source are denoted as and , respectively.is the corresponding matrix of the sound distribution mode on the holographic plane generated by the source and is a the vector with expansion coefficients corresponding to the sound radiation modes of the target source and the interference source, and .
By combining Eq. (3) and Eq. ( 4), we can obtain: where '+' denotes the pseudo-inverse operation.When the condition number of is comparatively large, regularization should be performed during the inversion process.Considering that the high-order radiation modes represent evanescent waves, the optimal order of the modes could be selected during the solution of the mode coefficient, so that the time consumption and workload in the calculation can be effectively reduced.The selection principle for the optimal order of the modes will be described in detail in the following sections.
By substituting Eq. ( 5) and Eq. ( 6) into Eq.(3) and Eq. ( 4), the pressure responses on the holographic plane generated by the target source and interference source can be expressed as:

Patch technique based on data interpolation and extrapolation
As described in Eq. ( 2), the pressure generated by the sound source on the holographic plane  Based on the position information of all spatial points on the holographic plane , we can obtain the corresponding field distribution mode: where denotes the transfer function matrix between the holographic plane and the target source plane, denotes the sound radiation modes on the target source surface.To obtain a cut-off of with the same mode's order and the expansion coefficient vector obtained in Eq. ( 5), the pressure matrix on the holographic plane containing more data can be estimated as: In order to retain the true value as much as possible and further reduce the estimation error, the pressure values in at the points whose spatial positions are identical to the measuring points are substituted by the values in , and finally, the holographic data for reconstructing the NAH sound field can be obtained.
To verify the sensibility of the errors of the proposed method, the error index, Error, can be defined as: where ‖•‖ denotes the matrix 2-norm.denotes the pressure response generated by the target source after the two processing steps, and denotes the theoretical free field radiation radiated by the target source.

Selection of the optimal cut-off order of modes
Assuming that ( × ) and ( × ) denote the measurement data at measuring points on the measurement planes and , respectively, ( × ) denotes the matrix of the sound field distribution modes corresponding to each source, and ( × ) denotes the expansion coefficient vectors of sound radiation modes.Due to the matrix inversions are always ill-posed in the GUO LIANG, ZHU HAICHAO, MAO RONGFU, SU JUNBO, SU CHANGWEI separation processes, some small measurement errors may be magnified to result a failure by the higher-order modes which are sensitive to errors.Therefore, there is no need to use all field distribution modes when performing sound field separation.In order to acquire a minimum separation error, the order of the sound field distribution modes, ( × ) , should be cut off, and only the first -order modes ( × ) are needed for further calculation.Based on the Helmholtz equation least squares method [29], the following steps are conducted in order to determine the optimal cut-off order.First, the sound field distribution modes are cut off, with a cut-off order equal to 1, i.e., = 1.Then, the obtained cut-off field-distribution modes denoted as ( × ) are substituted into Eq.( 5) and Eq. ( 6) to solve the expansion coefficients of the first order field distribution modes, ( × ) and ( × ) .Next, the corresponding reconstructed pressure on the measurement plane is reconstructed in Eq. ( 3) to calculate the relative error between the reconstructed value and the measurement data by using the following ‖ ‖ -norm method: Subsequently, the cut-off order of modes are increased one by one, i.e., = + 1.By repeating the steps described above, the expansion coefficients with + 1 order ( × ) and ( × ) are calculated again.Moreover, the relative error between the results on the reconstructed plane and the measurement plane can be calculated.
Finally, through the traversal on the cut-off order of modes , from 1 to , a set of relative errors which are correlated with could be obtained.The cut-off order corresponding to the minimum relative error was thus regarded as the optimal cut-off order .It should be noted that the calculated may be different given different numbers of measuring points at different frequencies.Generally, the higher the frequency value is, the larger the cut-off order is.

Construction of the model
To verify the effectiveness and accuracy of the proposed method, two rigid pulsating spheres as the target source and interference source have been adopted for numerical simulations.The centre of the rigid pulsating sphere as the target source was located at the origin of coordinates, whose radius and vibration velocity were 0.1 m and 0.08 m/s, respectively.The radius and vibration velocity of another rigid pulsating sphere located at (0, 0, 0.85 m) were 0.15 m and 0.15 m/s, respectively.The two measurement planes, and , have identical sizes (0.5 m×0.5 m), on which 64 measuring points (8×8) were distributed, centres of which were (0, 0, 0.2 m) and (0, 0, 0.35 m), respectively.In numerical calculations, the sound was transmitted through the air, and its velocity was set to 343 m/s.The theoretical values of the measurement planes were calculated according to the radiation induced by a rigid pulsating sphere [30].In order to obtain simulation results closer to the actual measurements, a random white noise with an intensity of 30 dB was added to the simulated an actual measurement.
When processing the data collected on measurement planes using the proposed method, we first calculated the sound-field distribution modes of two rigid pulsating spheres.As stated above, the sound radiation modes are different at different frequencies.In the present work, was set to 0.5, which means the frequency was set to 272.95 Hz.Using the target source as an example, the first nine-order sound radiation modes and the first nine field distribution modes on the measurement plane were calculated, with the results displayed in Fig. 3 and Fig. 4, respectively.
As described above, = 0.5.Subsequently, the data on the measurement plane were GUO LIANG, ZHU HAICHAO, MAO RONGFU, SU JUNBO, SU CHANGWEI calculated.Fig. 5(a) displays the amplitudes of pressure at the 64 measuring points (8×8).Then the pressure response generated by the target source was separated by using the proposed double-layer field separation technique, as displayed in Fig. 5(b).For comparison, the theoretical values of size 8×8 corresponding to the holographic data were calculated, as displayed in Fig. 5(c).Two-fold data interpolation and extrapolation were then conducted on to obtain the holographic data of size 17×17, in which the pressure values at the points whose spatial positions identical to the measuring point were substituted by the data , with the amplitudes displayed in Fig. 5(d).For further comparison, theoretical data of size 17×17 without interference was calculated, as displayed in Fig. 5(e).

Analysis of the results
The comparative analysis between Fig. 5 source exists, a larger deviation occurs in high-pressure regions on the measurement plane, and dislocation appears in the regions with high and low pressures.Moreover, several regions with pseudo pressure appear, leading to great losses of image details.According to Eq. ( 11), the calculated measurement error compared with the theoretical value is approximately 52.23 %.
After the two proposed steps, the results in the high-pressure regions returns to normal, and the pseudo pressure regions can be eliminated.The separation error with respect to the theoretical values for 8×8 measuring points is also calculated in Eq. ( 11) and reduced to only 0.86 % which means that the disturbance of the interference source can be effectively eliminated.holographic image and the theoretical image without interference and a favourable match between the estimated values and the separated values can be observed.According to Eq. ( 11), the error between the calculated pressure amplitudes in Fig. 5(d) and in Fig. 5(e) is only 0.89 %.Conclusively, from the results at the 64 measuring points (8×8) as shown in Fig. 5(a), 289 holographic points (17×17) can be obtained, as shown in Fig. 5(d).By using the proposed method, we can not only eliminate the image distortion induced by the interference source but also accurately estimate the lost sound field information caused by fewer measuring points.In actual applications, the sound field information, which should be measured using 289 measuring points (17×17), can now be measured using 64 measuring points (8×8).This means that the workload due to the measurements is reduced by 77.85 %, which further proves the superiority of the proposed method described in this article.

Applicability research
In order to further verify the universal applicability of the proposed method, the relative separation and Patch errors for different values of (from 0.1 to 3) which means the frequency was traversed within the range from 54.6 Hz to 1637.7 Hz were calculated.Fig. 6 shows the relationship between the relative errors and .The relative errors are always below 3.0 %, which demonstrates that the processed field can be used to replace the free radiation field since the influence of the disturbing effect can be neglected at these frequencies.GUO LIANG, ZHU HAICHAO, MAO RONGFU, SU JUNBO, SU CHANGWEI variation of the signal-to-noise ratio (SNR).Fig. 7 shows the relative errors with the variation of the SNR.It can be found that the relative errors between the theoretical free-field pressure and measurement data are larger than the separation and Patch errors.As shown in Fig. 7, the relative errors increased to 8 % at lower SNR, while the errors decreased slightly as the SNR increased, even below 1 %.These results again demonstrates that the interference pressure and the free-field radiated pressure can be separated effectively, and the interpolated and extrapolated data show good agreement with their matching theoretical values.When the frequency fluctuates across a wide range, satisfactory results can be obtained when the number of measuring points is selected between 36 and 196.Fig. 8 shows the separation errors with different numbers of measuring points.If a conventional separation method requires 9×9 points, 6×6 measuring points can be used to obtain the same results by using the proposed method, which means a saving of 40 % workload and time measurement cost.

Fig. 8. The separation errors with different numbers of measuring points
To further verify the generality of the proposed method, we also calculated the conditions when different kinds of sound sources were located on the two sides of the measurement planes.The target source was a 0.5 m×0.5 m rectangular panel, 0.008 m thickness, 7800 kg/m 3 density, 2×10 11 Pa Young modulus, 0.28 Poisson's ratio and centred at the origin of coordinates.The panel was driven at the position (0.125 m, 0.125 m, 0) by a point force of amplitude 10 N. Its radiation was calculated with a numerical approximation to Rayleigh's integral.The interference source was a rigid pulsating sphere placed at (0, 0, 0.12 m), 0.05 m radius and 0.08 m/s vibration velocity.The sound pressure was measured at 8×8 points in two measurement planes of identical dimensions 1 m×1 m with a distance of 0.05 m between them.The two measurement planes were 0.04 m and 0.09 m from the plate.As the panel mainly vibrates with an order of (2, 2) at a modal frequency 612 Hz and the vibration amplitude of the panel was the largest, 600 Hz was selected as the excitation frequency for obvious results.A random white noise with an intensity of 30 dB was added to the measured pressure on each measurement surface to simulate an actual measurement.
Fig. 9(a) shows the amplitudes of the measurement data of size 8×8 on the holographic plane .After the proposed separation step, the pressure responses generated by the target source independently on the measurement plane are displayed in Fig. 9(b).Fig. 9(c) compares the theoretical values of size 8×8 without the influence of the disturbing source.By performing the data interpolation and extrapolation step, holographic data of size 17×17 based on separated data of size 17×17 could be obtained, as depicted in Fig. 9(d).For comparison, theoretical data of size 17×17 without interference of the disturbing source was calculated, as displayed in Fig. 9(e).
As shown in Fig. 9, when the interference source exists, the high-sound-pressure regions obviously shift to greater peak values than the theoretical values, and the vibration mode of the panel presented in the holographic image approaches the (2, 2) mode.After processing the proposed two steps, the pressure distribution in the sound field can be accurately reconstructed, and the spatial resolution of the sound field can be effectively enhanced.According to Eq. ( 11)

Experiments
As shown in Fig. 10, the experiments were performed in an ordinary factory building.Two loudspeakers were then selected as sources for an in-depth study.The size of each loudspeaker was 0.23 m×0.14 m×0.13 m.One loudspeaker was adopted as the target sound source, centred as the origin of coordinates.The other loudspeaker driven coherently was used as disturbing source, placed opposite to the target source, at (0, -0.05 m, 0.6 m).The excitation signals of the two loudspeakers were generated by the same channel in the signal generator and then transmitted.The sound pressures were measured in a grid of 12×12 points with 0.05 m interspacing at the First, the loudspeaker as the target source was initiated, and the pressure distribution on the holographic plane was scanned and adopted as the theoretical value for further comparison, as displayed in Fig. 11(a).Subsequently, the two loudspeakers were turned on simultaneously and the sound pressure on two holographic planes, and , were acquired, as displayed in Fig. 11(b) and (c).By comparison, it can be observed that, due to the disturbance of the interference source, the pseudo-peaks appeared.The measurement errors were calculated to be 88.22 % and 95.78 %, respectively.Then, the data with a dimension of 6×6 were extracted and displayed in Fig. 11(e), (f) and (g), respectively.By performing the proposed separation step, the pressure responses generated by the target source independently on the measurement plane were separated, displayed in Fig. 11(h).According to Eq. (11), the separation error between the sound-field separation results and the theoretical values using 6×6 measuring points was calculated as 11.22 %.Then, the processed results via data interpolation and extrapolation could be obtained, depicted in Fig. 11(b).The Patch error between the processed results and the theoretical values using 12×12 measuring points was calculated as 12.38 %.To verify the effectiveness of the experimental results, experiments at other frequencies were also conducted and the relative errors were listed in Table 1.
To further verify the proposed method, a boundary constraint steel plate with a size of 0.5 m×0.5 m and the thickness of 0.001 m was selected as the target source, and a loudspeaker with a size of 0.23 m×0.14 m×0.13 m was selected as the interference source, as shown in Fig. 12.The plate was set as the coordinate plane and a rectangular coordinate system was constructed in which the interference source was located at (0, -0.05 m, 0.7 m).In total, 12 pressure sensors were included in the microphone array.The sound pressure was measured at 12×12 points in two measurement planes, and , with identical dimensions 0.66 m×0.66 m, located at (0, 0, 0.04 m) and (0, 0, 0.1 m), respectively.The excitation frequency was set to 292 Hz, the sampling frequency was set to 2048 Hz and the sampling time was 1 s.The complex pressure at each measuring point was calculated using the transfer function of a single reference source.First, the plate was excited, and the pressure distribution on the holographic plane was scanned and extracted as the theoretical value for further comparison and analysis, as displayed in Fig. 13(f) and 13(a).Subsequently, the plate and loudspeaker were turned on simultaneously, and the pressure distributions on the two holographic planes, and , could be obtained.By comparison, due to the disturbance of the interference source, the errors are comparatively large, calculated as 60.16 % and 83.91 %, respectively.To verify the accuracy of the proposed technique, the data with a dimension of 6×6 points were then extracted from the pressure distributions on the two holographic planes, and , depicted in Fig. 13(b) and 13(c).
By performing the field separation step, the separated result was obtained, as displayed in Fig. 13(d).According to Eq. ( 11), the separated error was 12.75 %.Then, via data interpolation and extrapolation on the separated results, the Patch results were obtained, as displayed in Fig. 13(e).Compared with the theoretical values in Fig. 13(f), the processed results were quite close to the theoretical results, and the error was calculated to be only 10.24 %.The results indicated that the effects induced by the interference noise could be eliminated, and moreover, the spatial resolution of the holographic images could be effectively enhanced.On the other hand, 6 pressure sensors were sufficient to replace the 12 pressure sensors.The workload just be 6×6/(12×12) = 0.25 = 25 % of the actual workload, i.e., the workload and time consumption in the measurement could be significantly reduced.
Guo Liang made a major contribution to the research.Zhu Haichao made guidance work to the research.Mao Rongfu gave consultative opinions to the paper.Su Junbo and Su Changwei made simulation calculations to the paper.

Conclusions
To solve the problems of NAH in actual applications (such as the required free field and aperture effect), a double layer Patch NAH based on sound radiation modes was established.First, the mathematical relationship between the pressure on the measurement plane and the field distribution modes was derived.The effects of the interference sources were then eliminated using the double layer sound field separation technique.Then, the pressure values at more discrete points using data interpolation and extrapolation technique were estimated, which could provide more accurate holographic data for reconstructing the sound field.When performing the sound field separation and Patch step, the field distribution modes with the same order were selected and therefore, a uniform study and analysis of the relationship between the error and various parameters could be studied.With comprehensive consideration of the effects induced by the proposed method, the optimization selection and control on each parameter could be realized.
As indicated by the simulated and experimental results, the proposed method exhibits favourable universal applicability, and holographic data that are closer to the actual results using fewer measuring points were obtained, indicating that the processed technique is simple, economic and convenient for engineering applications.

Fig. 1 .
Fig. 1.Illustration of the spatial arrangement of the target source, interference source and measurement planes.On the measurement plane , the pressure responses generated by the target and interference

d 1 d 2 d
1977.SOUND FIELD SEPARATION TECHNIQUE USING THE PRINCIPLE OF DOUBLE LAYER PATCH ACOUSTIC RADIATION MODES.GUO LIANG, ZHU HAICHAO, MAO RONGFU, SU JUNBO, SU CHANGWEI is related to the field distribution modes of the sound source on the measurement plane, while the sound field distribution modes have a close relationship with the corresponding spatial positions.Fig.2displays the spatial distribution of the measurement plane and the holographic plane after data interpolation and extrapolation.

Fig. 2 .
Fig. 2. Illustration of data interpolation and extrapolation when using the Patch technique

1 P
1977.SOUND FIELD SEPARATION TECHNIQUE USING THE PRINCIPLE OF DOUBLE LAYER PATCH ACOUSTIC RADIATION MODES.
1977.SOUND FIELD SEPARATION TECHNIQUE USING THE PRINCIPLE OF DOUBLE LAYER PATCH ACOUSTIC RADIATION MODES.

Fig. 3 .Fig. 4 .
Fig. 3.The first nine-order sound radiation modes of the rigid pulsating sphere (a), (b) and (c) indicate that, when the interference

5 .
a) Initial data of size 8×8 b) Separated data of size 8×8 c) Theoretical data of size 8×8 d) Processed data of size 17×17 e) Theoretical data of size 17×17 Fig. Distributions of the pressure amplitudes on the measurement plane when = 0.5 As shown in Fig. 5(b), (d) and (e), after two-fold data interpolation and extrapolation, the gradually-varied regions in the image become smoother.The sound field in the high-sound-pressure regions appears to be closer to the actual sound field generated by the rigid pulsating sphere, and some details in the images emerge.There is little difference between the xFIELD SEPARATION TECHNIQUE USING THE PRINCIPLE OF DOUBLE LAYER PATCH ACOUSTIC RADIATION MODES.GUO LIANG, ZHU HAICHAO, MAO RONGFU, SU JUNBO, SU CHANGWEI

Fig. 9 .
was calculated, and found to be 5.59 %.The results after performing the two-fold data interpolation and extrapolation on the separated data were compared to the theoretical values and the calculated Patch error was 6.46 %. a) Initial data of size 8×8 b) Separated data of size 8×8 c) Theoretical data of size 8×8 d) Processed data of size 17×17 e) Theoretical data of size 17×17 Distributions of the pressure amplitudes on the holographic plane when = 0.55 FIELD SEPARATION TECHNIQUE USING THE PRINCIPLE OF DOUBLE LAYER PATCH ACOUSTIC RADIATION MODES.GUO LIANG, ZHU HAICHAO, MAO RONGFU, SU JUNBO, SU CHANGWEI holographic planes = 0.1 m and = 0.3 m.

Fig. 10 .
Fig. 10.A picture of the experiment set-up

13 .
a) The extraction of the measured values of size 6×6 without interference on b) The extraction of the measured values of size 6×6 with interference on c) The extraction of the measured values of size 6×6 with interference on d) The separated values of size 6×6 on e) The processed results of size 12×12 using the double-layer sound field separation method on f) The measured values of size 12×12 without interference on Fig. Distributions of the pressure amplitude xFIELD SEPARATION TECHNIQUE USING THE PRINCIPLE OF DOUBLE LAYER PATCH ACOUSTIC RADIATION MODES.GUO LIANG, ZHU HAICHAO, MAO RONGFU, SU JUNBO, SU CHANGWEI

Table 1 .
The relative errors at different frequencies SOUND FIELD SEPARATION TECHNIQUE USING THE PRINCIPLE OF DOUBLE LAYER PATCH ACOUSTIC RADIATION MODES.GUO LIANG, ZHU HAICHAO, MAO RONGFU, SU JUNBO, SU CHANGWEI