Active noise control method, apparatus, device, and storage medium

Abstract

The application discloses an active noise control method and device, equipment and a storage medium. The sound field sound pressure of an external sound field is obtained based on a high-order sensor, the external sound field being a sound field outside a preset circular region of the high-order sensor; a first sound field coefficient of the external sound field is calculated according to the sound field sound pressure, and a target sound field intensity of a secondary sound field is determined based on a first preset relationship between a sound field intensity of the secondary sound field and a second sound field coefficient of the secondary sound field and a second preset relationship between the first sound field coefficient and the second sound field coefficient, the secondary sound field being a sound field within the preset circular region; and finally, the sound power of the secondary sound field is determined according to the target sound field intensity, the sound power being used to control the secondary sound source to generate the secondary sound field so as to offset the external sound field, thereby improving the noise reduction effect.

Description

Technical Field

[0001] This application relates to the field of noise reduction technology, and in particular to an active noise control method, apparatus, device and storage medium. Background Technology

[0002] In free-space active noise control (ANC), to achieve noise reduction across the entire space, the distance between the secondary and primary sound sources must be less than half a wavelength. Currently, ANC cost functions are constructed using sound field information picked up by error sensors to reduce the total radiated sound power. The close proximity of the error sensor and secondary sound source facilitates the implementation of compact systems, making it suitable for noise radiation control applications such as transformers. However, in traditional methods, the layout of the error sensor and the secondary sound source are coupled and influence each other, making it difficult to obtain the total radiated sound power and the optimal layout result, thus leading to poor noise reduction performance. Summary of the Invention

[0003] This application provides an active noise control method, apparatus, device, and storage medium to solve the technical problem of poor noise reduction effect in current noise control methods.

[0004] To address the aforementioned technical problems, in a first aspect, this application provides an active noise control method, comprising:

[0005] The sound pressure of the external sound field is obtained based on a high-order sensor. The external sound field is the sound field outside the preset circular area of ​​the high-order sensor.

[0006] Calculate the first sound field coefficient of the external sound field based on the sound pressure of the sound field;

[0007] Based on the first preset relationship between the sound field intensity of the secondary sound field and the second sound field coefficient of the secondary sound field, and the second preset relationship between the first sound field coefficient and the second sound field coefficient, the target sound field intensity of the secondary sound field is determined, and the secondary sound field is the sound field within a preset circular region.

[0008] Based on the target sound field intensity, the sound power of the secondary sound field is determined. The sound power is used to control the secondary sound source to generate a secondary sound field to cancel out the external sound field.

[0009] In some implementations, the first sound field coefficient of the external sound field is calculated based on the sound pressure level, including:

[0010] Using the least squares method, based on the sound pressure level of the sound field, the first sound field coefficients of the external sound field are calculated. The least squares method is as follows:

[0011] β s =(H H H) -1 H H p;

[0012] Where, β s H represents the first sound field coefficient, and H is the matrix form of the second type of column Hankel function. H denoted as the conjugate transpose in matrix form, where p is the sound pressure level of the sound field.

[0013] In some implementations, the target sound field intensity of the secondary sound field is determined based on a first preset relationship between the sound field intensity of the secondary sound field and the second sound field coefficients of the secondary sound field, and a second preset relationship between the first sound field coefficients and the second sound field coefficients, including:

[0014] Based on the first and second preset relationships, a cost function for the secondary sound field is established.

[0015] The target sound field intensity of the secondary sound field is determined using a cost function, the expression of which is:

[0016]

[0017] Among them, J WD (q s ) represents the output value of the cost function. Γ represents the conjugate transpose of the target sound field intensity, and Γ is the sound transmission impedance coefficient matrix. H Let q be the conjugate transpose of the acoustic transmission impedance coefficient matrix. s β represents the target sound field intensity. p Indicates the second sound field coefficient. This represents the conjugate transpose of the second sound field coefficient.

[0018] In some implementations, a cost function for the secondary sound field is established based on a first preset relationship and a second preset relationship, including:

