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 PDF

Info

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
Application number
CN201510112555.XA
Other languages
Chinese (zh)
Inventor
晏明峰
约夏·昆兹
曾世奇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shenzhen Acemile Electronic Co Ltd
Original Assignee
Shenzhen Acemile Electronic Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shenzhen Acemile Electronic Co Ltd filed Critical Shenzhen Acemile Electronic Co Ltd
Priority to CN201510112555.XA priority Critical patent/CN106034274A/en
Publication of CN106034274A publication Critical patent/CN106034274A/en
Pending legal-status Critical Current

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

3D PA-system based on the synthesis of sound field ripple and synthetic method thereof
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:
H ( z ) = a 7 ( z - D 1 1 - a 1 z - D 1 + z - D 2 1 - a 2 z - D 2 + z - D 3 1 - a 3 z - D 3 + z - D 4 1 - a 4 z - D 4 ) · a 5 + z - D 5 1 + a 5 z - D 5 · a 6 + z - D 6 1 + a 6 z - D 6 + 1
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.
CN201510112555.XA 2015-03-13 2015-03-13 3D sound device based on sound field wave synthesis and synthetic method Pending CN106034274A (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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)

Patent Citations (11)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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