CN105244040A - Audio signal consistency comparison method - Google Patents

Audio signal consistency comparison method Download PDF

Info

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
Application number
CN201510429221.5A
Other languages
Chinese (zh)
Other versions
CN105244040B (en
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.)
Hangzhou grandwell Polytron Technologies Inc
Original Assignee
HANGZHOU LINKER DIGITAL TECHNOLOGY 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 HANGZHOU LINKER DIGITAL TECHNOLOGY Co Ltd filed Critical HANGZHOU LINKER DIGITAL TECHNOLOGY Co Ltd
Priority to CN201510429221.5A priority Critical patent/CN105244040B/en
Publication of CN105244040A publication Critical patent/CN105244040A/en
Application granted granted Critical
Publication of CN105244040B publication Critical patent/CN105244040B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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

The control methods of a kind of sound signal consistance
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%.
X ( k ) = s u m x ( n ) * exp ( - j * 2 * π ( k - 1 ) * n - 1 n )
1≤K≤N
x ( n ) = ( 1 N ) s u m x ( k ) * exp ( j * 2 * π * ( k - 1 ) * n - 1 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
o f f s e t > N 2 ? o f f s e t = o f f s e t - N
r a t i o = r a t i o * N s u m
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
X [ j ] = max i ≤ k ≤ i f S a m p l e [ k ]
1≤i≤N;i+=96;j++
sum=sum+|fSample[i]|
1≤i≤N
m e a n X = s u m 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.
CN201510429221.5A 2015-07-20 2015-07-20 A kind of audio signal consistency control methods Active CN105244040B (en)

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)

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

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

Patent Citations (5)

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

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