GLOBAL ACTIVE NOISE CONTROL METHOD FOR ROTORCRAFT

Abstract

A global active noise control method for a rotorcraft, including: acquiring the acoustic pressure signal at a measuring point of the rotorcraft; predicting the holographic and global sound field of noise of the rotor; reconstructing the reverse sound field of the noise of the rotor; and performing adaptive sound field adjustment based on the optimal phase search.

Claims

1. A global active noise control method for a rotorcraft, comprising: measuring, by an acoustic measuring device array arranged on the rotorcraft, noise of the rotorcraft; inputting a noise pressure signal of the rotorcraft to acquire an acoustic mode expansion form of a global rotor noise by using an acoustic analysis method; online estimating optimal acoustic modal coefficients based on acoustic holography by using a measurement signal of the acoustic measuring device array to obtain an acoustic holographic global sound field of the rotor; based on the acoustic holographic global sound field, online generating a secondary sound field that achieves global sound-sound cancellation with the global rotor noise according to a sound field construction method by using a secondary acoustic source array; inputting the optimal acoustic modal coefficients to online calculate a real-time control signal of the secondary acoustic source array according to an acoustic modal orthogonal relationship; and online adjusting the real-time control signal of the secondary acoustic source array by using an adaptive method to realize global noise reduction of the rotor under different flight conditions.

2. The global active noise control method of claim 1, wherein the acoustic mode expansion form of the global rotor noise is obtained through steps of: for an aerodynamic noise of the rotor, a tip speed of which is less than speed of sound, setting a Ffowcs-Williams Hawking acoustic analogy equation as equation (1), wherein noise outside a rotor rotation area satisfies a passive homogeneous wave equation (2); and introducing a Fourier transform for derivation; shown as follows: $\begin{array}{cc}{\nabla }^{2}p-\frac{1}{c^2}\frac{{\partial }^{2}p}{\partial t^2}=-\frac{\partial }{\partial t}\left[{\rho }_0{v}_n.Math.\nabla f.Math.\delta \left(f\right)\right]+\frac{\partial }{\partial x_i}\left[l_i.Math.\nabla f.Math.\delta \left(f\right)\right];\mathrm{and}& \left(1\right)\\ {\nabla }^{2}p-\frac{1}{c^2}\frac{{\partial }^{2}p}{\partial t^2}=0;& \left(2\right)\end{array}$ wherein p is a sound pressure; c is the speed of sound; v.sub.n is a normal velocity of a blade surface; ρ.sub.0 is air density; l.sub.i is a load per unit area of a medium; f(x,t)=0 is a boundary of the blade surface; δ(f) indicates that thickness and load noise sources are only distributed on the blade surface, and are surface sound sources; r, θ, ϕ respectively represent a distance from an observation point to an origin, an elevation angle, and an azimuth angle; ω is a noise frequency; and k≙ω/c represents a wave number; expressing an arrangement position of the acoustic measuring device array in a spherical coordinate system as equation (3), wherein a frequency domain form of an acoustic wave equation of the spherical coordinate system is expressed by equation (4); shown as follows: $\begin{array}{cc}{r}_j=\left(r_j,{\theta }_j,{\varphi }_j\right),j=1\mathrm{.Math.}J;\mathrm{and}& \left(3\right)\\ {\nabla }^{2}p+k^{2}p=0;& \left(4\right)\end{array}$ based on a Fourier acoustic analysis method, a series expansion form of a rotor noise solution that meets Sommerfeld radiation condition is expressed as equation (5): $\begin{array}{cc}{p}^d\left(r,\theta ,\varphi ,k\right)=\underset{n=0}{\overset{\infty }{.Math.}}{h}_n^{\left(1\right)}\left(\mathrm{kr}\right)\underset{m=-n}{\overset{n}{.Math.}}{C}_{m,n}\left(k\right)Y_n^m\left(\theta ,\varphi \right);& \left(5\right)\end{array}$ wherein C.sub.m,n(k) is an acoustic modal coefficient; h.sub.n.sup.(1)(kr) represents a first-order Spherical Hankel function; and Y.sub.n.sup.m(θ,ϕ) represents a spherical harmonics function.

