CN106034274A - 3D sound device based on sound field wave synthesis and synthetic method - Google Patents
3D sound device based on sound field wave synthesis and synthetic method Download PDFInfo
- Publication number
- CN106034274A CN106034274A CN201510112555.XA CN201510112555A CN106034274A CN 106034274 A CN106034274 A CN 106034274A CN 201510112555 A CN201510112555 A CN 201510112555A CN 106034274 A CN106034274 A CN 106034274A
- Authority
- CN
- China
- Prior art keywords
- audio
- sound
- frequency
- signal
- digital
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Landscapes
- Stereophonic System (AREA)
Abstract
The invention discloses a 3D sound device based on sound field wave synthesis and a synthetic method. The 3D sound device comprises an analog-to-digital converter for performing analog-to-digital conversion on an inputted an audio signal and outputting a digital signal, a microcontroller which is connected to the analog-to-digital converter and performing pre-processing on the digital signal to form multiple paths of audio digital signals, a 3D digital signal processor which is connected to the microcontroller and is used for performing frequency division coding processing on each path of audio digital signals, performing sound effect equalization processing on digital signals after frequency division processing , forming a virtual space effect through calculation and analysis, and synthesizing and outputting multiple paths of radio frequency signals capable of forming surround sound after diffraction processing. Through above arrangement, the 3D sound device based on sound field wave synthesis and the synthetic method enable people to listen to high quality 3D sound in a wide space with a single sound device as a center, enable listening effects on various nodes on and under a plane surrounding the sound device to be uniform, do not need to put a listener to an appointed position and realize the surround sound effect.
Description
Technical field
The present invention relates to field of mobile equipment, especially relate to a kind of 3D based on the synthesis of sound field ripple
PA-system and synthetic method thereof.
Background technology
After the mankind enter mobile Internet society, media termination increasingly miniaturization, it is allowed to possess shifting
Moving and immanent characteristic, its sound effects content quality supported is more and more higher.In the case,
Traditional stereo sound field in double track monolayer face, can not meet the needs of people.People need
Small portable and can embody the acoustic product of high-quality audio realize more preferable third dimension by and
Space perception.
In daily life, we listen thing with two ears, obtain information from source of sound everywhere,
Location sound is carried out again by the calculating of human brain.The 3D audio of computer simulation human brain calculates, and passes through
Digital tone source plays back, and let us is felt oneself to place oneself in the virtual world.Traditional double track
Stereo speaker can realize certain direction feeling, but is difficulty with preferable Space.With many
The 3D audio that individual audio amplifier reaches, then need prohibitively expensive hardware and complicated loudspeaker position to set
Fixed.And, in this setting, people only sets focus (sweet spot) ability at some
Experience optimal 3D audio.Once leave focus too remote, then can be because of the interference between sound
And cause audio quality to fall sharply.
Summary of the invention
The technical problem that present invention mainly solves is to provide a kind of 3D sound equipment based on the synthesis of sound field ripple
Device and synthetic method thereof, it is possible in the space on a large scale centered by single stereo set, listen to height
The effectiveness of quality 3D sound, and uniform at the audibility of plane about and upper and lower each point,
Hearer need not be placed in appointment position, it is achieved that the effect of surround sound.
For solving above-mentioned technical problem, the technical scheme that the present invention uses is: provide a kind of base
In the 3D PA-system of sound field ripple synthesis, including analog-digital converter, for will the audio frequency of input
Signal carries out analog digital conversion, output digit signals;Microcontroller, is connected with analog-digital converter, uses
MCVF multichannel voice frequency digital signal is formed in digital signal being carried out pretreatment;3D digital signal processor,
It is connected with microcontroller, for each road audio digital signals being entered according to frequency and the wavelength of audio frequency
Row frequency dividing coded treatment, and carry out the digital signal after scaling down processing equalizing audio effect processing, pass through
Computational analysis forms Virtual space effect, synthesizes output and can be formed around three-dimensional after diffraction processes
The multipath audio signal of sound.
Wherein, 3D PA-system includes five loudspeaker, exports five tunnel audio signals, one of them
The audio signal of loudspeaker output below 100Hz, remaining two pairs of loudspeaker shares a sound chamber respectively, point
Do not export the audio signal of the left and right acoustic channels of full frequency band.
Wherein, 3D digital signal processor carries out dividing coded treatment, bag to audio digital signals
Include: digital signal samples analysis is divided into multiple frequency range;Value of calculation according to sampling and wave filter
Coefficient modifying algorithm is to adjust filter factor;Filter factor according to adjusting is entered by multiple wave filter
Row filtering calculates;3D digital signal processor is also by by after multiple IIR comb filter parallel connections
Cascade to do the computational analysis of Virtual space effect, wherein, the time delay of each wave filter and
Feedback oscillator is different.
Wherein, 3D digital signal processor carries out diffraction to the sound field forming Virtual space effect
Process, including: obtain coding parameter and initialize;The audio signal of input is carried out windowing
Process;The discrete cosine transform improveing the data block after windowing, obtains the frequency of audio signal
Spectrum parameter, and carry out fast Fourier transform analysis, and become according to frequency spectrum parameter and fast Fourier
Change the output overall situation masking curve of analysis;Calculate substrate curve according to overall situation masking curve, and use base
End curve carries out albefaction to frequency spectrum parameter, removes and the incoherent part of audition, obtains residual signals;
Residual signals is coupled, and quantization encoding is bit stream.
Wherein, 3D digital signal processor also combines modulation with band filter or low pass filter
Device or with block-based conversion to diffraction process after signal synthesize, obtain reconstruct letter
Number can form the audio signal of surround sound.
For solving above-mentioned technical problem, another technical solution used in the present invention is: provide one
The synthetic method of 3D PA-system based on the synthesis of sound field ripple, including: by the audio signal of input
Carry out analog digital conversion, output digit signals;Digital signal is carried out pretreatment and forms MCVF multichannel voice frequency number
Word signal;Each road audio digital signals is carried out dividing at coding by frequency and wavelength according to audio frequency
Reason, and carry out the digital signal after scaling down processing equalizing audio effect processing, formed by computational analysis
Virtual space effect, after diffraction processes, synthesis output can form the MCVF multichannel voice frequency of surround sound
Signal.
Wherein, the step that the digital signal after scaling down processing carries out equalizing audio effect processing includes: will
Digital signal samples analysis is divided into multiple frequency range;Value of calculation and filter coefficient according to sampling are repaiied
Change algorithm to adjust filter factor;Filter factor according to adjusting is filtered by multiple wave filter
Calculate.
Wherein, included by the step of computational analysis formation Virtual space effect: by by multiple
After IIR comb filter parallel connection, cascade forms Virtual space effect, prolonging of the most each wave filter
Time and feedback oscillator are different late.
Wherein, the sound field forming Virtual space effect is carried out diffraction process, including: obtain and compile
Code parameter also initializes;The audio signal of input is carried out windowing process;To the number after windowing
Carry out the discrete cosine transform improved according to block, obtain the frequency spectrum parameter of audio signal, and carry out fast
Speed Fourier transform analysis, and cover according to the output overall situation of frequency spectrum parameter and fast Fourier transform analysis
Cover curve;Calculate substrate curve according to overall situation masking curve, and with substrate curve, frequency spectrum parameter is entered
Row albefaction, removes and the incoherent part of audition, obtains residual signals;Residual signals is carried out coupling
Close, and quantization encoding is bit stream.
Wherein, the step of the multipath audio signal that synthesis output can form surround sound includes: use
Band filter or low pass filter combine manipulator or process diffraction with block-based conversion
After signal synthesize, the signal obtaining reconstruct can form the audio signal of surround sound.
The invention has the beneficial effects as follows: can be in the space on a large scale centered by single stereo set
Listen to the effectiveness of high-quality 3D sound, and at plane about and the audibility of upper and lower each point
Uniformly, hearer need not be placed in appointment position, it is achieved that the effect of surround sound.
Accompanying drawing explanation
Fig. 1 is the structural representation of the 3D PA-system based on the synthesis of sound field ripple of the embodiment of the present invention
Figure;
Fig. 2 is the schematic diagram that the 3D digital signal processor in Fig. 1 carries out dividing coded treatment;
Fig. 3 is the schematic diagram that the 3D digital signal processor in Fig. 1 carries out equalizing audio effect processing;
Fig. 4 is the computational analysis that the 3D digital signal processor in Fig. 1 does Virtual space effect
Schematic diagram;
Fig. 5 is the schematic diagram that the 3D digital signal processor in Fig. 1 carries out diffraction process to sound field;
Fig. 6 is that signal is synthesized by the block-based conversion of 3D digital signal processor in Fig. 1
The schematic diagram of output;
Fig. 7 is the stereochemical structure of the 3D PA-system based on the synthesis of sound field ripple of the embodiment of the present invention
Schematic diagram;
Fig. 8 is the side view of the 3D PA-system based on the synthesis of sound field ripple of the embodiment of the present invention;
Fig. 9 is the section signal of the 3D PA-system based on the synthesis of sound field ripple of the embodiment of the present invention
Figure;
Figure 10 is the synthetic method of the 3D PA-system based on the synthesis of sound field ripple of the embodiment of the present invention
Schematic flow sheet.
Detailed description of the invention
Refer to the 3D sound equipment dress based on the synthesis of sound field ripple that Fig. 1, Fig. 1 are the embodiment of the present invention
The structural representation put.As it is shown in figure 1,3D PA-system 10 based on the synthesis of sound field ripple includes
Analog-digital converter 11, microcontroller 12,3D digital signal processor 13, power amplifier 15,
Multiple wave filter 16, antenna 17, mike 18 and modem 19.From antenna 17 or
The sound that mike 18 receives is modulated after demodulator 19 processes and transmits to analog-digital converter 11.Mould
Number converter 11 is for carrying out analog digital conversion, output digit signals by the audio signal of input;Micro-
Controller 12 is connected with analog-digital converter 11, forms multichannel for digital signal carries out pretreatment
Audio digital signals;3D digital signal processor 13 is connected with microcontroller 12, for basis
The frequency of audio frequency and wavelength carry out dividing coded treatment to each road audio digital signals, and will frequency dividing
Digital signal after process carries out equalizing audio effect processing, forms Virtual space by computational analysis and imitates
Really, after diffraction processes, synthesis output can form the multipath audio signal of surround sound.Above-mentioned many
Road audio signal carries out power amplification through power amplifier 15, and wave filter 16 is filtered
Exported by loudspeaker 141,142,143,144,145 afterwards.
In embodiments of the present invention, 3D digital signal processor is first according to frequency and the ripple of audio frequency
Long carrying out divides coded treatment.The frequency of audio signal and wavelength table are as shown in table 1, audio signal
Frequency between 20Hz-20kHz, wherein below 100Hz is bass or subwoofer.At coding
Reason process, as a example by frame length 128, can also use other the most in other embodiments of the invention
Frame length encode.As in figure 2 it is shown, before sample code starts, first have to receive
128 time domain samples of the first frame, this is that so-called framing postpones, and its value is the length of frame, so
Afterwards 128 time domain samples of the first frame received are made linear prediction, obtain 128 predictive values,
Before being placed in original sampling point, form the data block of a length of 256, then this data block is improved
Discrete cosine transform (modifed discrete cosine transform, MDCT).Due to framing
Postpone the overlap needing adjacent data blocks to have 50%, from the character of MDCT, in order to recover
128 sampling points of the first frame, need to combine its adjacent data block, and therefore, encoder also needs to connect
Receive 128 sampling points of the second frame.
Table 1
In embodiments of the present invention, after frequency dividing coded treatment, as it is shown on figure 3,3D digital signal
Audio digital signals is carried out equalizing audio effect processing by processor, including: by digital signal samples analysis
It is divided into multiple frequency range;Value of calculation and filter coefficient according to sampling revise algorithm to adjust filtering
Coefficient;It is filtered calculating by multiple wave filter according to the filter factor adjusted.Specifically, as
Shown in Fig. 3, it is multiple frequency by the frequency range of audio frequency 20Hz-20kHz according to certain regular partition
Section.I.e. the signal of each frequency range is carried out Gain tuning so that it is reach certain audition balance.Total is
System function is:
H (z)=G1H1(z)+G2H2(z)+…+GNHN(z)
Wherein HNZ () is the filter system function of n-th frequency range, GNZ () is n-th frequency range
Signal gain conciliation amount, x (n) is input signal, and y (n) is output signal, and H (z) is total system
Gain.
The multiple wave filter carrying out equilibrium treatment constitute a wave filter covering whole tonal range
Group, when gain is 0dB, total system frequency response is the most smooth.With mega bass sound
As a example by effect, by low frequency frequency range divide in order to two sections: 120Hz-100Hz or below 100Hz low
Frequency range.Carry out equilibrium treatment with a high pass filter and a low pass filter, make up to the greatest extent
Possible is smooth, i.e. by low-band gain strengthens the effect controlling to reach mega bass.
Audio signal is after equilibrium treatment, and 3D digital signal processor is combed also by multiple IIR
Cascade to do the computational analysis of Virtual space effect specifically after shape wave filter parallel connection, such as Fig. 4 institute
Showing, the first order is made up of 4 IIR comb filter in parallel, and the time delay of wave filter is respectively
For ZD1、ZD2、ZD3、ZD4, and feedback oscillator is respectively a1, a2, a3, a4.Each filtering
The time delay of device and feedback oscillator are all different.And second, third level to be structure identical
All-pass filter.The time delay of wave filter is respectively ZD5、ZD6, feedback oscillator be respectively a5,
a6.Time delay and the feedback oscillator of these two wave filter are the most different.Total transfer function is:
In embodiments of the present invention, prolongation time and the feedback oscillator of each wave filter can be adjusted
Whole, suitably choose suitable parameter.Specifically, several groups of optional parameters can be formed, cause not
With effect, then lead to sound field and be synthetically formed Virtual space effect.
After forming Virtual space effect, 3D digital signal processor is to forming Virtual space effect
The sound field of fruit carries out diffraction process, as it is shown in figure 5, include:
1, obtain coding parameter and initialize.The main feature ginseng by input audio signal
Number and the parameter acquiring of user's input.Encoder need to initialize accordingly according to coding parameter,
The most also being sent to decoding end with the form of head bag by these coding parameters, decoder is according in head bag
Coding parameter carry out necessity initialization, set up decoding needed for all information.
2, the audio signal of input is carried out windowing process.In order to carry out time frequency analysis, it is right to need
The audio signal windowing process of input, is conveniently carried out the data block of subsequent treatment.Wherein added
Window be referred to as converting window.The selection of window type considers time domain and the resolution of frequency domain, meets simultaneously
MDCT, revise inverse discrete cosine transform (inverse modified discrete cosine transform,
IMDCT) reconstruction condition.In embodiments of the present invention, it is preferred to use sinusoidal windows.
3, the discrete cosine transform improveing the data block after windowing, obtains audio signal
Frequency spectrum parameter, and carry out fast Fourier transform (Fast Fourier Transformation, FFT)
Analyze, and the output overall situation masking curve analyzed according to frequency spectrum parameter and fast Fourier transform.
Specifically, the data block after windowing is analyzed through MDCT, obtains the frequency domain of audio signal
Parameter, i.e. frequency spectrum.In MDCT analyzes, can partly remove signal by time-frequency conversion
Statistical redundancy, the energy showing as audio signal in a frequency domain is concentrated mainly on the spectral line of minority,
Be conducive to quantifying and coding.Frequently the frequency domain parameter of signal is also referred to as the spectral coefficient of MDCT, also may be used
For psychoacoustic model analysis, to calculate masking by noise curve.Data block after windowing is made
Fft analysis and using analysis result as the input of psychoacoustic model analysis.Further, according to
Spectral coefficient and the fft analysis result of MDCT carry out psychoacoustic model analysis, obtain noise
Masking curve and tone mask curve.Jointly calculated entirely by masking by noise curve and tone mask curve
Office's masking curve.Overall situation masking curve is the reference quantifying MDCT spectral coefficient, is to allow to draw
The thresholding of the quantizing noise entered.
4, calculate substrate curve according to overall situation masking curve, and with substrate curve, frequency spectrum parameter is entered
Row albefaction, removes and the incoherent part of audition, obtains residual signals.Specifically, low prolonging is applied
High quality audio encryption algorithm (low-delay and high-quality audio coding late
Algorithm, LDX) combine the concept of critical band, the method approached with sectional broken line replaces entirely
Office's masking curve.Though compared to the result of overall situation masking curve energy accurate Characterization psychoacoustic analysis,
If but directly it quantified, encode, then the shortcoming needing a large amount of bit, substrate curve is in low-frequency range
There is less piecewise interval, it is ensured that higher frequency domain resolution, thus reduce low frequency signal
Coding distortion, and relatively big at the piecewise interval of high band, can preferably with the auditory properties phase of people
It coincide.Further, with substrate curve, the spectral coefficient of MDCT is carried out albefaction, remove and listen
Feel incoherent composition.So for stereo mode, enhance the dependency between sound channel, more have
It is beneficial to stereo coding.Residual signals is the MDCT spectral coefficient after albefaction, for follow-up
Quantify and coding.
5, residual signals is coupled, and quantization encoding is bit stream.For multichannel pattern,
Residual signals carries out coupling to improve further compression ratio, and wherein coupled modes include square pole
At least one during coordinate maps and sound channel interweaves.Carry out quantifying to compile by the residual signals after coupling again
Code, obtains coded bit stream, completes diffraction and processes.
Wherein, calculate substrate curve according to overall situation masking curve, and with substrate curve to frequency spectrum parameter
Carry out albefaction and obtain residual signals, and residual signals is coupled, and carry out quantization encoding genus
In rate distortion loop control, to obtain more preferable diffraction treatment effect.
After sound field diffraction has processed, the signal after diffraction is processed by 3D digital signal processor enters
Row synthesis output.Specifically, 3D digital signal processor band filter or low pass filter
In conjunction with manipulator or with block-based conversion, the signal after diffraction process is synthesized, obtain
The signal of reconstruct can form the audio signal of surround sound.Wherein by band filter or low pass filtered
Ripple device combines the method for manipulator and is easy to analyze the condition of signal reconstruction, the most each passage in a frequency domain
The frequency response of analysis/synthetic filtering device be superposed to a constant.Have been devised by based on this principle
The bank of filters of many classics.And the method for block-based conversion is as shown in Figure 6, at diffraction
After reason, the bit stream of output is divided into multiple input block 1,2,3, carries out the caching that splices.Tool
The input block quantity divided is not construed as limiting, and is determined on a case-by-case basis.Wherein N is transform block,
Namely the length of input block.Conventional conversion has discrete Fourier transform (DFT), discrete remaining
String conversion (DCT) and discrete sine transform (DST).After original time-domain signal windowing, direct transform becomes
Frequency domain representation, can explain the characteristic of signal the most in a frequency domain.During synthesis, above-mentioned signal warp
After inverse transformation formed time domain sequences, with synthesis window be multiplied, then with previous part data splicing adding,
Obtain the signal of reconstruct.Each currently processed data block and next number in embodiments of the present invention
According to block overlap 50%, carry out windowing process with suitable window function, make to become to the data block after windowing
Change.So, during synthesis, first carry out inverse transformation, the time domain sequences obtained be multiplied with window function,
Again by the right half part splicing adding of the left-half of the data block of gained Yu last data block, obtain
The signal of reconstruct.In Fig. 6 and represent the block-based conversion of adjacent two input blocks simply,
And need all of input block is carried out above-mentioned block-based conversion in embodiments of the present invention,
The signal finally reconstructed.After reconstruct, the signal of output is only in multipath audio signal herein one
Road, in embodiments of the present invention, 3D digital signal processor needs to enter multipath audio signal simultaneously
The process that row is similar, to obtain being formed the multipath audio signal of surround sound.So can be with
The effectiveness of high-quality 3D sound is listened in space on a large scale centered by single stereo set, and
Plane and the audibility of upper and lower each point about are uniform, and hearer need not be placed in appointment position,
Achieve the effect of surround sound.
Preferably, 3D PA-system includes five loudspeaker 141,142,143,144,145, defeated
Going out five tunnel audio signals, the audio signal of one of them loudspeaker output below 100Hz, remaining is two right
Loudspeaker share a sound chamber respectively, export the audio signal of the left and right acoustic channels of full frequency band respectively.3D
The concrete structure of PA-system is as Figure 7-9.Fig. 7 is axonometric chart, and Fig. 8 is side view, its
Middle figure b is front view, and figure a is left view, and figure c is rearview.Figure d is right view.Fig. 9
For profile.As seen from the figure, there are two loudspeaker 143,144 in left and right in the front of 3D PA-system 10,
Side is respectively arranged with loudspeaker 142,145, after have loudspeaker 141.Specifically, in rearview
Loudspeaker 141 export the audio signal of below 100Hz, the loudspeaker 144 in front view and left view
In loudspeaker 145 share sound chamber 1, export the audio signal of the left and right sound channels of full frequency band respectively,
Loudspeaker 143 in front view and the loudspeaker 142 in left view share sound chamber 2, export full range respectively
The audio signal of the left and right sound channels of section.
Figure 10 is the synthetic method of the 3D PA-system based on the synthesis of sound field ripple of the embodiment of the present invention
Schematic flow sheet.As shown in Figure 10, the synthesis side of 3D PA-system based on the synthesis of sound field ripple
Method includes:
Step S10: the audio signal of input is carried out analog digital conversion, output digit signals.From sky
The sound that line or mike receive is converted into the audio signal of input after being modulated demodulation process, to enter
Row analog digital conversion becomes digital signal.
Step S11: digital signal is carried out pretreatment and forms MCVF multichannel voice frequency digital signal.For being formed
The 3D effect of surround sound, needs to divide audio digital signals according to the isoparametric difference in orientation
Multipath audio signal is become to process to carry out follow-up sound field ripple synthesis.
Step S12: each road audio digital signals is divided by frequency and wavelength according to audio frequency
Coded treatment, and carry out the digital signal after scaling down processing equalizing audio effect processing, divide by calculating
Analysis forms Virtual space effect, and after diffraction processes, can to form surround sound many in synthesis output
Road audio signal.
In step s 12, first carry out dividing coded treatment frequency dividing according to the frequency of audio frequency and wavelength
After coded treatment, carry out audio digital signals equalizing audio effect processing, including: digital signal is adopted
Sample analysis is divided into multiple frequency range;Value of calculation and filter coefficient according to sampling revise algorithm to adjust
Whole filter factor;It is filtered calculating by multiple wave filter according to the filter factor adjusted.Specifically
Ground, is multiple frequency range by the frequency range of audio frequency 20Hz-20kHz according to certain regular partition.I.e.
The signal of each frequency range is carried out Gain tuning so that it is reach certain audition balance.Carry out equilibrium treatment
Multiple wave filter constitute one cover whole tonal range bank of filters, be 0dB in gain
Time, total system frequency response is the most smooth.As a example by mega bass audio, by low again and again
Section divides the low-frequency range for two sections: 120Hz-100Hz or below 100Hz.By a high pass
Wave filter and a low pass filter carry out equilibrium treatment, make up to the most smooth, the most logical
Cross and low-band gain is strengthened the effect controlling to reach mega bass.
Audio signal is after equilibrium treatment, by being cascaded by after multiple IIR comb filter parallel connections
Forming Virtual space effect, time delay and the feedback oscillator of the most each wave filter are different.
The prolongation time of each wave filter and feedback oscillator can be adjusted, and suitably choose suitable parameter.
Specifically, several groups of optional parameters can be formed, cause different effects, then lead to sound field synthesis shape
Become Virtual space effect.
After forming Virtual space effect, the sound field forming Virtual space effect is carried out at diffraction
Reason, including:
1, obtain coding parameter and initialize.The main feature ginseng by input audio signal
Number and the parameter acquiring of user's input.Encoder need to initialize accordingly according to coding parameter,
The most also being sent to decoding end with the form of head bag by these coding parameters, decoder is according in head bag
Coding parameter carry out necessity initialization, set up decoding needed for all information.
2, the audio signal of input is carried out windowing process.In order to carry out time frequency analysis, it is right to need
The audio signal windowing process of input, is conveniently carried out the data block of subsequent treatment.Wherein added
Window be referred to as converting window.The selection of window type considers time domain and the resolution of frequency domain, meets simultaneously
The reconstruction condition of MDCT, IMDCT.In embodiments of the present invention, it is preferred to use sinusoidal windows.
3, the discrete cosine transform improveing the data block after windowing, obtains audio signal
Frequency spectrum parameter, and carry out fft analysis, and according to frequency spectrum parameter and fast Fourier transform analysis
Output overall situation masking curve.
Specifically, the data block after windowing is analyzed through MDCT, obtains the frequency domain of audio signal
Parameter, i.e. frequency spectrum.In MDCT analyzes, can partly remove signal by time-frequency conversion
Statistical redundancy, the energy showing as audio signal in a frequency domain is concentrated mainly on the spectral line of minority,
Be conducive to quantifying and coding.Frequently the frequency domain parameter of signal is also referred to as the spectral coefficient of MDCT, also may be used
For psychoacoustic model analysis, to calculate masking by noise curve.Data block after windowing is made
Fft analysis and using analysis result as the input of psychoacoustic model analysis.Further, according to
Spectral coefficient and the fft analysis result of MDCT carry out psychoacoustic model analysis, obtain noise
Masking curve and tone mask curve.Jointly calculated entirely by masking by noise curve and tone mask curve
Office's masking curve.Overall situation masking curve is the reference quantifying MDCT spectral coefficient, is to allow to draw
The thresholding of the quantizing noise entered.
4, calculate substrate curve according to overall situation masking curve, and with substrate curve, frequency spectrum parameter is entered
Row albefaction, removes and the incoherent part of audition, obtains residual signals.Specifically, application LDX
In conjunction with the concept of critical band, the method approached with sectional broken line replaces overall situation masking curve.Phase
Though than in the result of overall situation masking curve energy accurate Characterization psychoacoustic analysis, if but directly to its amount
Change, encode, then the shortcoming needing a large amount of bit, substrate curve has less segmentation in low-frequency range
Interval, it is ensured that higher frequency domain resolution, thus reduce the coding distortion of low frequency signal, and
Relatively big at the piecewise interval of high band, can preferably match with the auditory properties of people.Further,
With substrate curve, the spectral coefficient of MDCT is carried out albefaction, remove and the incoherent composition of audition.
So for stereo mode, enhance the dependency between sound channel, be more beneficial for stereo coding.
Residual signals is the MDCT spectral coefficient after albefaction, for follow-up quantization and coding.
5, residual signals is coupled, and quantization encoding is bit stream.For multichannel pattern,
Residual signals carries out coupling to improve further compression ratio, and wherein coupled modes include square pole
At least one during coordinate maps and sound channel interweaves.Carry out quantifying to compile by the residual signals after coupling again
Code, obtains coded bit stream, completes diffraction and processes.
Wherein, calculate substrate curve according to overall situation masking curve, and with substrate curve to frequency spectrum parameter
Carry out albefaction and obtain residual signals, and residual signals is coupled, and carry out quantization encoding genus
In rate distortion loop control, to obtain more preferable diffraction treatment effect.
After sound field diffraction has processed, the signal after processing diffraction carries out synthesis output.Specifically,
With band filter or low pass filter combine manipulator or with block-based conversion to diffraction at
Signal after reason synthesizes, and the signal obtaining reconstruct can form the audio signal of surround sound.
Wherein combine the method for manipulator with band filter or low pass filter to be easy in a frequency domain point
The condition of analysis signal reconstruction, the frequency response of the most each multichannel analysis/composite filter to be superposed to one normal
Number.Many classical bank of filters are had been designed that based on this principle.Block-based conversion
Method includes: after being processed by diffraction, the bit stream of output is divided into multiple input block 1,2,3,
Carry out the caching that splices.The input block quantity that tool divides is not construed as limiting, and is determined on a case-by-case basis.
Conventional conversion has discrete Fourier transform (DFT), discrete cosine transform (DCT) and discrete sine to become
Change (DST).After original time-domain signal windowing, direct transform becomes frequency domain representation, can exist well
Frequency domain is explained the characteristic of signal.During synthesis, after above-mentioned signal is inverse transformed, form time domain sequences,
With synthesis window be multiplied, then with previous part data splicing adding, obtain reconstruct signal.At this
The data block that in inventive embodiments, each is currently processed overlapping with subsequent data chunk 50%, with suitably
Window function carry out windowing process, the data block after windowing is converted.So, during synthesis, first
Carry out inverse transformation, the time domain sequences obtained is multiplied with window function, then the left side by the data block of gained
Half part and the right half part splicing adding of last data block, obtain the signal of reconstruct.The present invention is real
Execute example and all of input block is carried out above-mentioned block-based conversion, the letter finally reconstructed
Number.In embodiments of the present invention, the process being simultaneously similar to multipath audio signal, to obtain
The multipath audio signal of surround sound can be formed.So can be centered by single stereo set
The effectiveness of high-quality 3D sound is listened on a large scale in space, and in plane about and the most each
The audibility of point is uniform, hearer need not be placed in appointment position, it is achieved that the effect of surround sound
Really.
In sum, the audio signal of input is carried out modulus by analog-digital converter and turns by the present invention
Change, output digit signals;Microcontroller carries out pretreatment and forms MCVF multichannel voice frequency numeral digital signal
Signal;Each road digital audio is believed by 3D digital signal processor according to frequency and the wavelength of audio frequency
Number carry out dividing coded treatment, and carry out the digital signal after scaling down processing equalizing audio effect processing,
Form Virtual space effect by computational analysis, after diffraction processes, synthesize output can form cincture
Stereosonic multipath audio signal, it is possible to receive in the space on a large scale centered by single stereo set
Listen the effectiveness of high-quality 3D sound, and equal at the audibility of plane about and upper and lower each point
Even, hearer need not be placed in appointment position, it is achieved that the effect of surround sound.
The foregoing is only embodiments of the invention, not thereby limit the scope of the claims of the present invention,
Every equivalent structure utilizing description of the invention and accompanying drawing content to be made or equivalence flow process conversion, or
Directly or indirectly being used in other relevant technical fields, the patent being the most in like manner included in the present invention is protected
In the range of protecting.
Claims (10)
1. a 3D PA-system based on the synthesis of sound field ripple, it is characterised in that described device bag
Include:
Analog-digital converter, for carrying out analog digital conversion, output digit signals by the audio signal of input;
Microcontroller, is connected with described analog-digital converter, for described digital signal is carried out pre-place
Reason forms MCVF multichannel voice frequency digital signal;
3D digital signal processor, is connected with described microcontroller, for the frequency according to audio frequency
And audio digital signals described in each road is carried out dividing coded treatment by wavelength, and by after scaling down processing
Described digital signal carry out equalize audio effect processing, by computational analysis formed Virtual space imitate
Really, after diffraction processes, synthesis output can form the multipath audio signal of surround sound.
3D PA-system the most according to claim 1, it is characterised in that described 3D sound
Ring device and also include five loudspeaker, export five tunnel audio signals, one of them loudspeaker output 100Hz
Following audio signal, remaining two pairs of loudspeaker shares a sound chamber respectively, exports full frequency band respectively
The audio signal of left and right acoustic channels.
3D PA-system the most according to claim 2, it is characterised in that
Described audio digital signals is carried out equalizing audio effect processing by described 3D digital signal processor,
Including: described digital signal samples analysis is divided into multiple frequency range;According to sampling value of calculation and
Filter coefficient amendment algorithm is to adjust filter factor;According to the described filter factor adjusted by many
Individual wave filter is filtered calculating;
Described 3D digital signal processor is also by cascading after multiple IIR comb filter parallel connections
To do the computational analysis of Virtual space effect, the time delay of the most each wave filter and feedback increase
Benefit is different.
3D PA-system the most according to claim 1, it is characterised in that described 3D number
Word signal processor carries out diffraction process to the sound field forming Virtual space effect, including:
Obtain coding parameter and initialize;
The audio signal of input is carried out windowing process;
The discrete cosine transform improveing the described data block after windowing, obtains described audio frequency letter
Number frequency spectrum parameter, and carry out fast Fourier transform analysis, and according to frequency spectrum parameter and quickly
The output overall situation masking curve that Fourier transform is analyzed;
Calculate substrate curve according to overall situation masking curve, and with described substrate curve, described frequency spectrum is joined
Number carries out albefaction, removes and the incoherent part of audition, obtains residual signals;
Described residual signals is coupled, and quantization encoding is bit stream.
3D PA-system the most according to claim 1, it is characterised in that described 3D number
Word signal processor also combines manipulator or with based on block with band filter or low pass filter
Conversion diffraction is processed after signal synthesize, the signal obtaining reconstruct can be formed around three-dimensional
The audio signal of sound.
6. the synthetic method of a 3D PA-system based on the synthesis of sound field ripple, it is characterised in that
Described method includes:
The audio signal of input is carried out analog digital conversion, output digit signals;
Described digital signal is carried out pretreatment and forms MCVF multichannel voice frequency digital signal;
Frequency and wavelength according to audio frequency carry out frequency dividing coding to audio digital signals described in each road
Process, and carry out the described digital signal after scaling down processing equalizing audio effect processing, divide by calculating
Analysis forms Virtual space effect, and after diffraction processes, can to form surround sound many in synthesis output
Road audio signal.
Method the most according to claim 6, it is characterised in that described by after scaling down processing
Digital signal carries out equalizing the step of audio effect processing and includes:
Described digital signal samples analysis is divided into multiple frequency range;
Value of calculation and filter coefficient according to sampling revise algorithm to adjust filter factor;
It is filtered calculating by multiple wave filter according to the described filter factor adjusted.
Method the most according to claim 6, it is characterised in that described by computational analysis shape
The step becoming Virtual space effect includes:
By being formed Virtual space effect, wherein by cascade after multiple IIR comb filter parallel connections
Time delay and the feedback oscillator of each wave filter are different.
Method the most according to claim 6, it is characterised in that described to forming virtual sky
Between the sound field of effect carry out diffraction process, including:
Obtain coding parameter and initialize;
The audio signal of input is carried out windowing process;
The discrete cosine transform improveing the described data block after windowing, obtains described audio frequency letter
Number frequency spectrum parameter, and carry out fast Fourier transform analysis, and according to frequency spectrum parameter and quickly
The output overall situation masking curve that Fourier transform is analyzed;
Calculate substrate curve according to overall situation masking curve, and with described substrate curve, described frequency spectrum is joined
Number carries out albefaction, removes and the incoherent part of audition, obtains residual signals;
Described residual signals is coupled, and quantization encoding is bit stream.
Method the most according to claim 6, it is characterised in that described synthesis output can shape
Cyclization includes around the step of stereosonic multipath audio signal:
Manipulator is combined or with block-based conversion to spreading out with band filter or low pass filter
Penetrating the signal after process to synthesize, the signal obtaining reconstruct can form the audio frequency letter of surround sound
Number.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510112555.XA CN106034274A (en) | 2015-03-13 | 2015-03-13 | 3D sound device based on sound field wave synthesis and synthetic method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510112555.XA CN106034274A (en) | 2015-03-13 | 2015-03-13 | 3D sound device based on sound field wave synthesis and synthetic method |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106034274A true CN106034274A (en) | 2016-10-19 |
Family
ID=57150673
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510112555.XA Pending CN106034274A (en) | 2015-03-13 | 2015-03-13 | 3D sound device based on sound field wave synthesis and synthetic method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106034274A (en) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110580914A (en) * | 2019-07-24 | 2019-12-17 | 安克创新科技股份有限公司 | Audio processing method and equipment and device with storage function |
CN110636407A (en) * | 2018-06-21 | 2019-12-31 | 刘云轩 | Full-digital loudspeaker system and working method thereof |
CN111556405A (en) * | 2020-04-09 | 2020-08-18 | 北京金茂绿建科技有限公司 | Power amplifier chip and electronic equipment |
WO2020177095A1 (en) * | 2019-03-06 | 2020-09-10 | Harman International Industries, Incorporated | Virtual height and surround effect in soundbar without up-firing and surround speakers |
CN116437268A (en) * | 2023-06-14 | 2023-07-14 | 武汉海微科技有限公司 | Adaptive frequency division surround sound upmixing method, device, equipment and storage medium |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101014107A (en) * | 2007-01-26 | 2007-08-08 | 深圳创维-Rgb电子有限公司 | Television voice processing apparatus and method |
CN101155439A (en) * | 2006-09-30 | 2008-04-02 | 安富科技股份有限公司 | Method for intensifying sound effect processing function of Bluetooth stereo device |
CN101193090A (en) * | 2006-11-27 | 2008-06-04 | 华为技术有限公司 | Signal processing method and its device |
CN101361124A (en) * | 2006-11-27 | 2009-02-04 | 索尼计算机娱乐公司 | Audio processing device and audio processing method |
CN101511047A (en) * | 2009-03-16 | 2009-08-19 | 东南大学 | Three-dimensional sound effect processing method for double track stereo based on loudspeaker box and earphone separately |
CN101997500A (en) * | 2009-08-26 | 2011-03-30 | 展讯通信(上海)有限公司 | Audio equalization treatment system and method thereof |
US20110280407A1 (en) * | 2008-11-14 | 2011-11-17 | Scott Skinner | Compressor Based Dynamic Bass Enhancement with EQ |
CN102332266A (en) * | 2010-07-13 | 2012-01-25 | 炬力集成电路设计有限公司 | Audio data encoding method and device |
CN202353798U (en) * | 2011-12-07 | 2012-07-25 | 广州声德电子有限公司 | Audio processor of digital cinema |
CN103780214A (en) * | 2012-10-24 | 2014-05-07 | 华为终端有限公司 | Method and device for adjusting audio equalizer |
CN103916727A (en) * | 2012-12-31 | 2014-07-09 | 广州励丰文化科技股份有限公司 | Active integrated sound box with multiple digital signal processors (DSPs) |
-
2015
- 2015-03-13 CN CN201510112555.XA patent/CN106034274A/en active Pending
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101155439A (en) * | 2006-09-30 | 2008-04-02 | 安富科技股份有限公司 | Method for intensifying sound effect processing function of Bluetooth stereo device |
CN101193090A (en) * | 2006-11-27 | 2008-06-04 | 华为技术有限公司 | Signal processing method and its device |
CN101361124A (en) * | 2006-11-27 | 2009-02-04 | 索尼计算机娱乐公司 | Audio processing device and audio processing method |
CN101014107A (en) * | 2007-01-26 | 2007-08-08 | 深圳创维-Rgb电子有限公司 | Television voice processing apparatus and method |
US20110280407A1 (en) * | 2008-11-14 | 2011-11-17 | Scott Skinner | Compressor Based Dynamic Bass Enhancement with EQ |
CN101511047A (en) * | 2009-03-16 | 2009-08-19 | 东南大学 | Three-dimensional sound effect processing method for double track stereo based on loudspeaker box and earphone separately |
CN101997500A (en) * | 2009-08-26 | 2011-03-30 | 展讯通信(上海)有限公司 | Audio equalization treatment system and method thereof |
CN102332266A (en) * | 2010-07-13 | 2012-01-25 | 炬力集成电路设计有限公司 | Audio data encoding method and device |
CN202353798U (en) * | 2011-12-07 | 2012-07-25 | 广州声德电子有限公司 | Audio processor of digital cinema |
CN103780214A (en) * | 2012-10-24 | 2014-05-07 | 华为终端有限公司 | Method and device for adjusting audio equalizer |
CN103916727A (en) * | 2012-12-31 | 2014-07-09 | 广州励丰文化科技股份有限公司 | Active integrated sound box with multiple digital signal processors (DSPs) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110636407A (en) * | 2018-06-21 | 2019-12-31 | 刘云轩 | Full-digital loudspeaker system and working method thereof |
CN110636407B (en) * | 2018-06-21 | 2024-06-07 | 刘云轩 | All-digital loudspeaker system and working method thereof |
WO2020177095A1 (en) * | 2019-03-06 | 2020-09-10 | Harman International Industries, Incorporated | Virtual height and surround effect in soundbar without up-firing and surround speakers |
CN110580914A (en) * | 2019-07-24 | 2019-12-17 | 安克创新科技股份有限公司 | Audio processing method and equipment and device with storage function |
CN111556405A (en) * | 2020-04-09 | 2020-08-18 | 北京金茂绿建科技有限公司 | Power amplifier chip and electronic equipment |
CN116437268A (en) * | 2023-06-14 | 2023-07-14 | 武汉海微科技有限公司 | Adaptive frequency division surround sound upmixing method, device, equipment and storage medium |
CN116437268B (en) * | 2023-06-14 | 2023-08-25 | 武汉海微科技有限公司 | Adaptive frequency division surround sound upmixing method, device, equipment and storage medium |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101390443B (en) | Audio encoding and decoding | |
JP4606507B2 (en) | Spatial downmix generation from parametric representations of multichannel signals | |
RU2381571C2 (en) | Synthesisation of monophonic sound signal based on encoded multichannel sound signal | |
TW412719B (en) | Method and apparatus for reproducing speech signals and method for transmitting same | |
US8917874B2 (en) | Method and apparatus for decoding an audio signal | |
EP2107833B1 (en) | Audio wave field encoding | |
CN104285390B (en) | The method and device that compression and decompression high-order ambisonics signal are represented | |
KR101589942B1 (en) | Cross product enhanced harmonic transposition | |
CN101067931B (en) | Efficient configurable frequency domain parameter stereo-sound and multi-sound channel coding and decoding method and system | |
RU2376726C2 (en) | Device and method for generating encoded stereo signal of audio part or stream of audio data | |
CN102158198B (en) | Filter generator, filter system and method for providing intermediate filters defined signal | |
Necciari et al. | The ERBlet transform: An auditory-based time-frequency representation with perfect reconstruction | |
CN101253556B (en) | Energy shaping device and energy shaping method | |
CN106034274A (en) | 3D sound device based on sound field wave synthesis and synthetic method | |
CN101401455A (en) | Binaural rendering using subband filters | |
CN108780649A (en) | Use the device and method of broadband alignment parameter and multiple narrowband alignment parameters coding or decoding multi-channel signal | |
CN102157156B (en) | Single-channel voice enhancement method and system | |
CN104854655A (en) | Method and apparatus for compressing and decompressing higher order ambisonics representation for sound field | |
MX2012010416A (en) | Apparatus and method for processing an audio signal using patch border alignment. | |
CN108600935A (en) | Acoustic signal processing method and equipment | |
CN104103276B (en) | Sound coding device, sound decoding device, sound coding method and sound decoding method | |
CN107005778A (en) | The audio signal processing apparatus and method rendered for ears | |
US9595267B2 (en) | Method and apparatus for decoding an audio signal | |
CN107610710B (en) | Audio coding and decoding method for multiple audio objects | |
CN109448741A (en) | A kind of 3D audio coding, coding/decoding method and device |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20161019 |