Spatial response compensation method and apparatus, electronic device, and storage medium
By decomposing the spatial transfer function into a minimum-phase system and an all-pass system, determining their inverse systems and constructing a compensation system, the problem of high complexity in the compensation model in the prior art is solved, and efficient and accurate spatial response compensation is achieved.
Patent Information
- Authority / Receiving Office
- CN · China
- Patent Type
- Patents(China)
- Current Assignee / Owner
- SHANGHAI AWINIC TECH CO LTD
- Filing Date
- 2023-05-23
- Publication Date
- 2026-05-22
Smart Images
Figure CN116582795B_ABST
Abstract
Description
Technical Field
[0001] This application relates to the field of digital signal processing technology, and in particular to a spatial response compensation method, apparatus, electronic device, and storage medium. Background Technology
[0002] No device responds ideally in space; reflections, standing waves, and other environmental effects exist. Pre-designed compensation filters can mitigate these effects. However, different environments or listening areas may have varying physical characteristics, generating positive or negative interference in sound waves. This can lead to the enhancement or weakening of various frequencies or acoustic information, diminishing the effectiveness of pre-designed compensation filters in reducing environmental impact. Therefore, to mitigate environmental effects, manual observation and adjustment are typically employed.
[0003] Currently, traditional manual observation and debugging methods require a lot of manpower to achieve the desired effect. Therefore, automatic compensation methods have been gradually developed, such as using pre-compensation models to compensate for spatial response, which greatly saves the cost of manual debugging.
[0004] However, existing automatic compensation methods compensate for the entire spatial response, resulting in a very high complexity of the compensation model. Summary of the Invention
[0005] This invention provides a spatial response compensation method, apparatus, electronic device, and storage medium to solve or alleviate technical problems in the prior art.
[0006] The technical solution adopted in this invention is as follows:
[0007] In a first aspect, embodiments of this application provide a spatial response compensation method, which includes:
[0008] Based on the system's spatial input information and corresponding spatial response information, a spatial transfer function is constructed, wherein the spatial transfer function indicates the mapping relationship between the system's input signal and output signal. The spatial transfer function is then decomposed into a minimum-phase system and an all-pass system. A first inverse system corresponding to the minimum-phase system and a second inverse system corresponding to the all-pass system are determined. A compensation system is constructed based on a preset target frequency response, the first inverse system, and the second inverse system. The spatial transfer function is compensated using the compensation system to obtain a compensated spatial transfer function. The spatial input information of the system is converted into target spatial response information using the compensated spatial transfer function, wherein the difference between the frequency response of the target spatial response information and the target frequency response is less than a preset threshold.
[0009] Optionally, in one embodiment, splitting the spatial transfer function into a minimum-phase system and an all-pass system includes: determining the all-pass system based on the external zeros of the spatial transfer function; and determining the minimum-phase system based on the internal zeros, external zeros, and first internal poles of the spatial transfer function.
[0010] Optionally, in one embodiment, determining the all-pass system based on the external zeros of the spatial transfer function includes: calculating the reciprocal of the external zeros to obtain a second internal pole; and determining the all-pass system based on the second internal pole.
[0011] Optionally, in one embodiment, determining the minimum phase system based on the first inner circle zero, the outer circle zero, and the first inner circle pole of the spatial transfer function includes: calculating the reciprocal of the outer circle zero to obtain the second inner circle zero; and determining the minimum phase system based on the first inner circle zero, the second inner circle zero, and the first inner circle pole.
[0012] Optionally, in one embodiment, determining the first inverse system corresponding to the minimum phase system and determining the second inverse system corresponding to the all-pass system includes: calculating the reciprocal of the minimum phase system to obtain the first inverse system corresponding to the minimum phase system, obtaining the finite impulse response system with delay corresponding to the all-pass system, and determining the finite impulse response system as the inverse system of the all-pass system to obtain the second inverse system.
[0013] Optionally, in one embodiment, obtaining the delayed finite impulse response system corresponding to the all-pass system includes: converting the all-pass system into a finite impulse response system; and performing a flipping process on the finite impulse response system to obtain the delayed finite impulse response system.
[0014] Optionally, in one embodiment, the compensation system is:
[0015]
[0016] Among them, H -1 (z) represents the compensation system, H aap (z) represents the all-pass system H ap (z) corresponds to the second inverse system, namely the finite impulse response system with time delay. For the minimum phase system H min (z) corresponds to the first inverse system.
[0017] Secondly, embodiments of this application provide a spatial response compensation device, comprising: a first construction module, configured to construct a spatial transfer function based on spatial input information and corresponding spatial response information of the system, wherein the spatial transfer function is used to indicate the mapping relationship between the input signal and the output signal of the system; a splitting module, configured to split the spatial transfer function into a minimum-phase system and an all-pass system; a determining module, configured to determine a first inverse system corresponding to the minimum-phase system and a second inverse system of the all-pass system function; a second construction module, configured to construct a compensation system based on a preset target frequency response, the first inverse system, and the second inverse system; a compensation module, configured to compensate the spatial transfer function through the compensation system to obtain a compensated spatial transfer function; and a conversion module, configured to convert the spatial input information of the system into target spatial response information through the compensated spatial transfer function, wherein the difference between the frequency response of the target spatial response information and the target frequency response is less than a preset threshold.
[0018] Thirdly, an electronic device is provided, comprising: a processor, a memory, a communication interface, and a communication bus, wherein the processor, the memory, and the communication interface communicate with each other via the communication bus; the memory is used to store at least one executable instruction, wherein the executable instruction causes the processor to perform an operation corresponding to the method described above.
[0019] Fourthly, a computer storage medium is provided, wherein a computer program is stored therein, which, when executed by a processor, implements the method described above.
[0020] Fifthly, a computer program product is provided, wherein the computer program product includes computer instructions that instruct a computing device to perform operations corresponding to the method described above.
[0021] As shown in the above scheme, by decomposing the spatial transfer function into a minimum-phase system and an all-pass system, and then obtaining their corresponding inverse systems, and finally combining the first and second inverse systems, the inverse system of any spatial response, i.e., the compensated system, is obtained. By decomposing any stable and causal linear time-invariant system into two systems with different properties, a minimum-phase system and an all-pass system, and then specifically solving for the inverse systems of the two systems, the complexity of the solution is greatly reduced, and the accuracy of the inverse system calculation is improved. Attached Figure Description
[0022] Figure 1A This is a schematic diagram of the impact response of a loudspeaker system in a room.
[0023] Figure 1B This is a schematic diagram of the frequency response of the loudspeaker system in the room;
[0024] Figure 2 This is a flowchart of a spatial response compensation method according to an embodiment of this application;
[0025] Figure 3A This is a schematic diagram of the amplitude-frequency response of the loudspeaker system in the room, including the minimum phase system and the first inverse system.
[0026] Figure 3B This is a schematic diagram of the phase frequency response of the loudspeaker system in the room, including the minimum phase system and the first inverse system.
[0027] Figure 3C This is a schematic diagram of the amplitude-frequency response of the loudspeaker system in the room's all-through system and the second inverse system;
[0028] Figure 3D This is a schematic diagram of the phase frequency response of the loudspeaker system in the room's all-through system and the second inverse system;
[0029] Figure 3E A schematic diagram showing the spatial transfer function of the loudspeaker system in the room and the frequency response of the compensation system;
[0030] Figure 3F This diagram illustrates the spatial transfer function of the loudspeaker system within the room and the group delay of the compensation system.
[0031] Figure 3G A schematic diagram showing the frequency response compensation results of the loudspeaker system in the room;
[0032] Figure 3H This is a schematic diagram showing the group delay compensation results of the speaker system in the room.
[0033] Figure 4 This is a schematic diagram of a spatial response compensation device according to an embodiment of this application;
[0034] Figure 5 This is a schematic diagram of an electronic device according to an embodiment of this application. Detailed Implementation
[0035] As used herein, "spatial response" can refer to the response of any system in space, and the response can be of any form, such as the impulse response of a loudspeaker system in a room. The spatial response compensation method described herein can be used to compensate for the spatial transfer function of any system. In some examples, the system can be a loudspeaker system in a room, using signals from multiple audio channels, or signals corresponding to multiple audio channels. In examples, the audio signals can include digital signals (containing information corresponding to multiple audio channels) and can include other information or metadata. The system can also be other systems such as an image system. Figure 1AThis is a schematic diagram of the impulse response of a loudspeaker system in a room, such as... Figure 1A As shown, the speaker system's response within the room is not ideal. Figure 1A In the diagram, "impulse response" refers to the impulse response, with the vertical axis representing amplitude and the horizontal axis representing time. The frequency response of this loudspeaker system in the room corresponding to its impulse response is shown below. Figure 1B As shown, it becomes more apparent that the amplitude and phase frequency responses of the loudspeaker system are not ideal. Figure 1B In this context, magnitude response (dB) and phase response represent the amplitude-frequency response and phase-frequency response. The vertical axis, magnitude, refers to the quantitative or qualitative magnitude or intensity, while the horizontal axis, frequency, refers to the frequency magnitude.
[0036] It should be noted that all systems mentioned in this article are stable and causal linear time-invariant systems.
[0037] The technical solution of the present invention will be described in detail below using a speaker system in a room as an example.
[0038] To make the technical problems, technical solutions and advantages of the present invention clearer, a detailed description will be given below in conjunction with the accompanying drawings and specific embodiments.
[0039] Spatial response compensation method
[0040] Figure 2 Here is a flowchart of a spatial response compensation method according to an embodiment of this application, as follows: Figure 2 As shown, the spatial response compensation method specifically includes the following steps:
[0041] Step 101: Construct a spatial transfer function based on the system's spatial input information and corresponding spatial response information.
[0042] In order to compensate the system's response in space to the target range, spatial input information is first input into the system, and then spatial response information is obtained. Based on the spatial input information and the corresponding spatial response information, a spatial transfer function can be constructed. The spatial transfer function is used to indicate the mapping relationship between the system's input signal and output signal.
[0043] For a loudspeaker system, the spatial input information can be a sound signal, such as an audio clip. Then, the spatial response information output by the loudspeaker system based on the input sound signal is collected. Since the spatial response information output by the loudspeaker system is also a sound signal, the output spatial response information can be collected using a sound signal collection device such as a microphone. Based on the sound signal input to the loudspeaker system and the sound signal output by the loudspeaker system, the spatial transfer function of the loudspeaker system can be constructed.
[0044] Step 102: Decompose the spatial transfer function into a minimum-phase system and an all-pass system.
[0045] After constructing the spatial transfer function of the system, since the system is a linear time-invariant system, and is stable and causal, the spatial transfer function can be decomposed into a minimum-phase system and an all-pass system.
[0046] A minimum-phase system is a system whose phase shift is minimal under certain amplitude-frequency characteristics.
[0047] A pass-through system is a system that allows all signals to pass through.
[0048] Step 103: Determine the first inverse system corresponding to the minimum phase system, and determine the second inverse system corresponding to the all-pass system.
[0049] After decomposing the minimum-phase system and the all-pass system, the first inverse system corresponding to the minimum-phase system and the second inverse system corresponding to the all-pass system are determined respectively.
[0050] An inverse system is a system that has the opposite input-output propagation relationship to a given system. The first inverse system is a system that has the opposite input-output propagation relationship to a minimum-phase system, and the second inverse system is a system that has the opposite input-output propagation relationship to an all-pass system.
[0051] For example, the frequency responses of the minimum phase system and the first inverse system of the above speaker system, decomposed from the spatial transfer function of the room, are as follows: Figure 3A and Figure 3B As shown, Figure 3A This is a schematic diagram of the amplitude-frequency response of the loudspeaker system in the room, including the minimum phase system and the first inverse system. Figure 3A In this context, magnitude response refers to the amplitude-frequency response. The vertical axis, magnitude, indicates the quantitative or qualitative magnitude or intensity, while the horizontal axis, frequency, indicates the frequency magnitude. Figure 3B This is a schematic diagram of the phase frequency response of the loudspeaker system in the room, including the minimum phase system and the first inverse system. Figure 3BIn the diagram, "phase response" represents the phase frequency response. The vertical axis "phase" (radians) refers to the phase shift value, and the horizontal axis "frequency" refers to the frequency magnitude. The frequency responses of the above loudspeaker system, decomposed from the spatial transfer function of the room, for the all-pass system and the second inverse system are as follows: Figure 3C and Figure 3D As shown, Figure 3C This is a schematic diagram of the amplitude-frequency response of the loudspeaker system in the all-through system and the second inverse system of the room. Figure 3C In this context, magnitude response refers to the amplitude-frequency response. The vertical axis, magnitude, indicates the quantitative or qualitative magnitude or intensity, while the horizontal axis, frequency, indicates the frequency magnitude. Figure 3D This is a schematic diagram of the phase frequency response of the loudspeaker system in the all-through system and the second inverse system of the room. Figure 3D In the diagram, "group delay" represents the group delay, the vertical axis "group delay (in samples)" represents the sampling frequency of the group delay, and the horizontal axis "Frequency" indicates the frequency magnitude. It can be seen from the above phase changes that the phase frequency response of the all-pass system and the second inverse system is symmetrical from top to bottom, with only a certain delay in between.
[0052] Step 104: Construct a compensation system based on the preset target frequency response, the first inverse system, and the second inverse system.
[0053] After determining the first and second inverse systems, a compensation system is constructed by combining the pre-set target frequency response.
[0054] For example, the preset target frequency response can be 0dB, which indicates that the target amplitude frequency response and target phase frequency response of the output signal corresponding to the input signal of the space transfer function after compensation by the compensation system are 0dB.
[0055] For example, the spatial transfer function of the aforementioned loudspeaker system in the room and the frequency response of the compensation system are as follows: Figure 3E As shown, Figure 3E In the graph, magnitude response (dB) and phase response represent the amplitude-frequency response and phase-frequency response, respectively. The vertical axis, magnitude, refers to the quantitative or qualitative magnitude or intensity, while the horizontal axis, frequency, refers to the frequency. It can be seen that the compensation system is constructed based on the preset target frequency response, the first inverse system, and the second inverse system. Figure 3F This diagram illustrates the spatial transfer function of a loudspeaker system within a room and the group delay of a compensation system. Group delay is a quantity that describes how quickly phase changes with frequency. Figure 3FIn the figure, group delay represents the group delay, the vertical axis group delay(insamples) represents the sampling frequency of the group delay, and the horizontal axis Frequency refers to the frequency magnitude.
[0056] It should be noted that the preset target frequency response can be a single value, in which case the target amplitude frequency response and the target phase frequency response are equal, or it can be two values, respectively indicating the target values corresponding to the target amplitude frequency response and the target phase frequency response.
[0057] Step 105: Compensate the spatial transfer function using a compensation system to obtain the compensated spatial transfer function.
[0058] After constructing the compensation system, the spatial transfer function is compensated using the compensation system to obtain the compensated spatial transfer function.
[0059] Step 106: Convert the spatial input information of the system into the target spatial response information through the compensated spatial transfer function.
[0060] After compensating for the spatial transfer function, the system receives spatial input information and converts it into target spatial response information. The difference between the frequency response of the target spatial response and the target frequency response is less than a preset threshold.
[0061] The preset threshold can be a single value, such as 0.01dB, in which case the difference between the frequency response of the target spatial response information and the target frequency response is less than 0.01dB. The preset threshold can also be two values, such as 0.01dB and 0.02dB, in which case the difference between the amplitude frequency response of the target spatial response information and the amplitude frequency response included in the target frequency response is less than 0.01dB, and the difference between the phase frequency response of the target spatial response information and the phase frequency response included in the target frequency response is less than 0.02dB.
[0062] For example, the frequency response compensation result of the loudspeaker system in the room can be as follows: Figure 3G As shown, Figure 3G In this context, magnitude response (dB) and phase response represent the amplitude-frequency response and phase-frequency response, respectively. The vertical axis, magnitude, refers to the quantitative or qualitative magnitude or intensity, while the horizontal axis, frequency, refers to the frequency range. The group delay compensation result of the loudspeaker system in the room can be expressed as follows: Figure 3H As shown, Figure 3HIn the figure, group delay represents the group delay, the vertical axis groupdelay(in samples) represents the sampling frequency of the group delay, and the horizontal axis Frequency refers to the frequency magnitude. It can be seen that the frequency response of the target spatial response information gradually approaches the target frequency response, and finally the difference between the two is less than the preset threshold.
[0063] In this embodiment, the spatial transfer function is decomposed into a minimum-phase system and an all-pass system, and then their corresponding inverse systems are obtained respectively. Finally, the first and second inverse systems are combined to obtain the inverse system of any spatial response, i.e., the compensation system. By decomposing any stable and causal linear time-invariant system into two systems with different properties, a minimum-phase system and an all-pass system, and specifically solving the inverse systems of the two systems, the complexity of the solution is greatly reduced and the accuracy of the inverse system calculation is improved.
[0064] In one possible implementation, the process of decomposing the spatial transfer function into a minimum-phase system and an all-pass system may also include:
[0065] The all-pass system is determined by the external zeros of the space transfer function, and the minimum-phase system is determined by the internal zeros, external zeros, and the first internal pole of the space transfer function.
[0066] Specifically, the spatial transfer function can be decomposed by dividing the zeros of the spatial transfer function into internal zeros, external zeros, and internal poles. Then, the all-pass system is determined based on the external zeros of the spatial transfer function, and the minimum-phase system is determined based on the internal zeros, external zeros, and the first internal pole of the spatial transfer function.
[0067] In the embodiments of this application, the all-pass system and the minimum-phase system are determined by the zeros inside the circle, the zeros outside the circle, and the poles inside the circle of the spatial transfer function, which can reduce the complexity of the solution.
[0068] In one possible implementation, the process of determining the all-pass system based on the external zeros of the spatial transfer function may also include:
[0069] Calculate the reciprocal of the zeros outside the circle to obtain the poles inside the second circle, and then determine the all-pass system based on the poles inside the second circle.
[0070] All zeros of an all-pass system are located outside the unit circle, making it a maximum-phase system. Furthermore, since all zero-pole pairs of an all-pass system are complex conjugates in the Z-plane, the all-pass system can be determined simply by identifying the poles inside the circle corresponding to the all-pass system.
[0071] In this embodiment of the application, by directly reciprocating the zeros outside the circle to obtain the poles inside the second circle, and determining the all-pass system based on the poles inside the second circle, the computational complexity of the all-pass system can be further reduced.
[0072] In one possible implementation, the process of determining the minimum-phase system based on the zeros inside and outside the first circle of the spatial transfer function, and the poles inside the first circle, may further include:
[0073] Calculate the reciprocal of the zeros outside the circle to obtain the zeros inside the second circle. Then, determine the minimum phase system based on the zeros inside the first circle, the zeros inside the second circle, and the poles inside the first circle.
[0074] The zeros and poles of the minimum-phase system are all inside the unit circle. Therefore, by reciprocating the zeros outside the circle, we can obtain the zeros inside the second circle. Then, the minimum-phase system can be determined based on the zeros of the spatial transfer function inside the unit circle and the zeros inside the second circle.
[0075] In this embodiment, by directly reciprocating the zeros outside the circle to obtain the zeros inside the second circle, and determining the minimum phase system based on the zeros inside the first circle, the zeros inside the second circle, and the poles inside the first circle, the complexity of the minimum phase system calculation can be further reduced.
[0076] In one possible implementation, the process of determining the first inverse system corresponding to the minimum-phase system and the second inverse system corresponding to the all-pass system may further include:
[0077] Calculate the reciprocal of the minimum-phase system to obtain the first inverse system corresponding to the minimum-phase system, and obtain the finite impulse response system with delay corresponding to the all-pass system. Then, determine the finite impulse response system as the inverse system of the all-pass system to obtain the second inverse system.
[0078] For a minimum-phase system, since it has no zeros outside the circle, its inverse system is easy to obtain. The first inverse system corresponding to the minimum-phase system can be obtained by taking the reciprocal. However, for an all-pass system, since it has zeros outside the circle, the inverse system cannot be obtained by taking the reciprocal directly. The all-pass system can be compensated by a finite impulse response (FIR) system with a delay. That is, there exists a finite impulse response system with a delay that is the inverse system of the all-pass system. Therefore, the finite impulse response system is determined as the inverse system of the all-pass system, and the second inverse system is obtained.
[0079] In this embodiment, the computational complexity is reduced by calculating the reciprocal of the minimum phase system to solve the inverse system. By determining the finite impulse response system with delay as the second inverse system of the all-pass system, the problem of the difficulty in solving the inverse system of the all-pass system can be solved.
[0080] In one possible implementation, the process of obtaining the time-delayed finite impulse response system corresponding to the all-pass system may also include:
[0081] The all-pass system is converted into a finite impulse response system, and the finite impulse response system is flipped to obtain a finite impulse response system with delay.
[0082] One specific way to solve for the second inverse system corresponding to an all-pass system is to first convert the all-pass system into a finite impulse response system. Since the all-pass system is an infinite impulse response system (IIR), it needs to be converted into a finite impulse response system first. Then, the obtained finite impulse response system is flipped to finally obtain a finite impulse response system (FIR) with delay.
[0083] In this embodiment, the computational complexity is reduced by first converting the all-pass system into a finite impulse response system and then solving the finite impulse response system with delay.
[0084] In one possible implementation, the compensation system can be:
[0085]
[0086] H -1 (z) is a compensation system, H aap (z) represents the all-pass system H ap (z) corresponds to the second inverse system, namely the finite impulse response system with time delay. For the minimum phase system H min The first inverse system corresponding to (z).
[0087] In the embodiments of this application, the relationship between the compensation system, the minimum phase system, and the all-pass system, and the first inverse system and the second inverse system can be understood more intuitively, thereby improving the accuracy of the compensation system calculation.
[0088] Spatial response compensation device
[0089] Figure 4 This is a schematic diagram of a spatial response compensation device according to an embodiment of this application. Figure 4 As shown, the spatial response compensation device 400 includes: a first construction module 401, a splitting module 402, a determining module 403, a second construction module 404, a compensation module 405, and a conversion module 406.
[0090] The first construction module 401 is used to construct a spatial transfer function based on the spatial input information and the corresponding spatial response information of the system. The spatial transfer function is used to indicate the mapping relationship between the input signal and the output signal of the system.
[0091] In order to compensate the system's response in space to the target range, spatial input information is first input into the system, and then spatial response information is obtained. The first construction module 401 can construct a spatial transfer function based on the spatial input information and the corresponding spatial response information. The spatial transfer function is used to indicate the mapping relationship between the system's input signal and output signal.
[0092] For a loudspeaker system, the spatial input information can be a sound signal, such as an audio clip. Then, the spatial response information output by the loudspeaker system based on the input sound signal is collected. Since the spatial response information output by the loudspeaker system is also a sound signal, the output spatial response information can be collected using a sound signal collection device such as a microphone. Based on the sound signal input to the loudspeaker system and the sound signal output by the loudspeaker system, the spatial transfer function of the loudspeaker system can be constructed.
[0093] The splitting module 402 is used to split the spatial transfer function into a minimum-phase system and an all-pass system.
[0094] After the first building module 401 constructs the spatial transfer function of the system, since the system is a linear time-invariant system and is stable and causal, the splitting module 402 can split the spatial transfer function into a minimum-phase system and an all-pass system.
[0095] A minimum-phase system is a system whose phase shift is minimal under certain amplitude-frequency characteristics.
[0096] A pass-through system is a system that allows all signals to pass through.
[0097] Module 403 is used to determine the first inverse system corresponding to the minimum phase system and the second inverse system of the all-pass system function.
[0098] After the splitting module 402 decomposes the minimum phase system and the all-pass system, the determining module 403 determines the first inverse system corresponding to the minimum phase system and the second inverse system corresponding to the all-pass system, respectively.
[0099] The second construction module 404 is used to construct a compensation system based on a preset target frequency response, a first inverse system, and a second inverse system.
[0100] After the determining module 403 determines the first inverse system and the second inverse system, the second construction module 404 constructs a compensation system by combining the pre-set target frequency response.
[0101] For example, the preset target frequency response can be 0dB, which indicates that the target amplitude frequency response and target phase frequency response of the output signal corresponding to the input signal of the space transfer function after compensation by the compensation system are 0dB.
[0102] It should be noted that the preset target frequency response can be a single value, in which case the target amplitude frequency response and the target phase frequency response are equal, or it can be two values, respectively indicating the target values corresponding to the target amplitude frequency response and the target phase frequency response.
[0103] The compensation module 405 is used to compensate the space transfer function through the compensation system to obtain the compensated space transfer function.
[0104] After the second construction module 404 constructs the compensation system, the compensation module 405 compensates the spatial transfer function through the compensation system to obtain the compensated spatial transfer function.
[0105] The conversion module 406 is used to convert the spatial input information of the system into the target spatial response information through the compensated spatial transfer function, wherein the difference between the frequency response of the target spatial response information and the target frequency response is less than a preset threshold.
[0106] After the compensation module 405 compensates for the spatial transfer function, the conversion module 406 inputs spatial input information into the system, and the system converts the spatial input information into target spatial response information. The difference between the frequency response of the target spatial response information and the target frequency response is less than a preset threshold.
[0107] The preset threshold can be a single value, such as 0.01dB, in which case the difference between the frequency response of the target spatial response information and the target frequency response is less than 0.01dB. The preset threshold can also be two values, such as 0.01dB and 0.02dB, in which case the difference between the amplitude frequency response of the target spatial response information and the amplitude frequency response included in the target frequency response is less than 0.01dB, and the difference between the phase frequency response of the target spatial response information and the phase frequency response included in the target frequency response is less than 0.02dB.
[0108] In this embodiment, the spatial transfer function is decomposed into a minimum-phase system and an all-pass system by the decomposition module 402. Then, the determination module 403 obtains the corresponding inverse systems. Finally, the second construction module 404 combines the first and second inverse systems to obtain the inverse system of arbitrary spatial response, i.e., the compensation system. By decomposing any stable and causal linear time-invariant system into two systems with different properties, a minimum-phase system and an all-pass system, and specifically solving the inverse systems of the two systems, the complexity of the solution is greatly reduced, and the accuracy of the inverse system calculation is improved.
[0109] electronic devices
[0110] Figure 5 This is a schematic diagram of an electronic device according to one embodiment of this application. The specific embodiments of this application do not limit the specific implementation of the electronic device. Figure 5 As shown, the electronic device 500 may include: a processor 501, a communications interface 502, a memory 503, and a communications bus 504. Wherein:
[0111] The processor 501, communication interface 502, and memory 503 communicate with each other through the communication bus 504.
[0112] Communication interface 502 is used for communication with other electronic devices or servers.
[0113] The processor 501 is used to execute program 505, which can specifically execute the relevant steps in any of the aforementioned method embodiments.
[0114] Specifically, program 505 may include program code that includes computer operation instructions.
[0115] The processor 501 may be a CPU, an application-specific integrated circuit (ASIC), or one or more integrated circuits configured to implement the embodiments of this application. The smart device includes one or more processors, which may be processors of the same type, such as one or more CPUs; or processors of different types, such as one or more CPUs and one or more ASICs.
[0116] Memory 503 is used to store program 505. Memory 503 may include high-speed RAM memory, and may also include non-volatile memory, such as at least one disk storage device.
[0117] Specifically, program 505 can be used to cause processor 501 to execute any of the methods in the foregoing embodiments.
[0118] The specific implementation of each step in program 505 can be found in the corresponding steps and units described in the aforementioned spatial response compensation method embodiments, and will not be repeated here. Those skilled in the art will understand that, for the sake of convenience and brevity, the specific working processes of the devices and modules described above can be referred to the corresponding process descriptions in the aforementioned method embodiments, and will not be repeated here.
[0119] The electronic device of this application decomposes the spatial transfer function into a minimum-phase system and an all-pass system, then obtains their corresponding inverse systems, and finally combines the first and second inverse systems to obtain an inverse system with arbitrary spatial response, i.e., a compensation system. By decomposing any stable and causal linear time-invariant system into two systems with different properties, a minimum-phase system and an all-pass system, and specifically solving the inverse systems of the two systems, the complexity of the solution is greatly reduced and the accuracy of the inverse system calculation is improved.
[0120] Computer storage media
[0121] This application also provides a computer-readable storage medium storing instructions for causing a machine to perform any of the methods described in the plurality of method embodiments herein. Specifically, a system or apparatus equipped with a storage medium storing software program code that implements the functions of any of the embodiments described above, and enabling the computer (or CPU or MPU) of the system or apparatus to read and execute the program code stored in the storage medium.
[0122] In this case, the program code read from the storage medium can itself implement the function of any of the above embodiments, and therefore the program code and the storage medium storing the program code constitute part of this application.
[0123] Examples of storage media used to provide program code include floppy disks, hard disks, magneto-optical disks, optical disks (such as CD-ROM, CD-R, CD-RW, DVD-ROM, DVD-RAM, DVD-RW, DVD+RW), magnetic tapes, non-volatile memory cards, and ROMs. Alternatively, program code can be downloaded from a server computer via a communication network.
[0124] Computer program products
[0125] This application also provides a computer program product, including computer instructions that instruct a computing device to perform any corresponding operation in the above-described plurality of method embodiments.
[0126] It should be noted that, depending on the implementation needs, the various components / steps described in the embodiments of this application can be broken down into more components / steps, or two or more components / steps or parts of the operation of components / steps can be combined into new components / steps to achieve the purpose of the embodiments of this application.
[0127] The methods described in the embodiments of this application can be implemented in hardware, firmware, or as software or computer code that can be stored in a recording medium (such as a CD-ROM, RAM, floppy disk, hard disk, or magneto-optical disk), or as computer code downloaded over a network that is originally stored in a remote recording medium or a non-transitory machine-readable medium and will be stored in a local recording medium. Thus, the methods described herein can be processed by software stored on a recording medium using a general-purpose computer, a dedicated processor, or programmable or dedicated hardware (such as an ASIC or FPGA). It is understood that the computer, processor, microprocessor controller, or programmable hardware includes storage components (e.g., RAM, ROM, flash memory, etc.) capable of storing or receiving software or computer code that, when accessed and executed by the computer, processor, or hardware, implements the methods described herein. Furthermore, when a general-purpose computer accesses code used to implement the methods shown herein, the execution of the code transforms the general-purpose computer into a dedicated computer for executing the methods shown herein.
[0128] It should be noted that not all steps and modules in the above processes and system structure diagrams are necessary; some steps or modules can be omitted as needed. The execution order of each step is not fixed and can be adjusted as required. The system structure described in the above embodiments can be a physical structure or a logical structure. That is, some modules may be implemented by the same physical entity, or some modules may be implemented by multiple physical entities, or they may be jointly implemented by certain components in multiple independent devices.
[0129] In the above embodiments, the hardware modules can be implemented mechanically or electrically. For example, a hardware module may include permanent dedicated circuitry or logic (such as a dedicated processor, FPGA, or ASIC) to perform the corresponding operation. The hardware module may also include programmable logic or circuitry (such as a general-purpose processor or other programmable processor), which can be temporarily configured by software to perform the corresponding operation. The specific implementation method (mechanical, dedicated permanent circuitry, or temporarily configured circuitry) can be determined based on cost and time considerations.
[0130] The present application has been shown and described in detail above with reference to the accompanying drawings and preferred embodiments. However, the present application is not limited to these disclosed embodiments. Based on the above multiple embodiments, those skilled in the art will know that more embodiments of the present application can be obtained by combining the code review methods in the different embodiments above. These embodiments are also within the protection scope of the present application.
Claims
1. A spatial response compensation method, characterized in that, include: Based on the spatial input information and the corresponding spatial response information of the system, a spatial transfer function is constructed, wherein the spatial transfer function is used to indicate the mapping relationship between the input signal and the output signal of the system; The spatial transfer function is decomposed into a minimum-phase system and an all-pass system; Determine the first inverse system corresponding to the minimum phase system, and determine the second inverse system corresponding to the all-pass system; A compensation system is constructed based on the preset target frequency response, the first inverse system, and the second inverse system. The space transfer function is compensated by the compensation system to obtain the compensated space transfer function; The spatial input information of the system is converted into target spatial response information through the compensated spatial transfer function, wherein the difference between the frequency response of the target spatial response information and the target frequency response is less than a preset threshold. The step of splitting the spatial transfer function into a minimum-phase system and an all-pass system includes: The all-pass system is determined based on the external zeros of the spatial transfer function; The minimum phase system is determined based on the first circle zeros, the outside zeros, and the first circle poles of the spatial transfer function; Determining the all-pass system based on the external zeros of the spatial transfer function includes: Calculate the reciprocal of the zero point outside the circle to obtain the pole inside the second circle; The all-pass system is determined based on the inner pole of the second circle; Determining the minimum phase system based on the first circle zeros, the outside zeros, and the first circle poles of the spatial transfer function includes: Calculate the reciprocal of the outer zero point to obtain the inner zero point of the second circle; The minimum phase system is determined based on the zero point inside the first circle, the zero point inside the second circle, and the pole inside the first circle.
2. The method according to claim 1, characterized in that, Determining the first inverse system corresponding to the minimum-phase system and determining the second inverse system corresponding to the all-pass system includes: Calculate the reciprocal of the minimum phase system to obtain the first inverse system corresponding to the minimum phase system; Obtain the finite impulse response system with delay corresponding to the all-pass system, and determine the finite impulse response system as the inverse system of the all-pass system to obtain the second inverse system.
3. The method according to claim 2, characterized in that, The method for obtaining the finite impulse response system with delay corresponding to the all-pass system includes: The all-pass system is converted into a finite impulse response system; The finite impulse response system is flipped to obtain the finite impulse response system with delay.
4. The method according to claim 2 or 3, characterized in that, The compensation system is as follows: ； in, For the compensation system, For the all-pass system The inverse system, namely the finite impulse response system with time delay, For the minimum phase system The inverse system.
5. A spatial response compensation device, characterized in that, include: The first construction module is used to construct a spatial transfer function based on the spatial input information and the corresponding spatial response information of the system, wherein the spatial transfer function is used to indicate the mapping relationship between the input signal and the output signal of the system; The splitting module is used to split the spatial transfer function into a minimum-phase system and an all-pass system; The determining module is used to determine the first inverse system corresponding to the minimum phase system and to determine the second inverse system of the all-pass system function; The second construction module is used to construct a compensation system based on the preset target frequency response, the first inverse system, and the second inverse system; The compensation module is used to compensate the spatial transfer function through the compensation system to obtain the compensated spatial transfer function; The conversion module is used to convert the spatial input information of the system into target spatial response information through the compensated spatial transfer function, wherein the difference between the frequency response of the target spatial response information and the target frequency response is less than a preset threshold. The splitting module is also used to determine the all-pass system based on the external zeros of the spatial transfer function, and to determine the minimum-phase system based on the first internal zeros, external zeros, and first internal poles of the spatial transfer function. The splitting module is also used to calculate the reciprocal of the zero point outside the circle to obtain the pole inside the second circle, and to determine the all-pass system based on the pole inside the second circle; The splitting module is also used to calculate the reciprocal of the outer zero point to obtain the second inner zero point, and to determine the minimum phase system based on the first inner zero point, the second inner zero point, and the first inner pole.
6. An electronic device, comprising: The processor, communication interface, memory, and communication bus communicate with each other through the communication bus. The memory is used to store at least one executable instruction that causes the processor to perform the operation corresponding to the space response compensation method as described in any one of claims 1-4.
7. A computer storage medium having a computer program stored thereon, which, when executed by a processor, implements the spatial response compensation method as described in any one of claims 1-4.
8. A computer program product comprising computer instructions that instruct a computing device to perform an operation corresponding to the spatial response compensation method as described in any one of claims 1-4.