3. The global active noise control method of claim 2, wherein the optimal acoustic modal coefficients are estimated online through steps of: in a specified basis function Ψ.sub.n,m.sup.(1), performing an optimal approximation on a noise pressure signal of a measuring point to estimate the optimal acoustic modal coefficients, wherein the acoustic modal coefficient and the noise pressure signal of the measuring point of the acoustic measuring device array meet the following equations: $\begin{array}{cc}{\left\{p^d\left(r_j,{\theta }_j,{\varphi }_j,k\right)\right\}}_{J\times 1}={\left[{\Psi }^{\left(1\right)}\left(r_j,{\theta }_j,{\varphi }_j,k\right)\right]}_{J\times \left({\left(N^{\prime }+1\right)}^2\right)}{\left\{C_{m,n}\left(k\right)\right\}}_{\left({\left(N^{\prime }+1\right)}^2\right)\times 1};\mathrm{and}& \left(6\right)\\ {\Psi }_{n,m}^{\left(1\right)}\left(r_j,{\theta }_j,{\varphi }_j,k\right)\equiv h_n^{\left(1\right)}\left({\mathrm{kr}}_j\right)Y_n^m\left({\theta }_j,{\varphi }_j\right);\mathrm{and}& \left(7\right)\end{array}$ solving the optimal acoustic modal coefficients by using a regularization method, expressed as: $\begin{array}{cc}\left\{C_{m,n}\left(k\right)\right\}={{\left({\left[{\Psi }^{\left(1\right)}\right]}^H\left[{\Psi }^{\left(1\right)}\right]\right)}^{-1}\left[{\Psi }^{\left(1\right)}\right]}^H\left\{p^d\right\}.& \left(8\right)\end{array}$

4. The global active noise control method of claim 3, wherein the secondary sound field is generated online through steps of: expressing an arrangement position of the secondary acoustic source array in the spherical coordinate system as r.sub.s=(r.sub.s, θ.sub.s, ϕ.sub.s), s=1 . . . S, wherein a sound field generated by the secondary acoustic source array is expressed as equation (9); and a reconstructed target sound field meets equation (10); shown as follows: $\begin{array}{cc}{p}^S\left(r,\theta ,\varphi ,k\right)=\underset{s=1}{\overset{S}{.Math.}}{p}^s\left(r,\theta ,\varphi ,k\right).Math.r.Math.\ge \underset{1\le s\le S}{\max }.Math.r_s.Math.=\underset{n=0}{\overset{\infty }{.Math.}}{h}_n^{\left(1\right)}\left(\mathrm{kr}\right)\underset{m=-n}{\overset{n}{.Math.}}\mathrm{ik}\left(\underset{s=1}{\overset{S}{.Math.}}{Q}_s{j}_n\left({\mathrm{kr}}_s\right){Y_n^m\left({\theta }_s,{\varphi }_s\right)}^{*}\right)Y_n^m\left(\theta ,\varphi \right);\mathrm{and}& \left(9\right)\\ p^S\left(r,\theta ,\varphi ,k\right)=-p^d\left(r,\theta ,\varphi ,k\right)=\underset{n=0}{\overset{\infty }{.Math.}}{h}_n^{\left(1\right)}\left(\mathrm{kr}\right)\underset{m=-n}{\overset{n}{.Math.}}-C_{m,n}\left(k\right)Y_n^m\left(\theta ,\varphi \right);& \left(10\right)\end{array}$ wherein Q.sub.S is a mass-source intensity of a secondary acoustic source, and Q.sub.S=−iωρ.sub.0q.sub.S; and q.sub.S is a volume-source intensity of the secondary acoustic source.

5. The global active noise control method of claim 4, wherein the real-time control signal of the secondary acoustic source array is calculated online through steps of: according to equation (9) and equation (10) and based on an orthogonality of the acoustic modal coefficient, adjusting the real-time control signal of the secondary acoustic source array until a source intensity meets equation (11) to generate a reverse sound field of a rotor noise, wherein matrix parameters in the equation (11) are expressed by equations (12)-(15): $\begin{array}{cc}\mathrm{ikJTQ}=-C;& \left(11\right)\\ Q={{\left[\begin{array}{ccc}{Q}_1& \mathrm{.Math.}& Q_S\end{array}\right]}^T}_{S\times 1};& \left(12\right)\\ T={\left\{{Y_n^m\left({\theta }_s,{\varphi }_s\right)}^{*}\right\}}_{{\left(N^{\prime }+1\right)}^2\times S};& \left(13\right)\\ J={\left\{j_n\left({\mathrm{kr}}_s\right)\right\}}_{{\left(N^{\prime }+1\right)}^2\times {\left(N^{\prime }+1\right)}^2};\mathrm{and}& \left(14\right)\\ C={\left\{C_{m,n}\left(k\right)\right\}}_{{\left(N^{\prime }+1\right)}^2\times 1};& \left(15\right)\end{array}$ wherein matrix Q represents a sound source intensity of each unit of the secondary acoustic source array; matrix T indicates that an independent vector set of an acoustic modal space generated by the secondary acoustic source array is determined by an azimuth angle ϕ.sub.s and an elevation angle θ.sub.s of the secondary acoustic source array; characteristics of the sound field generated by the secondary acoustic source array are determined by a radius r.sub.s of the secondary acoustic source array, and are reflected in low-pass characteristics of a function j.sub.n(kr.sub.s) of a diagonal matrix J with respect to order n; matrix T is not a square matrix; the real-time control signal of the secondary acoustic source array is calculated by regularization.

6. The global active noise control method of claim 1, wherein the sound field construction method is a high-order ambient stereo method, a wave field synthesis method, a spherical harmonic decomposition method, or a combination thereof.

7. The global active noise control method of claim 1, wherein the adaptive method is an exponential phase online search method.

8. The global active noise control method of claim 5, wherein during reconstruction of the reverse sound field of the rotor, when a phase change caused by a speed fluctuation or flight condition of the rotor exceeds a threshold, the adaptive method is used to update a phase and adjust the real-time control signal of the secondary acoustic source array online to realize adaptive reconstruction of the reverse sound field.

9. The global active noise control method of claim 1, wherein the acoustic measuring device array and the secondary acoustic source array are arranged in the rotorcraft; the acoustic measuring device array is configured to collect a noise pressure signal data at a measuring point; and the secondary acoustic source array is configured to online generate the secondary sound field that offsets the global noise of the rotor.

10. A computer-readable storage medium, wherein the computer-readable storage medium stores a computer program; and the computer program is configured to be executed by a processor to implement the global active noise control method of claim 1.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0040] FIG. 1 is a flow chart of a global active noise control method for a rotorcraft according to an embodiment of the present disclosure;

[0041] FIG. 2 is a flow chart of an adaptive online adjustment of a sound field according to an embodiment of the present disclosure; and

[0042] FIGS. 3A-3D schematically shows spherical acoustic pressure distributions before and after the global noise reduction according to an embodiment of the present disclosure, where 3A: original sound field of the rotor noise (r=0.7 m); 3B: sound field after the global noise control (r=0.7 m); 3C: original sound field of the rotor noise (r=1.4 m); and 3D: sound field after the global noise control (r=1.4 m).

DETAILED DESCRIPTION OF EMBODIMENTS

[0043] The technical solutions of the present disclosure will be clearly and completely described below with reference to the accompanying drawings and embodiments. Obviously, the described embodiments are only illustrative, and are not intended to limit the disclosure. Based on the embodiments of the present disclosure, all other embodiments obtained by those skilled in the art without paying creative efforts shall fall within the scope of the present disclosure.

