WO2023103468A1 - Method for determining transformer multi-sound-source noise equivalent model, terminal, and storage medium - Google Patents

Method for determining transformer multi-sound-source noise equivalent model, terminal, and storage medium Download PDF

Info

Publication number
WO2023103468A1
WO2023103468A1 PCT/CN2022/115386 CN2022115386W WO2023103468A1 WO 2023103468 A1 WO2023103468 A1 WO 2023103468A1 CN 2022115386 W CN2022115386 W CN 2022115386W WO 2023103468 A1 WO2023103468 A1 WO 2023103468A1
Authority
WO
WIPO (PCT)
Prior art keywords
preset
equivalent
octave band
sound
sound pressure
Prior art date
Application number
PCT/CN2022/115386
Other languages
French (fr)
Chinese (zh)
Inventor
吴鹏
胡源
刘长江
邢琳
王宁
张帅
何晓阳
段剑
邵华
李燕
赵彭辉
Original Assignee
国网河北省电力有限公司经济技术研究院
河北汇智电力工程设计有限公司
国网河北省电力有限公司
国家电网有限公司
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 国网河北省电力有限公司经济技术研究院, 河北汇智电力工程设计有限公司, 国网河北省电力有限公司, 国家电网有限公司 filed Critical 国网河北省电力有限公司经济技术研究院
Publication of WO2023103468A1 publication Critical patent/WO2023103468A1/en

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • HELECTRICITY
    • H01ELECTRIC ELEMENTS
    • H01FMAGNETS; INDUCTANCES; TRANSFORMERS; SELECTION OF MATERIALS FOR THEIR MAGNETIC PROPERTIES
    • H01F27/00Details of transformers or inductances, in general
    • H01F27/33Arrangements for noise damping
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/04Power grid distribution networks
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/10Noise analysis or noise optimisation

