WO2015050006A1 - 音場測定装置、音場測定方法および音場測定プログラム - Google Patents
音場測定装置、音場測定方法および音場測定プログラム Download PDFInfo
- Publication number
- WO2015050006A1 WO2015050006A1 PCT/JP2014/074848 JP2014074848W WO2015050006A1 WO 2015050006 A1 WO2015050006 A1 WO 2015050006A1 JP 2014074848 W JP2014074848 W JP 2014074848W WO 2015050006 A1 WO2015050006 A1 WO 2015050006A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- frequency
- signal
- measurement signal
- sound field
- measurement
- Prior art date
Links
- 238000000034 method Methods 0.000 title claims description 59
- 238000005259 measurement Methods 0.000 claims abstract description 286
- 230000007274 generation of a signal involved in cell-cell signaling Effects 0.000 claims abstract description 18
- 239000002131 composite material Substances 0.000 claims abstract description 8
- 238000001514 detection method Methods 0.000 claims description 32
- 238000004364 calculation method Methods 0.000 claims description 21
- 238000000605 extraction Methods 0.000 claims description 14
- 238000000691 measurement method Methods 0.000 claims description 14
- 239000000284 extract Substances 0.000 claims description 9
- 230000002194 synthesizing effect Effects 0.000 claims description 5
- 230000008569 process Effects 0.000 description 45
- 238000012935 Averaging Methods 0.000 description 35
- 230000006870 function Effects 0.000 description 35
- 238000001228 spectrum Methods 0.000 description 35
- 238000012545 processing Methods 0.000 description 32
- 230000015572 biosynthetic process Effects 0.000 description 6
- 238000003786 synthesis reaction Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 238000009434 installation Methods 0.000 description 5
- 230000003111 delayed effect Effects 0.000 description 4
- 101100328887 Caenorhabditis elegans col-34 gene Proteins 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- FFBHFFJDDLITSX-UHFFFAOYSA-N benzyl N-[2-hydroxy-4-(3-oxomorpholin-4-yl)phenyl]carbamate Chemical compound OC1=C(NC(=O)OCC2=CC=CC=C2)C=CC(=C1)N1CCOCC1=O FFBHFFJDDLITSX-UHFFFAOYSA-N 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000007429 general method Methods 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000008054 signal transmission Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000005236 sound signal Effects 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04R—LOUDSPEAKERS, MICROPHONES, GRAMOPHONE PICK-UPS OR LIKE ACOUSTIC ELECTROMECHANICAL TRANSDUCERS; DEAF-AID SETS; PUBLIC ADDRESS SYSTEMS
- H04R29/00—Monitoring arrangements; Testing arrangements
- H04R29/001—Monitoring arrangements; Testing arrangements for loudspeakers
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
- G10L21/00—Speech or voice signal processing techniques to produce another audible or non-audible signal, e.g. visual or tactile, in order to modify its quality or its intelligibility
- G10L21/06—Transformation of speech into a non-audible representation, e.g. speech visualisation or speech processing for tactile aids
- G10L21/10—Transforming into visible information
- G10L21/14—Transforming into visible information by displaying frequency domain information
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
- G10L21/00—Speech or voice signal processing techniques to produce another audible or non-audible signal, e.g. visual or tactile, in order to modify its quality or its intelligibility
- G10L21/06—Transformation of speech into a non-audible representation, e.g. speech visualisation or speech processing for tactile aids
- G10L21/18—Details of the transformation process
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H3/00—Measuring characteristics of vibrations by using a detector in a fluid
- G01H3/10—Amplitude; Power
- G01H3/12—Amplitude; Power by electric means
- G01H3/125—Amplitude; Power by electric means for representing acoustic field distribution
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04R—LOUDSPEAKERS, MICROPHONES, GRAMOPHONE PICK-UPS OR LIKE ACOUSTIC ELECTROMECHANICAL TRANSDUCERS; DEAF-AID SETS; PUBLIC ADDRESS SYSTEMS
- H04R29/00—Monitoring arrangements; Testing arrangements
- H04R29/008—Visual indication of individual signal levels
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04S—STEREOPHONIC SYSTEMS
- H04S1/00—Two-channel systems
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04S—STEREOPHONIC SYSTEMS
- H04S2400/00—Details of stereophonic systems covered by H04S but not provided for in its groups
- H04S2400/15—Aspects of sound capture and related signal processing for recording or reproduction
Definitions
- the present invention relates to a sound field measuring device, a sound field measuring method, and a sound field measuring program. More specifically, the present invention can measure frequency characteristics quickly and accurately in a sound field environment in which a pair of speakers are installed at a narrow interval.
- the present invention relates to a possible sound field measuring device, a sound field measuring method, and a sound field measuring program.
- frequency characteristics are measured in a sound field environment where speakers are installed, and the equalizer of the audio device is adjusted based on the measured frequency characteristics, or the output sound is corrected in advance according to the sound field.
- the equalizer of the audio device is adjusted based on the measured frequency characteristics, or the output sound is corrected in advance according to the sound field.
- an M-sequence code Maximum Length Sequence
- a TSP signal Time Stretched Pulse
- a method for obtaining the frequency characteristics of the sound field environment is used (see, for example, Patent Document 1 and Patent Document 2). It is also possible to obtain an impulse response by obtaining a cross-correlation characteristic with a measurement signal recorded by a microphone using the output measurement signal as a reference.
- the frequency characteristics of the sound field environment vary depending on the listening position with respect to the speaker. Furthermore, the measurement accuracy of frequency characteristics tends to vary depending on the distance between installed speakers and the like. For example, as shown in FIG. 14, in the case of a portable audio device 102 in which a music playback function unit 100 and a pair of right speaker 101a and left speaker 101b installed on the left and right sides are integrated, the left and right speakers 101a are generally used. , 101b tends to be narrower. Further, in the portable audio device 102, the device main body is placed near the listener because of its portability and the like, and the distance between the portable audio device 102 and the listener tends to be close.
- the measurement signal is output simultaneously from the right speaker 101a and the left speaker 101b using the same M-sequence code as the measurement signal, as shown in FIG. Then, frequency characteristics between the left and right speakers 101a and 101b and the microphone 7 are measured. Then, since the installation interval of the left and right speakers 101a and 101b is narrow and the distance to the microphone 7 is short, the delay time difference between channels (signal transmission path from one speaker to the microphone 7) tends to be short. Thus, when the delay time difference between channels becomes short, inter-channel interference occurs, and a large dip tends to occur at a specific frequency. Since the position of the dip largely fluctuates according to the delay time difference, there is a problem that the frequency characteristic changes depending on the measurement position, and the measurement accuracy may be greatly deteriorated.
- a measurement signal is output from only the left speaker 101b and the frequency characteristic is measured, and then the measurement signal is output only from the right speaker 101a.
- a method is known in which the frequency characteristics are measured by performing time division on the L channel and the R channel, respectively, by measuring the frequency characteristics. By measuring the frequency characteristics by performing time division in this way, interference between channels and interference between codes can be avoided.
- the measurement time increases and the measurement processing load increases.
- the present invention has been made in view of the above problems, and by outputting a measurement signal simultaneously from a pair of speakers installed at a narrow interval, a sound field that can accurately measure the frequency characteristics of the sound field environment. It is an object to provide a measuring device, a sound field measuring method, and a sound field measuring program.
- the sound field measuring apparatus collects output sound output from an audio device in which a pair of speakers are installed at a narrow interval, thereby improving the frequency characteristics of the sound field environment.
- a sound field measuring device to be obtained wherein a low-pass filter unit that extracts a low-frequency component of a first measurement signal, a high-pass filter unit that extracts a middle-high frequency component of a second measurement signal different from the first measurement signal, A synthesized signal generating unit that generates a synthesized signal by synthesizing the low frequency component of the first measurement signal extracted by the low pass filter unit and the middle and high frequency component of the second measurement signal extracted by the high pass filter unit; The first measurement signal before the low-frequency component is extracted by the low-pass filter unit is output from one of the pair of speakers and generated by the synthesized signal generation unit.
- An external output unit for outputting the first measurement signal and the synthesized signal to the audio device, and the pair of speakers for outputting the synthesized signal simultaneously from one of the pair of speakers.
- the sound field measuring method is a sound field measuring apparatus for obtaining frequency characteristics of a sound field environment by collecting output sound output from an audio device in which a pair of speakers are installed at a narrow interval.
- a low-frequency component extraction step in which a low-pass filter unit extracts a low-frequency component of a first measurement signal, and a high-pass filter unit that extracts a middle-high frequency component of a second measurement signal different from the first measurement signal.
- the middle and high frequency components extracted in the low frequency component extraction step, the low frequency components of the first measurement signal extracted in the low frequency component extraction step, and the middle and high frequency components of the second measurement signal extracted in the middle and high frequency component extraction step And a synthesized signal generating step in which a synthesized signal generating unit generates a synthesized signal, and the first measurement signal before the low frequency component is extracted in the low frequency component extracting step,
- an external output unit outputs to the audio device
- An external output step for outputting the first measurement signal and the combined signal; and a sound collection step for collecting the first measurement signal and the combined signal simultaneously output from the pair of speakers;
- the Fourier transform unit has a Fourier transform step for obtaining a frequency characteristic in the sound field environment by performing a Fourier transform on the signal collected in the sound collecting step.
- the sound field measurement program is a sound field measuring apparatus for obtaining frequency characteristics of a sound field environment by collecting output sound output from an audio device in which a pair of speakers are installed at a narrow interval.
- a sound field measurement program comprising: a low-pass filter function that causes a computer of the sound field measurement device to extract a low-frequency component of a first measurement signal; and a medium-high frequency component of a second measurement signal that is different from the first measurement signal.
- a Fourier transform function for obtaining a frequency characteristic in the sound field environment by performing a Fourier transform on the obtained signal.
- the first measurement signal is output from one of the speakers of the audio device, and the low frequency is the first measurement signal from the other speaker.
- a composite signal is output in which the middle and high frequencies are composed of the second measurement signal. For this reason, in the mid-high range, frequency characteristics can be measured using a stereo measurement signal, and the occurrence of dip can be suppressed.
- frequency characteristics can be measured using a monaural measurement signal in a low frequency range, so that the intersymbol interference of the measurement signal Can be suppressed.
- the sound field measurement device, the sound field measurement method, and the sound field measurement program according to the present invention since the measurement of the frequency characteristics can be performed using the stereo measurement signal, the occurrence of dip is suppressed. It becomes possible to do. In addition, since it is possible to measure a low frequency characteristic using a monaural measurement signal, intersymbol interference of the measurement signal can be suppressed. Furthermore, in the sound field measurement device, the sound field measurement method, and the sound field measurement program according to the present invention, the first measurement signal and the synthesized signal are simultaneously output from the pair of speakers and the frequency characteristics are measured. Compared with the case where measurement signals are alternately output and measurement is performed, the measurement burden can be reduced and the measurement speed can be improved.
- an M-sequence code is used as the first measurement signal
- an M-sequence code orthogonal to the M-sequence code of the first measurement signal is used as the second measurement signal.
- the average value of the signal level in a predetermined second frequency interval is calculated for each frequency shorter than the second frequency interval based on a frequency characteristic composed of a maximum value detector to be obtained and the maximum value detected by the maximum value detector. It may have an average value calculation part which calculates a frequency characteristic in the sound field environment by calculating while shifting.
- an M-sequence code is used as the first measurement signal
- an M-sequence code orthogonal to the M-sequence code of the first measurement signal is used as the second measurement signal.
- the maximum value detecting unit Based on the frequency characteristic obtained in the step, the maximum value detecting unit detects the maximum value of the signal level in the predetermined first frequency interval while shifting the frequency level by a frequency shorter than the first frequency interval, thereby A maximum value detecting step for obtaining a frequency characteristic composed of a value, and an average value calculating unit, based on the frequency characteristic composed of the maximum value obtained in the maximum value detecting step, calculates an average of signal levels in a predetermined second frequency interval
- the frequency characteristic in the sound field environment is calculated by shifting the value by a frequency shorter than the second frequency interval. It may have an average value calculation step of obtaining a.
- an M-sequence code is used as the first measurement signal
- an M-sequence code orthogonal to the M-sequence code of the first measurement signal is used as the second measurement signal.
- the maximum value of the signal level in the predetermined first frequency interval is detected while being shifted by a frequency that is shorter than the first frequency interval.
- the frequency characteristics in the sound field environment can be obtained by calculating while shifting by a shorter frequency. Average value calculating function and may be one to realize.
- a method using an M-sequence code is known as a measurement signal for measuring frequency characteristics.
- the M-sequence code is collected with a microphone as a measurement signal and subjected to Fourier transform processing, the sample length of the Fourier transform may become a non-integer multiple with respect to the code length of the M-sequence code, resulting in an asynchronous relationship. is there.
- a line spectrum with a small level is generated in the frequency characteristic (line spectrum) subjected to Fourier transform, and the measurement accuracy of the frequency characteristic may be deteriorated.
- the M-sequence code is used as the first measurement signal
- the M-sequence code of the first measurement signal is used as the second measurement signal.
- the orthogonal M-sequence code it is possible to realize a stereo measurement signal in the middle and high range using the M-sequence code, and a monaural measurement signal in the low range.
- the maximum value of the signal level in the predetermined first frequency interval is detected while being shifted by a frequency shorter than the first frequency interval.
- the frequency characteristic composed of the maximum value can be obtained.
- the average value of the signal level in the predetermined second frequency interval is calculated based on the frequency characteristic consisting of the maximum value.
- the frequency characteristic consisting of the maximum value is averaged to prevent the frequency characteristic from changing every Fourier transform. Therefore, it is possible to improve the detection accuracy of the frequency characteristics.
- the cutoff frequency set in the low-pass filter unit and the high-pass filter unit outputs the first measurement signal before the low-frequency component is extracted from the pair of speakers simultaneously.
- the frequency may be set lower than the frequency value of the dip that can occur in the frequency characteristics obtained by the Fourier transform unit.
- the cutoff frequency set in the low-frequency component extraction step and the middle-high frequency component extraction step is the first measurement signal before the low-frequency component is extracted
- the frequency may be set lower than the frequency value of the dip that can be generated in the frequency characteristic obtained in the Fourier transform step when simultaneously output from the speaker.
- the cutoff frequency set in the low-pass filter function and the high-pass filter function outputs the first measurement signal before the low-frequency component is extracted from the pair of speakers simultaneously.
- the frequency may be set lower than the frequency value of the dip that can occur in the frequency characteristics obtained by the Fourier transform function.
- the dip occurs in the frequency characteristics measured using the monaural measurement signal as the cutoff frequency of the low-pass filter unit and the high-pass filter unit. Set a frequency lower than the frequency value to be obtained. For this reason, it becomes possible to set the measurement signal in the middle and high range where dip can occur as a stereo measurement signal, and it is possible to improve the measurement accuracy of the frequency characteristics.
- the sound field measurement device, the sound field measurement method, and the sound field measurement program according to the present invention it is possible to measure the middle and high frequencies of frequency characteristics using a stereo measurement signal, thereby suppressing the occurrence of dip. Can do.
- a low frequency characteristic using a monaural measurement signal since it is possible to measure a low frequency characteristic using a monaural measurement signal, intersymbol interference of the measurement signal can be suppressed.
- the first measurement signal and the synthesized signal are simultaneously output from the pair of speakers and the frequency characteristics are measured. Compared with the case where measurement signals are alternately output and measurement is performed, the measurement burden can be reduced and the measurement speed can be improved.
- the upper part shows the frequency characteristic of the signal S1
- (a) the lower part shows the frequency characteristic of the signal S2
- (b) the upper part shows the frequency characteristic of the low-frequency extracted signal S1
- (b) the lower part Shows the frequency characteristic of the signal S2 extracted in the middle and high range
- (c) is a diagram showing the frequency characteristic of the synthesized signal obtained by synthesizing the signal S1 and the signal S2 shown in the upper and lower stages of (b). is there. It is a figure explaining the process which calculates
- FIG. 1 is a block diagram showing an example of a schematic hardware configuration of a sound field measuring apparatus according to the present invention.
- the sound field measuring apparatus 1 includes a CPU 2, a ROM (Read Only Memory) 3, a RAM (Random Access Memory) 4, a storage unit 5, an external output unit 6, a microphone 7, And a display unit 8.
- ROM Read Only Memory
- RAM Random Access Memory
- the ROM 3 stores a processing program executed by the sound field measuring apparatus 1 and the like. For example, when the sound field measuring apparatus 1 is activated or the CPU 2 reads a processing program in the ROM 3 in response to a user operation, various processes such as measurement of frequency characteristics can be performed.
- the RAM 4 is used as a work area for processing performed in the CPU 2.
- the storage unit 5 is a so-called auxiliary storage device, and is mainly composed of a hard disk, an SSD (Solid State Drive), a nonvolatile memory (flash ROM, flash memory, etc.), or the like.
- a removable memory card such as an SD card can be used as the storage unit 5.
- Various data and the like used for various processes in the CPU 2 are recorded in the storage unit 5.
- an information portable terminal such as a smartphone is used as the sound field measuring device 1
- an application program downloaded or the like can be recorded in the storage unit 5. Based on this application program, frequency characteristic measurement processing can be performed.
- the external output unit 6 includes an external output terminal for outputting a measurement signal (signal S1 and synthesized signal) described later to an external input terminal of the portable audio device.
- the measurement signal output via the external output unit 6 is input to the external input terminal of the portable audio device (audio device) 102 (see FIG. 2), whereby the left and right speakers 101a and 101b ( The measurement sound can be output from FIG.
- the external output unit 6 is not limited to a physical terminal structure. For example, a configuration may be adopted in which a measurement signal is input to the portable audio device 102 by using a wireless technology such as Bluetooth (registered trademark) or a wireless LAN.
- the microphone 7 has a role of collecting measurement sound output from the portable audio device 102 or the like.
- the measurement sound collected by the microphone 7 is recorded in the RAM 4 or the storage unit 5 and used for frequency characteristic measurement processing described later.
- the display unit 8 corresponds to a general liquid crystal display, a CRT display (CRT display), or the like.
- the display unit 8 has a role of displaying the sound field frequency characteristics (for example, frequency characteristics shown in FIGS. 7 to 12 described later) that are obtained by the frequency characteristic measurement process so as to be visible to the user.
- FIG. 2 is a block diagram illustrating a schematic configuration of a functional unit of the sound field measuring apparatus 1 when the CPU 2 performs frequency characteristic measurement processing based on a processing program or an application program.
- FIG. 3 is a flowchart showing the processing contents of the CPU 2 based on the processing program and the like.
- the sound field measuring device 1 includes a measurement signal generation unit 11, a frequency division synthesis unit 12, a Fourier transform unit 13, an averaging processing unit 14, an external output unit 6, a microphone 7, And a display unit 8. 2 shows the portable audio device 102 shown in FIG. Since the external output unit 6, the microphone 7, and the display unit 8 have already been described with reference to FIG. 1, description thereof is omitted here.
- the measurement signal generation unit 11 includes a first measurement signal generation unit 21a and a second measurement signal generation unit 21b.
- the measurement signal generation unit 11 generates M-sequence codes having different generation polynomials as measurement signals.
- the M sequence code generated by the first measurement signal generation unit 21a and the M sequence code generated by the second measurement signal generation unit 21b are orthogonal M sequence codes.
- the M sequence code generated by the first measurement signal generation unit 21a is referred to as a signal S1
- the M sequence code generated by the second measurement signal generation unit 21b is referred to as a signal S2.
- the M sequence is a pseudo-random number sequence.
- the M-sequence code is generated by performing feedback using a shift register having a predetermined length and an exclusive OR.
- the code period (code length) is 2 n ⁇ 1
- the feedback position of the shift register is obtained by a generator polynomial. Since the output of the M-sequence code is a binary value of 0 and 1 and contains a large amount of DC component, 0 is converted to -1 and output.
- the CPU 2 performs a process of generating the signal S1 and the signal S2 in the measurement signal generator 11 according to a processing program or the like (S1 in FIG. 3).
- the generated signals S1 and S2 are output to the frequency division synthesizer 12.
- the frequency division synthesis unit 12 includes an LPF unit (low-pass filter unit) 22a, an HPF unit (high-pass filter unit) 22b, an addition unit (synthetic signal generation unit) 22c, and a delay unit 22d.
- LPF unit low-pass filter unit
- HPF unit high-pass filter unit
- addition unit synthetic signal generation unit
- delay unit 22d a delay unit 22d.
- the LPF unit 22a is a low-pass filter that allows low-frequency signal passage.
- the HPF unit 22b is a high-pass filter that allows passage of signals in the mid-high range.
- the cut-off frequency of the LPF unit 22a and the cut-off frequency of the HPF unit 22b are set to the same value. Details of the set cutoff frequency will be described later.
- FIG. 4A shows the frequency characteristic of the signal S1, and the lower part shows the frequency characteristic of the signal S2.
- the signal S1 shown in the upper part of FIG. 4A is input to the LPF unit 22a and low-pass filtering is performed, a low-frequency component is extracted as shown in the upper part of FIG. 3 S2, low-frequency component extraction step, low-pass filter function).
- the signal S2 shown in the lower part of FIG. 4A is input to the HPF unit 22b and subjected to high-pass filter processing, as shown in the lower part of FIG. (S3 in FIG. 3, middle and high frequency component extraction step, high-pass filter function).
- the adding unit 22c has a role of combining the signal S1 filtered by the LPF unit 22a and the signal S2 filtered by the HPF unit 22b.
- a combined signal is generated in which the low frequency side is the component of the signal S1 and the middle high frequency side is the component of the signal S2 (S4 in FIG. 3, a combined signal generating step, a combined signal generating function).
- the delay unit 22d has a role of delaying the input signal S1. Specifically, the timing at which the signal S1 is output from the delay unit 22d is delayed according to the processing time of the filtering process and the addition process in the LPF unit 22a, the HPF unit 22b, and the addition unit 22c described above. By this delay processing, timing adjustment of the output of the signal S1 output from the delay unit 22d and the combined signal output from the adder 22c is performed (S5 in FIG. 3).
- the CPU 2 generates a combined signal in which the low frequency side is the component of the signal S1 and the middle and high frequency side is the component of the signal S2 based on the signal S1 and the signal S2, and the signal S1. (S2-S5 in FIG. 3, frequency division synthesis process). Further, the CPU 2 outputs the generated composite signal as an L channel measurement signal to the portable audio device 102 via the external output unit 6, and the delayed signal S1 is used as an R channel measurement signal. As shown in FIG. 3, the portable audio device 102 outputs the data to the portable audio device 102 via the external output unit 6 (S6 in FIG. 3, external output step, external output function).
- the portable audio device 102 outputs the input measurement signals for L channel and R channel (synthesized signal and signal S1) from the right speaker 101a and the left speaker 101b. In this case, the portable audio device 102 outputs an L channel measurement signal and an R channel measurement signal simultaneously. In this way, by simultaneously outputting the measurement signal for the L channel and the measurement signal for the R channel, the signal component on the low frequency side becomes the signal S1 that is the same M-sequence code in both the L channel and the R channel, Monaural measurement signal.
- the signal component on the middle and high frequency side is an L-channel signal S2 and the R-channel signal S1, and thus becomes an orthogonal M-sequence code signal and a stereo measurement signal.
- the CPU 2 collects the measurement sound output from the left and right speakers 101a and 101b of the portable audio device 102 using the microphone 7 (S7 in FIG. 3, sound collection step, sound collection function).
- the collected measurement sound signal (measurement signal) is output to the Fourier transform unit 13.
- the Fourier transform unit 13 has a role of performing Fourier transform processing (high-speed FFT processing) on the collected measurement signal.
- the CPU 2 performs a Fourier transform process after weighting the collected measurement signal with a window function in the Fourier transform unit 13.
- the CPU 2 transforms the measurement signal in the time domain into the frequency domain, and outputs a line spectrum for each Fourier transform (S8 in FIG. 3, Fourier transform step, Fourier transform function).
- the line spectrum is a power spectrum
- the number of line spectra is half the Fourier transform sample length.
- the measurement signal subjected to the Fourier transform processing is output to the averaging processing unit 14.
- the averaging processing unit 14 has a role of detecting a maximum value for each predetermined number of samples of the measurement signal subjected to the Fourier transform process, and further calculating an average value of the detected maximum values.
- the averaging processor 14 includes a maximum value detector 23a and an average value calculator 23b.
- the maximum value detection unit 23a shifts the line spectrum of the measurement signal subjected to the Fourier transform from the low frequency side to the high frequency side by one sample at a time, and calculates the maximum value within a predetermined number of samples. It has the role of detecting and obtaining the envelope of the line spectrum (S9 in FIG. 3, maximum value detection step, maximum value detection function).
- the average value calculation unit 23b similarly shifts the maximum value line spectrum obtained while shifting the sample by one sample by the maximum value detection unit 23a from the low frequency side to the high frequency side by shifting one sample at a time. It has a role of calculating an average value within the number of samples (S10 in FIG. 3, average value calculation step, average value calculation function).
- FIG. 6 shows the maximum number detection sample number width (maximum value detection width) and average value calculation sample number width (average width), which are set according to the frequency sample when the line spectrum is shifted by one sample at a time.
- the sample length of Fourier transform is 65,536, and the number of line spectra is 32,768.
- the frequency samples on the horizontal axis in FIG. 6 are numbers corresponding to the number of line spectra.
- the predetermined number of samples at the time of maximum value detection and the predetermined number of samples at the time of calculating the average value change according to the frequency sample, and the maximum value detection width and the average from the low frequency side to the high frequency side. It is set so that the conversion width increases.
- the CPU 2 performs a maximum value detection process of 1/18 octave by setting the maximum value detection width as shown in FIG. 6 (S9 in FIG. 3). Further, the CPU 2 performs an average value calculation process of 1/6 octave by setting the averaging width as shown in FIG. 6 (S10 in FIG. 3). It is known that auditory resolution is about 1/3 octave. For this reason, the setting of the maximum value detection width and the setting of the average width as shown in FIG. 6 are performed by the averaging processing unit 14 (the maximum value detection unit 23a and the average value calculation unit 23b), so that the resolution is sufficiently large. An averaging process (maximum value detection process and average value calculation process) can be performed.
- the code length of the M-sequence code is 2 n ⁇ 1 as described above.
- the sample length when the recorded measurement signal is Fourier-transformed generally does not have an integer multiple of the code length of the M-sequence code, and therefore tends to be asynchronous.
- the frequency characteristic is obtained by performing Fourier transform processing on the recorded measurement signal (M-sequence code) using the M-sequence code as the measurement signal, the level is obtained for each Fourier transform.
- a small fluctuating line spectrum occurs between the uniform line spectra. There is a possibility that the line spectrum having a small level changes as noise in the detection of the frequency characteristic.
- the maximum value detecting unit 23a extracts the line spectrum having the largest level within the predetermined number of samples and obtains the frequency characteristics, thereby suppressing the line spectrum having a small level from being detected as noise. it can. Furthermore, by performing average value calculation processing together in the average value calculation unit 23b, it is possible to prevent the frequency characteristics from changing for each Fourier transform.
- the maximum value detection is performed by the maximum value detection unit 23a of the averaging processing unit 14, and the frequency characteristic calculated by the average value calculation unit 23b is output to the display unit 8 in accordance with an instruction from the CPU 2.
- the CPU 2 can output the frequency characteristic obtained by the Fourier transform unit 13 to the display unit 8 as it is without performing the averaging process in the averaging processing unit 14.
- the frequency characteristic received in the display unit 8 is displayed on the display screen or the like of the display unit 8 in a state visible to the user based on an instruction from the CPU 2 (S11 in FIG. 3).
- FIGS. 10 to 12 show examples of frequency characteristics displayed on the display unit 8.
- the frequency characteristic obtained from the collected measurement signal is averaged in the averaging processing unit 14 (maximum value detection processing by the maximum value detection unit 23a and average value calculation unit 23b. An average value calculation process) is not shown.
- (b) of each figure shows the case where the frequency characteristic calculated
- the measurement conditions of the frequency characteristics shown in FIGS. 7 to 9 and FIGS. 10 to 12 use an M-sequence code as a measurement signal, a measurement signal sampling rate of 44.1 kHz, an M-sequence code length of 32,767, and frequency division synthesis.
- the low-pass and middle-high pass filters in the LPF unit 22a and the HPF unit 22b of the unit 12 are 512-tap FIR filters, and the cutoff frequency is 1 kHz.
- the sample length of the Fourier transform by the Fourier transform unit 13 is 65,536.
- the window function used in the Fourier transform unit 13 is a Hamming window
- the maximum value detection width in the maximum value detection unit 23a is 1/18 octave
- the average width in the average value calculation unit 23b is 1/6 octave.
- FIG. 7 and 10 show a case where the L channel and R channel measurement signals input to the portable audio device 102 are the same M-sequence code and are input as monaural measurement signals.
- 8 and 11 show a case where the L channel and R channel measurement signals input to the portable audio device 102 are orthogonal M-sequence codes and are input as stereo measurement signals.
- 9 and 12 as described above, the R channel measurement signal input to the portable audio device 102 is the signal S1 (see the upper part of FIG. 4A), and the L channel measurement signal is the signal S1.
- the signal S2 are combined signals (see FIG. 4C) combined with 1 kHz as a boundary, and the case where they are input as combined measurement signals is shown.
- FIGS. 7 to 9 show the frequency characteristics by loopback in order to show the influence of the inter-channel interference that occurs in the monaural measurement signal and the influence of the inter-symbol interference that occurs in the stereo measurement signal due to a slight delay time difference.
- the case of measurement is shown.
- the measurement by the loopback means that “the measurement signal output from the external output unit 6 is output from the right speaker 101 a and the left speaker 101 b of the portable audio device 102, collected by the microphone 7, and output to the Fourier transform unit 13. It means not a measurement method but a method of measuring frequency characteristics by “outputting the measurement signal output from the external output unit 6 to the Fourier transform unit 13 as it is”.
- FIG. 10 to FIG. 12 show that “the measurement signal output from the external output unit 6 is output from the right speaker 101a and the left speaker 101b of the portable audio device 102, collected by the microphone 7, and Fourier-transformed unit 13”.
- the frequency characteristic is measured by Accordingly, in FIGS. 10 to 12, the frequency characteristics of the sound field at the installation position of the microphone 7 are measured. Also in the case of FIGS. 10 to 12, the measurement signal is collected by the microphone 7 in a state where the R channel measurement signal is delayed by 0.00009 sec (3.1 cm) with respect to the L channel measurement signal. The environment is realized and frequency characteristics are measured.
- the position of the dip varies greatly depending on the delay time difference between the L channel measurement signal and the R channel measurement signal. For this reason, the delay time difference changes according to the measurement position of the measurement signal by the microphone 7, the frequency characteristics also change, and the measurement accuracy tends to deteriorate greatly.
- the maximum value of the line spectrum when the averaging process is not performed (see FIG. 8A) varies due to intersymbol interference, and the low frequency range It is prominent in. Also in the frequency characteristics when the averaging process is performed (see FIG. 8B), the signal level greatly fluctuates in a low region where the number of line spectra is small. As described above, when the stereo measurement signal is used, since the intersymbol interference is added to the frequency characteristic of the sound field environment, the measurement accuracy of the frequency characteristic in the low frequency tends to deteriorate.
- the measurement signal in the low band is a monaural measurement signal
- the measurement signal in the middle and high band is a stereo measurement signal.
- the maximum value of the line spectrum indicates a uniform signal level in the entire band. Yes. Further, even when the averaging process is performed (see FIG. 9B), the signal levels of all the bands are uniform.
- a monaural measurement signal, a stereo measurement signal, and a synthesized measurement signal are actually output from the left and right speakers 101a and 101b of the portable audio device 102 to generate sound.
- the same characteristics as the frequency characteristics shown in FIGS. 7 to 9 are exhibited.
- FIG. 10 as in FIG. 7, dip occurs in the middle and high range due to inter-channel interference, and in FIG. 11, the signal level varies greatly in the low range due to inter-symbol interference as in FIG.
- FIG. 12 shows a state in which the dip is suppressed in the mid-high range and the fluctuation of the signal level in the low range is suppressed as in FIG.
- the frequency characteristics of the sound field environment including the portable audio device 102 and the listening position where the measurement signal is collected by the microphone 7 are obtained.
- the sound field correction is performed by adjusting the equalizer of the portable audio device 102 based on the obtained frequency characteristic.
- the music output from the left and right speakers 101a and 101b of the portable audio device 102 is corrected in advance, and the corrected music is output from the left and right speakers 101a and 101b of the portable audio device 102.
- the code length of the M-sequence code is 2 n ⁇ 1.
- the sample length when the measurement signal is subjected to Fourier transform generally does not have an integer multiple of the code length of the M-sequence code, and thus tends to be asynchronous. Therefore, when the recorded measurement signal (M-sequence code) is subjected to Fourier transform processing by the Fourier transform unit 13 to obtain frequency characteristics, a line spectrum with a small level may be generated between uniform line spectra. There is a risk of noise in the detection of frequency characteristics.
- the maximum value detection unit 23a extracts a line spectrum of the largest level within a predetermined number of samples and obtains a frequency characteristic, thereby causing a small level fluctuation. By suppressing the detection of the line spectrum as noise and performing the average value calculation process together, it is possible to prevent the frequency characteristic from changing for each Fourier transform.
- the frequency interval of the dip in the inter-channel interference changes according to the propagation delay difference as shown in FIG.
- the propagation delay difference between the channels is about several centimeters, and the frequency interval of the dip is around 10 kHz.
- the relationship between the installation distance of the left and right speakers relative to the listening position is asymmetric because of the positional relationship between the left and right speaker installation positions and the seating position.
- the propagation delay difference between channels tends to be several tens of centimeters.
- the frequency interval at which dip occurs may be 1 kHz.
- the cut-off frequency in the LPF unit 22a and the HPF unit 22b is lower than the mid-high frequency where dip can occur due to inter-channel interference, and lower than the low frequency where inter-symbol interference may occur. It is desirable to set a high frequency. By setting the cut-off frequency in this way, it is possible to simultaneously suppress inter-channel interference and reduce inter-symbol interference, and to accurately measure the frequency characteristics of the sound field environment.
- the sound field measurement device, the sound field measurement method, and the sound field measurement program according to the present invention have been described in detail with reference to the drawings.
- the sound field measurement device, the sound field measurement method, and the sound field according to the present invention have been described above.
- the measurement program is not limited to the example described in the embodiment. It will be apparent to those skilled in the art that various changes and modifications can be made within the scope of the claims, and these are also within the technical scope of the present invention. .
- the CPU 2 functions as shown in FIG. 2 based on the processing program or application program recorded in the ROM 3 or the storage unit 5 as shown in FIG.
- the number of CPUs that realize the function of each functional unit is not limited to one.
- a configuration may be adopted in which dedicated processing units (for example, a plurality of CPUs or the like) for realizing a part of the functions of each functional unit are installed, and each dedicated processing unit realizes at least one or more functions.
- dedicated processing units for example, a plurality of CPUs or the like
- each dedicated processing unit realizes at least one or more functions.
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Health & Medical Sciences (AREA)
- Signal Processing (AREA)
- Quality & Reliability (AREA)
- Audiology, Speech & Language Pathology (AREA)
- Computational Linguistics (AREA)
- Data Mining & Analysis (AREA)
- Human Computer Interaction (AREA)
- Multimedia (AREA)
- General Health & Medical Sciences (AREA)
- Otolaryngology (AREA)
- Stereophonic System (AREA)
- Circuit For Audible Band Transducer (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
一対のスピーカ(101a),(101b)が狭間隔に設置されて出力された出力音を集音して周波数特性を求める音場測定装置(1)であって、第1測定信号の低域成分を抽出するローパスフィルタ部(22a)と、第1測定信号とは異なる第2測定信号の中高域成分を抽出するハイパスフィルタ部(22b)と、第1測定信号の低域成分と第2測定信号の中高域成分とを合成して合成信号を生成する合成信号生成部(22c)と、オーディオ機器(102)に対して第1測定信号と合成信号とを出力する外部出力部(6)と、一対のスピーカから同時に出力された第1測定信号と合成信号とを集音するマイクロフォン(7)と、集音された信号をフーリエ変換して、音場環境の周波数特性を求めるフーリエ変換部(13)とを有する。
Description
本発明は音場測定装置、音場測定方法および音場測定プログラムに関し、より詳細には、一対のスピーカが狭間隔に設置される音場環境において、周波数特性を迅速かつ精度良く測定することが可能な音場測定装置、音場測定方法および音場測定プログラムに関する。
従来より、スピーカ等が設置される音場環境において周波数特性を測定し、測定された周波数特性に基づいてオーディオ機器のイコライザを調節したり、出力音を音場に合わせて予め補正したりすることにより、音場環境に最適な音質の音楽を提供する方法が知られている。
周波数特性を測定するための測定信号として、M系列符号(Maximum Length Sequence)やTSP信号(Time Stretched Pulse)が知られている。これらの測定信号を用いて周波数特性を測定する場合には、例えば、スピーカより出力された測定信号を、聴取位置に設置されるマイクロフォンで収録したのちに、収録された信号をフーリエ変換して、音場環境の周波数特性を求める方法を用いる(例えば、特許文献1、特許文献2参照)。また、出力される測定信号をリファレンスとして用い、マイクロフォンで収録された測定信号との相互相関特性を求めることによって、インパルス応答を求めることも可能である。
音場環境の周波数特性は、スピーカに対する聴取位置によって変化する。さらに、周波数特性の測定精度は、設置されるスピーカの間隔等によって異なる傾向がある。例えば、図14に示すように、音楽再生機能部100と、その左右に設置される一対の右スピーカ101aおよび左スピーカ101bとが一体化されたポータブルオーディオ機器102の場合には、一般に左右スピーカ101a,101bの設置間隔が狭くなる傾向がある。さらに、ポータブルオーディオ機器102では、その携帯性等から聴取者の近くに機器本体が置かれ、ポータブルオーディオ機器102と聴取者との距離が近くなる傾向がある。
このような場合において、例えば、同一のM系列符号を測定信号として用いて、図14に示すように、右スピーカ101aと左スピーカ101bとから同時に測定信号を出力する。そして、左右のスピーカ101a,101bとマイクロフォン7との間の周波数特性を測定する。すると、左右スピーカ101a,101bの設置間隔が狭く、マイクロフォン7との距離が近いため、チャンネル(一方のスピーカからマイクロフォン7までの信号の伝達経路)間の遅延時間差が短くなる傾向がある。このように、チャンネル間の遅延時間差が短くなると、チャンネル間干渉が生じ、特定の周波数に大きなディップが発生する傾向が生じる。このディップの位置は遅延時間差に応じて大きく変動するため、測定位置によって周波数特性が変化してしまい、測定精度が大きく劣化するおそれがあるという問題があった。
一方で、直交したM系列符号を測定信号として用いて、右スピーカ101aと左スピーカ101bとから同時に測定信号を出力して周波数特性を測定する場合は、同一周波数の線スペクトルにおいて位相が異なるため干渉が生じてしまう。その結果、合成された各々の線スペクトルに振幅変動が発生するおそれがあった。このため、音場環境の周波数特性に対して符号間の干渉が加わってしまい、測定精度が劣化するという問題があった。また、測定時間やフーリエ変換処理におけるにメモリ使用量を低減させるために、周期の短い符号を測定信号として用いる場合には、線スペクトルの間隔が大きくなる。このため、線スペクトルの平均化処理を行っても符号間干渉の影響は軽減できず、測定精度が大きく低下してしまう傾向があった。
さらに、チャンネル間の干渉や、線スペクトルにおける符号間の干渉を回避するために、まず左スピーカ101bのみから測定信号を出力して周波数特性を測定した後に、右スピーカ101aのみから測定信号を出力して周波数特性を測定することにより、LチャンネルとRチャンネルとでそれぞれ時間分割を行って周波数特性を測定する方法が知られている。このように時間分割を行って周波数特性を測定することにより、チャンネル間の干渉や、符号間の干渉を避けることはできる。しかしながら、LチャンネルとRチャンネルとで2回も測定処理を行う必要があるため、測定時間が増加すると共に、測定の処理負担が増大するという問題があった。
本発明は上記問題に鑑みてなされたものであり、狭間隔に設置される一対のスピーカから同時に測定信号を出力することにより、精度良く音場環境の周波数特性を測定することが可能な音場測定装置、音場測定方法および音場測定プログラムを提供することを課題とする。
上記課題を解決するために、本発明に係る音場測定装置は、一対のスピーカが狭間隔に設置されたオーディオ機器から出力される出力音を集音することにより、音場環境の周波数特性を求める音場測定装置であって、第1測定信号の低域成分を抽出するローパスフィルタ部と、前記第1測定信号とは異なる第2測定信号の中高域成分を抽出するハイパスフィルタ部と、前記ローパスフィルタ部により抽出された前記第1測定信号の低域成分と、前記ハイパスフィルタ部により抽出された前記第2測定信号の中高域成分とを合成して合成信号を生成する合成信号生成部と、前記ローパスフィルタ部により低域成分が抽出される前の前記第1測定信号を、前記一対のスピーカの何れか一方より出力させると共に、前記合成信号生成部により生成された前記合成信号を前記一対のスピーカの何れか他方より同時に出力させるために、前記オーディオ機器に対して前記第1測定信号と前記合成信号とを出力する外部出力部と、前記一対のスピーカから同時に出力された前記第1測定信号と前記合成信号とを集音するためのマイクロフォンと、該マイクロフォンにより集音された信号をフーリエ変換することにより、前記音場環境における周波数特性を求めるフーリエ変換部とを有することを特徴とする。
また、本発明に係る音場測定方法は、一対のスピーカが狭間隔に設置されたオーディオ機器から出力される出力音を集音することにより、音場環境の周波数特性を求める音場測定装置の音場測定方法であって、第1測定信号の低域成分をローパスフィルタ部が抽出する低域成分抽出ステップと、前記第1測定信号とは異なる第2測定信号の中高域成分をハイパスフィルタ部が抽出する中高域成分抽出ステップと、前記低域成分抽出ステップにおいて抽出された前記第1測定信号の低域成分と、前記中高域成分抽出ステップにおいて抽出された前記第2測定信号の中高域成分とを合成して、合成信号生成部が合成信号を生成する合成信号生成ステップと、前記低域成分抽出ステップにおいて低域成分が抽出される前の前記第1測定信号を、前記一対のスピーカの何れか一方より出力させると共に、前記合成信号生成ステップにおいて生成された前記合成信号を前記一対のスピーカの何れか他方より同時に出力させるために、外部出力部が、前記オーディオ機器に対して前記第1測定信号と前記合成信号とを出力する外部出力ステップと、前記一対のスピーカから同時に出力された前記第1測定信号と前記合成信号とを、マイクロフォンが集音する集音ステップと、該集音ステップにおいて集音された信号をフーリエ変換することにより、フーリエ変換部が、前記音場環境における周波数特性を求めるフーリエ変換ステップとを有することを特徴とする。
さらに、本発明に係る音場測定プログラムは、一対のスピーカが狭間隔に設置されたオーディオ機器から出力される出力音を集音することにより、音場環境の周波数特性を求める音場測定装置の音場測定プログラムであって、前記音場測定装置のコンピュータに、第1測定信号の低域成分を抽出させるローパスフィルタ機能と、前記第1測定信号とは異なる第2測定信号の中高域成分を抽出させるハイパスフィルタ機能と、前記ローパスフィルタ機能により抽出された前記第1測定信号の低域成分と、前記ハイパスフィルタ機能により抽出された前記第2測定信号の中高域成分とを合成して合成信号を生成させる合成信号生成機能と、前記ローパスフィルタ機能により低域成分が抽出される前の前記第1測定信号を、前記一対のスピーカの何れか一方より出力させると共に、前記合成信号生成機能により生成された前記合成信号を前記一対のスピーカの何れか他方より同時に出力させるために、前記オーディオ機器に対して前記第1測定信号と前記合成信号とを出力させる外部出力機能と、前記一対のスピーカから同時に出力された前記第1測定信号と前記合成信号とをマイクロフォンを用いて集音させる集音機能と、該集音機能により集音された信号をフーリエ変換することにより、前記音場環境における周波数特性を求めさせるフーリエ変換機能とを実現させるための音場測定プログラムであることを特徴とする。
一対のスピーカの設置間隔が狭いオーディオ機器に対して正面位置にマイクロフォンが設置される場合には、一対のスピーカからマイクロフォンまでの距離が短くなる傾向がある。このような状況において一対のスピーカから同じ測定信号(モノラルの測定信号)を出力して周波数特性を測定すると、チャンネル間の出力信号の伝搬遅延差によって、中高域にディップが発生するおそれがある。したがって、精度良く音場環境の周波数特性を測定することが容易ではなかった。
一方で、一対のスピーカより異なる測定信号(ステレオの測定信号)を出力する場合には、中高域においてディップが発生しにくくなる。しかしながら、低域において測定信号の符号間干渉が生じるおそれがあり、精度良く音場環境の周波数特性を測定することが容易ではなかった。
本発明に係る音場測定装置、音場測定方法および音場測定プログラムでは、オーディオ機器の何れか一方のスピーカから第1測定信号が出力され、他方のスピーカから、低域が第1測定信号であって中高域が第2測定信号からなる合成信号が出力される。このため、中高域においては、ステレオの測定信号を用いた周波数特性の測定を行うことができ、ディップの発生を抑制することが可能となる。
さらに、本発明に係る音場測定装置、音場測定方法および音場測定プログラムでは、低域において、モノラルの測定信号を用いて周波数特性の測定を行うことができるため、測定信号の符号間干渉を抑制することが可能になる。
このように、本発明に係る音場測定装置、音場測定方法および音場測定プログラムでは、ステレオの測定信号を用いて周波数特性の中高域の測定を行うことができるため、ディップの発生を抑制することが可能になる。また、モノラルの測定信号を用いて周波数特性の低域の測定を行うことができるため、測定信号の符号間干渉を抑制することが可能になる。さらに、本発明に係る音場測定装置、音場測定方法および音場測定プログラムでは、一対のスピーカから第1測定信号と合成信号とを同時に出力して周波数特性を測定するので、一対のスピーカから交互に測定信号を出力して測定する場合に比べて、測定負担の軽減と測定スピードの向上を図ることが可能となる。
また、上述した音場測定装置において、前記第1測定信号としてM系列符号を用いると共に、前記第2測定信号として前記第1測定信号のM系列符号に直交したM系列符号を用い、前記フーリエ変換部により求められた周波数特性に基づいて、所定の第1周波数間隔における信号レベルの最大値を、当該第1周波数間隔よりも短い周波数ずつシフトさせながら検出することにより、最大値から成る周波数特性を求める最大値検出部と、該最大値検出部により検出された最大値から成る周波数特性に基づいて、所定の第2周波数間隔における信号レベルの平均値を、当該第2周波数間隔よりも短い周波数ずつシフトさせながら算出することにより、前記音場環境における周波数特性を求める平均値算出部とを有するものであってもよい。
また、上述した音場測定方法において、前記第1測定信号としてM系列符号を用いると共に、前記第2測定信号として前記第1測定信号のM系列符号に直交したM系列符号を用い、前記フーリエ変換ステップにおいて求められた周波数特性に基づいて、最大値検出部が、所定の第1周波数間隔における信号レベルの最大値を、当該第1周波数間隔よりも短い周波数ずつシフトさせながら検出することにより、最大値から成る周波数特性を求める最大値検出ステップと、該最大値検出ステップにおいて求められた、最大値から成る周波数特性に基づいて、平均値算出部が、所定の第2周波数間隔における信号レベルの平均値を、当該第2周波数間隔よりも短い周波数ずつシフトさせながら算出することにより、前記音場環境における周波数特性を求める平均値算出ステップとを有するものであってもよい。
さらに、上述した音場測定プログラムにおいて、前記第1測定信号としてM系列符号を用いると共に、前記第2測定信号として前記第1測定信号のM系列符号に直交したM系列符号を用い、前記コンピュータに、前記フーリエ変換機能により求められた周波数特性に基づいて、所定の第1周波数間隔における信号レベルの最大値を、当該第1周波数間隔よりも短い周波数ずつシフトさせながら検出させることにより、最大値から成る周波数特性を求めさせる最大値検出機能と、該最大値検出機能により求められた最大値から成る周波数特性に基づいて、所定の第2周波数間隔における信号レベルの平均値を、当該第2周波数間隔よりも短い周波数ずつシフトさせながら算出させることにより、前記音場環境における周波数特性を求めさせる平均値算出機能とを実現させるものであってもよい。
周波数特性を測定するための測定信号として、従来よりM系列符号を用いる方法が知られている。しかしながら、M系列符号を測定信号としてマイクロフォンで集音してフーリエ変換処理を行うと、M系列符号の符号長に対してフーリエ変換のサンプル長が非整数倍となり非同期の関係になってしまうおそれがある。このように非同期となる場合には、フーリエ変換された周波数特性(線スペクトル)において、レベルの小さな変動する線スペクトルが発生し、周波数特性の測定精度を劣化させてしまうおそれがある。
このため、本発明に係る音場測定装置、音場測定方法及び音場測定プログラムでは、まず、第1測定信号としてM系列符号を用い、第2測定信号として第1測定信号のM系列符号に直交したM系列符号を用いることにより、M系列符号を用いて中高域ではステレオの測定信号を実現すると共に、低域ではモノラルの測定信号を実現することができる。そして、集音された信号をフーリエ変換して求めた周波数特性に基づいて、所定の第1周波数間隔における信号レベルの最大値を、当該第1周波数間隔よりも短い周波数ずつシフトさせながら検出することにより、最大値から成る周波数特性を求めることができる。このようにして最大値からなる周波数特性を求めることにより、レベルの小さな変動する線スペクトルが、ノイズとして周波数特性に現れてしまうことを抑制する(マスク化する)ことができる。
さらに、本発明に係る音場測定装置、音場測定方法および音場測定プログラムでは、最大値から成る周波数特性に基づいて、所定の第2周波数間隔における信号レベルの平均値を、第2周波数間隔よりも短い周波数ずつシフトさせながら算出することにより、平均化処理された周波数特性を求めることができる。このように、最大値を検出することにより周波数特性を求めるだけでなく、さらに、最大値から成る周波数特性の平均化処理を行うことによって、フーリエ変換毎に周波数特性が変化してしまうことを抑制することができ、周波数特性の検出精度の向上を図ることが可能となる。
また、上述した音場測定装置において、前記ローパスフィルタ部および前記ハイパスフィルタ部において設定されるカットオフ周波数は、低域成分が抽出される前の前記第1測定信号を前記一対のスピーカより同時に出力した場合に、前記フーリエ変換部により求められる周波数特性に発生し得るディップの周波数値よりも低い周波数に設定されるものであってもよい。
さらに、上述した音場測定方法において、前記低域成分抽出ステップおよび前記中高域成分抽出ステップにおいて設定されるカットオフ周波数は、低域成分が抽出される前の前記第1測定信号を前記一対のスピーカより同時に出力した場合に、前記フーリエ変換ステップにおいて求められる周波数特性に発生し得るディップの周波数値よりも低い周波数に設定されるものであってもよい。
また、上述した音場測定プログラムにおいて、前記ローパスフィルタ機能および前記ハイパスフィルタ機能において設定されるカットオフ周波数は、低域成分が抽出される前の前記第1測定信号を前記一対のスピーカより同時に出力した場合に、前記フーリエ変換機能により求められる周波数特性に発生し得るディップの周波数値よりも低い周波数に設定されるものであってもよい。
本発明に係る音場測定装置、音場測定方法および音場測定プログラムでは、ローパスフィルタ部およびハイパスフィルタ部のカットオフ周波数を、モノラルの測定信号を用いて測定された周波数特性においてディップが発生し得る周波数値よりも低い周波数に設定する。このため、ディップが発生し得る中高域の測定信号をステレオの測定信号に設定することが可能となり、周波数特性の測定精度を向上させることが可能となる。
本発明に係る音場測定装置、音場測定方法および音場測定プログラムによれば、ステレオの測定信号を用いて周波数特性の中高域の測定を行うことができるため、ディップの発生を抑制することができる。また、モノラルの測定信号を用いて周波数特性の低域の測定を行うことができるため、測定信号の符号間干渉を抑制することが可能となる。さらに、本発明に係る音場測定装置、音場測定方法および音場測定プログラムでは、一対のスピーカから第1測定信号と合成信号とを同時に出力して周波数特性を測定するため、一対のスピーカから交互に測定信号を出力して測定する場合に比べて、測定負担の軽減と測定スピードの向上を図ることが可能となる。
以下、本発明に係る音場測定装置について、図面を用いて詳細に説明する。図1は、本発明に係る音場測定装置のハードウエア的な概略構成の一例を示したブロック図である。音場測定装置1は、図1に示すように、CPU2と、ROM(Read Only Memory)3と、RAM(Random Access Memory)4と、記憶部5と,外部出力部6と、マイクロフォン7と、表示部8とを有している。
ROM3には、音場測定装置1において実行される処理プログラム等が記憶されている。例えば、音場測定装置1の起動時、あるいは、ユーザによる操作に応じてCPU2がROM3の処理プログラム等を読み出すことによって、周波数特性の測定等の各種処理を行うことが可能になっている。RAM4は、CPU2において行われる処理のワークエリア等として用いられる。
記憶部5は、いわゆる補助記憶装置であって、主としてハードディスク、SSD(Solid State Drive)、あるいは不揮発性メモリ(フラッシュROM・フラッシュメモリ等)等により構成されている。また、SDカードのように取り外し可能なメモリーカードを記憶部5として利用することも可能である。記憶部5には、CPU2における各種処理に用いられる様々なデータ等が記録される。
なお、音場測定装置1としてスマートフォンなどの情報携帯端末を利用する場合には、記憶部5にダウンロード等されたアプリケーションプログラムを記録させることも可能である。このアプリケーションプログラムに基づいて、周波数特性の測定処理を行うことが可能になる。
外部出力部6は、後述する測定信号(信号S1および合成信号)を、ポータブルオーディオ機器の外部入力端子に対して出力するための外部出力端子等により構成される。この外部出力部6を介して出力された測定信号が、ポータブルオーディオ機器(オーディオ機器)102(図2参照)の外部入力端子に入力されることにより、ポータブルオーディオ機器102の左右スピーカ101a,101b(図2参照)から測定音を出力することが可能となる。なお、外部出力部6は、物理的な端子構造によるものには限定されない。例えば、Bluetooth(登録商標)や無線LANなどの無線技術を用いることにより測定信号をポータブルオーディオ機器102に入力させる構成であっても良い。
マイクロフォン7は、ポータブルオーディオ機器102等より出力される測定音を集音する役割を有している。マイクロフォン7により集音された測定音は、RAM4あるいは記憶部5に記録され、後述する周波数特性の測定処理に利用される。表示部8は、一般的な液晶ディスプレイやCRTディスプレイ(ブラウン管ディスプレイ)等が該当する。表示部8は、周波数特性の測定処理により求められた音場の周波数特性(例えば、後述する図7~図12に示す周波数特性)を、ユーザに視認可能に表示させる役割を有している。
CPU2は、ROM3に記録される処理プログラム、あるいは、記憶部5に記録される周波数特性の測定用アプリケーションプログラムに従って、ポータブルオーディオ機器102とマイクロフォン7との間の周波数特性の測定処理を実行する役割を有している。図2は、CPU2が処理プログラムあるいはアプリケーションプログラムに基づいて、周波数特性の測定処理を行う場合における音場測定装置1の機能部の概略構成を示したブロック図である。また、図3は、処理プログラム等に基づくCPU2の処理内容を示したフローチャートである。
音場測定装置1は、図2に示すように、測定信号生成部11と、周波数分割合成部12と、フーリエ変換部13と、平均化処理部14と、外部出力部6と、マイクロフォン7と表示部8とを有している。また、図2には、図14に示したポータブルオーディオ機器102が示されている。なお、外部出力部6とマイクロフォン7と表示部8とについては、既に図1において説明したので、ここでの説明は省略する。
測定信号生成部11は、図2に示すように、第1測定信号生成部21aと、第2測定信号生成部21bとを有している。測定信号生成部11では、測定信号として、それぞれ生成多項式の異なるM系列符号を生成する。具体的には、第1測定信号生成部21aで生成されるM系列符号と、第2測定信号生成部21bで生成されるM系列符号とは、それぞれ直交したM系列符号となる。本実施の形態では、便宜上、第1測定信号生成部21aで生成されるM系列符号を信号S1とし、第2測定信号生成部21bで生成されるM系列符号を信号S2とする。
ここで、M系列とは、擬似乱数列である。M系列符号は、所定の長さのシフトレジスタと排他的論理和とを用いてフィードバックすることにより生成される。シフトレジスタの長さがnの場合には、符号の周期(符号長)は2n-1となり、シフトレジスタのフィードバック位置は、生成多項式で求められる。M系列符号の出力は0と1との2値で直流成分が多く含まれた信号となるため、0を-1に変換して出力される。
CPU2は、処理プログラム等に従って、測定信号生成部11において、信号S1と信号S2とを生成する処理を行う(図3のS1)。生成された信号S1および信号S2は、周波数分割合成部12に出力される。
周波数分割合成部12は、図2に示すように、LPF部(ローパスフィルタ部)22aと、HPF部(ハイパスフィルタ部)22bと、加算部(合成信号生成部)22cと、遅延部22dとを有している。測定信号生成部11において生成された信号S1は、LPF部22aおよび遅延部22dに入力され、信号S2は、HPF部22bに入力される。
LPF部22aは、低域の信号通過を許容するローパスフィルタである。HPF部22bは、中高域の信号通過を許容するハイパスフィルタである。ここで、LPF部22aのカットオフ周波数と、HPF部22bのカットオフ周波数とは、同じ値に設定される。設定されるカットオフ周波数の詳細については後述する。
図4(a)の上段は、信号S1の周波数特性を示し、下段は信号S2の周波数特性を示している。図4(a)の上段に示す信号S1がLPF部22aに入力されてローパスフィルタ処理されると、図4(b)の上段に示すように、低域成分が抽出された状態となる(図3のS2、低域成分抽出ステップ、ローパスフィルタ機能)。一方で、図4(a)の下段に示す信号S2がHPF部22bに入力されてハイパスフィルタ処理されると、図4(b)の下段に示すように、中高域成分が抽出された状態となる(図3のS3、中高域成分抽出ステップ、ハイパスフィルタ機能)。
加算部22cは、LPF部22aでフィルタ処理された信号S1と、HPF部22bでフィルタ処理された信号S2とを合成する役割を有している。加算部22cにおいて、LPF部22aより出力された信号S1(図4(b)上段)と、HPF部22bより出力された信号S2(図4(b)下段)とを合成すると、図4(c)に示すように、低域側が信号S1の成分となり、中高域側が信号S2の成分となった合成信号が生成される(図3のS4、合成信号生成ステップ、合成信号生成機能)。
遅延部22dは、入力される信号S1の遅延を行う役割を有している。具体的には、上述したLPF部22a,HPF部22bおよび加算部22cにおけるフィルタ処理および加算処理の処理時間に応じて、信号S1が遅延部22dから出力されるタイミングを遅らせる。この遅延処理により、遅延部22dより出力される信号S1と加算部22cより出力される合成信号との出力のタイミング調整を行っている(図3のS5)。
このように、CPU2は、周波数分割合成部12において、信号S1及び信号S2に基づいて、低域側が信号S1の成分となり、中高域側が信号S2の成分となる合成信号を生成すると共に、信号S1の遅延処理を行う(図3のS2~S5、周波数分割合成処理)。また、CPU2は、生成された合成信号を、Lチャンネル用の測定信号として、外部出力部6を介して、ポータブルオーディオ機器102に出力させ、遅延処理された信号S1を、Rチャンネル用の測定信号として、外部出力部6を介して、ポータブルオーディオ機器102に出力させる(図3のS6、外部出力ステップ、外部出力機能)。
ポータブルオーディオ機器102では、入力されたLチャンネルおよびRチャンネル用の測定信号(合成信号および信号S1)を右スピーカ101aおよび左スピーカ101bから出力する。この場合、ポータブルオーディオ機器102では、Lチャンネル用の測定信号とRチャンネル用の測定信号とを同時に出力する。このように、Lチャンネル用の測定信号とRチャンネル用の測定信号とを同時に出力することにより、低域側の信号成分は、LチャンネルもRチャンネルも同一のM系列符号である信号S1となり、モノラルの測定信号となる。一方で、中高域側の信号成分は、Lチャンネルが信号S2となり、Rチャンネルが信号S1となるため、直交したM系列符号の信号となり、ステレオの測定信号となる。
CPU2は、ポータブルオーディオ機器102の左右スピーカ101a,101bより出力された測定音を、マイクロフォン7を用いて集音する(図3のS7、集音ステップ、集音機能)。集音された測定音の信号(測定信号)は、フーリエ変換部13へ出力される。
フーリエ変換部13は、集音された測定信号に対してフーリエ変換処理(高速FFT処理)を行う役割を有している。CPU2は、フーリエ変換部13において、集音した測定信号に対して窓関数で重みづけを行った後、フーリエ変換処理を行う。このフーリエ変換処理により、CPU2は、時間領域の測定信号を周波数領域に変換し、フーリエ変換毎に線スペクトルを出力する(図3のS8、フーリエ変換ステップ、フーリエ変換機能)。ここで、線スペクトルはパワースペクトルであり、線スペクトル数はフーリエ変換のサンプル長の半分となる。フーリエ変換処理された測定信号は、平均化処理部14へ出力される。
平均化処理部14は、フーリエ変換処理された測定信号の所定サンプル数毎の最大値を検出し、さらに検出された最大値の平均値を算出する役割を有している。平均化処理部14は、図2に示すように最大値検出部23aと、平均値算出部23bとを有している。最大値検出部23aは、図5に示すように、フーリエ変換処理された測定信号の線スペクトルを、低域側から高域側にわたって、1サンプルずつシフトさせながら、所定サンプル数内の最大値を検出して、線スペクトルの包絡線を求める役割を有している(図3のS9、最大値検出ステップ、最大値検出機能)。また、平均値算出部23bも同様に、最大値検出部23aによって1サンプルずつシフトさせながら求められた最大値の線スペクトルを、低域側から高域側にわたって、1サンプルずつシフトさせながら、所定サンプル数内の平均値を算出する役割を有している(図3のS10、平均値算出ステップ、平均値算出機能)。
図6は、線スペクトルを1サンプルずつシフトさせる時の周波数サンプルに応じて設定される、最大値検出のサンプル数幅(最大値検出幅)と、平均値算出のサンプル数幅(平均化幅)とを示した図である。図6において、フーリエ変換のサンプル長は65,536、線スペクトル数は32,768となっている。図6の横軸の周波数サンプルは、線スペクトル数に対応する数字となっている。図6に示すように、周波数サンプルに応じて、最大値検出時の所定サンプル数も平均値算出時の所定サンプル数も変化しており、低域側から高域側にわたって最大値検出幅と平均化幅が増加するようにして設定される。CPU2では、図6のような最大値検出幅の設定により、1/18オクターブの最大値検出処理を行う(図3のS9)。さらに、CPU2では、図6のような平均化幅の設定により、1/6オクターブの平均値算出処理を行う(図3のS10)。聴覚の分解能は1/3オクターブ程度であることが知られている。このため、図6に示すような最大値検出幅の設定および平均化幅の設定を平均化処理部14(最大値検出部23aおよび平均値算出部23b)で行うことにより、十分に大きな分解能で平均化処理(最大値検出処理および平均値算出処理)を行うことが可能となる。
なお、M系列符号の符号長は、上述したように2n-1となる。一方で、収録された測定信号をフーリエ変換する場合のサンプル長は、一般的にM系列符号の符号長の整数倍の関係とはならないため、非同期となってしまう傾向がある。このように非同期となってしまうことから、M系列符号を測定信号として用い、収録された測定信号(M系列符号)をフーリエ変換処理して周波数特性を求める場合には、フーリエ変換毎に、レベルの小さな変動する線スペクトルが、一様な線スペクトルの間に発生する。このレベルの小さな変動する線スペクトルは、周波数特性の検出においてノイズとなってしまうおそれがある。しかしながら、最大値検出部23aにおいて、所定サンプル数内で最も大きなレベルの線スペクトルを抽出して周波数特性を求めることにより、レベルの小さな変動する線スペクトルがノイズとして検出されることを抑制することができる。さらに、平均値算出部23bにおいて平均値算出処理を併せて行うことにより、フーリエ変換毎に周波数特性が変化しないようにすることが可能となる。
平均化処理部14の最大値検出部23aで最大値検出が行われ、平均値算出部23bで平均値算出された周波数特性は、CPU2の指示に従って、表示部8へと出力される。なお、CPU2では、平均化処理部14における平均化処理を行うことなく、フーリエ変換部13により求められた周波数特性をそのまま表示部8へと出力することも可能となっている。表示部8において受信された周波数特性は、CPU2の指示に基づいて、表示部8のディスプレイ画面等において、ユーザに視認可能な状態で表示される(図3のS11)。
図7~図9および図10~12は、表示部8に表示された周波数特性の一例を示している。各図の(a)には、集音された測定信号によって求められた周波数特性を、平均化処理部14において平均化処理(最大値検出部23aによる最大値検出処理および平均値算出部23bによる平均値算出処理)しない場合が示されている。また、各図の(b)には、集音された測定信号によって求められた周波数特性を、平均化処理部14において平均化処理した場合が示されている。
図7~図9および図10~12に示す周波数特性の測定条件は、測定信号としてM系列符号を用い、測定信号のサンプリング速度を44.1kHz、M系列符号長を32,767、周波数分割合成部12のLPF部22aおよびHPF部22bにおける低域および中高域通過フィルタを512タップのFIRフィルタとし、カットオフ周波数を1kHzとする。また、フーリエ変換部13によるフーリエ変換のサンプル長は65,536とする。フーリエ変換部13において使用する窓関数はハミング窓、最大値検出部23aにおける最大値検出幅は1/18オクターブ、平均値算出部23bにおける平均化幅は1/6オクターブとする。
また、図7及び図10は、ポータブルオーディオ機器102に入力されるLチャンネルおよびRチャンネルの測定信号が、同一のM系列符号であって、モノラル測定信号として入力される場合を示している。図8及び図11は、ポータブルオーディオ機器102に入力されるLチャンネルおよびRチャンネルの測定信号が、直交したM系列符号であって、ステレオ測定信号として入力される場合を示している。図9および図12は、既に説明したように、ポータブルオーディオ機器102に入力されるRチャンネルの測定信号が信号S1(図4(a)の上段参照)であり、Lチャンネルの測定信号が信号S1と信号S2とが1kHzを境にして合成された合成信号(図4(c)参照)であって、合成測定信号として入力される場合を示している。
さらに、図7~図9は、微少な遅延時間差によって、モノラル測定信号において生じるチャンネル間干渉の影響や、ステレオ測定信号において生じる符号間干渉の影響を顕著に示すために、ループバックにより周波数特性を測定した場合を示している。ループバックによる測定とは、「外部出力部6より出力された測定信号を、ポータブルオーディオ機器102の右スピーカ101aおよび左スピーカ101bより出力させて、マイクロフォン7で集音してフーリエ変換部13へ出力する」測定方法ではなく、「外部出力部6より出力された測定信号をそのまま、フーリエ変換部13へ出力する」ことにより、周波数特性を測定する方法を意味している。このようなループバック測定を用いることにより、理想的にはフラットな周波数特性を得ることができ、測定処理におけるノイズ等を容易に判断することが可能になる。また、図7~図9では、Lチャンネルの測定信号に対して、Rチャンネルの測定信号が0.00009sec(3.1cm)の伝搬遅延した状態における周波数特性を音場測定装置1で測定した場合を示している。
一方で、図10~図12は、「外部出力部6より出力された測定信号を、ポータブルオーディオ機器102の右スピーカ101aおよび左スピーカ101bより出力させ、マイクロフォン7で集音してフーリエ変換部13へ出力する」ことにより、周波数特性を測定した場合を示している。従って、図10~図12では、マイクロフォン7の設置位置における音場の周波数特性が測定されることになる。図10~図12の場合においても、Lチャンネルの測定信号に対して、Rチャンネルの測定信号が0.00009sec(3.1cm)だけ伝搬遅延した状態において、マイクロフォン7で測定信号が集音される環境を実現させて、周波数特性の測定を行っている。
図7に示すモノラル測定信号の周波数特性では、平均化処理を行った場合(図7(b)参照)であっても、あるいは、平均化処理を行わなかった場合(図7(a)参照)であっても、低域における線スペクトルの最大値が顕著に変動することなく一様な状態(信号レベルが一様に0dBを示している状態)となっている。一方で、中高域では、平均化処理を行った場合であっても、あるいは、平均化処理を行わなかった場合であっても、チャンネル間干渉によって大きなディップが発生しており、ディップの間隔は伝搬遅延時間の逆数である11kHz程度になっている(図7(a)(b)参照)。ディップの位置は、Lチャンネルの測定信号とRチャンネルの測定信号との遅延時間差に応じて大きく変動する。このため、マイクロフォン7による測定信号の測定位置に応じて遅延時間差が変化して周波数特性も変化し、測定精度が大きく劣化するという傾向がある。
一方で、図8に示すステレオ測定信号の周波数特性では、平均化処理を行っていない場合(図8(a)参照)における線スペクトルの最大値が、符号間干渉により変動しており、低域において顕著になっている。また、平均化処理を行った場合(図8(b)参照)における周波数特性においても、線スペクトル数が少ない低域において信号レベルが大きく変動している。このように、ステレオ測定信号を用いた場合には、音場環境の周波数特性に対して符号間の干渉が加わるため、低域における周波数特性の測定精度は劣化する傾向がある。
図9に示す合成信号を用いた周波数特性の場合では、低域における測定信号はモノラルの測定信号となり、中高域における測定信号はステレオの測定信号となる。このような測定信号に基づいて求められた周波数特性では、平均化処理を行わなかったときの線スペクトル(図9(a)参照)の最大値が、全帯域において一様な信号レベルを示している。また、平均化処理を行った場合(図9(b)参照)であっても、全帯域の信号レベルが一様になっている。
図7において説明したように、モノラル測定信号を用いて周波数特性を測定すると、Lチャンネルの測定信号とRチャンネルの測定信号との遅延時間差によって、中高域でディップを発生させるおそれが生じる。このため、図9に示すように、中高域においてステレオ測定信号を用いて周波数特性を測定することにより、中高域におけるディップの発生を抑制することができ、周波数特性の測定精度を向上させることが可能となる。さらに、図8において説明したように、ステレオ測定信号を用いて周波数特性を測定すると、符号間の干渉が発生してしまう傾向があり、特に、線スペクトル数が少ない低域では信号レベルが大きく変動してしまうおそれがある。このため、図9に示すように、低域においてモノラル測定信号を用いて周波数特性を測定することにより、符号間の干渉を低減させることができ、線スペクトル数が少ない低域であっても信号レベルの変動を抑制することが可能となる。従って、図9に示すように、低域における測定信号としてモノラルの測定信号を用い、中高域における測定信号としてステレオの測定信号を用いることにより、中高域におけるチャンネル間干渉の抑制と、低域における符号間干渉の低減とを同時に実現することができ、精度良く音場環境の周波数特性を測定することが可能となる。
図10~図12に示すように、実際にポータブルオーディオ機器102の左右のスピーカ101a,101bから、モノラルの測定信号、ステレオの測定信号および合成された測定信号(合成信号)を出力して、音場の周波数特性を測定した場合であっても、図7~図9に示す周波数特性と同様の特徴を示すことになる。図10では、図7と同様に、チャンネル間干渉により中高域でディップが発生し、図11では、図8と同様に、符号間干渉により低域で信号レベルが大きく変動している。さらに、図12においては、図9と同様に、中高域においてディップが抑制され、低域における信号レベルの変動が抑制された状態が示されている。このため、図12においても、図9の場合と同様に、低域における測定信号としてモノラルの測定信号を用い、中高域における測定信号としてステレオの測定信号を用いることによって、中高域におけるチャンネル間干渉の抑制と、低域における符号間干渉の低減とを同時に実現することができ、精度良く音場環境の周波数特性を測定することが可能となる。
従って、図12に示す平均化処理を行った後の周波数特性に基づいて、ポータブルオーディオ機器102とマイクロフォン7により測定信号の集音が行われた聴取位置とからなる音場環境の周波数特性を求めて、求められた周波数特性に基づいてポータブルオーディオ機器102のイコライザを調節して音場補正を行う。この音場補正により、聴取位置における音質の向上を図ることが可能となる。また、求められた周波数特性に基づいて、ポータブルオーディオ機器102の左右スピーカ101a,101bより出力する音楽に予め補正を施し、ポータブルオーディオ機器102の左右スピーカ101a,101bから補正された音楽を出力することによって、聴取位置において聴取される音質の向上を図ることが可能となる。
また、既に説明したように、測定信号としてM系列符号を用いる場合には、M系列符号の符号長が2n-1となる。一方で、測定信号をフーリエ変換する場合のサンプル長は、一般的にM系列符号の符号長の整数倍の関係とはならないため、非同期となってしまう傾向がある。従って、収録された測定信号(M系列符号)をフーリエ変換部13でフーリエ変換処理して周波数特性を求める場合に、レベルの小さな変動する線スペクトルが、一様な線スペクトルの間に発生するおそれがあり、周波数特性の検出においてノイズとなってしまうおそれがある。このため、本実施の形態に係る音場測定装置1では、最大値検出部23aにおいて、所定サンプル数内で最も大きなレベルの線スペクトルを抽出して周波数特性を求めることにより、レベルの小さな変動する線スペクトルがノイズとして検出されることを抑制し、さらに、平均値算出処理を併せて行うことにより、フーリエ変換毎に周波数特性が変化しないようにすることが可能になる。
なお、チャンネル間干渉におけるディップの周波数間隔は、図13に示すように、伝搬遅延差に応じて変化する。ポータブルオーディオ機器102のように左右のスピーカ101a,101bの間隔が狭く、さらに聴取者に対して比較的近い位置に設置されて、間接音よりも直接音の方がより顕著に聞こえる場合には、チャンネル間の伝搬遅延差は数センチメートル程度となり、ディップの周波数間隔が10kHz前後となる。このため周波数分割合成部12におけるLPF部22aおよびHPF部22bのフィルタカットオフ周波数は、それよりも十分に小さい周波数である1kHz前後に設定することが好ましい。1kHz前後にカットオフ周波数を設定することにより、ディップが発生しやすい10kHz前後の中高域の測定信号を、ディップが発生しにくいモノラルの測定信号にすることが可能となる。これにより、中高域における周波数特性の測定精度を向上させることが可能になる。
また、車室内に設置される車載用オーディオ機器の音場環境を考えると、左右のスピーカ設置位置と着座位置との位置関係から、聴取位置に対する左右のスピーカの設置距離の関係は非対称な関係となる。チャンネル間の伝搬遅延差は、数10センチメートルになる傾向がある。このように、伝搬遅延差が数10センチメートルになる場合には、ディップが発生する周波数間隔が1kHzとなる場合がある。このように伝搬遅延差が大きくなる場合には、伝搬遅延差により生じるディップの発生周波数に応じて、フィルタカットオフ周波数をより小さい周波数に設定することが望ましい。
さらに、カットオフ周波数をあまりにも小さい周波数に設定すると、符号間干渉により信号レベルが大きく変動する傾向ある低域を、ステレオ測定信号を用いて測定することができなくなってしまう。このため、LPF部22aおよびHPF部22bにおけるカットオフ周波数は、チャンネル間干渉によりディップが発生され得る中高域の周波数よりも低い周波数であって、符号間干渉が発生するおそれのある低域よりも高い周波数に設定することが望ましい。このようにカットオフ周波数を設定することにより、チャンネル間干渉の抑制と、符号間干渉の低減とを同時に実現することができ、精度良く音場環境の周波数特性を測定することが可能となる。
以上、本発明に係る音場測定装置、音場測定方法および音場測定プログラムについて、図面を用いて詳細に説明を行ったが、本発明に係る音場測定装置、音場測定方法および音場測定プログラムは、実施の形態で説明した例には限定されない。当業者であれば、請求の範囲に記載された範疇内において、各種の変更例または修正例に想到しうることは明らかであり、それらについても本発明の技術的範囲に属するものと了解される。
例えば、実施の形態に係る音場測定装置1では、図1に示すようにROM3や記憶部5に記録される処理プログラムあるいはアプリケーションプログラムに基づいて、CPU2が図2に示すような機能部の機能を実現する旨の説明を行ったが、各機能部の機能を実現するCPUは、1つには限定されない。各機能部の一部の機能を実現するための専用の処理部(例えば、複数のCPUなど)を設置し、各専用の処理部が少なくとも1以上の機能を実現する構成であってもよい。このように、複数の専用の処理部を設定した場合であっても、1つのCPUが処理プログラム等に基づいて音場測定処理を行う場合であっても、モノラルとステレオの測定信号を組み合わせることによって、チャンネル間干渉や符号間干渉の影響を低減し、チャンネル毎に時間分割することなく、短時間に精度良く音場環境の周波数特性を測定することが可能である。
1 …音場測定装置
2 …CPU
3 …ROM
4 …RAM
5 …記憶部
6 …外部出力部
7 …マイクロフォン
8 …表示部
11 …測定信号生成部
12 …周波数分割合成部
13 …フーリエ変換部
14 …平均化処理部
21a …第1測定信号生成部
21b …第2測定信号生成部
22a …LPF部(ローパスフィルタ部)
22b …HPF部(ハイパスフィルタ部)
22c …加算部(合成信号生成部)
22d …遅延部
23a …最大値検出部
23b …平均値算出部
100 …音楽再生機能部
101a …右スピーカ
101b …左スピーカ
102 …ポータブルオーディオ機器(オーディオ機器)
2 …CPU
3 …ROM
4 …RAM
5 …記憶部
6 …外部出力部
7 …マイクロフォン
8 …表示部
11 …測定信号生成部
12 …周波数分割合成部
13 …フーリエ変換部
14 …平均化処理部
21a …第1測定信号生成部
21b …第2測定信号生成部
22a …LPF部(ローパスフィルタ部)
22b …HPF部(ハイパスフィルタ部)
22c …加算部(合成信号生成部)
22d …遅延部
23a …最大値検出部
23b …平均値算出部
100 …音楽再生機能部
101a …右スピーカ
101b …左スピーカ
102 …ポータブルオーディオ機器(オーディオ機器)
Claims (9)
- 一対のスピーカが狭間隔に設置されたオーディオ機器から出力される出力音を集音することにより、音場環境の周波数特性を求める音場測定装置であって、
第1測定信号の低域成分を抽出するローパスフィルタ部と、
前記第1測定信号とは異なる第2測定信号の中高域成分を抽出するハイパスフィルタ部と、
前記ローパスフィルタ部により抽出された前記第1測定信号の低域成分と、前記ハイパスフィルタ部により抽出された前記第2測定信号の中高域成分とを合成して合成信号を生成する合成信号生成部と、
前記ローパスフィルタ部により低域成分が抽出される前の前記第1測定信号を、前記一対のスピーカの何れか一方より出力させると共に、前記合成信号生成部により生成された前記合成信号を前記一対のスピーカの何れか他方より同時に出力させるために、前記オーディオ機器に対して前記第1測定信号と前記合成信号とを出力する外部出力部と、
前記一対のスピーカから同時に出力された前記第1測定信号と前記合成信号とを集音するためのマイクロフォンと、
該マイクロフォンにより集音された信号をフーリエ変換することにより、前記音場環境における周波数特性を求めるフーリエ変換部と
を有することを特徴とする音場測定装置。 - 前記第1測定信号としてM系列符号を用いると共に、前記第2測定信号として前記第1測定信号のM系列符号に直交したM系列符号を用い、
前記フーリエ変換部により求められた周波数特性に基づいて、所定の第1周波数間隔における信号レベルの最大値を、当該第1周波数間隔よりも短い周波数ずつシフトさせながら検出することにより、最大値から成る周波数特性を求める最大値検出部と、
該最大値検出部により検出された最大値から成る周波数特性に基づいて、所定の第2周波数間隔における信号レベルの平均値を、当該第2周波数間隔よりも短い周波数ずつシフトさせながら算出することにより、前記音場環境における周波数特性を求める平均値算出部と
を有することを特徴とする請求項1に記載の音場測定装置。 - 前記ローパスフィルタ部および前記ハイパスフィルタ部において設定されるカットオフ周波数は、低域成分が抽出される前の前記第1測定信号を前記一対のスピーカより同時に出力した場合に、前記フーリエ変換部により求められる周波数特性に発生し得るディップの周波数値よりも低い周波数に設定されること
を特徴とする請求項1に記載の音場測定装置。 - 一対のスピーカが狭間隔に設置されたオーディオ機器から出力される出力音を集音することにより、音場環境の周波数特性を求める音場測定装置の音場測定方法であって、
第1測定信号の低域成分をローパスフィルタ部が抽出する低域成分抽出ステップと、
前記第1測定信号とは異なる第2測定信号の中高域成分をハイパスフィルタ部が抽出する中高域成分抽出ステップと、
前記低域成分抽出ステップにおいて抽出された前記第1測定信号の低域成分と、前記中高域成分抽出ステップにおいて抽出された前記第2測定信号の中高域成分とを合成して、合成信号生成部が合成信号を生成する合成信号生成ステップと、
前記低域成分抽出ステップにおいて低域成分が抽出される前の前記第1測定信号を、前記一対のスピーカの何れか一方より出力させると共に、前記合成信号生成ステップにおいて生成された前記合成信号を前記一対のスピーカの何れか他方より同時に出力させるために、外部出力部が、前記オーディオ機器に対して前記第1測定信号と前記合成信号とを出力する外部出力ステップと、
前記一対のスピーカから同時に出力された前記第1測定信号と前記合成信号とを、マイクロフォンが集音する集音ステップと、
該集音ステップにおいて集音された信号をフーリエ変換することにより、フーリエ変換部が、前記音場環境における周波数特性を求めるフーリエ変換ステップと
を有することを特徴とする音場測定装置の音場測定方法。 - 前記第1測定信号としてM系列符号を用いると共に、前記第2測定信号として前記第1測定信号のM系列符号に直交したM系列符号を用い、
前記フーリエ変換ステップにおいて求められた周波数特性に基づいて、最大値検出部が、所定の第1周波数間隔における信号レベルの最大値を、当該第1周波数間隔よりも短い周波数ずつシフトさせながら検出することにより、最大値から成る周波数特性を求める最大値検出ステップと、
該最大値検出ステップにおいて求められた、最大値から成る周波数特性に基づいて、平均値算出部が、所定の第2周波数間隔における信号レベルの平均値を、当該第2周波数間隔よりも短い周波数ずつシフトさせながら算出することにより、前記音場環境における周波数特性を求める平均値算出ステップと
を有することを特徴とする請求項4に記載の音場測定装置の音場測定方法。 - 前記低域成分抽出ステップおよび前記中高域成分抽出ステップにおいて設定されるカットオフ周波数は、低域成分が抽出される前の前記第1測定信号を前記一対のスピーカより同時に出力した場合に、前記フーリエ変換ステップにおいて求められる周波数特性に発生し得るディップの周波数値よりも低い周波数に設定されること
を特徴とする請求項4に記載の音場測定装置の音場測定方法。 - 一対のスピーカが狭間隔に設置されたオーディオ機器から出力される出力音を集音することにより、音場環境の周波数特性を求める音場測定装置の音場測定プログラムであって、
前記音場測定装置のコンピュータに、
第1測定信号の低域成分を抽出させるローパスフィルタ機能と、
前記第1測定信号とは異なる第2測定信号の中高域成分を抽出させるハイパスフィルタ機能と、
前記ローパスフィルタ機能により抽出された前記第1測定信号の低域成分と、前記ハイパスフィルタ機能により抽出された前記第2測定信号の中高域成分とを合成して合成信号を生成させる合成信号生成機能と、
前記ローパスフィルタ機能により低域成分が抽出される前の前記第1測定信号を、前記一対のスピーカの何れか一方より出力させると共に、前記合成信号生成機能により生成された前記合成信号を前記一対のスピーカの何れか他方より同時に出力させるために、前記オーディオ機器に対して前記第1測定信号と前記合成信号とを出力させる外部出力機能と、
前記一対のスピーカから同時に出力された前記第1測定信号と前記合成信号とをマイクロフォンを用いて集音させる集音機能と、
該集音機能により集音された信号をフーリエ変換することにより、前記音場環境における周波数特性を求めさせるフーリエ変換機能と
を実現させるための音場測定プログラム。 - 前記第1測定信号としてM系列符号を用いると共に、前記第2測定信号として前記第1測定信号のM系列符号に直交したM系列符号を用い、
前記コンピュータに、
前記フーリエ変換機能により求められた周波数特性に基づいて、所定の第1周波数間隔における信号レベルの最大値を、当該第1周波数間隔よりも短い周波数ずつシフトさせながら検出させることにより、最大値から成る周波数特性を求めさせる最大値検出機能と、
該最大値検出機能により求められた最大値から成る周波数特性に基づいて、所定の第2周波数間隔における信号レベルの平均値を、当該第2周波数間隔よりも短い周波数ずつシフトさせながら算出させることにより、前記音場環境における周波数特性を求めさせる平均値算出機能と
を実現させるための請求項7に記載の音場測定プログラム。 - 前記ローパスフィルタ機能および前記ハイパスフィルタ機能において設定されるカットオフ周波数は、低域成分が抽出される前の前記第1測定信号を前記一対のスピーカより同時に出力した場合に、前記フーリエ変換機能により求められる周波数特性に発生し得るディップの周波数値よりも低い周波数に設定されること
を特徴とする請求項7に記載の音場測定プログラム。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
EP14851266.8A EP3054707B1 (en) | 2013-10-01 | 2014-09-19 | Device, method, and program for measuring sound field |
CN201480053267.5A CN105659631B (zh) | 2013-10-01 | 2014-09-19 | 音场测量装置和音场测量方法 |
US14/916,567 US9883303B2 (en) | 2013-10-01 | 2014-09-19 | Sound field measuring device, method and program |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2013-206163 | 2013-10-01 | ||
JP2013206163A JP6138015B2 (ja) | 2013-10-01 | 2013-10-01 | 音場測定装置、音場測定方法および音場測定プログラム |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2015050006A1 true WO2015050006A1 (ja) | 2015-04-09 |
Family
ID=52778599
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/JP2014/074848 WO2015050006A1 (ja) | 2013-10-01 | 2014-09-19 | 音場測定装置、音場測定方法および音場測定プログラム |
Country Status (5)
Country | Link |
---|---|
US (1) | US9883303B2 (ja) |
EP (1) | EP3054707B1 (ja) |
JP (1) | JP6138015B2 (ja) |
CN (1) | CN105659631B (ja) |
WO (1) | WO2015050006A1 (ja) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113965859A (zh) * | 2020-07-20 | 2022-01-21 | Jvc建伍株式会社 | 头外定位滤波器确定系统、方法以及程序 |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6573657B2 (ja) * | 2015-02-27 | 2019-09-11 | パイオニア株式会社 | 音量制御装置、音量制御方法、及び、音量制御プログラム |
US11140483B2 (en) | 2019-03-05 | 2021-10-05 | Maxim Integrated Products, Inc. | Management of low frequency components of an audio signal at a mobile computing device |
WO2021120033A1 (zh) * | 2019-12-18 | 2021-06-24 | 瑞声声学科技(深圳)有限公司 | 振动信号的处理方法、装置、设备及可读介质 |
JP2021164109A (ja) * | 2020-04-02 | 2021-10-11 | アルプスアルパイン株式会社 | 音場補正方法、音場補正プログラムおよび音場補正システム |
JP7552089B2 (ja) * | 2020-06-18 | 2024-09-18 | ヤマハ株式会社 | 音響特性の補正方法および音響特性補正装置 |
CN111751683A (zh) * | 2020-07-03 | 2020-10-09 | 国网河南省电力公司检修公司 | Gis设备局部放电检测方法 |
CN114501283B (zh) * | 2022-04-15 | 2022-06-28 | 南京天悦电子科技有限公司 | 一种针对数字助听器的低复杂度双麦克风定向拾音方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0775190A (ja) | 1993-09-03 | 1995-03-17 | Fujitsu Ten Ltd | M系列信号を用いた車室内音響自動補正システム |
JP2006262015A (ja) * | 2005-03-16 | 2006-09-28 | Clarion Co Ltd | 音場特性測定方法、音場特性測定システム、増幅装置、及び音場特性測定装置 |
JP2007232492A (ja) | 2006-02-28 | 2007-09-13 | Kawai Musical Instr Mfg Co Ltd | 伝達特性測定方法および装置 |
JP2008177795A (ja) * | 2007-01-17 | 2008-07-31 | Alpine Electronics Inc | ノイズ出力装置、ノイズ出力方法及びノイズ出力プログラム |
WO2009008068A1 (ja) * | 2007-07-11 | 2009-01-15 | Pioneer Corporation | 自動音場補正装置 |
WO2009095965A1 (ja) * | 2008-01-31 | 2009-08-06 | Mitsubishi Electric Corporation | 帯域分割時間補正信号処理装置 |
JP2013527491A (ja) * | 2010-04-09 | 2013-06-27 | ディーティーエス・インコーポレイテッド | オーディオ再生のための適応的環境ノイズ補償 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4171675B2 (ja) * | 2003-07-15 | 2008-10-22 | パイオニア株式会社 | 音場制御システム、および音場制御方法 |
JP2005252467A (ja) | 2004-03-02 | 2005-09-15 | Sony Corp | 音響再生方法、音響再生装置および記録メディア |
JP2007142875A (ja) * | 2005-11-18 | 2007-06-07 | Sony Corp | 音響特性補正装置 |
WO2009095331A1 (en) * | 2008-01-29 | 2009-08-06 | Audioasics A/S | Signal conditioner with suppression of interfering signals |
US20110058684A1 (en) * | 2008-05-27 | 2011-03-10 | Pioneer Corporation | Audio reproduction device |
CN104689402B (zh) * | 2009-06-26 | 2017-06-13 | 甘布罗伦迪亚股份公司 | 透析机以及用于对信号进行处理的装置和方法 |
-
2013
- 2013-10-01 JP JP2013206163A patent/JP6138015B2/ja active Active
-
2014
- 2014-09-19 WO PCT/JP2014/074848 patent/WO2015050006A1/ja active Application Filing
- 2014-09-19 EP EP14851266.8A patent/EP3054707B1/en active Active
- 2014-09-19 US US14/916,567 patent/US9883303B2/en active Active
- 2014-09-19 CN CN201480053267.5A patent/CN105659631B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0775190A (ja) | 1993-09-03 | 1995-03-17 | Fujitsu Ten Ltd | M系列信号を用いた車室内音響自動補正システム |
JP2006262015A (ja) * | 2005-03-16 | 2006-09-28 | Clarion Co Ltd | 音場特性測定方法、音場特性測定システム、増幅装置、及び音場特性測定装置 |
JP2007232492A (ja) | 2006-02-28 | 2007-09-13 | Kawai Musical Instr Mfg Co Ltd | 伝達特性測定方法および装置 |
JP2008177795A (ja) * | 2007-01-17 | 2008-07-31 | Alpine Electronics Inc | ノイズ出力装置、ノイズ出力方法及びノイズ出力プログラム |
WO2009008068A1 (ja) * | 2007-07-11 | 2009-01-15 | Pioneer Corporation | 自動音場補正装置 |
WO2009095965A1 (ja) * | 2008-01-31 | 2009-08-06 | Mitsubishi Electric Corporation | 帯域分割時間補正信号処理装置 |
JP2013527491A (ja) * | 2010-04-09 | 2013-06-27 | ディーティーエス・インコーポレイテッド | オーディオ再生のための適応的環境ノイズ補償 |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113965859A (zh) * | 2020-07-20 | 2022-01-21 | Jvc建伍株式会社 | 头外定位滤波器确定系统、方法以及程序 |
Also Published As
Publication number | Publication date |
---|---|
EP3054707A1 (en) | 2016-08-10 |
JP2015070589A (ja) | 2015-04-13 |
US20160219385A1 (en) | 2016-07-28 |
CN105659631A (zh) | 2016-06-08 |
US9883303B2 (en) | 2018-01-30 |
JP6138015B2 (ja) | 2017-05-31 |
EP3054707B1 (en) | 2018-04-11 |
CN105659631B (zh) | 2017-06-09 |
EP3054707A4 (en) | 2017-05-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP6138015B2 (ja) | 音場測定装置、音場測定方法および音場測定プログラム | |
JP5898534B2 (ja) | 音響信号処理装置および音響信号処理方法 | |
EP2291002A1 (en) | Acoustic processing apparatus | |
EP3016412A1 (en) | Propagation delay correction apparatus and propagation delay correction method | |
JP6401521B2 (ja) | 信号処理装置及び信号処理方法 | |
US11032644B2 (en) | Subband spatial and crosstalk processing using spectrally orthogonal audio components | |
JP6533959B2 (ja) | 音声信号処理装置および音声信号処理方法 | |
EP2949133A1 (en) | Automatic loudspeaker polarity detection | |
WO2016140121A1 (ja) | フィルタ生成装置、フィルタ生成方法およびフィルタ生成プログラム | |
JP6096956B2 (ja) | 周波数に依存して入力信号の雑音を抑圧するための方法 | |
CN105324815A (zh) | 信号处理装置和信号处理方法 | |
JP6151619B2 (ja) | 音場測定装置、音場測定方法および音場測定プログラム | |
CN109791773B (zh) | 音频输出产生系统、音频通道输出方法和计算机可读介质 | |
US10524052B2 (en) | Dominant sub-band determination | |
JP6155132B2 (ja) | 低域補完装置および低域補完方法 | |
JP2016123059A (ja) | 信号生成装置、および端末装置 | |
US9571950B1 (en) | System and method for audio reproduction | |
JP6438786B2 (ja) | 信号中に含まれる正弦成分抽出装置、正弦成分抽出方法及びプログラム | |
JP2012100117A (ja) | 音響処理装置及び方法 | |
JP2011035573A (ja) | 音信号処理装置および音信号処理方法 | |
JP7286532B2 (ja) | 信号処理装置、音響装置、信号処理方法及び信号処理プログラム | |
JP5652515B2 (ja) | 信号処理装置 | |
JP6314803B2 (ja) | 信号処理装置、信号処理方法及びプログラム |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 14851266 Country of ref document: EP Kind code of ref document: A1 |
|
WWE | Wipo information: entry into national phase |
Ref document number: 14916567 Country of ref document: US |
|
REEP | Request for entry into the european phase |
Ref document number: 2014851266 Country of ref document: EP |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2014851266 Country of ref document: EP |
|
NENP | Non-entry into the national phase |
Ref country code: DE |