[0044] Referring to FIG. 1, provided is a global active noise control method based on acoustic holography and sound field reconstruction, which includes the acoustic pressure signal acquisition of noise at a measuring point by means of an acoustic measuring device array on the rotorcraft, the holographic and global sound field calculation of the rotor, the generation of a sound field that offsets the global noise of the rotor, and the online sound field adjustment based on an adaptive method. The global active noise control method specifically includes the following steps.

[0045] (S1) Noise Acoustic Pressure Signals Acquisition

[0046] Measuring points, the acoustic measuring device array, and the secondary acoustic source array are arranged on the rotorcraft based on the analysis of the shape and basic structure of the rotorcraft. The acoustic measuring device array is configured to collect an acoustic pressure signal data of noise at a measuring point. The secondary acoustic source array is configured to online generate the secondary sound field that offsets the global noise of the rotor.

[0047] The common sampling forms of the acoustic measuring device array and the secondary acoustic source array outside the rotating region include but are not limited to uniform sampling, Gaussian sampling, approximately uniform sampling, etc.

[0048] (S2) Holographic and Global Sound Field Calculation

[0049] An acoustic pressure signal of the noise of the rotorcraft is inputted to acquire an acoustic mode expansion form of a global noise of a rotor by using an acoustic analysis method. Then a measurement signal of the acoustic measuring device array is used to online estimate an optimal acoustic modal coefficient based on an acoustic holography method to obtain an acoustic holographic global sound field of the rotor.

[0050] For the aerodynamic noise of the rotor, the Ffowcs-Williams Hawking acoustic analogy equation is equation (1). As the Dirichlet function δ(f) is only meaningful on the object plane, the sound source term on the right side of equation (1) only appears in the bounded rotor rotation area. When noise outside a rotor rotation area satisfies a passive homogeneous wave equation (2), a Fourier transform is introduced for derivation, shown as follows:

[00008]$\begin{array}{cc}{\nabla }^{2}p-\frac{1}{c^2}\frac{{\partial }^{2}p}{\partial t^2}=-\frac{\partial }{\partial t}\left[{\rho }_0{v}_n.Math.\nabla f.Math.\delta \left(f\right)\right]+\frac{\partial }{\partial x_i}\left[l_i.Math.\nabla f.Math.\delta \left(f\right)\right]-\frac{{\partial }^2}{\partial x_i{\mathrm{\psi x}}_j}\left[T_{\mathrm{ij}}H\left(f\right)\right];\mathrm{and}& \left(1\right)\\ {\nabla }^{2}p-\frac{1}{c^2}\frac{{\partial }^{2}p}{\partial t^2}=0;& \left(2\right)\end{array}$

[0051] where p is a sound pressure; c is the speed of sound; v.sub.n is a normal velocity of a surface of a blade; ρ.sub.0 is an air density; l.sub.i is a load per unit area of a medium; f(x,t)=0 represents a surface motion equation; δ(f) indicates that thickness and load noise sources are only distributed on the surface of the blade, and are surface sound sources; r, θ, ϕ respectively represent a distance from an observation point to an origin, an elevation angle, and an azimuth angle; ω is a noise frequency; and k≙ω/c represents a wave number.

[0052] In a spherical coordinate system, an arrangement position of the acoustic measuring device array is expressed as equation (3), and a frequency domain form of an acoustic wave equation of the spherical coordinate system is expressed by equation (4), shown as follows:

[00009]$\begin{array}{cc}{r}_j=\left(r_j,{\theta }_j,{\varphi }_j\right),j=1\mathrm{.Math.}J;\mathrm{and}& \left(3\right)\\ {\nabla }^{2}p+k^{2}p=0.& \left(4\right)\end{array}$

[0053] In addition, the sound field solution represented by equation (2) should also satisfy the two boundary conditions (Somerfield radiation condition), namely, the sound pressure is continuous at the measurement point and the rotor noise sound pressure approaches 0 at infinity. Then based on a Fourier acoustic analysis method, a series expansion form of a rotor noise solution that meets a Sommerfeld radiation condition is expressed as equation (5):