Definitions

  • the invention relates to the technical field of noise equivalent models, in particular to a method for determining an equivalent model of noise from multiple sound sources of a transformer, a terminal and a storage medium.
  • the transformer is the largest single device in the substation, and it is also the main source of noise in the substation. Constructing an accurate multi-source noise equivalent model of a transformer based on the acoustic equivalent source theory is of great significance for predicting the noise model of a substation.
  • near-field acoustic holography is usually used to construct the equivalent model of transformer multi-source noise.
  • this method is complex and computationally intensive.
  • Embodiments of the present invention provide a method for determining an equivalent model of transformer multi-sound source noise, a terminal and a storage medium, so as to solve the problems of complex calculation and large amount of calculation in the prior art.
  • an embodiment of the present invention provides a method for determining an equivalent model of transformer multi-source noise, including:
  • a univariate linear regression model corresponding to the preset octave band is constructed , and solve the univariate linear regression model corresponding to the preset octave band to obtain the sound pressure levels of multiple equivalent sound sources in the preset octave band;
  • the transformer multi-source noise equivalent model corresponding to the preset octave band is obtained.
  • the preset The univariate linear regression model corresponding to the octave band, and solve the univariate linear regression model corresponding to the preset octave band, to obtain the sound pressure level of multiple equivalent sound sources in the preset octave band, including:
  • a univariate linear regression model corresponding to the preset octave band is constructed ;
  • the univariate linear regression model corresponding to the preset octave band is solved.
  • Layer optimization obtains the sound pressure levels of multiple equivalent sound sources in preset octave bands.
  • the unit corresponding to the preset octave band The variable linear regression model is solved, and the sound pressure levels of multiple equivalent sound sources in the preset octave bands are obtained through double-layer optimization, including:
  • the first batch of random variables are generated as the input variables for the first optimization of the univariate linear regression model corresponding to the preset octave band;
  • the second batch of random variables is solved to obtain the sound pressure level of each equivalent sound source in the preset octave band.
  • the MSE loss function is:
  • MSE is the error value calculated by the MSE loss function
  • m is the number of predicted detection points
  • Lm j is the actual sound pressure level of the jth preset detection point in the preset octave band
  • Lw j is the jth preset The predicted sound pressure level of the detection point in the preset octave band
  • Lw ij Lp i -20lg(d ij )-0.001* ⁇ *d ij -11
  • n is the number of equivalent sound sources
  • Lw ij is the i-th equivalent sound source generated at the j-th preset detection point
  • Lp i is the sound pressure level of the i-th equivalent sound source in the preset octave band, which is an unknown quantity
  • d ij is the i-th equivalent sound source and the j-th
  • is the atmospheric absorption attenuation coefficient of noise during
  • the cross-entropy loss function is:
  • L is the value calculated by the cross-entropy loss function
  • P(Lp i ) is the value probability of Lp i in the normal distribution function.
  • the number of equivalent sound sources is 24;
  • the distribution of equivalent sound sources is as follows: the walls of the long box are arranged at equal intervals according to the 4*2 specification, and the walls of the short box are arranged at equal intervals according to the 2*2 specification.
  • an embodiment of the present invention provides a device for determining an equivalent model of transformer multi-source noise, including:
  • the first acquisition module is used to acquire the spatial coordinates of multiple preset detection points around the transformer and the actual sound pressure level in the preset octave band;
  • the second acquisition module is used to acquire the number of equivalent sound sources of the transformer and the spatial coordinates of each equivalent sound source;
  • the solving module is used to construct the corresponding sound pressure level of the preset octave band according to the spatial coordinates of multiple preset detection points, the actual sound pressure level in the preset octave band, the number of equivalent sound sources, and the spatial coordinates of each equivalent sound source.
  • a univariate linear regression model and solve the univariate linear regression model corresponding to the preset octave band to obtain the sound pressure levels of multiple equivalent sound sources in the preset octave band;
  • the model determination module is used to obtain the transformer multi-source noise corresponding to the preset octave band according to the number of equivalent sound sources, the sound pressure level of each equivalent sound source in the preset octave band, and the spatial coordinates of each equivalent sound source effective model.
  • the solving module is specifically used for:
  • a univariate linear regression model corresponding to the preset octave band is constructed ;
  • the univariate linear regression model corresponding to the preset octave band is solved.
  • Layer optimization obtains the sound pressure levels of multiple equivalent sound sources in preset octave bands.
  • an embodiment of the present invention provides a terminal, including a memory, a processor, and a computer program stored in the memory and operable on the processor.
  • the processor executes the computer program, the The steps of the method for determining the equivalent model of transformer multi-source noise as described in the first aspect or any possible implementation manner of the first aspect.
  • an embodiment of the present invention provides a computer-readable storage medium, the computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, it implements any of the above first aspect or the first aspect. Steps in the method for determining an equivalent model of transformer multi-source noise described in a possible implementation manner.
  • An embodiment of the present invention provides a method for determining an equivalent model of a transformer multi-sound source noise, a terminal, and a storage medium, by obtaining the spatial coordinates of multiple preset detection points around the transformer and the actual sound pressure level in a preset octave band; Obtain the number of equivalent sound sources of the transformer and the spatial coordinates of each equivalent sound source; according to the spatial coordinates of multiple preset detection points and the actual sound pressure level in the preset octave band, the number of equivalent sound sources and each The spatial coordinates of the effective sound source, construct the univariate linear regression model corresponding to the preset octave band, and solve the univariate linear regression model corresponding to the preset octave band, and obtain the sound of multiple equivalent sound sources in the preset octave band pressure level; according to the number of equivalent sound sources, the sound pressure level of each equivalent sound source and the spatial coordinates of each equivalent sound source, the transformer multi-source noise equivalent model corresponding to the preset
  • Fig. 1 is the implementation flowchart of the method for determining the transformer multi-source noise equivalent model provided by the embodiment of the present invention
  • Fig. 2 is a schematic diagram of a transformer multi-source noise equivalent model provided by an embodiment of the present invention
  • Fig. 3 is a plan view of equivalent source location distribution provided by an embodiment of the present invention.
  • FIG. 4 is a schematic structural diagram of a device for determining an equivalent model of transformer multi-source noise provided by an embodiment of the present invention
  • Fig. 5 is a schematic diagram of a terminal provided by an embodiment of the present invention.
  • FIG. 1 shows a flow chart of a method for determining an equivalent model of transformer multi-source noise provided by an embodiment of the present invention.
  • the subject of execution of the method for determining the equivalent model of transformer multi-sound source noise may be the terminal.
  • the methods for determining the equivalent model of the above-mentioned transformer multi-source noise include:
  • the actual typical substation can be selected, the planar distribution map and three-dimensional space model of the substation can be obtained, and a certain point on the ground of the substation is used as the coordinate origin to establish a spatial rectangular coordinate system; several preset detection points in the free space around the substation (marked as m) to measure the sound pressure level, and record the spatial coordinates (x Rj , y Rj , z Rj ) of the preset detection point and the actual sound pressure level Lm j corresponding to the preset octave band.
  • the detection point may also be referred to as a field point.
  • the coordinate origin is used as a reference point
  • the east direction is used as the positive direction of the X-axis
  • the north direction is used as the positive direction of the Y-axis
  • the upward direction is used as the positive direction of the Z-axis.
  • a real-time noise signal analyzer can be used to measure the actual sound pressure level at each preset detection point.
  • Noise real-time signal analyzer is a pocket real-time analyzer with digital signal processing technology, which can analyze the frequency spectrum and amplitude of noise, vibration or other electrical signals.
  • each sampling point When measuring the sound pressure level of the detection point, choose a sunny and windless weather, and record the temperature, air humidity and atmospheric pressure of the day.
  • each sampling point When using a noise real-time signal analyzer to measure the sound pressure level at the detection point, each sampling point should be continuously sampled for 30 to 60 seconds, 5 times in a row, and the average value should be taken. Spectrum analysis is performed on the sampling results of the sampling points, and the sound pressure levels of the eight octave bands of 63HZ, 125HZ, 250HZ, 500HZ, 1000HZ, 2000HZ, 4000HZ and 8000HZ are respectively extracted.
  • each preset detection point can be selected according to actual needs, and the preset octave band can be any one of the above eight octave bands.
  • the equivalent sound source may also be called an equivalent point sound source.
  • This embodiment can be based on the theory of multi-point equivalent sources, and according to the obtained spatial position parameters and geometric parameters of the transformer, use a correlation algorithm to obtain the geometric space information of the transformer equivalent point sound source, including the number of equivalent point sound sources and etc. Based on the spatial distribution of sound sources at effective points, etc., a transformer multi-point equivalent source noise radiation model (i.e. transformer multi-source noise equivalent model) is established, and the spatial position coordinates of n equivalent sources are marked as (x NSi , y NSi , zNSi ).
  • transformer multi-point equivalent source noise radiation model i.e. transformer multi-source noise equivalent model
  • the transformer may be a 500KV three-phase main transformer, and the geometric parameters of the transformer are 16.0m in length, 5.0m in width, and 5.0m in height.
  • the number of the above-mentioned equivalent sound sources is 24;
  • the distribution of equivalent sound sources is as follows: the walls of the long box are arranged at equal intervals according to the 4*2 specification, and the walls of the short box are arranged at equal intervals according to the 2*2 specification.
  • the equivalent model of transformer multi-sound source noise when constructing the equivalent model of transformer multi-sound source noise, the number of equivalent sound sources is set to 24, and the arrangement method is 4x2 equivalent sound sources arranged at equal intervals on the wall of the long box, etc. effect, the equivalent sound sources arranged at 2x2 equal intervals on the wall of the short box are equivalent, as shown in Figure 2, the cuboid in Figure 2 is the virtual model of the transformer, the black dots are the equivalent sound sources, and the equivalent sound sources are distributed in the There are no equivalent sound sources for the four sides, the upper and lower sides. It should be noted that, in order to clearly show the equivalent sound source in Fig. 2, Fig. 2 only draws the equivalent sound source of a long box wall and a short box wall, but in the actual model, the two opposite long box walls are There are equivalent sound sources, and the two opposite short box walls have equivalent sound sources.
  • the spatial position coordinates of the equivalent sound source are denoted as (x NSi , y NSi , z NSi ).
  • the sound pressure levels of multiple equivalent sound sources in the preset octave band can be obtained .
  • the sound pressure levels of the equivalent sound source of the source noise equivalent model may include sound pressure levels of eight octave bands. That is to say, different univariate linear regression models need to be constructed for different octave bands and solved to obtain the sound pressure level of the equivalent sound source in the octave band.
  • the above S103 may include:
  • a univariate linear regression model corresponding to the preset octave band is constructed ;
  • the univariate linear regression model corresponding to the preset octave band is solved.
  • Layer optimization obtains the sound pressure levels of multiple equivalent sound sources in preset octave bands.
  • the univariate linear regression model corresponding to the preset octave band is solved through double-layer optimization, which can improve the accuracy.
  • the univariate linearity corresponding to the preset octave band is solved, and the sound pressure levels of multiple equivalent sound sources in the preset octave bands are obtained through double-layer optimization, including:
  • the first batch of random variables are generated as the input variables for the first optimization of the univariate linear regression model corresponding to the preset octave band;
  • is the mathematical expectation value of the random variable with the preset number in front
  • ⁇ 2 is the variance value of the random variable with the preset number in the front.
  • the second batch of random variables is solved, and the second optimization is carried out, and finally the sound pressure level of each equivalent sound source in the preset octave band is obtained through the solution.
  • the Uniform function can be used to generate the first batch of random variables with uniform distribution.
  • Each random variable contains the sound pressure levels of 24 equivalent sound sources in preset octave bands, and each random variable can be different.
  • the global error between the predicted value and the actual measured value is calculated according to the MSE loss function, and a preset number of random variables with smaller errors are selected, and according to the value of the preset number of random variables with smaller errors
  • the mathematical expectation value and variance value are used to generate the second batch of random variables with a normal distribution.
  • the cross-entropy loss function is used to finally solve each equivalent sound source from the second batch of random variables in the preset Sound pressure level in octave bands. The sound pressure level of each equivalent sound source in the preset octave band is obtained through the second optimization solution.
  • the above MSE loss function is:
  • MSE is the error value calculated by the MSE loss function
  • m is the number of predicted detection points
  • Lm j is the actual sound pressure level of the jth preset detection point in the preset octave band
  • Lw j is the jth preset The predicted sound pressure level of the detection point in the preset octave band
  • Lw ij Lp i -20lg(d ij )-0.001* ⁇ *d ij -11
  • n is the number of equivalent sound sources
  • Lw ij is the i-th equivalent sound source generated at the j-th preset detection point
  • Lp i is the sound pressure level of the i-th equivalent sound source in the preset octave band, which is an unknown quantity
  • d ij is the i-th equivalent sound source and the j-th
  • is the atmospheric absorption attenuation coefficient of noise during
  • Figure 3 shows 4 field points (field point 1, field point 2, field point 3 and field point 4) and an equivalent source plane in which there are 4 equivalent acoustic source, d in Figure 3 represents the distance between one of the field points (detection points) and one of the equivalent sound sources.
  • the MSE loss function is a globally sensitive loss function, which is used to locate the candidate interval of Lp i in the second optimization.
  • the sound pressure level of the equivalent sound source in the space is the quantity to be sought, which is set as Lp i first.
  • the propagation distance d ii determines the geometric divergence attenuation and atmospheric absorption attenuation in the propagation process.
  • the above-mentioned cross-entropy loss function is:
  • L is the value calculated by the cross-entropy loss function
  • P(Lp i ) is the value probability of Lp i in the normal distribution function.
  • the above requirements are only the sound pressure level of the single octave band of the equivalent sound source. If you want to obtain the sound pressure level of the eight octave bands of 63HZ, 125HZ, 250HZ, 500HZ, 1000HZ, 2000HZ, 4000HZ and 8000HZ, you need to build an 8-time model.
  • the sound pressure level of the 8 octave bands of the equivalent sound source can be obtained by using the double-layer optimization method to perform 8 optimal solutions.
  • the equivalent model of transformer multi-source noise corresponding to the preset octave band can be obtained.
  • the final multi-source noise equivalent model of the transformer After solving the sound pressure levels of each equivalent sound source in eight octave bands, the final multi-source noise equivalent model of the transformer can be obtained.
  • This embodiment obtains the spatial coordinates of multiple preset detection points around the transformer and the actual sound pressure level in the preset octave band; obtains the number of equivalent sound sources of the transformer and the spatial coordinates of each equivalent sound source; according to multiple The spatial coordinates of a preset detection point and the actual sound pressure level in the preset octave band, the number of equivalent sound sources and the spatial coordinates of each equivalent sound source, construct a univariate linear regression model corresponding to the preset octave band, and Solve the univariate linear regression model corresponding to the preset octave band to obtain the sound pressure levels of multiple equivalent sound sources in the preset octave band; according to the number of equivalent sound sources, the sound pressure levels of each equivalent sound source and The spatial coordinates of each equivalent sound source obtain the equivalent model of transformer multi-source noise corresponding to the preset octave band, which can overcome the problems of complex calculation and large amount of calculation in near-field acoustic holography
  • a univariate linear regression model is constructed based on the noise radiation attenuation characteristics, using MSE and cross-entropy as loss functions, and the sound pressure level of the multi-sound source equivalent model of the transformer is obtained through double-layer optimization.
  • the equivalent process is convenient and simple, and the calculation amount is small, and the equivalent result is more accurate through double-layer optimization.
  • Fig. 4 shows a schematic structural diagram of a device for determining the equivalent model of transformer multi-source noise provided by an embodiment of the present invention.
  • Fig. 4 shows a schematic structural diagram of a device for determining the equivalent model of transformer multi-source noise provided by an embodiment of the present invention.
  • the details are as follows:
  • the device 30 for determining the equivalent model of transformer multi-source noise includes: a first acquisition module 31 , a second acquisition module 32 , a solution module 33 and a model determination module 34 .
  • the first acquisition module 31 is used to acquire the spatial coordinates of a plurality of preset detection points around the transformer and the actual sound pressure level in the preset octave band;
  • the second acquisition module 32 is used to acquire the quantity of equivalent sound sources of the transformer and the spatial coordinates of each equivalent sound source;
  • the solution module 33 is used to construct a preset octave band corresponding and solve the univariate linear regression model corresponding to the preset octave band to obtain the sound pressure levels of multiple equivalent sound sources in the preset octave band;
  • the model determination module 34 is used to obtain the transformer multi-source noise corresponding to the preset octave band according to the number of equivalent sound sources, the sound pressure level of each equivalent sound source in the preset octave band and the spatial coordinates of each equivalent sound source equivalent model.
  • the solving module 33 is specifically used for:
  • a univariate linear regression model corresponding to the preset octave band is constructed ;
  • the univariate linear regression model corresponding to the preset octave band is solved.
  • Layer optimization obtains the sound pressure levels of multiple equivalent sound sources in preset octave bands.
  • the solving module 33 is specifically used for:
  • the first batch of random variables are generated as the input variables for the first optimization of the univariate linear regression model corresponding to the preset octave band;
  • the first optimization is carried out, and the error value calculated by the MSE loss function is selected from the first batch of random variables in order from small to large, and the first preset number of random variables are ranked;
  • the second batch of random variables is generated by using the normal distribution method, as the second optimization of the univariate linear regression model corresponding to the preset octave band input variable;
  • the second optimization is carried out, and the sound pressure level of each equivalent sound source in the preset octave band is finally obtained from the second batch of random variables.
  • the MSE loss function is:
  • MSE is the error value calculated by the MSE loss function
  • m is the number of predicted detection points
  • Lm j is the actual sound pressure level of the jth preset detection point in the preset octave band
  • Lw j is the jth preset The predicted sound pressure level of the detection point in the preset octave band
  • Lw ij Lp i -20lg(d ij )-0.001* ⁇ *d ij -11
  • n is the number of equivalent sound sources
  • Lw ij is the i-th equivalent sound source generated at the j-th preset detection point
  • Lp i is the sound pressure level of the i-th equivalent sound source in the preset octave band, which is an unknown quantity
  • d ij is the i-th equivalent sound source and the j-th
  • is the atmospheric absorption attenuation coefficient of noise during
  • the cross-entropy loss function is:
  • L is the value calculated by the cross-entropy loss function
  • P(Lp i ) is the value probability of Lp i in the normal distribution function.
  • the number of equivalent sound sources is 24;
  • the distribution of equivalent sound sources is as follows: the walls of the long box are arranged at equal intervals according to the 4*2 specification, and the walls of the short box are arranged at equal intervals according to the 2*2 specification.
  • the above-mentioned device for determining the multi-source noise equivalent model of a transformer includes: a processor, wherein the processor is used to execute the above-mentioned program modules stored in the memory, including: a first acquisition module 31, a second acquisition module 32. A solution module 33 and a model determination module 34.
  • Fig. 5 is a schematic diagram of a terminal provided by an embodiment of the present invention.
  • the terminal 4 of this embodiment includes: a processor 40 , a memory 41 and a computer program 42 stored in the memory 41 and operable on the processor 40 .
  • the processor 40 executes the computer program 42, it implements the steps in the embodiment of the method for determining the equivalent model of each transformer multi-source noise, such as S101 to S104 shown in FIG. 1 .
  • the processor 40 executes the computer program 42, it realizes the functions of the modules/units in the above-mentioned device embodiments, such as the functions of the modules/units 31 to 34 shown in FIG. 4 .
  • the computer program 42 can be divided into one or more modules/units, and the one or more modules/units are stored in the memory 41 and executed by the processor 40 to complete this invention.
  • the one or more modules/units may be a series of computer program instruction segments capable of accomplishing specific functions, and the instruction segments are used to describe the execution process of the computer program 42 in the terminal 4 .
  • the computer program 42 may be divided into the modules/units 31 to 34 shown in FIG. 4 .
  • the terminal 4 can be a computing device such as a desktop computer, a notebook, a palmtop computer, or a cloud server.
  • the terminal 4 may include, but not limited to, a processor 40 and a memory 41 .
  • FIG. 5 is only an example of the terminal 4 and does not constitute a limitation to the terminal 4. It may include more or less components than those shown in the figure, or combine certain components, or different components, such as
  • the terminal may also include an input and output device, a network access device, a bus, and the like.
  • the so-called processor 40 can be a central processing unit (Central Processing Unit, CPU), and can also be other general-purpose processors, digital signal processors (Digital Signal Processor, DSP), application specific integrated circuits (Application Specific Integrated Circuit, ASIC), Field-Programmable Gate Array (Field-Programmable Gate Array, FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc.
  • a general-purpose processor may be a microprocessor, or the processor may be any conventional processor, or the like.
  • the memory 41 may be an internal storage unit of the terminal 4 , such as a hard disk or memory of the terminal 4 .
  • the memory 41 can also be an external storage device of the terminal 4, such as a plug-in hard disk equipped on the terminal 4, a smart memory card (Smart Media Card, SMC), a secure digital (Secure Digital, SD) card, Flash card (Flash Card), etc.
  • the memory 41 may also include both an internal storage unit of the terminal 4 and an external storage device.
  • the memory 41 is used to store the computer program and other programs and data required by the terminal.
  • the memory 41 can also be used to temporarily store data that has been output or will be output.
  • the disclosed device/terminal and method may be implemented in other ways.
  • the device/terminal embodiments described above are only illustrative.
  • the division of the modules or units is only a logical function division.
  • the mutual coupling or direct coupling or communication connection shown or discussed may be through some interfaces, and the indirect coupling or communication connection of devices or units may be in electrical, mechanical or other forms.
  • the units described as separate components may or may not be physically separated, and the components shown as units may or may not be physical units, that is, they may be located in one place, or may be distributed to multiple network units. Part or all of the units can be selected according to actual needs to achieve the purpose of the solution of this embodiment.
  • each functional unit in each embodiment of the present invention may be integrated into one processing unit, each unit may exist separately physically, or two or more units may be integrated into one unit.
  • the above-mentioned integrated units can be implemented in the form of hardware or in the form of software functional units.
  • the integrated module/unit is realized in the form of a software function unit and sold or used as an independent product, it can be stored in a computer-readable storage medium.
  • the present invention realizes all or part of the processes in the methods of the above embodiments, and can also be completed by instructing related hardware through a computer program.
  • the computer program can be stored in a computer-readable storage medium, and the computer When the program is executed by the processor, the steps in the above embodiments of the method for determining the equivalent model of the transformer multi-source noise can be realized.
  • the computer program includes computer program code, and the computer program code may be in the form of source code, object code, executable file or some intermediate form.
  • the computer-readable medium may include: any entity or device capable of carrying the computer program code, a recording medium, a USB flash drive, a removable hard disk, a magnetic disk, an optical disk, a computer memory, and a read-only memory (Read-Only Memory, ROM) , random access memory (Random Access Memory, RAM), electric carrier signal, telecommunication signal and software distribution medium, etc. It should be noted that the content contained in the computer-readable medium may be appropriately increased or decreased according to the requirements of legislation and patent practice in the jurisdiction. For example, in some jurisdictions, computer-readable media Excluding electrical carrier signals and telecommunication signals.

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Power Engineering (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

The present invention provides a method for determining a transformer multi-sound-source noise equivalent model, a terminal, and a storage medium. The method comprises: acquiring spatial coordinates of multiple preset detection points around a transformer and actual sound pressure levels in a preset octave band; acquiring the number of equivalent sound sources of the transformer and spatial coordinates of each equivalent sound source; according to the spatial coordinates of the multiple preset detection points and the actual sound pressure levels in the preset octave band, the number of the equivalent sound sources, and the spatial coordinates of each equivalent sound source, constructing a univariate linear regression model corresponding to the preset octave band, and solving the univariate linear regression model corresponding to the preset octave band to obtain sound pressure levels of multiple equivalent sound sources in the preset octave band; and according to the number of the equivalent sound sources, the sound pressure level of each equivalent sound source, and the spatial coordinates of each equivalent sound source, obtaining a transformer multi-sound-source noise equivalent model corresponding to the preset octave band. The process of determining a transformer multi-sound-source noise equivalent model of the present invention is convenient and simple, and the amount of calculation is small.

Description

变压器多声源噪声等效模型的确定方法、终端及存储介质Determination method, terminal and storage medium of transformer multi-source noise equivalent model 技术领域technical field
本发明涉及噪声等效模型技术领域,尤其涉及一种变压器多声源噪声等效模型的确定方法、终端及存储介质。The invention relates to the technical field of noise equivalent models, in particular to a method for determining an equivalent model of noise from multiple sound sources of a transformer, a terminal and a storage medium.
背景技术Background technique
随着变电站选址越来越接近居民区,变电站带来的噪声污染的问题也越来越受到人们的重视。变压器是变电站内最大的单体设备,也是变电站内最主要的噪声来源。根据声学等效源理论构建准确的变压器多声源噪声等效模型,对于进行变电站噪声模型预测具有十分重要的意义。As the location of substations is getting closer to residential areas, the problem of noise pollution caused by substations has attracted more and more attention. The transformer is the largest single device in the substation, and it is also the main source of noise in the substation. Constructing an accurate multi-source noise equivalent model of a transformer based on the acoustic equivalent source theory is of great significance for predicting the noise model of a substation.
目前,通常采用近场声全息技术来构建变压器多声源噪声等效模型,然而,这种方法计算复杂,计算量大。At present, near-field acoustic holography is usually used to construct the equivalent model of transformer multi-source noise. However, this method is complex and computationally intensive.
发明内容Contents of the invention
本发明实施例提供了一种变压器多声源噪声等效模型的确定方法、终端及存储介质,以解决现有技术计算复杂,计算量大的问题。Embodiments of the present invention provide a method for determining an equivalent model of transformer multi-sound source noise, a terminal and a storage medium, so as to solve the problems of complex calculation and large amount of calculation in the prior art.
第一方面,本发明实施例提供了一种变压器多声源噪声等效模型的确定方法,包括:In the first aspect, an embodiment of the present invention provides a method for determining an equivalent model of transformer multi-source noise, including:
获取变压器周围的多个预设检测点的空间坐标与在预设倍频带的实际声压级;Obtain the spatial coordinates of multiple preset detection points around the transformer and the actual sound pressure level in the preset octave band;
获取变压器的等效声源的数量和各个等效声源的空间坐标;Obtain the number of equivalent sound sources of the transformer and the spatial coordinates of each equivalent sound source;
根据多个预设检测点的空间坐标与在预设倍频带的实际声压级、等效声源的数量以及各个等效声源的空间坐标,构建预设倍频带对应的单变量线性回归模型,并对预设倍频带对应的单变量线性回归模型进行求解,得到多个等效声源在预设倍频带的声压级;According to the spatial coordinates of multiple preset detection points and the actual sound pressure level in the preset octave band, the number of equivalent sound sources and the spatial coordinates of each equivalent sound source, a univariate linear regression model corresponding to the preset octave band is constructed , and solve the univariate linear regression model corresponding to the preset octave band to obtain the sound pressure levels of multiple equivalent sound sources in the preset octave band;
根据等效声源的数量、各个等效声源在预设倍频带的声压级和各个等效声源的空间坐标得到预设倍频带对应的变压器多声源噪声等效模型。According to the number of equivalent sound sources, the sound pressure level of each equivalent sound source in the preset octave band and the spatial coordinates of each equivalent sound source, the transformer multi-source noise equivalent model corresponding to the preset octave band is obtained.
在一种可能的实现方式中,根据多个预设检测点的空间坐标与在预设倍频带的实际声压级、等效声源的数量以及各个等效声源的空间坐标,构建预设倍频带对应的单变量线性回归模型,并对预设倍频带对应的单变量线性回归模型进行求解,得到多个等效声源在预设倍 频带的声压级,包括:In a possible implementation manner, the preset The univariate linear regression model corresponding to the octave band, and solve the univariate linear regression model corresponding to the preset octave band, to obtain the sound pressure level of multiple equivalent sound sources in the preset octave band, including:
根据多个预设检测点的空间坐标与在预设倍频带的实际声压级、等效声源的数量以及各个等效声源的空间坐标,构建预设倍频带对应的单变量线性回归模型;According to the spatial coordinates of multiple preset detection points and the actual sound pressure level in the preset octave band, the number of equivalent sound sources and the spatial coordinates of each equivalent sound source, a univariate linear regression model corresponding to the preset octave band is constructed ;
根据各个预设检测点在预设倍频带的预测声压级以及各个预设检测点在预设倍频带的实际声压级,对预设倍频带对应的单变量线性回归模型进行求解,通过双层优化得到多个等效声源在预设倍频带的声压级。According to the predicted sound pressure level of each preset detection point in the preset octave band and the actual sound pressure level of each preset detection point in the preset octave band, the univariate linear regression model corresponding to the preset octave band is solved. Layer optimization obtains the sound pressure levels of multiple equivalent sound sources in preset octave bands.
在一种可能的实现方式中,根据各个预设检测点在预设倍频带的预测声压级以及各个预设检测点在预设倍频带的实际声压级,对预设倍频带对应的单变量线性回归模型进行求解,通过双层优化得到多个等效声源在预设倍频带的声压级,包括:In a possible implementation, according to the predicted sound pressure level of each preset detection point in the preset octave band and the actual sound pressure level of each preset detection point in the preset octave band, the unit corresponding to the preset octave band The variable linear regression model is solved, and the sound pressure levels of multiple equivalent sound sources in the preset octave bands are obtained through double-layer optimization, including:
采用平均分布的方法,生成第一批随机变量,作为预设倍频带对应的单变量线性回归模型第一次优化时的输入变量;Using the method of average distribution, the first batch of random variables are generated as the input variables for the first optimization of the univariate linear regression model corresponding to the preset octave band;
根据MSE损失函数计算第一批随机变量中各个随机变量的误差值,从按照误差值由小到大排序的所述第一批随机变量中按排序选取预设数量的随机变量;Calculate the error value of each random variable in the first batch of random variables according to the MSE loss function, and select a preset number of random variables in order from the first batch of random variables sorted according to the error value from small to large;
采用正态分布的方法对每一选取到的随机变量的数学期望值和方差值进行计算,生成第二批随机变量,作为预设倍频带对应的单变量线性回归模型第二次优化时的输入变量;Calculate the mathematical expectation and variance of each selected random variable using the normal distribution method to generate the second batch of random variables as the input for the second optimization of the univariate linear regression model corresponding to the preset octave band variable;
根据交叉熵损失函数对第二批随机变量进行求解,得到各个等效声源在预设倍频带的声压级。According to the cross-entropy loss function, the second batch of random variables is solved to obtain the sound pressure level of each equivalent sound source in the preset octave band.
在一种可能的实现方式中,MSE损失函数为:In one possible implementation, the MSE loss function is:
Figure PCTCN2022115386-appb-000001
Figure PCTCN2022115386-appb-000001
其中,MSE为MSE损失函数计算得到的误差值;m为预测检测点的数量;Lm j为第j个预设检测点在预设倍频带的实际声压级;Lw j为第j个预设检测点在预设倍频带的预测声压级,
Figure PCTCN2022115386-appb-000002
Lw ij=Lp i-20lg(d ij)-0.001*α*d ij-11;n为等效声源的数量;Lw ij为第i个等效声源在第j个预设检测点产生的在预设倍频带的声压级;Lp i为第i个等效声源在预设倍频带待求的声压级,为未知量;d ij为第i个等效声源与第j个预设检测点之间的距离;α为噪声在传播过程中的大气吸收衰减系数。
Among them, MSE is the error value calculated by the MSE loss function; m is the number of predicted detection points; Lm j is the actual sound pressure level of the jth preset detection point in the preset octave band; Lw j is the jth preset The predicted sound pressure level of the detection point in the preset octave band,
Figure PCTCN2022115386-appb-000002
Lw ij =Lp i -20lg(d ij )-0.001*α*d ij -11; n is the number of equivalent sound sources; Lw ij is the i-th equivalent sound source generated at the j-th preset detection point The sound pressure level in the preset octave band; Lp i is the sound pressure level of the i-th equivalent sound source in the preset octave band, which is an unknown quantity; d ij is the i-th equivalent sound source and the j-th The distance between preset detection points; α is the atmospheric absorption attenuation coefficient of noise during propagation.
在一种可能的实现方式中,交叉熵损失函数为:In one possible implementation, the cross-entropy loss function is:
Figure PCTCN2022115386-appb-000003
Figure PCTCN2022115386-appb-000003
其中,L为交叉熵损失函数计算得到的值;P(Lp i)为Lp i在正态分布函数中的取值概率。 Among them, L is the value calculated by the cross-entropy loss function; P(Lp i ) is the value probability of Lp i in the normal distribution function.
在一种可能的实现方式中,等效声源的数量为24个;In a possible implementation manner, the number of equivalent sound sources is 24;
等效声源的分布方式为:长箱壁面按照4*2规格等间距布置,短箱壁面按照2*2规格等间距布置。The distribution of equivalent sound sources is as follows: the walls of the long box are arranged at equal intervals according to the 4*2 specification, and the walls of the short box are arranged at equal intervals according to the 2*2 specification.
第二方面,本发明实施例提供了变压器多声源噪声等效模型的确定装置,包括:In the second aspect, an embodiment of the present invention provides a device for determining an equivalent model of transformer multi-source noise, including:
第一获取模块,用于获取变压器周围的多个预设检测点的空间坐标与在预设倍频带的实际声压级;The first acquisition module is used to acquire the spatial coordinates of multiple preset detection points around the transformer and the actual sound pressure level in the preset octave band;
第二获取模块,用于获取变压器的等效声源的数量和各个等效声源的空间坐标;The second acquisition module is used to acquire the number of equivalent sound sources of the transformer and the spatial coordinates of each equivalent sound source;
求解模块,用于根据多个预设检测点的空间坐标与在预设倍频带的实际声压级、等效声源的数量以及各个等效声源的空间坐标,构建预设倍频带对应的单变量线性回归模型,并对预设倍频带对应的单变量线性回归模型进行求解,得到多个等效声源在预设倍频带的声压级;The solving module is used to construct the corresponding sound pressure level of the preset octave band according to the spatial coordinates of multiple preset detection points, the actual sound pressure level in the preset octave band, the number of equivalent sound sources, and the spatial coordinates of each equivalent sound source. A univariate linear regression model, and solve the univariate linear regression model corresponding to the preset octave band to obtain the sound pressure levels of multiple equivalent sound sources in the preset octave band;
模型确定模块,用于根据等效声源的数量、各个等效声源在预设倍频带的声压级和各个等效声源的空间坐标得到预设倍频带对应的变压器多声源噪声等效模型。The model determination module is used to obtain the transformer multi-source noise corresponding to the preset octave band according to the number of equivalent sound sources, the sound pressure level of each equivalent sound source in the preset octave band, and the spatial coordinates of each equivalent sound source effective model.
在一种可能的实现方式中,求解模块具体用于:In a possible implementation, the solving module is specifically used for:
根据多个预设检测点的空间坐标与在预设倍频带的实际声压级、等效声源的数量以及各个等效声源的空间坐标,构建预设倍频带对应的单变量线性回归模型;According to the spatial coordinates of multiple preset detection points and the actual sound pressure level in the preset octave band, the number of equivalent sound sources and the spatial coordinates of each equivalent sound source, a univariate linear regression model corresponding to the preset octave band is constructed ;
根据各个预设检测点在预设倍频带的预测声压级以及各个预设检测点在预设倍频带的实际声压级,对预设倍频带对应的单变量线性回归模型进行求解,通过双层优化得到多个等效声源在预设倍频带的声压级。According to the predicted sound pressure level of each preset detection point in the preset octave band and the actual sound pressure level of each preset detection point in the preset octave band, the univariate linear regression model corresponding to the preset octave band is solved. Layer optimization obtains the sound pressure levels of multiple equivalent sound sources in preset octave bands.
第三方面,本发明实施例提供了一种终端,包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现如上第一方面或第一方面的任一种可能的实现方式所述的变压器多声源噪声等效模型的确定方法的步骤。In a third aspect, an embodiment of the present invention provides a terminal, including a memory, a processor, and a computer program stored in the memory and operable on the processor. When the processor executes the computer program, the The steps of the method for determining the equivalent model of transformer multi-source noise as described in the first aspect or any possible implementation manner of the first aspect.
第四方面,本发明实施例提供了一种计算机可读存储介质,所述计算机可读存储介质存储有计算机程序,所述计算机程序被处理器执行时实现如上第一方面或第一方面的任一种可能的实现方式所述的变压器多声源噪声等效模型的确定方法的步骤。In a fourth aspect, an embodiment of the present invention provides a computer-readable storage medium, the computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, it implements any of the above first aspect or the first aspect. Steps in the method for determining an equivalent model of transformer multi-source noise described in a possible implementation manner.
本发明实施例提供一种变压器多声源噪声等效模型的确定方法、终端及存储介质,通过获取变压器周围的多个预设检测点的空间坐标与在预设倍频带的实际声压级;获取变压器 的等效声源的数量和各个等效声源的空间坐标;根据多个预设检测点的空间坐标与在预设倍频带的实际声压级、等效声源的数量以及各个等效声源的空间坐标,构建预设倍频带对应的单变量线性回归模型,并对预设倍频带对应的单变量线性回归模型进行求解,得到多个等效声源在预设倍频带的声压级;根据等效声源的数量、各个等效声源的声压级和各个等效声源的空间坐标得到预设倍频带对应的变压器多声源噪声等效模型,能够克服近场声全息技术中的计算复杂和计算量大的问题,只需测量变压器附近少量检测点的声压级,并基于单变量线性回归的方法获得变压器多声源噪声等效模型,等效过程方便简单,运算量小。An embodiment of the present invention provides a method for determining an equivalent model of a transformer multi-sound source noise, a terminal, and a storage medium, by obtaining the spatial coordinates of multiple preset detection points around the transformer and the actual sound pressure level in a preset octave band; Obtain the number of equivalent sound sources of the transformer and the spatial coordinates of each equivalent sound source; according to the spatial coordinates of multiple preset detection points and the actual sound pressure level in the preset octave band, the number of equivalent sound sources and each The spatial coordinates of the effective sound source, construct the univariate linear regression model corresponding to the preset octave band, and solve the univariate linear regression model corresponding to the preset octave band, and obtain the sound of multiple equivalent sound sources in the preset octave band pressure level; according to the number of equivalent sound sources, the sound pressure level of each equivalent sound source and the spatial coordinates of each equivalent sound source, the transformer multi-source noise equivalent model corresponding to the preset octave band can be obtained, which can overcome the near-field sound For the problems of complex calculation and large amount of calculation in holographic technology, it is only necessary to measure the sound pressure level of a small number of detection points near the transformer, and obtain the equivalent model of transformer multi-sound source noise based on the method of univariate linear regression. The equivalent process is convenient and simple. The amount of calculation is small.
附图说明Description of drawings
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the following will briefly introduce the accompanying drawings that need to be used in the descriptions of the embodiments or the prior art. Obviously, the accompanying drawings in the following description are only of the present invention. For some embodiments, those of ordinary skill in the art can also obtain other drawings based on these drawings without paying creative efforts.
图1是本发明实施例提供的变压器多声源噪声等效模型的确定方法的实现流程图;Fig. 1 is the implementation flowchart of the method for determining the transformer multi-source noise equivalent model provided by the embodiment of the present invention;
图2是本发明实施例提供的变压器多声源噪声等效模型的示意图;Fig. 2 is a schematic diagram of a transformer multi-source noise equivalent model provided by an embodiment of the present invention;
图3是本发明实施例提供的等效源位置分布平面图;Fig. 3 is a plan view of equivalent source location distribution provided by an embodiment of the present invention;
图4是本发明实施例提供的变压器多声源噪声等效模型的确定装置的结构示意图;4 is a schematic structural diagram of a device for determining an equivalent model of transformer multi-source noise provided by an embodiment of the present invention;
图5是本发明实施例提供的终端的示意图。Fig. 5 is a schematic diagram of a terminal provided by an embodiment of the present invention.
具体实施方式Detailed ways
以下描述中,为了说明而不是为了限定,提出了诸如特定系统结构、技术之类的具体细节,以便透彻理解本发明实施例。然而,本领域的技术人员应当清楚,在没有这些具体细节的其它实施例中也可以实现本发明。在其它情况中,省略对众所周知的系统、装置、电路以及方法的详细说明,以免不必要的细节妨碍本发明的描述。In the following description, specific details such as specific system structures and technologies are presented for the purpose of illustration rather than limitation, so as to thoroughly understand the embodiments of the present invention. It will be apparent, however, to one skilled in the art that the invention may be practiced in other embodiments without these specific details. In other instances, detailed descriptions of well-known systems, devices, circuits, and methods are omitted so as not to obscure the description of the present invention with unnecessary detail.
为使本发明的目的、技术方案和优点更加清楚,下面将结合附图通过具体实施例来进行说明。In order to make the purpose, technical solution and advantages of the present invention clearer, specific embodiments will be described below in conjunction with the accompanying drawings.
参见图1,其示出了本发明实施例提供的变压器多声源噪声等效模型的确定方法的实现流程图。其中,变压器多声源噪声等效模型的确定方法的执行主体可以是终端。Referring to FIG. 1 , it shows a flow chart of a method for determining an equivalent model of transformer multi-source noise provided by an embodiment of the present invention. Wherein, the subject of execution of the method for determining the equivalent model of transformer multi-sound source noise may be the terminal.
参见图1,上述变压器多声源噪声等效模型的确定方法包括:Referring to Figure 1, the methods for determining the equivalent model of the above-mentioned transformer multi-source noise include:
在S101中,获取变压器周围的多个预设检测点的空间坐标与在预设倍频带的实际声压 级。In S101, the spatial coordinates of multiple preset detection points around the transformer and the actual sound pressure level in the preset octave band are obtained.
在本实施例中,可以选择实际典型变电站,获取变电站的平面分布图与三维空间模型,以变电站地面上某一点为坐标原点,建立空间直角坐标系;对变电站周围自由空间中若干预设检测点(记为m个)进行声压级测量,并记录预设检测点的空间坐标(x Rj,y Rj,z Rj)与预设倍频带对应的实际声压级大小Lm j。其中,检测点也可以称为场点。 In this embodiment, the actual typical substation can be selected, the planar distribution map and three-dimensional space model of the substation can be obtained, and a certain point on the ground of the substation is used as the coordinate origin to establish a spatial rectangular coordinate system; several preset detection points in the free space around the substation (marked as m) to measure the sound pressure level, and record the spatial coordinates (x Rj , y Rj , z Rj ) of the preset detection point and the actual sound pressure level Lm j corresponding to the preset octave band. Wherein, the detection point may also be referred to as a field point.
在一种可能的实现方式中,以坐标原点作为参考点,以东向作为X轴的正方向,以北向作为Y轴的正方向,向上为Z轴的正方向。根据检测点与坐标原点之间的空间位置关系,建立检测点的空间位置坐标,记为(x Rj,y Rj,z Rj),单位采用“米”。 In a possible implementation manner, the coordinate origin is used as a reference point, the east direction is used as the positive direction of the X-axis, the north direction is used as the positive direction of the Y-axis, and the upward direction is used as the positive direction of the Z-axis. According to the spatial position relationship between the detection point and the coordinate origin, the spatial position coordinates of the detection point are established, which are recorded as (x Rj , y Rj , z Rj ), and the unit is "meter".
本实施例可以采用噪声实时信号分析仪对各个预设检测点处的实际声压级进行测量。噪声实时信号分析仪是一种数字信号处理技术的袖珍式实时分析仪,它可以对噪声、振动或其它电信号进行频谱及幅值分析。In this embodiment, a real-time noise signal analyzer can be used to measure the actual sound pressure level at each preset detection point. Noise real-time signal analyzer is a pocket real-time analyzer with digital signal processing technology, which can analyze the frequency spectrum and amplitude of noise, vibration or other electrical signals.
在测量检测点声压级时应选取晴朗、无风的天气,并记录下当天的温度、空气湿度以及大气压强。用噪声实时信号分析仪对检测点处的声压级进行测量时,每个采样点应连续采样30~60秒,连续采5次,取其平均值。对采样点的采样结果进行频谱分析,分别提取63HZ、125HZ、250HZ、500HZ、1000HZ、2000HZ、4000HZ和8000HZ八个倍频带的声压级。When measuring the sound pressure level of the detection point, choose a sunny and windless weather, and record the temperature, air humidity and atmospheric pressure of the day. When using a noise real-time signal analyzer to measure the sound pressure level at the detection point, each sampling point should be continuously sampled for 30 to 60 seconds, 5 times in a row, and the average value should be taken. Spectrum analysis is performed on the sampling results of the sampling points, and the sound pressure levels of the eight octave bands of 63HZ, 125HZ, 250HZ, 500HZ, 1000HZ, 2000HZ, 4000HZ and 8000HZ are respectively extracted.
各个预设检测点的位置可以根据实际需求选取,预设倍频带可以是上述八个倍频带中的任意一个。The positions of each preset detection point can be selected according to actual needs, and the preset octave band can be any one of the above eight octave bands.
在S102中,获取变压器的等效声源的数量和各个等效声源的空间坐标。In S102, the number of equivalent sound sources of the transformer and the spatial coordinates of each equivalent sound source are acquired.
其中,等效声源也可以称为等效点声源。Wherein, the equivalent sound source may also be called an equivalent point sound source.
本实施例可以基于多点等效源理论,根据获得的变压器的空间位置参数与几何参数,利用相关算法获得变压器等效点声源的几何空间信息,包括等效点声源的个数以及等效点声源的空间分布情况等,据此建立变压器多点等效源噪声辐射模型(即变压器多声源噪声等效模型),n个等效源的空间位置坐标记 为(x NSi,y NSi,z NSi)。 This embodiment can be based on the theory of multi-point equivalent sources, and according to the obtained spatial position parameters and geometric parameters of the transformer, use a correlation algorithm to obtain the geometric space information of the transformer equivalent point sound source, including the number of equivalent point sound sources and etc. Based on the spatial distribution of sound sources at effective points, etc., a transformer multi-point equivalent source noise radiation model (i.e. transformer multi-source noise equivalent model) is established, and the spatial position coordinates of n equivalent sources are marked as (x NSi , y NSi , zNSi ).
其中,变压器可以是500KV三相主变压器,变压器的几何参数为长16.0m,宽5.0m,高5.0m。Wherein, the transformer may be a 500KV three-phase main transformer, and the geometric parameters of the transformer are 16.0m in length, 5.0m in width, and 5.0m in height.
在一些实施例中,上述等效声源的数量为24个;In some embodiments, the number of the above-mentioned equivalent sound sources is 24;
等效声源的分布方式为:长箱壁面按照4*2规格等间距布置,短箱壁面按照2*2规格等间距布置。The distribution of equivalent sound sources is as follows: the walls of the long box are arranged at equal intervals according to the 4*2 specification, and the walls of the short box are arranged at equal intervals according to the 2*2 specification.
基于上述变压器的几何参数信息,在构建变压器多声源噪声等效模型时,等效声源的个数设置为24个,布置方式为长箱壁面按4x2个等间距布置的等效声源等效,短箱壁面按2x2个等间距布置的等效声源等效,如图2所示,图2中长方体为变压器虚拟模型,黑点为等效声源,等效声源分布在长方体的四个侧面,上下两个面是没有等效声源的。需要说明的是,为了使图2清楚表明等效声源,图2仅画出了一个长箱壁面和一个短箱壁面的等效声源,而实际模型中,两个相对的长箱壁面均具有等效声源,两个相对的短箱壁面均具有等效声源。Based on the above-mentioned geometric parameter information of the transformer, when constructing the equivalent model of transformer multi-sound source noise, the number of equivalent sound sources is set to 24, and the arrangement method is 4x2 equivalent sound sources arranged at equal intervals on the wall of the long box, etc. effect, the equivalent sound sources arranged at 2x2 equal intervals on the wall of the short box are equivalent, as shown in Figure 2, the cuboid in Figure 2 is the virtual model of the transformer, the black dots are the equivalent sound sources, and the equivalent sound sources are distributed in the There are no equivalent sound sources for the four sides, the upper and lower sides. It should be noted that, in order to clearly show the equivalent sound source in Fig. 2, Fig. 2 only draws the equivalent sound source of a long box wall and a short box wall, but in the actual model, the two opposite long box walls are There are equivalent sound sources, and the two opposite short box walls have equivalent sound sources.
根据等效声源的空间位置,等效声源的空间位置坐标,记为(x NSi,y NSi,z NSi)。 According to the spatial position of the equivalent sound source, the spatial position coordinates of the equivalent sound source are denoted as (x NSi , y NSi , z NSi ).
在S103中,根据多个预设检测点的空间坐标与在预设倍频带的实际声压级、等效声源的数量以及各个等效声源的空间坐标,构建预设倍频带对应的单变量线性回归模型,并对预设倍频带对应的单变量线性回归模型进行求解,得到多个等效声源在预设倍频带的声压级。In S103, according to the spatial coordinates of a plurality of preset detection points and the actual sound pressure level in the preset octave band, the number of equivalent sound sources and the spatial coordinates of each equivalent sound source, a single corresponding to the preset octave band is constructed. Variable linear regression model, and solve the univariate linear regression model corresponding to the preset octave band to obtain the sound pressure levels of multiple equivalent sound sources in the preset octave band.
本实施例通过构建预设倍频带对应的单变量线性回归模型,并对预设倍频带对应的单变量线性回归模型进行求解,可以得到多个等效声源在预设倍频带的声压级。In this embodiment, by constructing the univariate linear regression model corresponding to the preset octave band and solving the univariate linear regression model corresponding to the preset octave band, the sound pressure levels of multiple equivalent sound sources in the preset octave band can be obtained .
在本实施例中,若想得到等效声源的八个倍频带的声压级,可以将八个倍频带分别作为预设倍频带,执行上述S101-S103,最后再通过S104得到的变压器多声源噪声等效模型的等效声源的声压级可以包括八个倍频带的声压级。也就是说,不同倍频带需构建不同的单变量线性回归模型,并进行求解,得到等 效声源在该倍频带的声压级。In this embodiment, if you want to obtain the sound pressure levels of the eight octave bands of the equivalent sound source, you can use the eight octave bands as the preset octave bands, perform the above S101-S103, and finally obtain the transformer multi-acoustic sound through S104. The sound pressure levels of the equivalent sound source of the source noise equivalent model may include sound pressure levels of eight octave bands. That is to say, different univariate linear regression models need to be constructed for different octave bands and solved to obtain the sound pressure level of the equivalent sound source in the octave band.
在一些实施例中,上述S103可以包括:In some embodiments, the above S103 may include:
根据多个预设检测点的空间坐标与在预设倍频带的实际声压级、等效声源的数量以及各个等效声源的空间坐标,构建预设倍频带对应的单变量线性回归模型;According to the spatial coordinates of multiple preset detection points and the actual sound pressure level in the preset octave band, the number of equivalent sound sources and the spatial coordinates of each equivalent sound source, a univariate linear regression model corresponding to the preset octave band is constructed ;
根据各个预设检测点在预设倍频带的预测声压级以及各个预设检测点在预设倍频带的实际声压级,对预设倍频带对应的单变量线性回归模型进行求解,通过双层优化得到多个等效声源在预设倍频带的声压级。According to the predicted sound pressure level of each preset detection point in the preset octave band and the actual sound pressure level of each preset detection point in the preset octave band, the univariate linear regression model corresponding to the preset octave band is solved. Layer optimization obtains the sound pressure levels of multiple equivalent sound sources in preset octave bands.
本实施例通过双层优化对预设倍频带对应的单变量线性回归模型进行求解,能够提高准确性。In this embodiment, the univariate linear regression model corresponding to the preset octave band is solved through double-layer optimization, which can improve the accuracy.
在一些实施例中,上述根据各个预设检测点在预设倍频带的预测声压级以及各个预设检测点在预设倍频带的实际声压级,对预设倍频带对应的单变量线性回归模型进行求解,通过双层优化得到多个等效声源在预设倍频带的声压级,包括:In some embodiments, according to the predicted sound pressure level of each preset detection point in the preset octave band and the actual sound pressure level of each preset detection point in the preset octave band, the univariate linearity corresponding to the preset octave band The regression model is solved, and the sound pressure levels of multiple equivalent sound sources in the preset octave bands are obtained through double-layer optimization, including:
采用平均分布的方法,生成第一批随机变量,作为预设倍频带对应的单变量线性回归模型第一次优化时的输入变量;Using the method of average distribution, the first batch of random variables are generated as the input variables for the first optimization of the univariate linear regression model corresponding to the preset octave band;
根据MSE损失函数计算第一批随机变量中各个随机变量的误差值,进行第一次优化,从按照误差值由小到大排序的第一批随机变量中按顺序选取,排在前预设数量的随机变量;Calculate the error value of each random variable in the first batch of random variables according to the MSE loss function, perform the first optimization, select in order from the first batch of random variables sorted according to the error value from small to large, and rank the first preset number the random variable;
采用正态分布的方法对每一选取到的随机变量的数学期望值和方差值进行计算,生成第二批随机变量,作为预设倍频带对应的单变量线性回归模型第二次优化时的输入变量;Calculate the mathematical expectation and variance of each selected random variable using the normal distribution method to generate the second batch of random variables as the input for the second optimization of the univariate linear regression model corresponding to the preset octave band variable;
正态分布的公式如下:The formula for the normal distribution is as follows:
Figure PCTCN2022115386-appb-000004
Figure PCTCN2022115386-appb-000004
其中,μ为排在前预设数量的随机变量的数学期望值,σ 2为排在前预设数量的随机变量的方差值。 Among them, μ is the mathematical expectation value of the random variable with the preset number in front, and σ 2 is the variance value of the random variable with the preset number in the front.
根据交叉熵损失函数对第二批随机变量进行求解,进行第二次优化,最终求解得到各个等效声源在预设倍频带的声压级。According to the cross-entropy loss function, the second batch of random variables is solved, and the second optimization is carried out, and finally the sound pressure level of each equivalent sound source in the preset octave band is obtained through the solution.
其中,可以采用Uniform函数产生平均分布的第一批随机变量。每个随机变量均包含24个等效声源在预设倍频带的声压级,各个随机变量可以不同。Among them, the Uniform function can be used to generate the first batch of random variables with uniform distribution. Each random variable contains the sound pressure levels of 24 equivalent sound sources in preset octave bands, and each random variable can be different.
本实施例在第一次优化时,根据MSE损失函数计算预测值与实际测量值的全局误差,选取误差较小的预设数量的随机变量,根据该误差较小的预设数 量的随机变量的数学期望值和方差值,以正态分布生成第二批随机变量,在第二次优化时,采用交叉熵损失函数,从第二批随机变量中,最终求解得到各个等效声源在预设倍频带的声压级。通过第二次优化求解得到各个等效声源在预设倍频带的声压级。In the first optimization of this embodiment, the global error between the predicted value and the actual measured value is calculated according to the MSE loss function, and a preset number of random variables with smaller errors are selected, and according to the value of the preset number of random variables with smaller errors The mathematical expectation value and variance value are used to generate the second batch of random variables with a normal distribution. In the second optimization, the cross-entropy loss function is used to finally solve each equivalent sound source from the second batch of random variables in the preset Sound pressure level in octave bands. The sound pressure level of each equivalent sound source in the preset octave band is obtained through the second optimization solution.
在一些实施例中,上述MSE损失函数为:In some embodiments, the above MSE loss function is:
Figure PCTCN2022115386-appb-000005
Figure PCTCN2022115386-appb-000005
其中,MSE为MSE损失函数计算得到的误差值;m为预测检测点的数量;Lm j为第j个预设检测点在预设倍频带的实际声压级;Lw j为第j个预设检测点在预设倍频带的预测声压级,
Figure PCTCN2022115386-appb-000006
Lw ij=Lp i-20lg(d ij)-0.001*α*d ij-11;n为等效声源的数量;Lw ij为第i个等效声源在第j个预设检测点产生的在预设倍频带的声压级;Lp i为第i个等效声源在预设倍频带待求的声压级,为未知量;d ij为第i个等效声源与第j个预设检测点之间的距离;α为噪声在传播过程中的大气吸收衰减系数。
Among them, MSE is the error value calculated by the MSE loss function; m is the number of predicted detection points; Lm j is the actual sound pressure level of the jth preset detection point in the preset octave band; Lw j is the jth preset The predicted sound pressure level of the detection point in the preset octave band,
Figure PCTCN2022115386-appb-000006
Lw ij =Lp i -20lg(d ij )-0.001*α*d ij -11; n is the number of equivalent sound sources; Lw ij is the i-th equivalent sound source generated at the j-th preset detection point The sound pressure level in the preset octave band; Lp i is the sound pressure level of the i-th equivalent sound source in the preset octave band, which is an unknown quantity; d ij is the i-th equivalent sound source and the j-th The distance between preset detection points; α is the atmospheric absorption attenuation coefficient of noise during propagation.
其中,
Figure PCTCN2022115386-appb-000007
如图3所示,图3示出了4个场点(场点1、场点2、场点3和场点4)以及一个等效源平面该等效源平面中有4个等效声源,图3中的d表示的是其中一个场点(检测点)与其中一个等效声源之间的距离。
in,
Figure PCTCN2022115386-appb-000007
As shown in Figure 3, Figure 3 shows 4 field points (field point 1, field point 2, field point 3 and field point 4) and an equivalent source plane in which there are 4 equivalent acoustic source, d in Figure 3 represents the distance between one of the field points (detection points) and one of the equivalent sound sources.
MSE损失函数是一种全局敏感的损失函数,用于定位第二次优化时Lp i的候选区间。 The MSE loss function is a globally sensitive loss function, which is used to locate the candidate interval of Lp i in the second optimization.
在本实施例中,空间中等效声源的声压级为待求量,先设为Lp i。根据噪声在自由场传播中的衰减公式,每个等效声源在变压器到检测点的传播过程中,噪声传播只经历了几何发散衰减与大气吸收衰减,检测点处产生的声压级大小为:Lw ij=Lp i-20lg(d ij)-0.001*α*d ij-11,其中,α可以通过查询表1得到,在测量预设检测点的声压级时,会记录温度以及相对湿度。传播距离d ii决定了传播过程中的几何发散衰减量与大气吸收衰减量。 In this embodiment, the sound pressure level of the equivalent sound source in the space is the quantity to be sought, which is set as Lp i first. According to the attenuation formula of noise in free field propagation, during the propagation of each equivalent sound source from the transformer to the detection point, the noise propagation only undergoes geometric divergence attenuation and atmospheric absorption attenuation, and the sound pressure level generated at the detection point is : Lw ij =Lp i -20lg(d ij )-0.001*α*d ij -11, wherein, α can be obtained by looking up Table 1, when measuring the sound pressure level of the preset detection point, the temperature and relative humidity will be recorded . The propagation distance d ii determines the geometric divergence attenuation and atmospheric absorption attenuation in the propagation process.
将所有的等效声源在检测点处产生的声压级大小进行叠加,即为检测点处预测声压级大小,叠加公式如下:
Figure PCTCN2022115386-appb-000008
Superimpose the sound pressure levels generated by all equivalent sound sources at the detection point, which is the predicted sound pressure level at the detection point. The superposition formula is as follows:
Figure PCTCN2022115386-appb-000008
表1大气吸收衰减系数表Table 1 Atmospheric absorption attenuation coefficient table
Figure PCTCN2022115386-appb-000009
Figure PCTCN2022115386-appb-000009
在一些实施例中,上述交叉熵损失函数为:In some embodiments, the above-mentioned cross-entropy loss function is:
Figure PCTCN2022115386-appb-000010
Figure PCTCN2022115386-appb-000010
其中,L为交叉熵损失函数计算得到的值;P(Lp i)为Lp i在正态分布函数中的取值概率。 Among them, L is the value calculated by the cross-entropy loss function; P(Lp i ) is the value probability of Lp i in the normal distribution function.
交叉熵作为损失函数,同时为迭代收敛过程中的梯度下降提供数据,通过小批量随机梯度下降,收敛Lp i的取值区间,最后求解出Lp i的取值,即为所求的变压器多声源噪声等效模型的等效声源在预设倍频带的声压级大小。 As a loss function, cross entropy provides data for the gradient descent in the iterative convergence process. Through small-batch stochastic gradient descent, the value interval of Lp i is converged, and finally the value of Lp i is solved, which is the required transformer polyphony The sound pressure level of the equivalent sound source of the source noise equivalent model in the preset octave band.
以上所求只是等效声源单倍频带的声压级大小,若想获得63HZ、125HZ、250HZ、500HZ、1000HZ、2000HZ、4000HZ和8000HZ八个倍频带的声压级,需要构建8次模型,运用双层优化的方法进行8次最优求解,即可求得等效声源8个倍频带的声压级大小。The above requirements are only the sound pressure level of the single octave band of the equivalent sound source. If you want to obtain the sound pressure level of the eight octave bands of 63HZ, 125HZ, 250HZ, 500HZ, 1000HZ, 2000HZ, 4000HZ and 8000HZ, you need to build an 8-time model. The sound pressure level of the 8 octave bands of the equivalent sound source can be obtained by using the double-layer optimization method to perform 8 optimal solutions.
在S104中,根据等效声源的数量、各个等效声源在预设倍频带的声压级和各个等效声源的空间坐标得到预设倍频带对应的变压器多声源噪声等效模型。In S104, according to the number of equivalent sound sources, the sound pressure level of each equivalent sound source in the preset octave band and the spatial coordinates of each equivalent sound source, the transformer multi-source noise equivalent model corresponding to the preset octave band is obtained .
求解得到各个等效声源在预设倍频带的声压级之后,根据等效声源的数量、各个等效声源在预设倍频带的声压级和各个等效声源的空间坐标,即可得到预 设倍频带对应的变压器多声源噪声等效模型。After solving the sound pressure level of each equivalent sound source in the preset octave band, according to the number of equivalent sound sources, the sound pressure level of each equivalent sound source in the preset octave band and the spatial coordinates of each equivalent sound source, The equivalent model of transformer multi-source noise corresponding to the preset octave band can be obtained.
若求解得到各个等效声源在八个倍频带的声压级之后,则可以得到最终的变压器多声源噪声等效模型。After solving the sound pressure levels of each equivalent sound source in eight octave bands, the final multi-source noise equivalent model of the transformer can be obtained.
本实施例通过获取变压器周围的多个预设检测点的空间坐标与在预设倍频带的实际声压级;获取变压器的等效声源的数量和各个等效声源的空间坐标;根据多个预设检测点的空间坐标与在预设倍频带的实际声压级、等效声源的数量以及各个等效声源的空间坐标,构建预设倍频带对应的单变量线性回归模型,并对预设倍频带对应的单变量线性回归模型进行求解,得到多个等效声源在预设倍频带的声压级;根据等效声源的数量、各个等效声源的声压级和各个等效声源的空间坐标得到预设倍频带对应的变压器多声源噪声等效模型,能够克服近场声全息技术中的计算复杂和计算量大的问题,只需测量变压器附近少量检测点的声压级,并基于单变量线性回归的方法获得变压器多声源噪声等效模型,等效过程方便简单,运算量小。This embodiment obtains the spatial coordinates of multiple preset detection points around the transformer and the actual sound pressure level in the preset octave band; obtains the number of equivalent sound sources of the transformer and the spatial coordinates of each equivalent sound source; according to multiple The spatial coordinates of a preset detection point and the actual sound pressure level in the preset octave band, the number of equivalent sound sources and the spatial coordinates of each equivalent sound source, construct a univariate linear regression model corresponding to the preset octave band, and Solve the univariate linear regression model corresponding to the preset octave band to obtain the sound pressure levels of multiple equivalent sound sources in the preset octave band; according to the number of equivalent sound sources, the sound pressure levels of each equivalent sound source and The spatial coordinates of each equivalent sound source obtain the equivalent model of transformer multi-source noise corresponding to the preset octave band, which can overcome the problems of complex calculation and large amount of calculation in near-field acoustic holography technology, and only need to measure a small number of detection points near the transformer The sound pressure level of the transformer is obtained based on the method of univariate linear regression to obtain the equivalent model of transformer multi-sound source noise. The equivalent process is convenient and simple, and the amount of calculation is small.
本实施例通过噪声辐射衰减特性构建单变量线性回归模型,以MSE和交叉熵作为损失函数,通过双层优化获得变压器多声源等效模型的声压级。等效过程方便简单,运算量小,通过双层优化使得等效结果更加精确。In this embodiment, a univariate linear regression model is constructed based on the noise radiation attenuation characteristics, using MSE and cross-entropy as loss functions, and the sound pressure level of the multi-sound source equivalent model of the transformer is obtained through double-layer optimization. The equivalent process is convenient and simple, and the calculation amount is small, and the equivalent result is more accurate through double-layer optimization.
应理解,上述实施例中各步骤的序号的大小并不意味着执行顺序的先后,各过程的执行顺序应以其功能和内在逻辑确定,而不应对本发明实施例的实施过程构成任何限定。It should be understood that the sequence numbers of the steps in the above embodiments do not mean the order of execution, and the execution order of each process should be determined by its functions and internal logic, and should not constitute any limitation to the implementation process of the embodiment of the present invention.
以下为本发明的装置实施例,对于其中未详尽描述的细节,可以参考上述对应的方法实施例。The following are device embodiments of the present invention. For details that are not exhaustively described therein, reference may be made to the corresponding method embodiments above.
图4示出了本发明实施例提供的变压器多声源噪声等效模型的确定装置的结构示意图,为了便于说明,仅示出了与本发明实施例相关的部分,详述如下:Fig. 4 shows a schematic structural diagram of a device for determining the equivalent model of transformer multi-source noise provided by an embodiment of the present invention. For the convenience of description, only the parts related to the embodiment of the present invention are shown, and the details are as follows:
如图4所示,变压器多声源噪声等效模型的确定装置30包括:第一获取模块31、第二获取模块32、求解模块33和模型确定模块34。As shown in FIG. 4 , the device 30 for determining the equivalent model of transformer multi-source noise includes: a first acquisition module 31 , a second acquisition module 32 , a solution module 33 and a model determination module 34 .
第一获取模块31,用于获取变压器周围的多个预设检测点的空间坐标与在 预设倍频带的实际声压级;The first acquisition module 31 is used to acquire the spatial coordinates of a plurality of preset detection points around the transformer and the actual sound pressure level in the preset octave band;
第二获取模块32,用于获取变压器的等效声源的数量和各个等效声源的空间坐标;The second acquisition module 32 is used to acquire the quantity of equivalent sound sources of the transformer and the spatial coordinates of each equivalent sound source;
求解模块33,用于根据多个预设检测点的空间坐标与在预设倍频带的实际声压级、等效声源的数量以及各个等效声源的空间坐标,构建预设倍频带对应的单变量线性回归模型,并对预设倍频带对应的单变量线性回归模型进行求解,得到多个等效声源在预设倍频带的声压级;The solution module 33 is used to construct a preset octave band corresponding and solve the univariate linear regression model corresponding to the preset octave band to obtain the sound pressure levels of multiple equivalent sound sources in the preset octave band;
模型确定模块34,用于根据等效声源的数量、各个等效声源在预设倍频带的声压级和各个等效声源的空间坐标得到预设倍频带对应的变压器多声源噪声等效模型。The model determination module 34 is used to obtain the transformer multi-source noise corresponding to the preset octave band according to the number of equivalent sound sources, the sound pressure level of each equivalent sound source in the preset octave band and the spatial coordinates of each equivalent sound source equivalent model.
在一种可能的实现方式中,求解模块33具体用于:In a possible implementation, the solving module 33 is specifically used for:
根据多个预设检测点的空间坐标与在预设倍频带的实际声压级、等效声源的数量以及各个等效声源的空间坐标,构建预设倍频带对应的单变量线性回归模型;According to the spatial coordinates of multiple preset detection points and the actual sound pressure level in the preset octave band, the number of equivalent sound sources and the spatial coordinates of each equivalent sound source, a univariate linear regression model corresponding to the preset octave band is constructed ;
根据各个预设检测点在预设倍频带的预测声压级以及各个预设检测点在预设倍频带的实际声压级,对预设倍频带对应的单变量线性回归模型进行求解,通过双层优化得到多个等效声源在预设倍频带的声压级。According to the predicted sound pressure level of each preset detection point in the preset octave band and the actual sound pressure level of each preset detection point in the preset octave band, the univariate linear regression model corresponding to the preset octave band is solved. Layer optimization obtains the sound pressure levels of multiple equivalent sound sources in preset octave bands.
在一种可能的实现方式中,求解模块33具体用于:In a possible implementation, the solving module 33 is specifically used for:
采用平均分布的方法,生成第一批随机变量,作为预设倍频带对应的单变量线性回归模型第一次优化时的输入变量;Using the method of average distribution, the first batch of random variables are generated as the input variables for the first optimization of the univariate linear regression model corresponding to the preset octave band;
根据MSE损失函数,进行第一次优化,从第一批随机变量中选取MSE损失函数计算得到的误差值按照从小到大的顺序,排在前预设数量的随机变量;According to the MSE loss function, the first optimization is carried out, and the error value calculated by the MSE loss function is selected from the first batch of random variables in order from small to large, and the first preset number of random variables are ranked;
根据排在前预设数量的随机变量的数学期望值和方差值,采用正态分布的方法,生成第二批随机变量,作为预设倍频带对应的单变量线性回归模型第二次优化时的输入变量;According to the mathematical expectation and variance of the preset number of random variables, the second batch of random variables is generated by using the normal distribution method, as the second optimization of the univariate linear regression model corresponding to the preset octave band input variable;
根据交叉熵损失函数,进行第二次优化,从第二批随机变量中,最终求解 得到各个等效声源在预设倍频带的声压级。According to the cross-entropy loss function, the second optimization is carried out, and the sound pressure level of each equivalent sound source in the preset octave band is finally obtained from the second batch of random variables.
在一种可能的实现方式中,MSE损失函数为:In one possible implementation, the MSE loss function is:
Figure PCTCN2022115386-appb-000011
Figure PCTCN2022115386-appb-000011
其中,MSE为MSE损失函数计算得到的误差值;m为预测检测点的数量;Lm j为第j个预设检测点在预设倍频带的实际声压级;Lw j为第j个预设检测点在预设倍频带的预测声压级,
Figure PCTCN2022115386-appb-000012
Lw ij=Lp i-20lg(d ij)-0.001*α*d ij-11;n为等效声源的数量;Lw ij为第i个等效声源在第j个预设检测点产生的在预设倍频带的声压级;Lp i为第i个等效声源在预设倍频带待求的声压级,为未知量;d ij为第i个等效声源与第j个预设检测点之间的距离;α为噪声在传播过程中的大气吸收衰减系数。
Among them, MSE is the error value calculated by the MSE loss function; m is the number of predicted detection points; Lm j is the actual sound pressure level of the jth preset detection point in the preset octave band; Lw j is the jth preset The predicted sound pressure level of the detection point in the preset octave band,
Figure PCTCN2022115386-appb-000012
Lw ij =Lp i -20lg(d ij )-0.001*α*d ij -11; n is the number of equivalent sound sources; Lw ij is the i-th equivalent sound source generated at the j-th preset detection point The sound pressure level in the preset octave band; Lp i is the sound pressure level of the i-th equivalent sound source in the preset octave band, which is an unknown quantity; d ij is the i-th equivalent sound source and the j-th The distance between preset detection points; α is the atmospheric absorption attenuation coefficient of noise during propagation.
在一种可能的实现方式中,交叉熵损失函数为:In one possible implementation, the cross-entropy loss function is:
Figure PCTCN2022115386-appb-000013
Figure PCTCN2022115386-appb-000013
其中,L为交叉熵损失函数计算得到的值;P(Lp i)为Lp i在正太分布函数中的取值概率。 Among them, L is the value calculated by the cross-entropy loss function; P(Lp i ) is the value probability of Lp i in the normal distribution function.
在一种可能的实现方式中,等效声源的数量为24个;In a possible implementation manner, the number of equivalent sound sources is 24;
等效声源的分布方式为:长箱壁面按照4*2规格等间距布置,短箱壁面按照2*2规格等间距布置。The distribution of equivalent sound sources is as follows: the walls of the long box are arranged at equal intervals according to the 4*2 specification, and the walls of the short box are arranged at equal intervals according to the 2*2 specification.
在另一个实施示例中,上述变压器多声源噪声等效模型的确定装置包括:处理器,其中所述处理器用于执行存在存储器的上述程序模块,包括:第一获取模块31、第二获取模块32、求解模块33和模型确定模块34。In another implementation example, the above-mentioned device for determining the multi-source noise equivalent model of a transformer includes: a processor, wherein the processor is used to execute the above-mentioned program modules stored in the memory, including: a first acquisition module 31, a second acquisition module 32. A solution module 33 and a model determination module 34.
图5是本发明实施例提供的终端的示意图。如图5所示,该实施例的终端4包括:处理器40、存储器41以及存储在所述存储器41中并可在所述处理器40上运行的计算机程序42。所述处理器40执行所述计算机程序42时实现上述各个变压器多声源噪声等效模型的确定方法实施例中的步骤,例如图1所示的S101至S104。或者,所述处理器40执行所述计算机程序42时实现上述各装 置实施例中各模块/单元的功能,例如图4所示模块/单元31至34的功能。Fig. 5 is a schematic diagram of a terminal provided by an embodiment of the present invention. As shown in FIG. 5 , the terminal 4 of this embodiment includes: a processor 40 , a memory 41 and a computer program 42 stored in the memory 41 and operable on the processor 40 . When the processor 40 executes the computer program 42, it implements the steps in the embodiment of the method for determining the equivalent model of each transformer multi-source noise, such as S101 to S104 shown in FIG. 1 . Or, when the processor 40 executes the computer program 42, it realizes the functions of the modules/units in the above-mentioned device embodiments, such as the functions of the modules/units 31 to 34 shown in FIG. 4 .
示例性的,所述计算机程序42可以被分割成一个或多个模块/单元,所述一个或者多个模块/单元被存储在所述存储器41中,并由所述处理器40执行,以完成本发明。所述一个或多个模块/单元可以是能够完成特定功能的一系列计算机程序指令段,该指令段用于描述所述计算机程序42在所述终端4中的执行过程。例如,所述计算机程序42可以被分割成图4所示的模块/单元31至34。Exemplarily, the computer program 42 can be divided into one or more modules/units, and the one or more modules/units are stored in the memory 41 and executed by the processor 40 to complete this invention. The one or more modules/units may be a series of computer program instruction segments capable of accomplishing specific functions, and the instruction segments are used to describe the execution process of the computer program 42 in the terminal 4 . For example, the computer program 42 may be divided into the modules/units 31 to 34 shown in FIG. 4 .
所述终端4可以桌上型计算机、笔记本、掌上电脑及云端服务器等计算设备。所述终端4可包括,但不仅限于,处理器40、存储器41。本领域技术人员可以理解,图5仅仅是终端4的示例,并不构成对终端4的限定,可以包括比图示更多或更少的部件,或者组合某些部件,或者不同的部件,例如所述终端还可以包括输入输出设备、网络接入设备、总线等。The terminal 4 can be a computing device such as a desktop computer, a notebook, a palmtop computer, or a cloud server. The terminal 4 may include, but not limited to, a processor 40 and a memory 41 . Those skilled in the art can understand that FIG. 5 is only an example of the terminal 4 and does not constitute a limitation to the terminal 4. It may include more or less components than those shown in the figure, or combine certain components, or different components, such as The terminal may also include an input and output device, a network access device, a bus, and the like.
所称处理器40可以是中央处理单元(Central Processing Unit,CPU),还可以是其他通用处理器、数字信号处理器(Digital Signal Processor,DSP)、专用集成电路(Application Specific Integrated Circuit,ASIC)、现场可编程门阵列(Field-Programmable Gate Array,FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。The so-called processor 40 can be a central processing unit (Central Processing Unit, CPU), and can also be other general-purpose processors, digital signal processors (Digital Signal Processor, DSP), application specific integrated circuits (Application Specific Integrated Circuit, ASIC), Field-Programmable Gate Array (Field-Programmable Gate Array, FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. A general-purpose processor may be a microprocessor, or the processor may be any conventional processor, or the like.
所述存储器41可以是所述终端4的内部存储单元,例如终端4的硬盘或内存。所述存储器41也可以是所述终端4的外部存储设备,例如所述终端4上配备的插接式硬盘,智能存储卡(Smart Media Card,SMC),安全数字(Secure Digital,SD)卡,闪存卡(Flash Card)等。进一步地,所述存储器41还可以既包括所述终端4的内部存储单元也包括外部存储设备。所述存储器41用于存储所述计算机程序以及所述终端所需的其他程序和数据。所述存储器41还可以用于暂时地存储已经输出或者将要输出的数据。The memory 41 may be an internal storage unit of the terminal 4 , such as a hard disk or memory of the terminal 4 . The memory 41 can also be an external storage device of the terminal 4, such as a plug-in hard disk equipped on the terminal 4, a smart memory card (Smart Media Card, SMC), a secure digital (Secure Digital, SD) card, Flash card (Flash Card), etc. Further, the memory 41 may also include both an internal storage unit of the terminal 4 and an external storage device. The memory 41 is used to store the computer program and other programs and data required by the terminal. The memory 41 can also be used to temporarily store data that has been output or will be output.
所属领域的技术人员可以清楚地了解到,为了描述的方便和简洁,仅以上述各功能单元、模块的划分进行举例说明,实际应用中,可以根据需要而将上 述功能分配由不同的功能单元、模块完成,即将所述装置的内部结构划分成不同的功能单元或模块,以完成以上描述的全部或者部分功能。实施例中的各功能单元、模块可以集成在一个处理单元中,也可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个单元中,上述集成的单元既可以采用硬件的形式实现,也可以采用软件功能单元的形式实现。另外,各功能单元、模块的具体名称也只是为了便于相互区分,并不用于限制本申请的保护范围。上述系统中单元、模块的具体工作过程,可以参考前述方法实施例中的对应过程,在此不再赘述。Those skilled in the art can clearly understand that for the convenience and brevity of description, only the division of the above-mentioned functional units and modules is used for illustration. In practical applications, the above-mentioned functions can be assigned to different functional units, Completion of modules means that the internal structure of the device is divided into different functional units or modules to complete all or part of the functions described above. Each functional unit and module in the embodiment may be integrated into one processing unit, or each unit may exist separately physically, or two or more units may be integrated into one unit, and the above-mentioned integrated units may adopt hardware It can also be implemented in the form of software functional units. In addition, the specific names of the functional units and modules are only for the convenience of distinguishing each other, and are not used to limit the protection scope of the present application. For the specific working processes of the units and modules in the above system, reference may be made to the corresponding processes in the aforementioned method embodiments, and details will not be repeated here.
在上述实施例中,对各个实施例的描述都各有侧重,某个实施例中没有详述或记载的部分,可以参见其它实施例的相关描述。In the above-mentioned embodiments, the descriptions of each embodiment have their own emphases, and for parts that are not detailed or recorded in a certain embodiment, refer to the relevant descriptions of other embodiments.
本领域普通技术人员可以意识到,结合本文中所公开的实施例描述的各示例的单元及算法步骤,能够以电子硬件、或者计算机软件和电子硬件的结合来实现。这些功能究竟以硬件还是软件方式来执行,取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能,但是这种实现不应认为超出本发明的范围。Those skilled in the art can appreciate that the units and algorithm steps of the examples described in conjunction with the embodiments disclosed herein can be implemented by electronic hardware, or a combination of computer software and electronic hardware. Whether these functions are executed by hardware or software depends on the specific application and design constraints of the technical solution. Skilled artisans may use different methods to implement the described functions for each specific application, but such implementation should not be regarded as exceeding the scope of the present invention.
在本发明所提供的实施例中,应该理解到,所揭露的装置/终端和方法,可以通过其它的方式实现。例如,以上所描述的装置/终端实施例仅仅是示意性的,例如,所述模块或单元的划分,仅仅为一种逻辑功能划分,实际实现时可以有另外的划分方式,例如多个单元或组件可以结合或者可以集成到另一个系统,或一些特征可以忽略,或不执行。另一点,所显示或讨论的相互之间的耦合或直接耦合或通讯连接可以是通过一些接口,装置或单元的间接耦合或通讯连接,可以是电性,机械或其它的形式。In the embodiments provided in the present invention, it should be understood that the disclosed device/terminal and method may be implemented in other ways. For example, the device/terminal embodiments described above are only illustrative. For example, the division of the modules or units is only a logical function division. In actual implementation, there may be other division methods, such as multiple units or Components may be combined or integrated into another system, or some features may be omitted, or not implemented. In another point, the mutual coupling or direct coupling or communication connection shown or discussed may be through some interfaces, and the indirect coupling or communication connection of devices or units may be in electrical, mechanical or other forms.
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。The units described as separate components may or may not be physically separated, and the components shown as units may or may not be physical units, that is, they may be located in one place, or may be distributed to multiple network units. Part or all of the units can be selected according to actual needs to achieve the purpose of the solution of this embodiment.
另外,在本发明各个实施例中的各功能单元可以集成在一个处理单元中,也可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现,也可以采用软件功能单元的形式实现。In addition, each functional unit in each embodiment of the present invention may be integrated into one processing unit, each unit may exist separately physically, or two or more units may be integrated into one unit. The above-mentioned integrated units can be implemented in the form of hardware or in the form of software functional units.
所述集成的模块/单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。基于这样的理解,本发明实现上述实施例方法中的全部或部分流程,也可以通过计算机程序来指令相关的硬件来完成,所述的计算机程序可存储于一计算机可读存储介质中,该计算机程序在被处理器执行时,可实现上述各个变压器多声源噪声等效模型的确定方法实施例的步骤。其中,所述计算机程序包括计算机程序代码,所述计算机程序代码可以为源代码形式、对象代码形式、可执行文件或某些中间形式等。所述计算机可读介质可以包括:能够携带所述计算机程序代码的任何实体或装置、记录介质、U盘、移动硬盘、磁碟、光盘、计算机存储器、只读存储器(Read-Only Memory,ROM)、随机存取存储器(Random Access Memory,RAM)、电载波信号、电信信号以及软件分发介质等。需要说明的是,所述计算机可读介质包含的内容可以根据司法管辖区内立法和专利实践的要求进行适当的增减,例如在某些司法管辖区,根据立法和专利实践,计算机可读介质不包括是电载波信号和电信信号。If the integrated module/unit is realized in the form of a software function unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, the present invention realizes all or part of the processes in the methods of the above embodiments, and can also be completed by instructing related hardware through a computer program. The computer program can be stored in a computer-readable storage medium, and the computer When the program is executed by the processor, the steps in the above embodiments of the method for determining the equivalent model of the transformer multi-source noise can be realized. Wherein, the computer program includes computer program code, and the computer program code may be in the form of source code, object code, executable file or some intermediate form. The computer-readable medium may include: any entity or device capable of carrying the computer program code, a recording medium, a USB flash drive, a removable hard disk, a magnetic disk, an optical disk, a computer memory, and a read-only memory (Read-Only Memory, ROM) , random access memory (Random Access Memory, RAM), electric carrier signal, telecommunication signal and software distribution medium, etc. It should be noted that the content contained in the computer-readable medium may be appropriately increased or decreased according to the requirements of legislation and patent practice in the jurisdiction. For example, in some jurisdictions, computer-readable media Excluding electrical carrier signals and telecommunication signals.
以上所述实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围,均应包含在本发明的保护范围之内。The above-described embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit them; although the present invention has been described in detail with reference to the foregoing embodiments, those of ordinary skill in the art should understand that: it can still carry out the foregoing embodiments Modifications to the technical solutions recorded in the examples, or equivalent replacement of some of the technical features; and these modifications or replacements do not make the essence of the corresponding technical solutions deviate from the spirit and scope of the technical solutions of the various embodiments of the present invention, and should be included in within the protection scope of the present invention.

Claims (10)

  1. 一种变压器多声源噪声等效模型的确定方法,其特征在于,包括:A method for determining an equivalent model of transformer multi-source noise, characterized in that it includes:
    获取变压器周围的多个预设检测点的空间坐标与在预设倍频带的实际声压级;Obtain the spatial coordinates of multiple preset detection points around the transformer and the actual sound pressure level in the preset octave band;
    获取变压器的等效声源的数量和各个等效声源的空间坐标;Obtain the number of equivalent sound sources of the transformer and the spatial coordinates of each equivalent sound source;
    根据多个预设检测点的空间坐标与在预设倍频带的实际声压级、所述等效声源的数量以及各个等效声源的空间坐标,构建所述预设倍频带对应的单变量线性回归模型,并对所述预设倍频带对应的单变量线性回归模型进行求解,得到多个等效声源在所述预设倍频带的声压级;According to the spatial coordinates of a plurality of preset detection points and the actual sound pressure level in the preset octave band, the number of the equivalent sound sources and the spatial coordinates of each equivalent sound source, construct the unit corresponding to the preset octave band A variable linear regression model, and solving the univariate linear regression model corresponding to the preset octave band to obtain the sound pressure levels of multiple equivalent sound sources in the preset octave band;
    根据等效声源的数量、各个等效声源在所述预设倍频带的声压级和各个等效声源的空间坐标得到所述预设倍频带对应的变压器多声源噪声等效模型。According to the number of equivalent sound sources, the sound pressure level of each equivalent sound source in the preset octave band and the spatial coordinates of each equivalent sound source, the transformer multi-source noise equivalent model corresponding to the preset octave band is obtained .
  2. 根据权利要求1所述的变压器多声源噪声等效模型的确定方法,其特征在于,所述根据多个预设检测点的空间坐标与在预设倍频带的实际声压级、所述等效声源的数量以及各个等效声源的空间坐标,构建所述预设倍频带对应的单变量线性回归模型,并对所述预设倍频带对应的单变量线性回归模型进行求解,得到多个等效声源在所述预设倍频带的声压级,包括:The method for determining the equivalent model of transformer multi-sound source noise according to claim 1, wherein the spatial coordinates of multiple preset detection points and the actual sound pressure level in the preset octave band, the etc. The number of effective sound sources and the spatial coordinates of each equivalent sound source, construct the univariate linear regression model corresponding to the preset octave band, and solve the univariate linear regression model corresponding to the preset octave band, and obtain multiple The sound pressure level of an equivalent sound source in the preset octave band, including:
    根据多个预设检测点的空间坐标与在预设倍频带的实际声压级、所述等效声源的数量以及各个等效声源的空间坐标,构建所述预设倍频带对应的单变量线性回归模型;According to the spatial coordinates of a plurality of preset detection points and the actual sound pressure level in the preset octave band, the number of the equivalent sound sources and the spatial coordinates of each equivalent sound source, construct the unit corresponding to the preset octave band variable linear regression model;
    根据各个预设检测点在预设倍频带的预测声压级以及各个预设检测点在预设倍频带的实际声压级,对所述预设倍频带对应的单变量线性回归模型进行求解,通过双层优化得到多个等效声源在所述预设倍频带的声压级。According to the predicted sound pressure level of each preset detection point in the preset octave band and the actual sound pressure level of each preset detection point in the preset octave band, the univariate linear regression model corresponding to the preset octave band is solved, The sound pressure levels of multiple equivalent sound sources in the preset octave bands are obtained through double-layer optimization.
  3. 根据权利要求2所述的变压器多声源噪声等效模型的确定方法,其特征在于,所述根据各个预设检测点在预设倍频带的预测声压级以及各个预设检测点在预设倍频带的实际声压级,对所述预设倍频带对应的单变量线性回归模型进行求解,通过双层优化得到多个等效声源在所述预设倍频带的声压级,包括:The method for determining the equivalent model of transformer multi-sound source noise according to claim 2, wherein the predicted sound pressure level of each preset detection point in the preset octave band and each preset detection point in the preset For the actual sound pressure level of the octave band, the univariate linear regression model corresponding to the preset octave band is solved, and the sound pressure levels of multiple equivalent sound sources in the preset octave band are obtained through double-layer optimization, including:
    采用平均分布的方法,生成第一批随机变量,作为所述预设倍频带对应的单变量线性回归模型第一次优化时的输入变量;Using the method of average distribution, generate the first batch of random variables, as the input variable when the univariate linear regression model corresponding to the preset octave band is optimized for the first time;
    根据MSE损失函数计算第一批随机变量中各个随机变量的误差值,从按照误差值由小至大排序的所述第一批随机变量中按排序选取预设数量的随机变量;Calculate the error value of each random variable in the first batch of random variables according to the MSE loss function, and select a preset number of random variables in order from the first batch of random variables sorted from small to large according to the error value;
    采用正态分布的方法对每一选取到的随机变量的数学期望值和方差值进行计算,生成第二批随机变量,作为所述预设倍频带对应的单变量线性回归模型第二次优化时的输入变量;The mathematical expectation value and variance value of each selected random variable are calculated using the method of normal distribution, and the second batch of random variables are generated as the second optimization of the univariate linear regression model corresponding to the preset octave band the input variable;
    根据交叉熵损失函数对所述第二批随机变量进行求解,得到各个等效声源在所述预设倍频带的声压级。The second batch of random variables are solved according to the cross-entropy loss function to obtain the sound pressure level of each equivalent sound source in the preset octave band.
  4. 根据权利要求3所述的变压器多声源噪声等效模型的确定方法,其特征在于,所述MSE损失函数为:The method for determining the equivalent model of transformer multi-sound source noise according to claim 3, wherein the MSE loss function is:
    Figure PCTCN2022115386-appb-100001
    Figure PCTCN2022115386-appb-100001
    其中,MSE为MSE损失函数计算得到的误差值;m为预测检测点的数量;Lm j为第j个预设检测点在预设倍频带的实际声压级;Lw j为第j个预设检测点在预设倍频带的预测声压级,
    Figure PCTCN2022115386-appb-100002
    Lw ij=Lp i-20lg(d ij)-0.001*α*d ij-11;n为等效声源的数量;Lw ij为第i个等效声源在第j个预设检测点产生的在预设倍频带的声压级;Lp i为第i个等效声源在预设倍频带待求的声压级,为未知量;d ij为第i个等效声源与第j个预设检测点之间的距离;α为噪声在传播过程中的大气吸收衰减系数。
    Among them, MSE is the error value calculated by the MSE loss function; m is the number of predicted detection points; Lm j is the actual sound pressure level of the jth preset detection point in the preset octave band; Lw j is the jth preset The predicted sound pressure level of the detection point in the preset octave band,
    Figure PCTCN2022115386-appb-100002
    Lw ij =Lp i -20lg(d ij )-0.001*α*d ij -11; n is the number of equivalent sound sources; Lw ij is the i-th equivalent sound source generated at the j-th preset detection point The sound pressure level in the preset octave band; Lp i is the sound pressure level of the i-th equivalent sound source in the preset octave band, which is an unknown quantity; d ij is the i-th equivalent sound source and the j-th The distance between preset detection points; α is the atmospheric absorption attenuation coefficient of noise during propagation.
  5. 根据权利要求4所述的变压器多声源噪声等效模型的确定方法,其特征在于,所述交叉熵损失函数为:The method for determining the equivalent model of transformer multi-source noise according to claim 4, wherein the cross-entropy loss function is:
    Figure PCTCN2022115386-appb-100003
    Figure PCTCN2022115386-appb-100003
    其中,L为交叉熵损失函数计算得到的值;P(Lp i)为Lp i在正太分布函数中的取值概率。 Among them, L is the value calculated by the cross-entropy loss function; P(Lp i ) is the value probability of Lp i in the normal distribution function.
  6. 根据权利要求1至5任一项所述的变压器多声源噪声等效模型的确定方法,其特征在于,所述等效声源的数量为24个;The method for determining the equivalent model of transformer multi-sound source noise according to any one of claims 1 to 5, wherein the number of said equivalent sound sources is 24;
    所述等效声源的分布方式为:长箱壁面按照4*2规格等间距布置,短箱壁面按照2*2规格等间距布置。The distribution of the equivalent sound sources is as follows: the wall surface of the long box is arranged at equal intervals according to the 4*2 specification, and the wall surface of the short box is arranged at equal intervals according to the 2*2 specification.
  7. 一种变压器多声源噪声等效模型的确定装置,其特征在于,包括:A device for determining an equivalent model of transformer multi-source noise, characterized in that it includes:
    第一获取模块,用于获取变压器周围的多个预设检测点的空间坐标与在预设倍频带的实际声压级;The first acquisition module is used to acquire the spatial coordinates of multiple preset detection points around the transformer and the actual sound pressure level in the preset octave band;
    第二获取模块,用于获取变压器的等效声源的数量和各个等效声源的空间坐标;The second acquisition module is used to acquire the number of equivalent sound sources of the transformer and the spatial coordinates of each equivalent sound source;
    求解模块,用于根据多个预设检测点的空间坐标与在预设倍频带的实际声压级、所述等效声源的数量以及各个等效声源的空间坐标,构建所述预设倍频带对应的单变量线性回归模型,并对所述预设倍频带对应的单变量线性回归模型进行求解,得到多个等效声源在所述预设倍频带的声压级;A solution module, configured to construct the preset according to the spatial coordinates of multiple preset detection points, the actual sound pressure level in the preset octave band, the number of the equivalent sound sources, and the spatial coordinates of each equivalent sound source A univariate linear regression model corresponding to the octave band, and solving the univariate linear regression model corresponding to the preset octave band to obtain the sound pressure levels of multiple equivalent sound sources in the preset octave band;
    模型确定模块,用于根据等效声源的数量、各个等效声源在所述预设倍频带的声压级和各个等效声源的空间坐标得到所述预设倍频带对应的变压器多声源噪声等效模型。The model determination module is used to obtain the number of transformers corresponding to the preset octave band according to the number of equivalent sound sources, the sound pressure level of each equivalent sound source in the preset octave band and the spatial coordinates of each equivalent sound source Equivalent model of sound source noise.
  8. 根据权利要求7所述的变压器多声源噪声等效模型的确定装置,其特征在于,所述求解模块具体用于:The device for determining the equivalent model of transformer multi-sound source noise according to claim 7, wherein the solving module is specifically used for:
    根据多个预设检测点的空间坐标与在预设倍频带的实际声压级、所述等效声源的数量以及各个等效声源的空间坐标,构建所述预设倍频带对应的单变量线性回归模型;According to the spatial coordinates of a plurality of preset detection points and the actual sound pressure level in the preset octave band, the number of the equivalent sound sources and the spatial coordinates of each equivalent sound source, construct the unit corresponding to the preset octave band variable linear regression model;
    根据各个预设检测点在预设倍频带的预测声压级以及各个预设检测点在预设倍频带的实际声压级,对所述预设倍频带对应的单变量线性回归模型进行求解,通过双层优化得到多个等效声源在所述预设倍频带的声压级。According to the predicted sound pressure level of each preset detection point in the preset octave band and the actual sound pressure level of each preset detection point in the preset octave band, the univariate linear regression model corresponding to the preset octave band is solved, The sound pressure levels of multiple equivalent sound sources in the preset octave bands are obtained through double-layer optimization.
  9. 一种终端,包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现如上的权利要求1至6中任一项所述变压器多声源噪声等效模型的确定方法的步骤。A terminal, comprising a memory, a processor, and a computer program stored in the memory and operable on the processor, characterized in that, when the processor executes the computer program, the above claims 1 to 1 are implemented. Steps in the determination method of the multi-source noise equivalent model of the transformer described in any one of 6.
  10. 一种计算机可读存储介质,所述计算机可读存储介质存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实现如上的权利要求1至6中任一项所述变压器多声源噪声等效模型的确定方法的步骤。A computer-readable storage medium, the computer-readable storage medium stores a computer program, characterized in that, when the computer program is executed by a processor, it realizes the multi-sound transformer described in any one of claims 1 to 6 above. The steps of the determination method of the source noise equivalent model.
PCT/CN2022/115386 2021-12-10 2022-08-29 Method for determining transformer multi-sound-source noise equivalent model, terminal, and storage medium WO2023103468A1 (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN202111510147.1 2021-12-10
CN202111510147.1A CN114201875A (en) 2021-12-10 2021-12-10 Method for determining multi-sound-source noise equivalent model of transformer, terminal and storage medium

Publications (1)

Publication Number Publication Date
WO2023103468A1 true WO2023103468A1 (en) 2023-06-15

Family

ID=80652421

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2022/115386 WO2023103468A1 (en) 2021-12-10 2022-08-29 Method for determining transformer multi-sound-source noise equivalent model, terminal, and storage medium

Country Status (2)

Country Link
CN (1) CN114201875A (en)
WO (1) WO2023103468A1 (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114201875A (en) * 2021-12-10 2022-03-18 国网河北省电力有限公司经济技术研究院 Method for determining multi-sound-source noise equivalent model of transformer, terminal and storage medium

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001159559A (en) * 1999-12-03 2001-06-12 Chubu Electric Power Co Inc Noise measuring method
CN111159928A (en) * 2019-11-26 2020-05-15 中国电力科学研究院有限公司 Transformer noise calculation method and system based on multi-line sound source model
CN111261188A (en) * 2020-01-20 2020-06-09 中国电力科学研究院有限公司 High-voltage transformer noise spectrum determination method and device
CN114201875A (en) * 2021-12-10 2022-03-18 国网河北省电力有限公司经济技术研究院 Method for determining multi-sound-source noise equivalent model of transformer, terminal and storage medium

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001159559A (en) * 1999-12-03 2001-06-12 Chubu Electric Power Co Inc Noise measuring method
CN111159928A (en) * 2019-11-26 2020-05-15 中国电力科学研究院有限公司 Transformer noise calculation method and system based on multi-line sound source model
CN111261188A (en) * 2020-01-20 2020-06-09 中国电力科学研究院有限公司 High-voltage transformer noise spectrum determination method and device
CN114201875A (en) * 2021-12-10 2022-03-18 国网河北省电力有限公司经济技术研究院 Method for determining multi-sound-source noise equivalent model of transformer, terminal and storage medium

Also Published As

Publication number Publication date
CN114201875A (en) 2022-03-18

Similar Documents

Publication Publication Date Title
US8423328B2 (en) Method of distributing a random variable using statistically correct spatial interpolation continuously with spatially inhomogeneous statistical correlation versus distance, standard deviation, and mean
Marin Treatment of singularities in the method of fundamental solutions for two-dimensional Helmholtz-type equations
WO2023103468A1 (en) Method for determining transformer multi-sound-source noise equivalent model, terminal, and storage medium
CN112288254B (en) Construction quality detection method, terminal and storage medium
CN115878963B (en) Method, system, terminal and medium for predicting metering error of capacitive voltage transformer
JP2009098891A (en) Simulation device, simulation program, recording medium in which simulation program is stored and simulation method
CN112990330A (en) User energy abnormal data detection method and device
CN112308969B (en) Three-dimensional map generation method and device
CN114897032B (en) Current transformer fault diagnosis method and device based on width learning
Zhang et al. Short-term load forecasting method based on EWT and IDBSCAN
CN115239105A (en) Method and device for evaluating wind resources of in-service wind power plant
CN117557681A (en) High-precision topographic map generation method and device based on multi-source mapping data
CN104635206B (en) A kind of method and device of wireless location
KR102318309B1 (en) A method for determining a doping concentration of a three-dimensional structure and a method of manufacturing a semiconductor device using the same
CN117194872A (en) Wind power plant wind speed correction method, system, computer equipment and storage medium
CN116860562A (en) Method and system for monitoring data quality of data center
Sundström et al. Optimal transport based impulse response interpolation in the presence of calibration errors
Wojcikowski et al. A surrogate-assisted measurement correction method for accurate and low-cost monitoring of particulate matter pollutants
CN115983053A (en) Electromagnetic simulation method of antenna with dielectric substrate based on moment method
Liang et al. Fractal Design of Indoor and Outdoor Forms of Architectural Space Based on a Three‐Dimensional Box Dimension Algorithm
CN111787608B (en) Anchor node laying method and device
CN114519267A (en) Data updating method of underground cable model
CN112182739A (en) Aircraft structure non-probability credible reliability topological optimization design method
CN116713808B (en) Method and system for accurately detecting hole level in shielding case
TWI779561B (en) Three-dimensional sound localization method, electronic device and computer readable storage

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 22902900

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE