CN105244040A - Audio signal consistency comparison method - Google Patents
Audio signal consistency comparison method Download PDFInfo
- Publication number
- CN105244040A CN105244040A CN201510429221.5A CN201510429221A CN105244040A CN 105244040 A CN105244040 A CN 105244040A CN 201510429221 A CN201510429221 A CN 201510429221A CN 105244040 A CN105244040 A CN 105244040A
- Authority
- CN
- China
- Prior art keywords
- data
- sampled data
- similarity
- tunnel
- sound signal
- 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.)
- Granted
Links
Landscapes
- Stereo-Broadcasting Methods (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
The present invention discloses an audio signal consistency comparison method which comprises the steps of obtaining sampling data, expanding a length, calculating a fast correlation similarity, calculating a similarity and an offset, judging that two paths of data are consistent if the similarity is larger than or equal to a threshold value, judging whether the sampling data is subjected to approximated envelope line processing or not if the similarity is smaller than the threshold value, judging that the two paths of data are inconsistent if the sampling data is subjected to approximated envelope line processing, and carrying out approximated envelope line processing on the sampling data if the sampling data is not subjected to approximated envelope line processing and then repeating the above steps. According to the scheme of the invention, the consistency of two paths of audio signals can be accurately judged, and the method is suitable for the monitoring of a broadcast signal by a radio station.
Description
Technical field
The present invention relates to audio data analysis process field, especially relate to the control methods of a kind of sound signal consistance.
Background technology
In whole transmission in broadcasting station, monitoring link, be all carry out detecting for the physical index of audio frequency, such as: whether break, level is on the low side or higher, anti-phase etc.But whether correctly detection means is lacked for the content of audio frequency, whether the signal that such as master control finally exports is the signal that studio's sound console exports, whether air-launched signal is that the problems such as the signal that exports of master control final stage do not have suitable detection means at present, to only have etc. to there occurs accident and just can detect.
State Intellectual Property Office of the People's Republic of China discloses on Dec 01st, 2010 patent documentation that publication number is CN101902677A, and title is audio detection device and method, and this device comprises: MIC electric circuit inspection module and phonation circuit detection module; According to amplitude and frequency, MIC electric circuit inspection module, for obtaining amplitude and the frequency of the sine wave signal after MIC circuit sampling, determines that whether MIC circuit is qualified, and sine wave output signal; Phonation circuit detection module comprises: change-over panel and the processor with MIC circuit; Change-over panel, for the simulating signal that phonation circuit of sampling exports, and is converted to digital signal by simulating signal, is sent to processor; Processor, for carrying out analyzing and processing to the amplitude of digital signal and frequency, judges that whether phonation circuit is qualified.This scheme cannot carry out conforming judgement to sound signal.
Summary of the invention
The present invention mainly solves the technical matters can not carrying out accurately judgement to sound signal consistance existing for prior art, a kind of characteristic signal analysis can passed through audio frequency is provided, detect the sound signal consistance control methods whether two paths of signals content is consistent, be convenient to broadcasting station Timeliness coverage and broadcast fault.
The present invention is directed to that above-mentioned technical matters mainly solved by following technical proposals: the control methods of a kind of sound signal consistance, comprises the following steps:
S01, from first via sound signal, obtain N number of sampled data as first via sampled data, from the second tunnel sound signal, obtain N number of sampled data as the second tunnel sampled data;
S02, by mending 0 after first via sampled data, to extend to length be 2N, and by mending 0 after the second tunnel sampled data, to extend to length be 2N;
S03, the first via sampled data after expansion is carried out to Fast Fourier Transform (FFT) and then asked conjugation to obtain first via data, then carry out Fast Fourier Transform (FFT) to the second tunnel sampled data after expansion asks conjugation to obtain the second circuit-switched data, ask the cross correlation function of first via data and the second circuit-switched data, carrying out Fast Fourier Transform (FFT) to cross correlation function, obtain similarity data;
S04, choose the maximal value of absolute value in similarity data, then calculate Similarity value and maximal value side-play amount;
If S05 Similarity value is more than or equal to consistance threshold value, then think that first via sampled data is consistent with the second tunnel sampled data, contrast flow process terminates; If similarity is less than consistance threshold value, then enter step S06;
S06, judge that whether first via sampled data and the second tunnel sampled data are through approximate envelop process, if do not had, then carry out approximate envelop process to first via sampled data and the second tunnel sampled data, then repeat step S02 to S05; If through approximate envelop process, then judge first via sampled data and the second tunnel sampled data inconsistent.
Said process more adequately can judge the similarity of first via sampled signal and the second tunnel sampled signal.According to the side-play amount calculated to after sound signal skew, the consistance of sound signal more adequately can be calculated.
As preferably, Similarity value ratio computing formula is as follows:
ratio=fSample
max*M/sum
In formula, fSample
maxfor the maximal value of absolute value in similarity data, M is the number of similarity data, sum be the absolute value of all similarity datas and.
As preferably, maximal value side-play amount offset determines according to following methods: fSample
maxfor the data of i-th in similarity data, if i is less than or equal to M/2, then offset=i; If i is greater than M/2, then offset=i-M/2; FSample
maxfor the maximal value of absolute value in similarity data, M is the number of similarity data.
As preferably, approximate envelop process is specially: first via sampled data is divided into the window that length is 2 milliseconds, gets maximal value and form first window maximum data * X1 from each window; Second tunnel sampled data is divided into the window that length is 2 milliseconds, from each window, gets maximal value form Second Window maximum data * Y1; Calculate the average meanX1 of the absolute value of first via sampled data; Calculate the average meanY1 of the absolute value of the second tunnel sampled data; Each data in * X1 are deducted meanX1, and the data obtained are as first via sampled data; Each data in * Y1 are deducted meanY1, and the data obtained are as the second tunnel sampled data.
For volume smaller and make an uproar in the end obvious audio content, the similarity-rough set calculated directly carrying out step S02-S05 is little, generally can be less than consistance threshold value, can eliminate small-signal interference, improve differentiation rate after now carrying out approximate envelop process.
As preferably, the span of described consistance threshold value is 20-30.
By consistance threshold value accurately can be oriented after a certain amount of sampling analysis.
As preferably, when continuous three times judge that first via sampled data and the inconsistent then identification first via sound signal of the second tunnel sampled data and the second tunnel sound signal are inconsistent.
When there is first via sound signal and the inconsistent situation of the second tunnel sound signal, can carry out early warning, prompting staff carries out revising or safeguarding.
This programme may be used for following occasion:
Radio station master control:
For judging the consistance of master control final stage output signal and sound console output signal;
Master control final stage signal and empty consistance of to collect mail number;
Transmitting station:
Empty to collect mail number and the consistance of optical transmitter and receiver output signal;
Monitoring center:
User judges that whether broadcasting signal is disturbed or invade;
The substantial effect that the present invention brings is, reliably can calculate similarity and the side-play amount of two-way audio signal, eliminates the impact of sound signal decay, has high accuracy of judgement degree.
Accompanying drawing explanation
Fig. 1 is a kind of process flow diagram of the present invention.
Embodiment
Below by embodiment, and by reference to the accompanying drawings, technical scheme of the present invention is described in further detail.
Embodiment: a kind of sound signal consistance control methods of the present embodiment, as shown in Figure 1, comprises the following steps:
Record the voice data of certain sampled point number, first extension length is twice, is 0 by assignment below.
floatfSample[N*2]
fSample[i]=0N≤i≤N*2
Then sampled data is carried out FFT conversion, fft algorithm is described below:
Strange according to discrete fourier transform, even, empty, the characteristics such as reality, improvement is carried out to the algorithm of discrete Fourier transform (DFT) and obtains fast Fourier transformation algorithm FFT, if x (n) is the sequence of complex numbers of N item, converted by DFT, the calculating of arbitrary X (m) all needs N complex multiplication and N-1 complex addition, and a complex multiplication equals four real multiplications and twice real addition, one time complex addition equals twice real addition, even if a complex multiplication and a complex addition are defined as once " computing " (four real multiplications and four real additions), so obtain the X (m) of N item sequence of complex numbers, namely N point DFT conversion approximately just needs N^2 computing.In FFT, utilize periodicity and the symmetry of WN, a N item sequence (is established N=2k, k is positive integer), be divided into the subsequence of two N/2 items, each N/2 point DFT converts needs (N/2) 2 computings, then the DFT of two N/2 points conversion is combined into the DFT conversion of a N point with N computing.After such conversion, total operation times just becomes N+2* (N/2) ^2=N+ (N^2)/2, saves the operand of about 50%.
1≤K≤N
1≤n≤N
By the two-way sampled data after FFT budget, ask conjugation, then carry out a FFT computing, process is called calculating fast correlation similarity.
Then the data of Fourier transform are got maximal value, calculate maximum average value by number of samples, and calculate maximal value side-play amount by the start bit of maximal value, so just can calculate the sampled point side-play amount of Similarity value and two-way audio, computing method are:
Sampled point number: N
Side-play amount: offset
Similarity: ratio
Data after Fourier transform: fSample, because sound signal is sinusoidal wave, have positive and negative values, so absolute value should be used to carry out size contrast when calculating.
Data volume summation: sum
ratio<|fSample[i]|?ratio=|fSample[i]|offset=i
sum=sum+|fSample[i]|
1≤i≤N
After similarity (ratio) value is greater than certain value range (this value is commonly defined as 20-30 after long-term test), just two paths of signals can be regarded as consistent, and then according to the side-play amount calculated (offset), carry out the skew of two paths of signals, migration algorithm is:
fSample[i]=fSample[i+offset]
1≤i≤(N-offset)
Follow-uply calculating is carried out to shifted data more will press close to contrasting region by Similarity value, improve accuracy.
Smaller for volume, the obvious audio content and make an uproar in the end, the approximate value drawn in current algorithm is smaller, generally can be less than predetermined scope, this value fluctuation situation is obvious, and Similarity value is generally at about 10-15, during this situation, adopt the mode of getting the approximate envelop of voice band line to carry out secondary calculating, small-signal interference can be eliminated and improve differentiation rate.Get the maximal value of 2 milliseconds of windows at every turn, obtain * X1, * Y1, and the average of two paths of signals:
Setting sampling rate is 48000, and so the window sample number of 2 milliseconds is 96
1≤i≤N;i+=96;j++
sum=sum+|fSample[i]|
1≤i≤N
Average in * X1, * Y1 signal is removed:
rX=*X1-meanX1
rY=*Y1-meanY1
Then to data rX, the rY of secondary calculating again before step, calculate fast correlation similarity, the data calculated carried out again to the calculating of similarity and side-play amount, obtain new comparing result, improve accuracy further.The frequency occurred in actual environment due to this class audio frequency is not high, and twice calculating can promote CPU burden, therefore secondary authentication policy when being lower as Similarity value exists.
The basis of fast correlation similarity adds determination strategy again, and continuous 3 decision threshold are less than setting value (being generally 20-30), and think that two paths of signals is inconsistent during side-play amount fluctuation range comparatively large (being greater than the sampling of 1s).
After fast correlation Similarity measures, similarity and the side-play amount of two-way audio can be calculated reliably, again by after a certain amount of sampling analysis, just can orient decision threshold accurately, just can effectively detect the consistance of two paths of signals according to this threshold values, early warning, and owing to having carried out twice FFT computing, get the calculating that similarity is the signal curve carried out, the decay of sound signal does not affect result of calculation, and accuracy of judgement degree can reach more than 99.99%.
Specific embodiment described herein is only to the explanation for example of the present invention's spirit.Those skilled in the art can make various amendment or supplement or adopt similar mode to substitute to described specific embodiment, but can't depart from spirit of the present invention or surmount the scope that appended claims defines.
Although more employ the term such as sampled data, Fast Fourier Transform (FFT) herein, do not get rid of the possibility using other term.These terms are used to be only used to describe and explain essence of the present invention more easily; The restriction that they are construed to any one additional is all contrary with spirit of the present invention.
Claims (6)
1. the control methods of sound signal consistance, is characterized in that, comprises the following steps:
S01, from first via sound signal, obtain N number of sampled data as first via sampled data, from the second tunnel sound signal, obtain N number of sampled data as the second tunnel sampled data;
S02, by mending 0 after first via sampled data, to extend to length be 2N, and by mending 0 after the second tunnel sampled data, to extend to length be 2N;
S03, the first via sampled data after expansion is carried out to Fast Fourier Transform (FFT) and then asked conjugation to obtain first via data, then carry out Fast Fourier Transform (FFT) to the second tunnel sampled data after expansion asks conjugation to obtain the second circuit-switched data, ask the cross correlation function of first via data and the second circuit-switched data, carrying out Fast Fourier Transform (FFT) to cross correlation function, obtain similarity data;
S04, choose the maximal value of absolute value in similarity data, then calculate Similarity value and maximal value side-play amount;
If S05 Similarity value is more than or equal to consistance threshold value, then think that first via sampled data is consistent with the second tunnel sampled data, contrast flow process terminates; If similarity is less than consistance threshold value, then enter step S06;
S06, judge that whether first via sampled data and the second tunnel sampled data are through approximate envelop process, if do not had, then carry out approximate envelop process to first via sampled data and the second tunnel sampled data, then repeat step S02 to S05; If through approximate envelop process, then judge first via sampled data and the second tunnel sampled data inconsistent.
2. a kind of sound signal consistance according to claim 1 control methods, is characterized in that, Similarity value ratio computing formula is as follows:
ratio=fSample
max*M/sum
In formula, fSample
maxfor the maximal value of absolute value in similarity data, M is the number of similarity data, sum be the absolute value of all similarity datas and.
3. a kind of sound signal consistance according to claim 1 control methods, is characterized in that, maximal value side-play amount offset determines according to following methods: fSample
maxfor the data of i-th in similarity data, if i is less than or equal to M/2, then offset=i; If i is greater than M/2, then offset=i-M/2; FSample
maxfor the maximal value of absolute value in similarity data, M is the number of similarity data.
4. a kind of sound signal consistance control methods according to Claims 2 or 3, it is characterized in that, approximate envelop process is specially: first via sampled data is divided into the window that length is 2 milliseconds, gets maximal value and form first window maximum data * X1 from each window; Second tunnel sampled data is divided into the window that length is 2 milliseconds, from each window, gets maximal value form Second Window maximum data * Y1; Calculate the average meanX1 of the absolute value of first via sampled data; Calculate the average meanY1 of the absolute value of the second tunnel sampled data; Each data in * X1 are deducted meanX1, and the data obtained are as first via sampled data; Each data in * Y1 are deducted meanY1, and the data obtained are as the second tunnel sampled data.
5. a kind of sound signal consistance according to claim 1 control methods, is characterized in that, the span of described consistance threshold value is 20-30.
6. a kind of sound signal consistance according to claim 1 control methods, is characterized in that, when continuous three times judge that first via sampled data and the inconsistent then identification first via sound signal of the second tunnel sampled data and the second tunnel sound signal are inconsistent.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510429221.5A CN105244040B (en) | 2015-07-20 | 2015-07-20 | A kind of audio signal consistency control methods |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510429221.5A CN105244040B (en) | 2015-07-20 | 2015-07-20 | A kind of audio signal consistency control methods |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105244040A true CN105244040A (en) | 2016-01-13 |
CN105244040B CN105244040B (en) | 2018-12-11 |
Family
ID=55041659
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510429221.5A Active CN105244040B (en) | 2015-07-20 | 2015-07-20 | A kind of audio signal consistency control methods |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105244040B (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107274911A (en) * | 2017-05-03 | 2017-10-20 | 昆明理工大学 | A kind of similarity analysis method based on sound characteristic |
CN108259104A (en) * | 2017-12-08 | 2018-07-06 | 中国航空工业集团公司成都飞机设计研究所 | A kind of double remaining transmitting station means of voting |
CN109658941A (en) * | 2018-12-03 | 2019-04-19 | 四川虹美智能科技有限公司 | A kind of intelligent refrigerator information flow method for pushing and device based on Application on Voiceprint Recognition |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3919479A (en) * | 1972-09-21 | 1975-11-11 | First National Bank Of Boston | Broadcast signal identification system |
CN102142257A (en) * | 2010-12-28 | 2011-08-03 | 北大方正集团有限公司 | Audio signal processing method and device |
CN103426439A (en) * | 2013-05-08 | 2013-12-04 | 杭州联汇数字科技有限公司 | Method for detecting consistency of broadcast television audio signal content |
CN104468697A (en) * | 2014-10-10 | 2015-03-25 | 浙江广播电视集团 | Radio station data transmission load distribution method |
CN104505101A (en) * | 2014-12-24 | 2015-04-08 | 北京巴越赤石科技有限公司 | Real-time audio comparison method |
-
2015
- 2015-07-20 CN CN201510429221.5A patent/CN105244040B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3919479A (en) * | 1972-09-21 | 1975-11-11 | First National Bank Of Boston | Broadcast signal identification system |
CN102142257A (en) * | 2010-12-28 | 2011-08-03 | 北大方正集团有限公司 | Audio signal processing method and device |
CN103426439A (en) * | 2013-05-08 | 2013-12-04 | 杭州联汇数字科技有限公司 | Method for detecting consistency of broadcast television audio signal content |
CN104468697A (en) * | 2014-10-10 | 2015-03-25 | 浙江广播电视集团 | Radio station data transmission load distribution method |
CN104505101A (en) * | 2014-12-24 | 2015-04-08 | 北京巴越赤石科技有限公司 | Real-time audio comparison method |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107274911A (en) * | 2017-05-03 | 2017-10-20 | 昆明理工大学 | A kind of similarity analysis method based on sound characteristic |
CN108259104A (en) * | 2017-12-08 | 2018-07-06 | 中国航空工业集团公司成都飞机设计研究所 | A kind of double remaining transmitting station means of voting |
CN109658941A (en) * | 2018-12-03 | 2019-04-19 | 四川虹美智能科技有限公司 | A kind of intelligent refrigerator information flow method for pushing and device based on Application on Voiceprint Recognition |
CN109658941B (en) * | 2018-12-03 | 2021-02-19 | 四川虹美智能科技有限公司 | Intelligent refrigerator information flow pushing method and device based on voiceprint recognition |
Also Published As
Publication number | Publication date |
---|---|
CN105244040B (en) | 2018-12-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11909563B2 (en) | Method and apparatus for modulation recognition of signals based on cyclic residual network | |
US20200044899A1 (en) | Method for automatically identifying modulation mode for digital communication signal | |
CN101494508A (en) | Frequency spectrum detection method based on characteristic cyclic frequency | |
CN111832462B (en) | Frequency hopping signal detection and parameter estimation method based on deep neural network | |
CN105244040A (en) | Audio signal consistency comparison method | |
CN106357575A (en) | Multi-parameter jointly-estimated interference type identification method | |
CN112949387A (en) | Intelligent anti-interference target detection method based on transfer learning | |
CN103809198A (en) | Method for measuring particle accelerator beam position | |
CN106197650A (en) | A kind of method judging that audio signal is noiseless | |
CN105429719A (en) | Strong interference signal detection method based on power spectrum and multiple dimensioned wavelet transformation analysis | |
CN104392086B (en) | A kind of signal deteching circuit and method based on Pearson came rand variate coefficient correlation | |
CN102914325A (en) | Dissipation synchronization-based detection method of small signal under chaos background | |
CN112764003A (en) | Radar radiation source signal time-frequency feature identification method and device and storage medium | |
CN113300986B (en) | Unmanned aerial vehicle image transmission signal and hotspot signal identification method, medium and computer equipment | |
CN101588191A (en) | Method and device for radio signal recognition | |
CN102075202B (en) | Characteristic value-based passive channel interference detection method | |
Gravelle et al. | Deep learning-enabled real-time recognition of wireless signals | |
CN104868962A (en) | Spectrum detection method and device based on compressed sensing | |
CN104734793A (en) | Energy detection method based on p-power wireless cooperative spectrum sensing | |
CN108039182A (en) | A kind of voice-activation detecting method | |
CN106771928A (en) | A kind of online pick-up method of partial discharge pulse's initial time | |
CN104333423B (en) | Signal detection circuit and method based on Gini Correlation | |
Wu et al. | Interference detection and recognition based on signal reconstruction using recurrent neural network | |
JP5282914B2 (en) | Radio identification device | |
CN107341519A (en) | A kind of SVMs identification optimization method based on multiresolution analysis |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CP03 | Change of name, title or address | ||
CP03 | Change of name, title or address |
Address after: 310000 Fourth and Five Floors of Building 399, Qiuyi Road, Changhe Street, Binjiang District, Hangzhou City, Zhejiang Province Patentee after: Hangzhou grandwell Polytron Technologies Inc Address before: 310052 Fourth Floor, Building 399 Qiuyi Road, Changhe Street, Binjiang District, Hangzhou City, Zhejiang Province Patentee before: Hangzhou Linker Digital Technology Co., Ltd. |