[00010]$\begin{array}{cc}{p}^d\left(r,\theta ,\varphi ,k\right)=\underset{n=0}{\overset{\infty }{.Math.}}{h}_n^{\left(1\right)}\left(\mathrm{kr}\right)\underset{m=-n}{\overset{n}{.Math.}}{C}_{m,n}\left(k\right)Y_n^m\left(\theta ,\varphi \right);& \left(5\right)\end{array}$

[0054] where C.sub.m,n(k) is an acoustic modal coefficient, which is merely related to the acoustic mode order and wavenumber; the acoustic mode distribution of the rotor noise is closely related to the number of rotor blades; h.sub.n.sup.(1)(kr) represents a first-order Spherical Hankel function, which describes the changing law of the acoustic mode in the radius; and Y.sub.m,n(θ,ϕ) represents a spherical harmonics function, which can describe the changing law of the acoustic mode in the azimuth and elevation.

[0055] Moreover, considering that there are unavoidable errors in the installation of the acoustic measuring device array, which will affect the measurement signal of the noise of the rotor. The commonly used method of calculating the acoustic mode based on the weighting coefficient does not consider the effect of those errors. The HELS method (expressed by equation (6)) developed by S. F. Wu et al. is employed, through which an optimal approximation is performed on an acoustic pressure signal of noise of a measuring point in a specified basis function Ψ.sub.n,m.sup.(1) to estimate the optimal acoustic modal coefficient, where the acoustic modal coefficient and the acoustic pressure signal of the measuring point of the acoustic measuring device array meet the following equations:

[00011]$\begin{array}{cc}{\left\{p^d\left(r_j,{\theta }_j,{\varphi }_j,k\right)\right\}}_{J\times 1}={\left[{\Psi }^{\left(1\right)}\left(r_j,{\theta }_j,{\varphi }_j,k\right)\right]}_{J\times \left({\left(N^{\prime }+1\right)}^2\right)}{\left\{C_{m,n}\left(k\right)\right\}}_{\left({\left(N^{\prime }+1\right)}^2\right)\times 1};\mathrm{and}& \left(6\right)\\ {\Psi }_{n,m}^{\left(1\right)}\left(r_j,{\theta }_j,{\varphi }_j,k\right)\equiv h_n^{\left(1\right)}\left({\mathrm{kr}}_j\right)Y_n^m\left({\theta }_j,{\varphi }_j\right).& \left(7\right)\end{array}$

[0056] Since the number of measurement points is generally more than the truncation term, the optimal acoustic modal coefficient can be solved by using a regularization method, expressed as:

[00012]$\begin{array}{cc}\left\{C_{m,n}\left(k\right)\right\}={{\left({\left[{\Psi }^{\left(1\right)}\right]}^H\left[{\Psi }^{\left(1\right)}\right]\right)}^{-1}\left[{\Psi }^{\left(1\right)}\right]}^H\left\{p^d\right\}.& \left(8\right)\end{array}$

[0057] (S3) Generation of a Sound Field that Offsets the Global Noise of the Rotor

[0058] Based on the holographic and global sound field obtained in step (S2), the target sound field to be reconstructed is analyzed according to a sound field construction method (i.e., high-order ambient stereo, wave field synthesis, and spherical harmonic decomposition) by using a secondary acoustic source array. Based on this, the control signal of the monopole sound source group is extracted based on the acoustic modal orthogonal relationship and the matching relationship between sound fields.

[0059] In an embodiment, the monopole sound source group is generated by the secondary acoustic source array, and the high-order ambient stereo method is used to realize the reverse sound field reconstruction. This method can be unified with the Helmholtz equation least square method in step (S2) to facilitate modeling calculations. First, an arrangement position of the secondary acoustic source array is denoted as r.sub.s=(r.sub.s, θ.sub.s, ϕ.sub.s), s=1 . . . S. Further, a sound field generated by the secondary acoustic source array is expressed as equation (9), which indicates that the acoustic modal coefficient generated by the secondary acoustic source array is uniquely determined by the source intensity of the secondary acoustic source array. Through adjusting the control signal of the secondary acoustic source array, any target sound fields can be generated, and the target sound field reconstructed based on acoustic cancellation is the reverse sound field of the rotor noise to achieve global noise reduction. Therefore, the reconstructed target sound field meets equation (10), shown as follows:

[00013]$\begin{array}{cc}{p}^S\left(r,\theta ,\varphi ,k\right)=\underset{s=1}{\overset{S}{.Math.}}{p}^s\left(r,\theta ,\varphi ,k\right).Math.r.Math.\ge \underset{1\le s\le S}{\max }.Math.r_s.Math.=\underset{n=0}{\overset{\infty }{.Math.}}{h}_n^{\left(1\right)}\left(\mathrm{kr}\right)\underset{m=-n}{\overset{n}{.Math.}}\mathrm{ik}\left(\underset{s=1}{\overset{S}{.Math.}}{Q}_s{j}_n\left({\mathrm{kr}}_s\right){Y_n^m\left({\theta }_s,{\varphi }_s\right)}^{*}\right)Y_n^m\left(\theta ,\varphi \right);\mathrm{and}& \left(9\right)\\ p^S\left(r,\theta ,\varphi ,k\right)=-p^d\left(r,\theta ,\varphi ,k\right)=\underset{n=0}{\overset{\infty }{.Math.}}{h}_n^{\left(1\right)}\left(\mathrm{kr}\right)\underset{m=-n}{\overset{n}{.Math.}}-C_{m,n}\left(k\right)Y_n^m\left(\theta ,\varphi \right);& \left(10\right)\end{array}$

[0060] where Q.sub.S is a mass-source intensity of the secondary acoustic source array; and Q.sub.S=−iωρ.sub.0q.sub.S, and q.sub.S are a volume-source intensity of the secondary acoustic source array.

[0061] according to equation (9) and equation (10) and based on an orthogonality of the acoustic modal coefficient, adjusting the real-time control signal of the secondary acoustic source array until a source intensity meets equation (11) to generate a reverse sound field of a rotor noise, wherein matrix parameters in the equation (11) are expressed by equations (12)-(15):

[00014]$\begin{array}{cc}\mathrm{ikJTQ}=-C;& \left(11\right)\\ Q={{\left[\begin{array}{ccc}{Q}_1& \mathrm{.Math.}& Q_S\end{array}\right]}^T}_{S\times 1};& \left(12\right)\\ T={\left\{{Y_n^m\left({\theta }_s,{\varphi }_s\right)}^{*}\right\}}_{{\left(N^{\prime }+1\right)}^2\times S};& \left(13\right)\\ J={\left\{j_n\left({\mathrm{kr}}_s\right)\right\}}_{{\left(N^{\prime }+1\right)}^2\times {\left(N^{\prime }+1\right)}^2};\mathrm{and}& \left(14\right)\\ C={\left\{C_{m,n}\left(k\right)\right\}}_{{\left(N^{\prime }+1\right)}^2\times 1};& \left(15\right)\end{array}$

[0062] where matrix Q represents a sound source intensity of each unit of the secondary acoustic source array; matrix T indicates that an independent vector set of an acoustic modal space generated by the secondary acoustic source array is determined by an azimuth angle ϕ.sub.s and an elevation angle θ.sub.s of the secondary acoustic source array; characteristics of the sound field generated by the secondary acoustic source array are determined by a radius r.sub.s of the secondary acoustic source array, and are reflected in low-pass characteristics of a function j.sub.n(kr.sub.s) of a diagonal matrix J with respect to order n; matrix T is not a square matrix; the real-time control signal of the secondary acoustic source array is calculated by regularization.

[0063] (S4) On-Line Sound Field Adjustment Based on an Adaptive Method

[0064] The adaptive method is employed to realize the online sound field adjustment, which can overcome the adverse effects of the rotation speed fluctuation or flight states to realize the global noise reduction of the rotor under different flight states.

[0065] Ideally, the rotor noise is stable, and the control signal of the secondary acoustic source array obtained based on equation (11) can make the secondary acoustic source array accurately reconstruct the reverse sound field of noise of the rotor, which can realize the global noise reduction of the rotor noise. However, in actual situations, the inevitable rotation speed fluctuation of the rotor will change the phase of noise of the rotor, which will seriously affect the acoustic cancellation effect. Therefore, to guarantee the global noise reduction effect in the actual work of the rotorcraft, it is necessary to use an adaptive control technology to perform real-time adjustment of the reconstructed sound field. In this example, the control signal of the secondary acoustic source array is employed to online adjust the optimal phase and based on the optimal phase search method to suppress the adverse effects of the rotor speed fluctuation on the noise reduction effect.

[0066] As shown in FIG. 2, during reconstruction of a reverse sound field of the rotor, when a phase change caused by a speed fluctuation of the rotor exceeds a threshold, the adaptive method is used to update a phase and adjust the real-time control signal of the secondary acoustic source array online to realize adaptive reconstruction of the reverse sound field and further improve the practical value of the present disclosure.

[0067] The global noise reduction of the rotor based on acoustic holography and sound field reconstruction utilizes the constitutive relationship of the acoustic wave equation. On one hand, the complex noise control system with multiple inputs and outputs can be reduced to the optimal phase search, which greatly reduces the amount of calculation, facilitating to realize the online active control; on the other hand, it can wholly reduce the noise of the rotor, and realize the global noise reduction of the rotor noise.

[0068] The noise reduction simulation result of the rotor based on acoustic holography and acoustic reconstruction of the present disclosure shows that when the number of the secondary acoustic source array reaches eight, 22.70 dB noise suppression can be achieved at the measuring radius of the acoustic measuring device. As shown in FIGS. 3A-D, based on the simulation results, the effect of global noise reduction on the spherical sound pressure distribution is illustrated. It can be seen from FIGS. 3A-3B that the in-plane noise is the largest, and the out-of-plane noise attenuates faster with the increase in radius than the in-plane noise. FIGS. 3B-3D shows that the method provided herein can achieve in-plane noise reduction while having a significant noise reduction effect on out-of-plane noise. The test results show that the method provided herein can achieve an overall 17.1 dB noise attenuation of the rotor noise, and meanwhile, the average noise attenuation of the secondary acoustic source array at 0.7 m is 15.8 dB, indicating that the rotor noise suppression method of the present disclosure has a good noise reduction effect.

[0069] In addition, an embodiment of the present disclosure also provides a computer-readable storage medium, which can store a program. The program is executed by a processor to implement any part or all steps of the global active noise control method described in the above embodiments.

[0070] In some embodiments, the functional units can be integrated into one processing unit, or independent, or two or more units may be integrated into one unit. The above-mentioned integrated unit can be implemented in the form of a hardware or a software functional unit.

[0071] If the integrated unit is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a computer-readable memory. Based on this, the technical solutions of the present disclosure essentially or the part that contributes to the existing technology or all or part of the technical solutions can be embodied in the form of a software product, and the computer software product is stored in a memory, including a number of instructions to enable a computer device (or a personal computer, a server, or a network device, etc.) to perform all or part of the steps of the method described in each embodiment of the present disclosure. The aforementioned memory includes a U disk, a read-only memory (ROM), a random access memory (RAM), a mobile hard disk, a magnetic disk or an optical disk and other media that can store program codes.

[0072] It should be understood by those skilled in the art that all or part of the steps in the method of the above-mentioned embodiments can be implemented by relevant hardware instructed by a program. The program can be stored in a computer-readable memory, including a flash disk, a store media of a controller, a RAM, a magnetic disk, or an optical disc.

[0073] The above-mentioned embodiments are merely illustrative of the present disclosure, and are not intended to limit the disclosure. It should be noted that any modifications, changes and replacements made by those skilled in the art without departing from the spirit of the disclosure should fall within the scope of the disclosure defined by the appended claims.