[0019] Based on the first and second preset relationships, the sum of squares of the sound field coefficients of the total sound field is minimized to generate the cost function. The sum of squares of the sound field coefficients of the total sound field is: β represents the conjugate transpose of the total sound field coefficients. t This represents the total sound field coefficient.

[0020] In some implementations, the first preset relation is: β s =Γq s .

[0021] In some implementations, the second preset relationship is: β t =β p +β s .

[0022] In some implementations, the acoustic power of the secondary sound field is determined based on the target sound field intensity, including:

[0023] Using a preset sound power calculation formula, the sound power of the secondary sound field is calculated based on the target sound field intensity. The sound power calculation formula is as follows:

[0024]

[0025] Among them, W t The sound power to be controlled for the secondary sound source. q represents the conjugate transpose of the target sound field intensity. s The target sound field intensity is represented by A and b, which are preset sound field parameters. W0 is the sound power of the secondary sound source before control.

[0026] Secondly, this application provides an active noise control device, comprising:

[0027] The acquisition module is used to acquire the sound pressure of the external sound field based on the high-order sensor. The external sound field is the sound field outside the preset circular area of ​​the high-order sensor.

[0028] The calculation module is used to calculate the first sound field coefficient of the external sound field based on the sound pressure of the sound field.

[0029] The first determining module is used to determine the target sound field intensity of the secondary sound field based on the first preset relationship between the sound field intensity of the secondary sound field and the second sound field coefficient of the secondary sound field, and the second preset relationship between the first sound field coefficient and the second sound field coefficient. The secondary sound field is a sound field within a preset circular region.

[0030] The second determining module is used to determine the sound power of the secondary sound field based on the target sound field intensity. The sound power is used to control the secondary sound source to generate a secondary sound field to cancel out the external sound field.

[0031] Thirdly, this application provides a computer device including a processor and a memory, the memory being used to store a computer program, which, when executed by the processor, implements the active noise control method as described in the first aspect.

[0032] Fourthly, this application provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements the active noise control method as described in the first aspect.

[0033] Compared with existing technologies, it has at least the following beneficial effects:

[0034] By acquiring the sound pressure of the external sound field using a high-order sensor (the sound field outside a preset circular region of the high-order sensor), the noise reduction effect can be improved by utilizing the high-order sensor as an error sensor. Based on the sound pressure, the first sound field coefficient of the external sound field is calculated, along with a first preset relationship between the sound field intensity and the second sound field coefficient of the secondary sound field, and a second preset relationship between the first and second sound field coefficients. This determines the target sound field intensity of the secondary sound field (the sound field within a preset circular region), thereby accurately estimating the sound field coefficients without needing to consider the secondary sound source layout. This decouples the error sensor layout from the secondary sound source layout, allowing for separate optimization of the error sensor and secondary sound source layouts, thus improving the active noise control optimization effect. Finally, based on the target sound field intensity, the sound power of the secondary sound field is determined. This sound power is used to control the secondary sound source to generate a secondary sound field to cancel out the external sound field, thereby improving the noise reduction effect. Furthermore, this application can also optimize the layout of high-order sensors by combining the noise reduction effect with the preset circular area of ​​the high-order sensor, thereby further improving the noise reduction effect. Attached Figure Description

[0035] Figure 1 This is a schematic flowchart illustrating an active noise control method according to an embodiment of this application;

[0036] Figure 2 This is a schematic diagram of the structure of an active noise control system shown in an embodiment of this application;

[0037] Figure 3 This is a schematic diagram of a pre-defined circular region of a high-order sensor as shown in an embodiment of this application;

[0038] Figure 4 This is a schematic diagram of the structure of an active noise control device shown in an embodiment of this application;

[0039] Figure 5 This is a schematic diagram of the structure of a computer device shown in an embodiment of this application. Detailed Implementation

[0040] The technical solutions of the embodiments of this application will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of this application, and not all embodiments. Based on the embodiments of this application, all other embodiments obtained by those of ordinary skill in the art without creative effort are within the scope of protection of this application.

[0041] Please refer to Figure 1 , Figure 1This is a flowchart illustrating an active noise control method provided in an embodiment of this application. The active noise control method of this application can be applied to computer devices, including but not limited to smartphones, laptops, tablets, desktop computers, physical servers, and cloud servers. Figure 1 As shown, the active noise control method of this embodiment includes steps S101 to S104, which are described in detail below:

[0042] Step S101: Obtain the sound pressure of the external sound field based on the high-order sensor, wherein the external sound field is the sound field outside the preset circular area of ​​the high-order sensor.

[0043] In this step, the external sound field is the primary sound field, and the higher-order sensor is a higher-order microphone. In practical applications, they can be used in array form to collect the sound pressure of the external sound field. For example... Figure 2 A schematic diagram of the active noise control system is shown, as follows: Figure 3 A schematic diagram of a preset circular region is shown. The sound field outside the preset circular region is the external sound field, and the sound field inside the preset circular region is the secondary sound field.

[0044] Step S102: Calculate the first sound field coefficient of the external sound field based on the sound pressure of the sound field.

[0045] In this step, the sound pressure of the sound field at any position within the preset circular region can be represented by the cylindrical harmonic function. Then, by combining the Graf addition theorem of the cylindrical Hankel function, the relationship between the sound pressure of the sound field and the sound field coefficient can be established, thereby calculating the first sound field coefficient of the external sound field.

[0046] It should be noted that, based on the above principles, layout optimization can be performed during the layout stage of high-order sensors and secondary sound sources, thereby obtaining the optimal layout result for application in actual noise reduction environments and improving the actual noise reduction effect.

[0047] In some embodiments, step S102 includes:

[0048] Using the least squares method, based on the sound pressure of the sound field, the first sound field coefficient of the external sound field is calculated. The least squares method is as follows:

[0049] β s =(H H H) -1 H H p;

[0050] Where, β s H represents the first sound field coefficient, and H is the matrix form of the second type of column Hankel function. H Let p represent the conjugate transpose of the matrix form, where p is the sound pressure level of the sound field.

[0051] In this embodiment, for the local circular region Ω of the j-th (j=1,2,…,J) higher-order microphone... (j) Establish a local polar coordinate system, such as Figure 3 As shown. Since the sound sources are all located in region Ω (j) Therefore, it can be treated as an internal sound field problem. The position x of any microphone on a high-order microphone... (j) The sound pressure at (R, φ) can be expressed using cylindrical harmonic functions as follows:

[0052]

[0053] Where k is the wave number, e represents the natural constant, and J n (·) denotes the nth-order cylindrical Bessel function. It is relative to the origin O of the local coordinate system j =(R j ,θ j The local sound field coefficients and cutoff order of )

[0054] To avoid spatial aliasing, the number of microphones is set to M = 2N + 1. The above formula can be written in matrix form as follows:

[0055] p j =Y j α j (2);

[0056] Where, p j Let α be the M×1 pickup sound pressure vector. j Y is the local sound field coefficient vector of (2N+1)×1. j It is an M×(2N+1) matrix (where M=2N+1, hence a square matrix), with elements Where μ = M, ν = n + N + 1.

[0057] Since the sound sources are all located in region Ω in Therefore, the region Ω is within. in The external sound field is considered as the external sound field. The sound pressure at any point, x = (r, θ), can be expressed as (time factor e) iωt ):

[0058]

[0059] in Let β represent the m-th order second-order column Hankel function. m (k) is the external sound field coefficient, and the cutoff order is . Expressed in matrix form:

[0060] p = Hβ (4);

[0061] Where p is the JM×1 pickup sound pressure vector, and β is (2N s The external sound field coefficient vector is +1)×1.

[0062] The local sound field coefficients of group J were obtained using a high-order microphone array. Using the Graf addition theorem of the Column Hankel function, it can be expressed as:

[0063]

[0064] The local sound field coefficients α of the j-th group can be obtained. j The relationship between the external sound field coefficient β and the external sound field coefficient β:

[0065]

[0066] (6) can be written in matrix form as follows:

[0067] α j =T j β (7);

[0068] Where T j (2N+1)×(2N) s A matrix with +1 elements Where μ = n + N + 1, v = m + N s +1.

[0069] By combining the local sound field coefficients of group J, and combining equations (2), (4) and (7), we can obtain:

[0070]

[0071] It should be noted that in order to ensure that equation (4) has a unique solution, the system of equations must be overdetermined, i.e., JM≥2N. s +1.

[0072] Therefore, by using a high-order microphone array to pick up the sound pressure p of the sound field, the sound field coefficients can be solved using the least squares method:

[0073] β=(H H H) -1 H H p (9);

[0074] If H is a non-singular matrix with a condition number greater than a preset number, then regularization needs to be applied.

[0075] Step S103: Based on the first preset relationship between the sound field intensity of the secondary sound field and the second sound field coefficient of the secondary sound field, and the second preset relationship between the first sound field coefficient and the second sound field coefficient, determine the target sound field intensity of the secondary sound field, wherein the secondary sound field is the sound field within the preset circular region.

[0076] In this step, optionally, the first preset relationship is: β s =Γq s Optionally, the second preset relationship is: β t =β p +β s .

[0077] In some embodiments, step S103 includes:

[0078] Based on the first preset relationship and the second preset relationship, a cost function for the secondary sound field is established;

[0079] The target sound field intensity of the secondary sound field is determined using the cost function, which is expressed as follows:

[0080]

[0081] Among them, J WD (q s ) represents the output value of the cost function. Γ represents the conjugate transpose of the target sound field intensity, and Γ is the sound transmission impedance coefficient matrix. H Let q be the conjugate transpose of the acoustic transmission impedance coefficient matrix. s β represents the target sound field intensity. p Indicates the second sound field coefficient. This represents the conjugate transpose of the second sound field coefficient.

[0082] In this embodiment, optionally, based on the first preset relationship and the second preset relationship, the cost function is generated by minimizing the sum of squares of the sound field coefficients of the total sound field, where the sum of squares of the sound field coefficients of the total sound field is... β represents the conjugate transpose of the total sound field coefficients. t This represents the total sound field coefficient.

[0083] Specifically, assume the intensity of the l-th secondary sound source is q. l (k), then the secondary sound pressure generated at any point x = (r, θ) can be expressed as:

[0084]

[0085] The superscript (s) indicates the secondary sound field, z l (r,θ,k) represents the acoustic transmission impedance between the l-th secondary sound source and the observation point x = (r,θ).

[0086] For a two-dimensional free field, we have:

[0087]

[0088] Where z0 = ρck / 4, x l =(r l ,θ l ) represents the location of the l-th secondary sound source.

[0089] In the frequency domain, the acoustic transmission impedance between two points is numerically equal to the sound pressure produced at one point when a unit sound source intensity (q = 1 m² / s) is applied at the other point. Therefore, here z l (r,θ,k) can be expressed using cylindrical harmonic functions as follows:

[0090]

[0091] Where γ ml (k) is the acoustic transmission impedance coefficient, which is independent of the location of the receiving point.

[0092] Apply the following addition theorem:

[0093]

[0094] Combining equations (11) and (12), we get:

[0095]

[0096] Combining equations (3), (12), and (14), the sound field coefficient β of the secondary sound field can be obtained. s and secondary sound source intensity q s The relationship between them:

[0097]

[0098] Written in matrix form:

[0099] β s =Γq s (16);

[0100] Where Γ is (2N) s The acoustic transmission impedance coefficient matrix is ​​+1)×L, with elements... Where μ=m+N s +1, ν=l.

[0101] After expanding the sound field in the cylindrical harmonic domain, minimize the sum of squares of the sound field coefficients of the total sound field after control. Substitute β t =β p +β s Combined with equation (16), the cost function can be expressed as:

[0102]

[0103] Due to ΓH If Γ is a symmetric positive definite matrix, then the unique optimal secondary sound source intensity vector can be found as follows:

[0104] q s =-(Γ H Γ) -1 Γ H β p (18).

[0105] Step S104: Determine the sound power of the secondary sound field based on the target sound field intensity. The sound power is used to control the secondary sound source to generate the secondary sound field to cancel the external sound field.

[0106] In this step, the sound power of the secondary sound field is calculated based on the target sound field intensity using a preset sound power calculation formula.

[0107] Optionally, the formula for calculating the sound power is:

[0108]

[0109] Among them, W t The sound power to be controlled for the secondary sound source. q represents the conjugate transpose of the target sound field intensity. s The target sound field intensity is represented by A and b, which are both preset sound field parameters, and W0 is the sound power of the secondary sound source before control.

[0110] Furthermore, the noise reduction amount is:

[0111]

[0112] Optionally, the optimal layout result can be further optimized based on the noise reduction amount, thereby further optimizing the noise reduction effect.

[0113] To implement the active noise control method corresponding to the above method embodiments, and to achieve the corresponding functions and technical effects. See also Figure 4 , Figure 4 This diagram illustrates a structural block diagram of an active noise control device according to an embodiment of this application. For ease of explanation, only the parts relevant to this embodiment are shown. The active noise control device provided in this embodiment includes:

[0114] The acquisition module 401 is used to acquire the sound pressure of the external sound field based on the high-order sensor, wherein the external sound field is the sound field outside the preset circular area of ​​the high-order sensor.

[0115] Calculation module 402 is used to calculate the first sound field coefficient of the external sound field based on the sound pressure of the sound field;

[0116] The first determining module 403 is used to determine the target sound field intensity of the secondary sound field based on a first preset relationship between the sound field intensity of the secondary sound field and the second sound field coefficient of the secondary sound field, and a second preset relationship between the first sound field coefficient and the second sound field coefficient, wherein the secondary sound field is a sound field within the preset circular region.

[0117] The second determining module 404 is used to determine the sound power of the secondary sound field based on the target sound field intensity. The sound power is used to control the secondary sound source to generate the secondary sound field to cancel the external sound field.

[0118] In some embodiments, the computing module 402 is configured to:

[0119] Using the least squares method, based on the sound pressure of the sound field, the first sound field coefficient of the external sound field is calculated. The least squares method is as follows:

[0120] β s =(H H H) -1 H H p;

[0121] Where, β s H represents the first sound field coefficient, and H is the matrix form of the second type of column Hankel function. H Let p represent the conjugate transpose of the matrix form, where p is the sound pressure level of the sound field.

[0122] In some embodiments, the first determining module 403 includes:

[0123] A unit is established to establish a cost function for the secondary sound field based on the first preset relationship and the second preset relationship.

[0124] A determining unit is used to determine the target sound field intensity of the secondary sound field using the cost function, wherein the expression of the cost function is:

[0125]

[0126] Among them, J WD (q s ) represents the output value of the cost function. Γ represents the conjugate transpose of the target sound field intensity, and Γ is the sound transmission impedance coefficient matrix. H Let q be the conjugate transpose of the acoustic transmission impedance coefficient matrix. s β represents the target sound field intensity. p Indicates the second sound field coefficient. This represents the conjugate transpose of the second sound field coefficient.

[0127] In some embodiments, the establishing unit is configured to:

[0128] Based on the first preset relationship and the second preset relationship, the cost function is generated by minimizing the sum of squares of the sound field coefficients of the total sound field, where the sum of squares of the sound field coefficients of the total sound field is: β represents the conjugate transpose of the total sound field coefficients. t This represents the total sound field coefficient.

[0129] In some embodiments, the first preset relationship is: β s =Γq s .

[0130] In some embodiments, the second preset relationship is: β t =β p +β s .

[0131] In some embodiments, the second determining module 404 is configured to:

[0132] Using a preset sound power calculation formula, the sound power of the secondary sound field is calculated based on the target sound field intensity. The sound power calculation formula is as follows:

[0133]

[0134] Among them, W t The sound power to be controlled for the secondary sound source. q represents the conjugate transpose of the target sound field intensity. s The target sound field intensity is represented by A and b, which are both preset sound field parameters, and W0 is the sound power of the secondary sound source before control.

[0135] The active noise control device described above can implement the active noise control method of the above method embodiments. The options in the above method embodiments are also applicable to this embodiment, and will not be detailed here. The remaining content of this application's embodiments can be referred to the content of the above method embodiments, and will not be repeated in this embodiment.

