CN103826194B - Method and device for rebuilding sound source direction and distance in multichannel system - Google Patents

Abstract

The invention discloses a method and device for rebuilding a sound source direction and distance in a multichannel system. The particle speed and sound pressure of a sound source received by a sound listening point are calculated according to known sound source signals, the distance between the sound listening point and the sound source is (i)r(/i), four loudspeakers in the direction are selected according to the particle speed, the sound source signals are multiplied by four weight factors respectively and then are distributed to the four selected loudspeakers, then particle speed and sound pressure of a sound image inducted at the sound listening point are calculated at a playback end after the loudspeakers in a rebuilt sound field send signals, finally, direction and distance equivalent models are built according to the particle speed and sound pressure of the original sound source and the particle speed and sound pressure of the sound image in the rebuilt sound field, and the signals of the loudspeakers are distributed through the weight factors acquired according to a solution model. Compared with the prior art, the direction and distance information of the sound source in the original sound source space can be accurately restored, the operation is simple, calculation efficiency is high, and stability is good.

Description

The method and apparatus that in a kind of multi-channel system, Sounnd source direction and distance are rebuild

Technical field

The invention belongs to field of multimedia signal processing, to particularly relate in Audio Signal Processing direction the directions such as sound source reconstruction, sound field rebuilding, sound source dimensional orientation and range recovery, be specifically related to the method and apparatus that in a kind of multi-channel system, Sounnd source direction and distance are rebuild.

Background technology

Although the tonequality of stereophony and sound field effect are better than monophony greatly, it still has obvious limitation.Dual-channel stereo system can only reproduce the sound field in the sector region of one, people front, and the sound reset can not be allowed to give Sensurround.So multichannel technology also starts to grow up.Multi-sound channel digital audio system, by the expansion of number of channels, can recover multiple ambient sound, brings better feeling of immersion and higher-quality audio frequency to enjoy, surmounted the effect of traditional monophony and stereophonic sound system far away to audience.Along with the development of the storage mediums such as DVD, SACD, multichannel audio has progressively been entered multiple fields of people's life by initial exclusive cinema.

Typical multi-channel system has Dolby Digital AC-3 multichannel surround sound sound system, THX Surround EX system and 22.2 multichannel stereo systems.

Dolby digital audio compression standard (AC-3) is the international standard proposed, and it solves digital sound and simulated sound on a film and the problem of depositing.Dolby Digital AC-3 multichannel surround sound sound system by six independently sound channel form.Because the Hz-KHz of the first five separate channels in these six separate channels is all audio frequency Whole frequency band and 20Hz ~ 20kHz, and subwoofer sound channel Hz-KHz only has 15Hz ~ 150Hz, only account for 1/10th of whole frequency spectrum, therefore Dolby Digital AC-3 multichannel surround sound sound system is also called 5.l sound channel ambiophonic system.

THX Surround EX system has strictly been worked out cinema and to be correlated with the standard of audio-visual equipment and environment, as long as meet THX standard and through certification, just can have suitable level.As long as consumer selects the movie theatre with THX certification like this, just have excellent audio-visual enjoyment.THX was transplanted to home theater afterwards, for the seehears of certification high-quality, and had unique requirement for the difference of home environment.But THX is not Dolby Digital and DTS is a kind of audio format like that, but a kind of audio post-processing pattern, object obtains best seeing and hearing enjoyment.When the Dolby Digital EX of 6.1 sound channels and DTS ES out after, it is evolved into THX Surround EX system by THX further.In order to the former Bidirectional sounding of compatibility side sound channel and reinforcing ring is around sound effect Sensurround once again, so sound channel turn increase two on the basis of former side sound channel after, this just constitutes 7.1 sound channels.

22.2 multichannel stereo systems are developed as " ultra high-definition " image sound equipment by NHK.By the loud speaker of upper strata 9 sound channel, middle level 10 sound channel, lower floor 3 sound channel three layers configuration, and dual track low frequency audio (LFE) loud speaker, the sound that left and right and above-below direction are forwards, backwards propagated can be reproduced comparatively truly.The sound technique of 22.2 sound channels needs, with the position of audience's ear level, lay 10,9 and 3 loudspeaker respectively above and below the position of audience's ear, to be furnished with 2 subwoofer be made up of 36 little bass units in addition.Whole ambiophonic system forms by three layers, and the sound field of the bottom is made up of three front channels and two LFE sound channels.Middle level sound field is made up of surround channel after four front channels, two side surround channels and three.Upper strata sound field is by three front channels, two side surround channels, surround channel is added a top sound channel and formed after four.In addition, in main sound channel, in order to strengthen effect, the dynamics sense that two trumpet arrays be made up of 36 toy trumpets export to guarantee main sound channel can also be selected.

5.1 ambiophonic systems, 22.2 multi-channel systems can bring good Sensurround to audience, to the reconstruction of sound all based on the process of signal itself, gratifying effect is achieved in tonequality fidelity, but the recovery lacked the orientation of sound source, cannot rebuild the positional information of holding and recovering sound source.

Summary of the invention

The deficiency that Sounnd source direction and range information exist is being rebuild in order to overcome existing multi-channel system, the invention provides the method and apparatus that in a kind of multi-channel system, Sounnd source direction and distance are rebuild, can accurately recover Sounnd source direction and distance in multi-channel system.

The technical scheme that method of the present invention adopts is: the method and apparatus that in a kind of multi-channel system, Sounnd source direction and distance are rebuild, is characterized in that, comprise the following steps:

Step 1: according to time-domain signal s (t) of known sound source, to listen the point of articulation for after initial point sets up Descartes's rectangular coordinate system in original sound field, the particle rapidity pv of the sound source of listening point of articulation place to receive that to calculate with sound source distance be r ₀with acoustic pressure p ₀;

Step 2: according to the particle rapidity pv of the sound source calculated in step 1 ₀, select to comprise particle rapidity pv ₀4 loud speaker L in direction ₁, L ₂, L ₃, L ₄;

Step 3: by the loud speaker L in sound-source signal s (t) and step 2 ₁, L ₂, L ₃, L ₄, s (t) is multiplied by respectively 4 weight w ₁, w ₂, w ₃, w ₄after be assigned on four selected loud speakers;

Step 4: at playback end, calculates and rebuilds loud speaker L in sound field ₁, L ₂, L ₃, L ₄particle rapidity pv and the acoustic pressure p of the acoustic image of point of articulation place perception is listened after sending signal;