[0136] Figure 5 This is a schematic diagram of the structure of a computer device provided in an embodiment of this application. Figure 5 As shown, the computer device 5 of this embodiment includes: at least one processor 50 ( Figure 5 (Only one is shown) a processor, a memory 51, and a computer program 52 stored in the memory 51 and executable on the at least one processor 50, wherein the processor 50 executes the computer program 52 to implement the steps in any of the above method embodiments.

[0137] The computer device 5 may be a smartphone, tablet, desktop computer, or cloud server, among other computing devices. This computer device may include, but is not limited to, a processor 50 and a memory 51. Those skilled in the art will understand that... Figure 5 The computer device 5 is merely an example and does not constitute a limitation on the computer device 5. It may include more or fewer components than shown in the figure, or combine certain components, or different components, such as input / output devices, network access devices, etc.

[0138] The processor 50 may be a Central Processing Unit (CPU), or it may be other general-purpose processors, digital signal processors (DSPs), application-specific integrated circuits (ASICs), field-programmable gate arrays (FPGAs), or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. A general-purpose processor may be a microprocessor or any conventional processor.

[0139] In some embodiments, the memory 51 may be an internal storage unit of the computer device 5, such as a hard disk or memory of the computer device 5. In other embodiments, the memory 51 may be an external storage device of the computer device 5, such as a plug-in hard disk, smart media card (SMC), secure digital (SD) card, flash card, etc., equipped on the computer device 5. Furthermore, the memory 51 may include both internal and external storage units of the computer device 5. The memory 51 is used to store the operating system, applications, bootloader, data, and other programs, such as the program code of the computer program. The memory 51 can also be used to temporarily store data that has been output or will be output.

[0140] In addition, embodiments of this application also provide a computer-readable storage medium storing a computer program, which, when executed by a processor, implements the steps in any of the above method embodiments.

[0141] This application provides a computer program product that, when run on a computer device, enables the computer device to execute the steps described in the various method embodiments above.

[0142] In the several embodiments provided in this application, it will be understood that each block in the flowchart or block diagram may represent a module, segment, or portion of code containing one or more executable instructions for implementing a specified logical function. It should also be noted that in some alternative implementations, the functions marked in the blocks may occur in a different order than those shown in the figures. For example, two consecutive blocks may actually be executed substantially in parallel, and they may sometimes be executed in reverse order, depending on the functions involved.

[0143] If the aforementioned functions are implemented as software functional modules and sold or used as independent products, they can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of this application, in essence, or the part that contributes to the prior art, or a portion of the technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions to cause a computer device to execute all or part of the steps of the methods described in the various embodiments of this application. The aforementioned storage medium includes various media capable of storing program code, such as USB flash drives, portable hard drives, read-only memory (ROM), random access memory (RAM), magnetic disks, or optical disks.

[0144] The specific embodiments described above further illustrate the purpose, technical solution, and beneficial effects of this application. It should be understood that the above descriptions are merely specific embodiments of this application and are not intended to limit the scope of protection of this application. In particular, it should be noted that any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of this application should be included within the scope of protection of this application for those skilled in the art.

Claims

1. An active noise control method, characterized in that, include: The sound pressure of the external sound field is obtained based on a high-order sensor, wherein the external sound field is the sound field outside the preset circular region of the high-order sensor. Calculate the first sound field coefficient of the external sound field based on the sound pressure of the sound field; Based on a first preset relationship between the sound field intensity of the secondary sound field and the second sound field coefficient of the secondary sound field, and a second preset relationship between the first sound field coefficient and the second sound field coefficient, the target sound field intensity of the secondary sound field is determined, wherein the secondary sound field is the sound field within the preset circular region. Based on the target sound field intensity, the sound power of the secondary sound field is determined, and the sound power is used to control the secondary sound source to generate the secondary sound field to cancel the external sound field; The determination of the target sound field intensity of the secondary sound field, based on the first preset relationship between the sound field intensity of the secondary sound field and the second sound field coefficient of the secondary sound field, and the second preset relationship between the first sound field coefficient and the second sound field coefficient, includes: Based on the first preset relationship and the second preset relationship, a cost function for the secondary sound field is established; The target sound field intensity of the secondary sound field is determined using the cost function, which is expressed as follows: ； in, This represents the output value of the cost function. The conjugate transpose of the target sound field intensity. The acoustic transmission impedance coefficient matrix, This is the conjugate transpose of the acoustic transmission impedance coefficient matrix. Indicates the target sound field intensity. Indicates the second sound field coefficient. This represents the conjugate transpose of the second sound field coefficient; the first preset relationship is: The second preset relationship is: ; This represents the total sound field coefficient. This is the first sound field coefficient; The step of establishing the cost function for the secondary sound field based on the first preset relationship and the second preset relationship includes: Based on the first preset relationship and the second preset relationship, the cost function is generated by minimizing the sum of squares of the sound field coefficients of the total sound field, where the sum of squares of the sound field coefficients of the total sound field is: , This represents the conjugate transpose of the total sound field coefficients; The step of calculating the first sound field coefficient of the external sound field based on the sound pressure of the sound field includes: Using the least squares method, based on the sound pressure of the sound field, the first sound field coefficient of the external sound field is calculated. The least squares method is as follows: ； in, This is the matrix form of the second type of column Hankel function. This represents the conjugate transpose of the matrix form. This refers to the sound pressure level of the sound field.

2. The active noise control method as described in claim 1, characterized in that, Determining the acoustic power of the secondary sound field based on the target sound field intensity includes: Using a preset sound power calculation formula, the sound power of the secondary sound field is calculated based on the target sound field intensity. The sound power calculation formula is as follows: ； in, The sound power to be controlled for the secondary sound source. The conjugate transpose of the target sound field intensity. Indicates the target sound field intensity. and All are preset sound field parameters. The sound power of the secondary sound source before control.

3. An active noise control device, characterized in that, include: The acquisition module is used to acquire the sound pressure of the external sound field based on a high-order sensor, wherein the external sound field is the sound field outside the preset circular area of ​​the high-order sensor. The calculation module is used to calculate the first sound field coefficient of the external sound field based on the sound pressure of the sound field; The first determining module is used to determine the target sound field intensity of the secondary sound field based on a first preset relationship between the sound field intensity of the secondary sound field and the second sound field coefficient of the secondary sound field, and a second preset relationship between the first sound field coefficient and the second sound field coefficient, wherein the secondary sound field is a sound field within the preset circular region. The second determining module is used to determine the sound power of the secondary sound field based on the target sound field intensity. The sound power is used to control the secondary sound source to generate the secondary sound field to cancel the external sound field. The first determining module includes: A unit is established to establish a cost function for the secondary sound field based on the first preset relationship and the second preset relationship. A determining unit is used to determine the target sound field intensity of the secondary sound field using the cost function, wherein the expression of the cost function is: ； in, This represents the output value of the cost function. The conjugate transpose of the target sound field intensity. The acoustic transmission impedance coefficient matrix, This is the conjugate transpose of the acoustic transmission impedance coefficient matrix. Indicates the target sound field intensity. Indicates the second sound field coefficient. This represents the conjugate transpose of the second sound field coefficient; the first preset relationship is: The second preset relationship is: ; This represents the total sound field coefficient. This is the first sound field coefficient; The establishment unit is used for: Based on the first preset relationship and the second preset relationship, the cost function is generated by minimizing the sum of squares of the sound field coefficients of the total sound field, where the sum of squares of the sound field coefficients of the total sound field is: , This represents the conjugate transpose of the total sound field coefficients. Indicates the total sound field coefficient; The computing module is used for: Using the least squares method, based on the sound pressure of the sound field, the first sound field coefficient of the external sound field is calculated. The least squares method is as follows: ； in, The first sound field coefficient, This is the matrix form of the second type of column Hankel function. This represents the conjugate transpose of the matrix form. This refers to the sound pressure level of the sound field.

4. A computer device, characterized in that, It includes a processor and a memory, the memory being used to store a computer program that, when executed by the processor, implements the active noise control method as described in any one of claims 1 or 2.

5. A computer-readable storage medium, characterized in that, It stores a computer program that, when executed by a processor, implements the active noise control method as described in any one of claims 1 or 2.

Images