Step 5: according to the particle rapidity pv of the original sound source that step 1 calculates ₀with acoustic pressure p ₀, the particle rapidity pv of acoustic image and acoustic pressure p in the reconstruction sound field that calculates in step 4, set up direction, distance equivalence model;

Step 6: in the direction of step 5, apart from equivalence model, with the weight w arranged in step 3 ₁, w ₂, w ₃, w ₄for unknown quantity, solve the value of direction, distance equivalence model acquisition weight;

Step 7: utilize the weight w solved in step 6 ₁, w ₂, w ₃, w ₄, the signal carrying out loud speaker distributes, and makes the audio direction of reconstruction and distance with the direction of acoustic source with apart from consistent.

As preferably, the particle rapidity pv of the sound source of listening point of articulation place to receive that the calculating described in step 1 and sound source distance are r ₀with acoustic pressure p _0,its specific implementation comprises following sub-step:

Step 1.1: utilize Fourier transform that time-domain signal s (t) of sound source is transformed into frequency domain by time domain, obtains frequency-region signal s (ω)

$s\left(\mathrm{\ω}\right)=\underset{-\∞}{\overset{+\∞}{\∫}}s\left(t\right)e^{-\mathrm{iwt}}\mathrm{dt}$

Step 1.2: by frequency-region signal s (ω) and the sound source position vector ε in sound field=(ε _x, ε _y, ε _z), listen point of articulation position coordinates r=r (r _x, r _y, r _z), according to the definition of sound physical properties attribute particle rapidity, calculate the particle rapidity pv listening point of articulation place to receive ₀:

${\mathrm{pv}}_0(r,\mathrm{\ω})=G\frac{e^{-\mathrm{ik}|r-\mathrm{\ϵ}|}}{|r-\mathrm{\ϵ}|}(1+\frac{1}{\mathrm{ik}|r-\mathrm{\ϵ}|})\×\frac{1}{|r-\mathrm{\ϵ}|}\left(\begin{array}{c}{r}_x-{\mathrm{\ϵ}}_x\\ r_y-{\mathrm{\ϵ}}_y\\ r_z-{\mathrm{\ϵ}}_z\end{array}\right)s\left(\mathrm{\ω}\right);$

Calculate the acoustic pressure p listening point of articulation place to receive ₀:

$p_0(T_0,\mathrm{\ω})=G\frac{e^{-\mathrm{ik}|r-\mathrm{\ϵ}|}}{|r-\mathrm{\ϵ}|}s\left(\mathrm{\ω}\right);$

Wherein k represents wave number, relevant with the velocity of sound with the frequency of voice signal, and e is constant, and i is imaginary unit, and ω is the angular frequency of signal s (t), | r-ε | represent sound source and listen the distance between the point of articulation, G (ω) represents source strength.

As preferably, the loud speaker L of 4 described in step 2 ₁, L ₂, L ₃, L ₄, meet following condition: point of articulation r (x, y, z) and L will be listened respectively ₁, L ₂, L ₃, L ₄connect, in the polyhedron formed, the particle rapidity pv of sound source ₀unit vector to be positioned at this polyhedron inner.

As preferably, loud speaker L in sound field is rebuild in the calculating described in step 4 ₁, L ₂, L ₃, L ₄listen particle rapidity pv and the acoustic pressure p of the acoustic image of point of articulation place perception after sending signal, its specific implementation comprises following sub-step:

Step 4.1: at playback end, with the center of the ambiophonic system of multichannel composition for initial point, sets up cartesian coordinate system, loud speaker L _icoordinate be designated as L _i(L _ix, L _iy, L _iz), listen the coordinate of the point of articulation to be r (r _x, r _y, r _z);

Step 4.2: by loud speaker L ₁, L ₂, L ₃, L ₄the signal q sent ₁(t), q ₂(t), q ₃(t), q ₄t () is converted to frequency domain:

$q\left(\mathrm{\ω}\right)=\left(\begin{array}{c}{q}_1\left(\mathrm{\ω}\right)\\ q_2\left(\mathrm{\ω}\right)\\ q_3\left(\mathrm{\ω}\right)\\ q_4\left(\mathrm{\ω}\right)\end{array}\right),q_i\left(\mathrm{\ω}\right)=\underset{-\∞}{\overset{+\∞}{\∫}}{q}_i\left(t\right)e^{-\mathrm{iwt}}\mathrm{dt}$

Step 4.3: by loud speaker L ₁, L ₂, L ₃, L ₄send signal q ₁(t), q ₂(t), q ₃(t), q ₄t () calculates the frequency domain representation pv and the acoustic pressure p that listen the particle rapidity at point of articulation place:

$\begin{array}{c}{p}_v(r,\mathrm{\ω})=\underset{i=1}{\overset{3}{\mathrm{\Σ}}}G\frac{e^{-\mathrm{ik}|r-L_i|}}{|r-L_i|}(1+\frac{1}{\mathrm{ik}|r-L_i|})\×\frac{1}{|r-L_i|}\left(\begin{array}{c}{r}_x-L_{\mathrm{ix}}\\ r_y-L_{\mathrm{iy}}\\ r_z-L_{\mathrm{iz}}\end{array}\right)q\left(\mathrm{\ω}\right)\\ =\underset{i=1}{\overset{3}{\mathrm{\Σ}}}G\frac{e^{-\mathrm{ik}|r-L_i|}}{{|r-L_i|}^2}(1+\frac{1}{\mathrm{ik}|r-L_i|})\×\left(\begin{array}{c}{r}_x-L_{\mathrm{ix}}\\ r_y-L_{\mathrm{iy}}\\ r_z-L_{\mathrm{iz}}\end{array}\right)q\left(\mathrm{\ω}\right)\end{array}$

$p=G\underset{j=1}{\overset{n}{\mathrm{\Σ}}}\frac{e^{-\mathrm{ik}|r-L_i|}}{|r-L_i|}q\left(\mathrm{\ω}\right)$

Wherein k represents wave number, relevant with the velocity of sound with the frequency of voice signal, and e is constant, and i is imaginary unit, and ω is the angular frequency of signal s (t), | r-L _i| represent loud speaker and listen the distance between the point of articulation, G (ω) represents source strength.

As preferably, the particle rapidity pv of the original sound source calculated according to step 1 described in step 5 ₀with acoustic pressure p ₀, the particle rapidity pv of acoustic image and acoustic pressure p in the reconstruction sound field that calculates in step 4, set up direction, distance equivalence model, its specific implementation comprises following sub-step:

Step 5.1: according to the particle rapidity p of the original sound source that step 1 calculates _v0particle rapidity pv with acoustic image in the reconstruction sound field calculated in step 4, sets up direction equivalence relation as follows:

pv＝pv ₀

Namely

$\underset{i=1}{\overset{3}{\mathrm{\Σ}}}G\frac{e^{-\mathrm{ik}|r-L_i|}}{{|r-L_i|}^2}(1+\frac{1}{\mathrm{ik}|r-L_i|})\×\left(\begin{array}{c}{r}_{\mathrm{ix}}-L_{\mathrm{ix}}\\ r_{\mathrm{iy}}-L_{\mathrm{iy}}\\ r_{\mathrm{iz}}-L_{\mathrm{iz}}\end{array}\right)q\left(\mathrm{\ω}\right)=G\frac{e^{-\mathrm{ik}|r-\mathrm{\ϵ}|}}{|r-\mathrm{\ϵ}|}(1+\frac{1}{\mathrm{ik}|r-\mathrm{\ϵ}|})\×\frac{1}{|r-\mathrm{\ϵ}|}\left(\begin{array}{c}{r}_x-{\mathrm{\ϵ}}_x\\ r_y-{\mathrm{\ϵ}}_x\\ r_z-{\mathrm{\ϵ}}_z\end{array}\right)s\left(\mathrm{\ω}\right)$

Wherein k represents wave number, relevant with the velocity of sound with the frequency of voice signal, and e is constant, and i is imaginary unit, and ω is the angular frequency of signal s (t), | r-L _i| represent loud speaker and listen the distance between the point of articulation, | r-ε | represent sound source and listen the distance between the point of articulation, G (ω) represents source strength;

Step 5.2: according to the acoustic pressure p of the original sound source that step 1 calculates ₀acoustic pressure p with acoustic image in the reconstruction sound field calculated in step 4, sets up distance equivalence relation as follows:

p＝p ₀

Namely

$G\underset{j=1}{\overset{n}{\mathrm{\Σ}}}\frac{e^{-\mathrm{ik}|r-L_i|}}{|r-L_i|}q\left(\mathrm{\ω}\right)=G\frac{e^{-\mathrm{ik}|r-\mathrm{\ϵ}|}}{|r-\mathrm{\ϵ}|}s\left(\mathrm{\ω}\right)$

Wherein k represents wave number, relevant with the velocity of sound with the frequency of voice signal, and e is constant, and i is imaginary unit, and ω is the angular frequency of signal s (t), | r-L _i| represent loud speaker and listen the distance between the point of articulation, | r-ε | represent sound source and listen the distance between the point of articulation, G (ω) represents source strength;

Step 5.3: according to step 5.1 and 5.2 direction equivalence relation and distance equivalence relation, set up direction, distance equivalence model:

$\left\{\begin{array}{c}\mathrm{pv}={\mathrm{pv}}_0\\ p=p_0\end{array}\right.$

The technical scheme that device of the present invention adopts is: a kind of method that in multi-channel system utilized described in claim 1, Sounnd source direction and distance are rebuild carries out the device that in multi-channel system, Sounnd source direction and distance are rebuild, it is characterized in that, comprise: point of articulation directional information computing module (1) is listened in sound source space, point of articulation range information computing module (2) is listened in sound source space, module (3) selected by loud speaker, signal forward allocator module (4), rebuild sound field Sounnd source direction information computational module (5), rebuild sound field sound source range information computing module (6), model building module (7), model solution module (8), signal distribution module (9),

The direction of the sound source that described sound source space listens point of articulation directional information computing module (1) to perceive for the pleasant to the ear point of articulation in sound source space calculating acoustic source place, signal s (t) according to sound source calculates the particle rapidity pv listening the point of articulation ₀, and by pv ₀export model building module (7) to;

The distance of the sound source that described sound source space listens point of articulation range information computing module (2) to perceive for the pleasant to the ear point of articulation in sound source space calculating acoustic source place, signal s (t) according to sound source calculates the acoustic pressure p listening the point of articulation ₀, and by p ₀export model building module (7) to;

Described loud speaker select module (3) for set up reconstruction multi-channel system in the selection principle of loud speaker, determine the loud speaker L selected in multi-channel system ₁, L ₂, L ₃, L ₄, and export selected loud speaker to signal forward allocator module (4);

Described signal forward allocator module (4) is for being multiplied by weight w by original audio signal s (t) ₁, w ₂, w ₃, w ₄after distribute to the loud speaker L that loud speaker selects to select in module (3) ₁, L ₂, L ₃, L ₄, then the signal after distribution is exported to and rebuilds sound field Sounnd source direction information computational module (5) and rebuild sound field sound source range information computing module (6);

Described reconstruction sound field acoustic image directional information computing module (5) receives the direction of acoustic image for calculating the pleasant to the ear point of articulation place of multi-channel system, according in signal forward allocator module (4) to the preallocated signal of loud speaker, the pleasant to the ear point of articulation place of multi-channel system utilizing acoustic theory to calculate reconstruction receives the particle rapidity pv of acoustic image, and pv is exported to model building module (7);

Described reconstruction sound field audio-visual distance information computing module (6) receives the distance of acoustic image for calculating the pleasant to the ear point of articulation place of multi-channel system, according in signal forward allocator module (4) to the preallocated signal of loud speaker, the pleasant to the ear point of articulation place of multi-channel system utilizing acoustic theory to calculate reconstruction receives the acoustic pressure p of acoustic image, and p is exported to model building module (7);

Described model building module (7) for set up listen point of articulation prescription to perceived distance consistency model, the sound source particle rapidity p listening point of articulation directional information computing module (1) to export by sound source space _v0direction equivalence relation is set up, the sound source acoustic pressure p listening point of articulation range information computing module (2) to export by sound source space with the acoustic image particle rapidity pv rebuilding output in sound field Sounnd source direction information computational module (4) ₀distance equivalence relation is set up with the acoustic image acoustic pressure p rebuilding output in sound field Sounnd source direction information computational module (6), set up direction, distance equivalence model according to direction equivalence relation and distance equivalence relation, and export this model to model solution module (8);

Described model solution module (8) is set up the model of foundation in module (7) for solving model and then is obtained four weight w ₁, w ₂, w ₃, w ₄value, finally weights are exported to signal distribution module (9);

The weight w that described signal distribution module (9) solves for utilizing model solution module (8) ₁, w ₂, w ₃, w ₄, the signal carrying out loud speaker distributes, and makes the audio direction of reconstruction and distance with the direction of acoustic source with apart from consistent.

The present invention, relative to prior art, can recover direction and the range information of the sound source in acoustic source space accurately, and simple to operate, and computational efficiency is high, good stability.

Accompanying drawing explanation

Fig. 1: the device workflow diagram of the embodiment of the present invention.

Embodiment

By reference to the accompanying drawings technical scheme of the present invention and system are described further with specific embodiment below.

The technical scheme that method of the present invention adopts is: Sounnd source direction and distance reconstructing device and method in a kind of multi-channel system, comprise the following steps:

Step 1: according to time-domain signal s (t) of known sound source, to listen the point of articulation for after initial point sets up Descartes's rectangular coordinate system in original sound field, the particle rapidity pv of the sound source of listening point of articulation place to receive that to calculate with sound source distance be r ₀with acoustic pressure p ₀;

Step 1.1: utilize Fourier transform that time-domain signal s (t) of sound source is transformed into frequency domain by time domain, obtain frequency-region signal s (ω), in this example, time-domain signal s (t) adopts sinusoidal signal, and the sample rate of signal is 48000Hz.Computational process is as follows:

$s\left(\mathrm{\ω}\right)=\underset{-\∞}{\overset{+\∞}{\∫}}s\left(t\right)e^{-\mathrm{iwt}}\mathrm{dt}$

Step 1.2: this example, to listen the point of articulation for initial point, namely listens the coordinate r (x, y, z) of the point of articulation to be (0,0,0), by frequency-region signal s (ω) and the sound source position vector ε in sound field=(ε _x, ε _y, ε _z), listen point of articulation position coordinates (0,0,0), according to the definition of sound physical properties attribute particle rapidity, calculate listen point of articulation place to receive particle rapidity pv ₀, computational methods are as follows:

${\mathrm{pv}}_0(r,\mathrm{\ω})=G\frac{e^{-\mathrm{ik}\left|\mathrm{\ϵ}\right|}}{\left|\mathrm{\ϵ}\right|}(1+\frac{1}{\mathrm{ik}\left|\mathrm{\ϵ}\right|})\×\frac{1}{\left|\mathrm{\ϵ}\right|}\left(\begin{array}{c}-{\mathrm{\ϵ}}_x\\ -{\mathrm{\ϵ}}_y\\ -{\mathrm{\ϵ}}_z\end{array}\right)s\left(\mathrm{\ω}\right)$

Calculate the acoustic pressure p listening point of articulation place to receive simultaneously ₀:

$p_0(T_0,\mathrm{\ω})=G\frac{e^{-\mathrm{ik}\left|\mathrm{\ϵ}\right|}}{\left|\mathrm{\ϵ}\right|}s\left(\mathrm{\ω}\right)$

Wherein k represents wave number, relevant with the velocity of sound with the frequency of voice signal, and e is constant, and i is imaginary unit, and ω is the angular frequency of signal s (t), | ε | represent sound source and listen the distance between the point of articulation, G (ω) represents source strength.

Step 2: according to the particle rapidity pv of the sound source calculated in step 1 ₀, select to comprise particle rapidity pv ₀4 loud speaker L in direction ₁, L ₂, L ₃, L ₄, wherein L ₁, L ₂, L ₃, L ₄meet following condition: point of articulation r (x, y, z) and L will be listened respectively ₁, L ₂, L ₃, L ₄connect, in the polyhedron formed, the particle rapidity pv of sound source ₀unit vector to be positioned at this polyhedron inner.

Step 3: by the loud speaker L in sound-source signal s (t) and step 2 ₁, L ₂, L ₃, L ₄, s (t) is multiplied by respectively 4 weight w ₁, w ₂, w ₃, w ₄after be assigned on four selected loud speakers, obtain four preallocated signal q of loud speaker ₁(t), q ₂(t), q ₃(t), q ₄(t);

$q\left(t\right)=\left(\begin{array}{c}{q}_1\left(t\right)\\ q_2\left(t\right)\\ q_3\left(t\right)\\ q_4\left(t\right)\end{array}\right)=W\·s\left(t\right)=\left(\begin{array}{c}{w}_1\\ w_2\\ w_3\\ w_4\end{array}\right)s\left(t\right)$

Wherein W represents weight vector, and its value is:

$W=\left(\begin{array}{c}{w}_1\\ w_2\\ w_3\\ w_4\end{array}\right)$

Q (t)=(q ₁(t), q ₂(t), q ₃(t), q ₄(t)) T represents loud speaker L ₁, L ₂, L ₃, L ₄the signal phasor distributed, T represents transpose operation.

Step 4: calculate and rebuild in sound field by loud speaker L ₁, L ₂, L ₃, L ₄listen particle rapidity pv and the acoustic pressure p of the acoustic image of point of articulation place perception after sending signal, detailed process comprises following 3 sub-steps:

Step 4.1: at playback end, with the center of the ambiophonic system of multichannel composition for initial point, sets up cartesian coordinate system.The coordinate of loud speaker Li is designated as L _i(L _ix, L _iy, L _iz), in this example, listen the position of the point of articulation and the position consistency of central point, namely listen the coordinate of the point of articulation to be r (r _x, r _y, r _z)=(0,0,0);

Step 4.2: by loud speaker L ₁, L ₂, L ₃, L ₄the signal q sent ₁(t), q ₂(t), q ₃(t), q ₄t () is converted to frequency domain:

$q\left(\mathrm{\ω}\right)=\left(\begin{array}{c}{q}_1\left(\mathrm{\ω}\right)\\ q_2\left(\mathrm{\ω}\right)\\ q_3\left(\mathrm{\ω}\right)\\ q_4\left(\mathrm{\ω}\right)\end{array}\right),q_i\left(\mathrm{\ω}\right)=\underset{-\∞}{\overset{+\∞}{\∫}}{q}_i\left(t\right)e^{-\mathrm{iwt}}\mathrm{dt}$

Step 4.3: by loud speaker L ₁, L ₂, L ₃, L ₄send signal q ₁(t), q ₂(t), q ₃(t), q ₄t () calculates the frequency domain representation pv and the acoustic pressure p that listen the particle rapidity at point of articulation place:

$\begin{array}{c}{p}_v(r,w)=\underset{i=1}{\overset{3}{\mathrm{\Σ}}}G\frac{e^{-\mathrm{ik}\left|L_i\right|}}{\left|L_i\right|}(1+\frac{1}{\mathrm{ik}\left|L_i\right|})\×\frac{1}{\left|L_i\right|}\left(\begin{array}{c}-L_{\mathrm{ix}}\\ -L_{\mathrm{iy}}\\ -L_{\mathrm{iz}}\end{array}\right)q\left(\mathrm{\ω}\right)\\ =\underset{i=1}{\overset{3}{\mathrm{\Σ}}}G\frac{e^{-\mathrm{ik}\left|L_i\right|}}{{\left|L_i\right|}^2}(1+\frac{1}{\mathrm{ik}\left|L_i\right|})\×\left(\begin{array}{c}-L_{\mathrm{ix}}\\ -L_{\mathrm{iy}}\\ -L_{\mathrm{iz}}\end{array}\right)q\left(\mathrm{\ω}\right)\end{array}$

$p=G\underset{j=1}{\overset{n}{\mathrm{\Σ}}}\frac{e^{-\mathrm{ik}\left|L_i\right|}}{\left|L_i\right|}q\left(\mathrm{\ω}\right)$

Wherein k represents wave number, relevant with the velocity of sound with the frequency of voice signal, and e is constant, and i is imaginary unit, and ω is the angular frequency of signal s (t), | L _i| represent loud speaker and listen the distance between the point of articulation, G (ω) represents source strength;

Step 5: according to the particle rapidity pv of the original sound source that step 1 calculates ₀with acoustic pressure p ₀, the particle rapidity pv of acoustic image and acoustic pressure p in the reconstruction sound field that calculates in step 4, set up direction, distance equivalence model, comprise following 3 sub-steps:

Step 5.1: according to the particle rapidity pv of the original sound source that step 1 calculates ₀particle rapidity pv with acoustic image in the reconstruction sound field calculated in step 4, sets up direction equivalence relation as follows:

pv＝pv ₀

Namely

$\underset{i=1}{\overset{3}{\mathrm{\Σ}}}G\frac{e^{-\mathrm{ik}\left|L_i\right|}}{{\left|L_i\right|}^2}(1+\frac{1}{\mathrm{ik}\left|L_i\right|})\×\left(\begin{array}{c}-L_{\mathrm{ix}}\\ -L_{\mathrm{iy}}\\ -L_{\mathrm{iz}}\end{array}\right)q\left(\mathrm{\ω}\right)=G\frac{e^{-\mathrm{ik}\left|\mathrm{\ϵ}\right|\mathrm{\ϵ}}}{\left|\mathrm{\ϵ}\right|}(1+\frac{1}{\mathrm{ik}\left|\mathrm{\ϵ}\right|})\×\frac{1}{\left|\mathrm{\ϵ}\right|}\left(\begin{array}{c}{-\mathrm{\ϵ}}_x\\ -{\mathrm{\ϵ}}_x\\ -{\mathrm{\ϵ}}_z\end{array}\right)s\left(\mathrm{\ω}\right)$

Wherein k represents wave number, relevant with the velocity of sound with the frequency of voice signal, and e is constant, and i is imaginary unit, and ω is the angular frequency of signal s (t), | L _i| represent loud speaker and listen the distance between the point of articulation, | ε | represent sound source and listen the distance between the point of articulation, G (ω) represents source strength.

Step 5.2: according to the acoustic pressure p of the original sound source that step 1 calculates ₀acoustic pressure p with acoustic image in the reconstruction sound field calculated in step 4, sets up distance equivalence relation as follows:

p＝p ₀

Namely

$G\underset{j=1}{\overset{n}{\mathrm{\Σ}}}\frac{e^{-\mathrm{ik}\left|L_i\right|}}{\left|L_i\right|}q\left(\mathrm{\ω}\right)=G\frac{e^{-\mathrm{ik}\left|\mathrm{\ϵ}\right|}}{\left|\mathrm{\ϵ}\right|}s\left(\mathrm{\ω}\right)$

Wherein k represents wave number, relevant with the velocity of sound with the frequency of voice signal, and e is constant, and i is imaginary unit, and ω is the angular frequency of signal s (t), | L _i| represent loud speaker and listen the distance between the point of articulation, | ε | represent sound source and listen the distance between the point of articulation, G (ω) represents source strength.

Step 5.3: according to step 5.1 and 5.2 direction equivalence relation and distance equivalence relation, set up direction, distance equivalence model:

$\left\{\begin{array}{c}\mathrm{pv}={\mathrm{pv}}_0\\ p=p_0\end{array}\right.$

Step 6: in the direction of step 5, distance equivalence model, with the weight w arranged in step 3 ₁, w ₂, w ₃, w ₄for unknown quantity, solve the value of direction, distance equivalence model acquisition weight; Model is launched to obtain system of linear equations, separates this system of linear equations, can utilize but be not limited to the tool software such as Matlab, Mathematical, or other instrument, and manually calculate, obtain this solution of equations (w ₁, w ₂, w ₃, w ₄).

Step 7: utilize the weight w solved in step 6 ₁, w ₂, w ₃, w ₄, the signal carrying out loud speaker distributes, and makes the audio direction of reconstruction and distance with the direction of acoustic source with apart from consistent, in this example, all loud speakers are positioned on same sphere, namely each loud speaker is equal with listening the distance between the point of articulation, therefore in signal assigning process without the need to postponing, loud speaker L ₁, L ₂, L ₃, L ₄distribute signal as shown in the formula:

$\left\{\begin{array}{c}{q}_1\left(t\right)=w_{1}s\left(t\right)\\ q_2\left(t\right)=w_{2}s\left(t\right)\\ q_3\left(t\right)=w_{3}s\left(t\right)\\ q_4\left(t\right)=w_{4}s\left(t\right)\end{array}\right.$

Ask for an interview Fig. 1, the technical scheme that device of the present invention adopts is: the device recovering Sounnd source direction and range information in a kind of multi-channel system, it is characterized in that, comprising: point of articulation directional information computing module 1 is listened in sound source space, sound source space listens point of articulation range information computing module 2, loud speaker to select module 3, signal forward allocator module 4, rebuild sound field Sounnd source direction information computational module 5, rebuild sound field sound source range information computing module 6, model building module 7, model solution module 8, signal distribution module 9;

The direction of the sound source that described sound source space listens point of articulation directional information computing module 1 to perceive for the pleasant to the ear point of articulation in sound source space calculating acoustic source place, signal s (t) according to sound source calculates the particle rapidity pv listening the point of articulation ₀, and by pv ₀export model building module 7 to;

The distance of the sound source that described sound source space listens point of articulation range information computing module 2 to perceive for the pleasant to the ear point of articulation in sound source space calculating acoustic source place, signal s (t) according to sound source calculates the acoustic pressure p listening the point of articulation ₀, and by p ₀export model building module 7 to;

Described loud speaker select module 3 for set up reconstruction multi-channel system in the selection principle of loud speaker, determine the loud speaker L selected in multi-channel system ₁, L ₂, L ₃, L ₄, and export selected loud speaker to signal forward allocator module 4;

Described signal forward allocator module 4 is for being multiplied by weight w by original audio signal s (t) ₁, w ₂, w ₃, w ₄after distribute to the loud speaker L that loud speaker selects to select in module 3 ₁, L ₂, L ₃, L ₄, then the signal after distribution is exported to and rebuilds sound field Sounnd source direction information computational module 5 and rebuild sound field sound source range information computing module 6;

Described reconstruction sound field acoustic image directional information computing module 5 receives the direction of acoustic image for calculating the pleasant to the ear point of articulation place of multi-channel system, according in signal forward allocator module 4 to the preallocated signal of loud speaker, the pleasant to the ear point of articulation place of multi-channel system utilizing acoustic theory to calculate reconstruction receives the particle rapidity pv of acoustic image, and exports pv to model building module 7;

Described reconstruction sound field audio-visual distance information computing module 6 receives the distance of acoustic image for calculating the pleasant to the ear point of articulation place of multi-channel system, according in signal forward allocator module 4 to the preallocated signal of loud speaker, the pleasant to the ear point of articulation place of multi-channel system utilizing acoustic theory to calculate reconstruction receives the acoustic pressure p of acoustic image, and exports p to model building module 7;

Described model building module 7 for set up listen point of articulation prescription to perceived distance consistency model, the sound source particle rapidity pv listening point of articulation directional information computing module 1 to export by sound source space ₀direction equivalence relation is set up, the sound source acoustic pressure p listening point of articulation range information computing module 2 to export by sound source space with the acoustic image particle rapidity pv rebuilding output in sound field Sounnd source direction information computational module 4 ₀set up distance equivalence relation with the acoustic image acoustic pressure p rebuilding output in sound field Sounnd source direction information computational module 6, set up direction, distance equivalence model according to direction equivalence relation and distance equivalence relation, and export this model to model solution module 8;

Described model solution module 8 is set up the model of foundation in module 7 for solving model and then is obtained four weight w ₁, w ₂, w ₃, w ₄value, finally export weights to signal distribution module 9;

The weight w that described signal distribution module 9 solves for utilizing model solution module 8 ₁, w ₂, w ₃, w ₄, the signal carrying out loud speaker distributes, and makes the audio direction of reconstruction and distance with the direction of acoustic source with apart from consistent.

These are only preferred embodiment of the present invention, be not intended to limit protection scope of the present invention, therefore, all any amendments done within the spirit and principles in the present invention, equivalent replacement, improvement etc., all should be included within protection scope of the present invention.

Claims (5)

1. the method that in multi-channel system, Sounnd source direction and distance are rebuild, is characterized in that, comprise the following steps:

Step 1: according to time-domain signal s (t) of known sound source, to listen the point of articulation for after initial point sets up Descartes's rectangular coordinate system in original sound field, the particle rapidity pv of the sound source of listening point of articulation place to receive that to calculate with sound source distance be r ₀with acoustic pressure p ₀;

Step 2: according to the particle rapidity pv of the sound source calculated in step 1 ₀, select to comprise particle rapidity pv ₀4 loud speaker L in direction ₁, L ₂, L ₃, L ₄;

Step 3: by the loud speaker L in sound-source signal s (t) and step 2 ₁, L ₂, L ₃, L ₄, s (t) is multiplied by respectively 4 weight w ₁, w ₂, w ₃, w ₄after be assigned on four selected loud speakers;

Step 4: at playback end, calculates and rebuilds loud speaker L in sound field ₁, L ₂, L ₃, L ₄particle rapidity pv and the acoustic pressure p of the acoustic image of point of articulation place perception is listened after sending signal;

Step 5: according to the particle rapidity pv of the original sound source that step 1 calculates ₀with acoustic pressure p ₀, the particle rapidity pv of acoustic image and acoustic pressure p in the reconstruction sound field that calculates in step 4, set up direction, distance equivalence model; Its specific implementation comprises following sub-step:

Step 5.1: according to the particle rapidity p of the original sound source that step 1 calculates _v0particle rapidity pv with acoustic image in the reconstruction sound field calculated in step 4, sets up direction equivalence relation as follows:

pv＝pv ₀

Namely

$\underset{i=1}{\overset{3}{\mathrm{\Σ}}}G\frac{e^{-\mathrm{ik}|r-L_i|}}{{|r-L_i|}^2}(1+\frac{1}{\mathrm{ik}|r-L_i|})\×\left(\begin{array}{c}{r}_{\mathrm{ix}}-L_{\mathrm{ix}}\\ r_{\mathrm{iy}}-L_{\mathrm{iy}}\\ r_{\mathrm{iz}}-L_{\mathrm{iz}}\end{array}\right)q\left(\mathrm{\ω}\right)=G\frac{e^{-\mathrm{ik}|r-\mathrm{\ϵ}|}}{|r-\mathrm{\ϵ}|}(1+\frac{1}{\mathrm{ik}|r-\mathrm{\ϵ}|})\×\frac{1}{|r-\mathrm{\ϵ}|}\left(\begin{array}{c}{r}_x-{\mathrm{\ϵ}}_x\\ r_y-{\mathrm{\ϵ}}_x\\ r_z-{\mathrm{\ϵ}}_x\end{array}\right)s\left(\mathrm{\ω}\right)$

Wherein k represents wave number, relevant with the velocity of sound with the frequency of voice signal, and e is constant, and i is imaginary unit, and ω is the angular frequency of signal s (t), | r-L _i| represent loud speaker and listen the distance between the point of articulation, | r-ε | represent sound source and listen the distance between the point of articulation, G (ω) represents source strength;

Step 5.2: according to the acoustic pressure p of the original sound source that step 1 calculates ₀acoustic pressure p with acoustic image in the reconstruction sound field calculated in step 4, sets up distance equivalence relation as follows:

p＝p ₀

Namely

$G\underset{j=1}{\overset{n}{\mathrm{\Σ}}}\frac{e^{-\mathrm{ik}|r-L_i|}}{|r-L_i|}q\left(\mathrm{\ω}\right)=G\frac{e^{-\mathrm{ik}|r-\mathrm{\ϵ}|}}{|r-\mathrm{\ϵ}|}s\left(\mathrm{\ω}\right)$

Wherein k represents wave number, relevant with the velocity of sound with the frequency of voice signal, and e is constant, and i is imaginary unit, and ω is the angular frequency of signal s (t), | r-L _i| represent loud speaker and listen the distance between the point of articulation, | r-ε | represent sound source and listen the distance between the point of articulation, G (ω) represents source strength;

Step 5.3: according to step 5.1 and 5.2 direction equivalence relation and distance equivalence relation, set up direction, distance equivalence model:

$\left\{\begin{array}{c}\mathrm{pv}={\mathrm{pv}}_0\\ p=p_0\end{array}\right.;$

Step 6: in the direction of step 5, apart from equivalence model, with the weight w arranged in step 3 ₁, w ₂, w ₃, w ₄for unknown quantity, solve the value of direction, distance equivalence model acquisition weight;

Step 7: utilize the weight w solved in step 6 ₁, w ₂, w ₃, w ₄, the signal carrying out loud speaker distributes, and makes the audio direction of reconstruction and distance with the direction of acoustic source with apart from consistent.

2. the method that in multi-channel system according to claim 1, Sounnd source direction and distance are rebuild, is characterized in that: the calculating described in step 1 and sound source distance are the particle rapidity pv of the sound source of listening point of articulation place to receive of r ₀with acoustic pressure p ₀, its specific implementation comprises following sub-step:

Step 1.1: utilize Fourier transform that time-domain signal s (t) of sound source is transformed into frequency domain by time domain, obtains frequency-region signal s (ω)

$s\left(\mathrm{\ω}\right)=\underset{-\∞}{\overset{+\∞}{\∫}}s\left(t\right)e^{-\mathrm{iwt}}\mathrm{dt}$

Step 1.2: by frequency-region signal s (ω) and the sound source position vector ε in sound field=(ε _x, ε _y, ε _z), listen point of articulation position coordinates r=r (r _x, r _y, r _z), according to the definition of sound physical properties attribute particle rapidity, calculate the particle rapidity pv listening point of articulation place to receive ₀:

${\mathrm{pv}}_0(r,\mathrm{\ω})=G\frac{e^{-\mathrm{ik}|r-\mathrm{\ϵ}|}}{|r-\mathrm{\ϵ}|}(1+\frac{1}{\mathrm{ik}|r-\mathrm{\ϵ}|})\×\frac{1}{|r-\mathrm{\ϵ}|}\left(\begin{array}{c}{r}_x-{\mathrm{\ϵ}}_x\\ r_y-{\mathrm{\ϵ}}_y\\ r_z-{\mathrm{\ϵ}}_z\end{array}\right)s\left(\mathrm{\ω}\right);$

Calculate the acoustic pressure p listening point of articulation place to receive ₀:

$p_0(T_0,\mathrm{\ω})=G\frac{e^{-\mathrm{ik}|r-\mathrm{\ϵ}|}}{|r-\mathrm{\ϵ}|}s\left(\mathrm{\ω}\right);$

Wherein k represents wave number, relevant with the velocity of sound with the frequency of voice signal, and e is constant, and i is imaginary unit, and ω is the angular frequency of signal s (t), | r-ε | represent sound source and listen the distance between the point of articulation, G (ω) represents source strength.

3. the method that in multi-channel system according to claim 1, Sounnd source direction and distance are rebuild, is characterized in that: the loud speaker L of 4 described in step 2 ₁, L ₂, L ₃, L ₄, meet following condition: point of articulation r (x, y, z) and L will be listened respectively ₁, L ₂, L ₃, L ₄connect, in the polyhedron formed, the particle rapidity pv of sound source ₀unit vector to be positioned at this polyhedron inner.

4. the method that in multi-channel system according to claim 1, Sounnd source direction and distance are rebuild, is characterized in that: loud speaker L in sound field is rebuild in the calculating described in step 4 ₁, L ₂, L ₃, L ₄listen particle rapidity pv and the acoustic pressure p of the acoustic image of point of articulation place perception after sending signal, its specific implementation comprises following sub-step:

Step 4.1: at playback end, with the center of the ambiophonic system of multichannel composition for initial point, sets up cartesian coordinate system, loud speaker L _icoordinate be designated as L _i(L _ix, L _iy, L _iz), listen the coordinate of the point of articulation to be r (r _x, r _y, r _z);

Step 4.2: by loud speaker L ₁, L ₂, L ₃, L ₄the signal q sent ₁(t), q ₂(t), q ₃(t), q ₄t () is converted to frequency domain:

$q\left(\mathrm{\ω}\right)=\left(\begin{array}{c}{q}_1\left(\mathrm{\ω}\right)\\ q_2\left(\mathrm{\ω}\right)\\ q_3\left(\mathrm{\ω}\right)\\ q_4\left(\mathrm{\ω}\right)\end{array}\right),q_i\left(\mathrm{\ω}\right)=\underset{-\∞}{\overset{+\∞}{\∫}}{q}_i\left(t\right)e^{-\mathrm{iwt}}\mathrm{dt}$

Step 4.3: by loud speaker L ₁, L ₂, L ₃, L ₄send signal q ₁(t), q ₂(t), q ₃(t), q ₄t () calculates the frequency domain representation pv and the acoustic pressure p that listen the particle rapidity at point of articulation place:

$\begin{array}{c}{p}_v(r,\mathrm{\ω})=\underset{i=1}{\overset{3}{\mathrm{\Σ}}}G\frac{e^{-\mathrm{ik}|r-L_i|}}{|r-L_i|}(1+\frac{1}{\mathrm{ik}|r-L_i|})\×\frac{1}{|r-L_i|}\left(\begin{array}{c}{r}_x-L_{\mathrm{ix}}\\ r_y-L_{\mathrm{iy}}\\ r_z-L_{\mathrm{iz}}\end{array}\right)q\left(\mathrm{\ω}\right)\\ =\underset{i=1}{\overset{3}{\mathrm{\Σ}}}G\frac{e^{-\mathrm{ik}|r-L_i|}}{{|r-L_i|}^2}(1+\frac{1}{\mathrm{ik}|r-L_i|})\×\left(\begin{array}{c}{r}_x-L_{\mathrm{ix}}\\ r_y-L_{\mathrm{iy}}\\ r_z-L_{\mathrm{iz}}\end{array}\right)q\left(\mathrm{\ω}\right)\end{array}$

$p=G\underset{j=1}{\overset{n}{\mathrm{\Σ}}}\frac{e^{-\mathrm{ik}|r-L_i|}}{|r-L_i|}q\left(\mathrm{\ω}\right)$

Wherein k represents wave number, relevant with the velocity of sound with the frequency of voice signal, and e is constant, and i is imaginary unit, and ω is the angular frequency of signal s (t), | r-L _i| represent loud speaker and listen the distance between the point of articulation, G (ω) represents source strength.

5. one kind utilizes the method that in the multi-channel system described in claim 1, Sounnd source direction and distance are rebuild to carry out the device that in multi-channel system, Sounnd source direction and distance are rebuild, it is characterized in that, comprise: point of articulation directional information computing module (1) is listened in sound source space, point of articulation range information computing module (2) is listened in sound source space, module (3) selected by loud speaker, signal forward allocator module (4), rebuild sound field Sounnd source direction information computational module (5), rebuild sound field sound source range information computing module (6), model building module (7), model solution module (8), signal distribution module (9),

The direction of the sound source that described sound source space listens point of articulation directional information computing module (1) to perceive for the pleasant to the ear point of articulation in sound source space calculating acoustic source place, signal s (t) according to sound source calculates the particle rapidity pv listening the point of articulation ₀, and by pv ₀export model building module (7) to;

The distance of the sound source that described sound source space listens point of articulation range information computing module (2) to perceive for the pleasant to the ear point of articulation in sound source space calculating acoustic source place, signal s (t) according to sound source calculates the acoustic pressure p listening the point of articulation ₀, and by p ₀export model building module (7) to;

Described loud speaker select module (3) for set up reconstruction multi-channel system in the selection principle of loud speaker, determine the loud speaker L selected in multi-channel system ₁, L ₂, L ₃, L ₄, and export selected loud speaker to signal forward allocator module (4);

Described signal forward allocator module (4) is for being multiplied by weight w by original audio signal s (t) ₁, w ₂, w ₃, w ₄after distribute to the loud speaker L that loud speaker selects to select in module (3) ₁, L ₂, L ₃, L ₄, then the signal after distribution is exported to and rebuilds sound field Sounnd source direction information computational module (5) and rebuild sound field sound source range information computing module (6);

Described reconstruction sound field acoustic image directional information computing module (5) receives the direction of acoustic image for calculating the pleasant to the ear point of articulation place of multi-channel system, according in signal forward allocator module (4) to the preallocated signal of loud speaker, the pleasant to the ear point of articulation place of multi-channel system utilizing acoustic theory to calculate reconstruction receives the particle rapidity pv of acoustic image, and pv is exported to model building module (7);

Described reconstruction sound field audio-visual distance information computing module (6) receives the distance of acoustic image for calculating the pleasant to the ear point of articulation place of multi-channel system, according in signal forward allocator module (4) to the preallocated signal of loud speaker, the pleasant to the ear point of articulation place of multi-channel system utilizing acoustic theory to calculate reconstruction receives the acoustic pressure p of acoustic image, and p is exported to model building module (7);

Described model building module (7) for set up listen point of articulation prescription to perceived distance consistency model, the sound source particle rapidity p listening point of articulation directional information computing module (1) to export by sound source space _v0direction equivalence relation is set up, the sound source acoustic pressure p listening point of articulation range information computing module (2) to export by sound source space with the acoustic image particle rapidity pv rebuilding output in sound field Sounnd source direction information computational module (4) ₀distance equivalence relation is set up with the acoustic image acoustic pressure p rebuilding output in sound field Sounnd source direction information computational module (6), set up direction, distance equivalence model according to direction equivalence relation and distance equivalence relation, and export this model to model solution module (8);

Set up direction, distance equivalence model according to direction equivalence relation and distance equivalence relation, its specific implementation comprises following sub-step: step 5.1: according to the particle rapidity p of the original sound source that step 1 calculates _v0particle rapidity pv with acoustic image in the reconstruction sound field calculated in step 4, sets up direction equivalence relation as follows:

pv＝pv ₀

Namely

$\underset{i=1}{\overset{3}{\mathrm{\Σ}}}G\frac{e^{-\mathrm{ik}|r-L_i|}}{{|r-L_i|}^2}(1+\frac{1}{\mathrm{ik}|r-L_i|})\×\left(\begin{array}{c}{r}_{\mathrm{ix}}-L_{\mathrm{ix}}\\ r_{\mathrm{iy}}-L_{\mathrm{iy}}\\ r_{\mathrm{iz}}-L_{\mathrm{iz}}\end{array}\right)q\left(\mathrm{\ω}\right)=G\frac{e^{-\mathrm{ik}|r-\mathrm{\ϵ}|}}{|r-\mathrm{\ϵ}|}(1+\frac{1}{\mathrm{ik}|r-\mathrm{\ϵ}|})\×\frac{1}{|r-\mathrm{\ϵ}|}\left(\begin{array}{c}{r}_x-{\mathrm{\ϵ}}_x\\ r_y-{\mathrm{\ϵ}}_x\\ r_z-{\mathrm{\ϵ}}_x\end{array}\right)s\left(\mathrm{\ω}\right)$

Wherein k represents wave number, relevant with the velocity of sound with the frequency of voice signal, and e is constant, and i is imaginary unit, and ω is the angular frequency of signal s (t), | r-L _i| represent loud speaker and listen the distance between the point of articulation, | r-ε | represent sound source and listen the distance between the point of articulation, G (ω) represents source strength;

Step 5.2: according to the acoustic pressure p of the original sound source that step 1 calculates ₀acoustic pressure p with acoustic image in the reconstruction sound field calculated in step 4, sets up distance equivalence relation as follows:

p＝p ₀

Namely

$G\underset{j=1}{\overset{n}{\mathrm{\Σ}}}\frac{e^{-\mathrm{ik}|r-L_i|}}{|r-L_i|}q\left(\mathrm{\ω}\right)=G\frac{e^{-\mathrm{ik}|r-\mathrm{\ϵ}|}}{|r-\mathrm{\ϵ}|}s\left(\mathrm{\ω}\right)$

Wherein k represents wave number, relevant with the velocity of sound with the frequency of voice signal, and e is constant, and i is imaginary unit, and ω is the angular frequency of signal s (t), | r-L _i| represent loud speaker and listen the distance between the point of articulation, | r-ε | represent sound source and listen the distance between the point of articulation, G (ω) represents source strength;

Step 5.3: according to step 5.1 and 5.2 direction equivalence relation and distance equivalence relation, set up direction, distance equivalence model:

$\left\{\begin{array}{c}\mathrm{pv}=p{v}_0\\ p=p_0\end{array}\right.;$

Described model solution module (8) is set up the model of foundation in module (7) for solving model and then is obtained four weight w ₁, w ₂, w ₃, w ₄value, finally weights are exported to signal distribution module (9);

The weight w that described signal distribution module (9) solves for utilizing model solution module (8) ₁, w ₂, w ₃, w ₄, the signal carrying out loud speaker distributes, and makes the audio direction of reconstruction and distance with the direction of acoustic source with apart from consistent.