WO2015162736A1 - 検査装置 - Google Patents

検査装置 Download PDF

Info

Publication number
WO2015162736A1
WO2015162736A1 PCT/JP2014/061479 JP2014061479W WO2015162736A1 WO 2015162736 A1 WO2015162736 A1 WO 2015162736A1 JP 2014061479 W JP2014061479 W JP 2014061479W WO 2015162736 A1 WO2015162736 A1 WO 2015162736A1
Authority
WO
WIPO (PCT)
Prior art keywords
filter
level difference
signal
filters
inspection
Prior art date
Application number
PCT/JP2014/061479
Other languages
English (en)
French (fr)
Inventor
捷吉郎 日野
耕一 土屋
聡彦 庄司
友彦 遠藤
英二 樋渡
Original Assignee
エタニ電機株式会社
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 エタニ電機株式会社 filed Critical エタニ電機株式会社
Priority to PCT/JP2014/061479 priority Critical patent/WO2015162736A1/ja
Priority to JP2014550576A priority patent/JP5896434B1/ja
Publication of WO2015162736A1 publication Critical patent/WO2015162736A1/ja

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04RLOUDSPEAKERS, MICROPHONES, GRAMOPHONE PICK-UPS OR LIKE ACOUSTIC ELECTROMECHANICAL TRANSDUCERS; DEAF-AID SETS; PUBLIC ADDRESS SYSTEMS
    • H04R29/00Monitoring arrangements; Testing arrangements

Definitions

  • the present application relates to an inspection method for inspecting a device that generates sound or vibration, an inspection device, a filter determination method for determining a filter used in the inspection method, a program, and a filter determination device.
  • Patent Document 1 describes a detection device that detects a sticking sound generated by a speaker and determines whether the speaker is good or bad.
  • the detection device according to Patent Literature 1 detects a cluttering sound by extracting a signal component of the cluttering sound from a signal based on the sound generated by the speaker using a filter optimized in advance.
  • Patent Document 1 does not describe any specific description regarding “pre-optimized filter”.
  • the sound or vibration generated by an abnormal device varies depending on the specifications of the device. Therefore, in order to accurately determine the quality of various devices, a filter optimized for each specification is required. However, it is practically difficult to optimize the filter in advance for each device specification.
  • the present application has been made in view of such circumstances.
  • the purpose of the present application is to use an inspection method, an inspection apparatus, and an inspection method capable of inspecting the apparatus by determining an optimum inspection filter from a plurality of filters for each specification of the apparatus that generates sound or vibration. It is an object to provide a filter determination method, a program, and a filter determination device that can determine a filter. Furthermore, an object of the present application is to provide an inspection method and an inspection apparatus capable of inspecting a device that generates sound or vibration whose frequency changes with time, using a filter whose characteristics can be changed with time. is there.
  • the inspection method processes a signal related to a sound or vibration whose frequency changes over time generated by a device to be inspected using a filter capable of changing characteristics over time, and inspects the device. To do.
  • a signal relating to sound or vibration generated by one or a plurality of non-defective products and one or a plurality of defective devices prepared in advance is provided for each of a plurality of filters.
  • a test filter is determined, and using the determined test filter, a signal related to sound or vibration generated by the device to be tested is processed, and the processed signal Based on the above, the quality of the device to be inspected is determined.
  • the inspection method calculates a level difference between signals related to sound or vibration generated by the non-defective product and the defective product, processes a signal corresponding to the calculated level difference with each of a plurality of filters, and each of the plurality of filters.
  • the filter corresponding to the maximum level difference is determined as the inspection filter.
  • a level difference of a signal related to sound or vibration generated by each of the non-defective product and the plurality of defective products is calculated for each of the plurality of defective products, and a signal corresponding to the calculated level difference is a plurality of filters.
  • the maximum level difference corresponding to the processed signal is calculated for each of the plurality of defective products, and among each maximum value of the level difference calculated for each of the plurality of defective products for each of the plurality of filters.
  • the minimum value of the level difference is selected, and the filter corresponding to the maximum level difference among the minimum values of the level difference selected for each of the plurality of filters is determined as the inspection filter.
  • the inspection method is based on a signal based on time, frequency, and level information respectively calculated from a plurality of signals related to non-defective products prepared in advance, and the maximum for each time and frequency defined at a predetermined interval.
  • Each level is extracted, each extracted maximum level is set to the level of the signal related to the sound or vibration generated by the non-defective product, and the maximum value of the level difference is determined based on the set level of the signal related to the non-defective product. It is calculated for each of a plurality of defective products.
  • the selection of the minimum value of the level difference is finished, and the filter corresponding to the minimum value of the level difference exceeding the predetermined value Is determined as a filter for inspection.
  • the inspection apparatus processes a signal relating to a sound or vibration whose frequency changes with time, which is generated by the apparatus to be inspected, using a filter capable of changing characteristics over time, and inspects the apparatus. To do.
  • the inspection apparatus is based on a filter output unit that outputs a plurality of filters and a signal relating to sound or vibration generated by one or more non-defective products and one or more defective products prepared in advance.
  • a signal relating to sound or vibration generated by a device to be inspected is determined using a determination unit that determines a filter for inspection from a plurality of filters output by the filter output unit, and the filter for inspection determined by the determination unit.
  • a processing unit for processing and a determination unit for determining the quality of the device to be inspected based on a signal processed by the processing unit are provided.
  • the inspection apparatus calculates a level difference between a non-defective signal storage unit that stores a signal related to the non-defective product, a signal related to the non-defective product stored in the non-defective product signal storage unit, and a signal related to the defective product prepared in advance.
  • a level difference calculation unit and a signal processing unit that processes a signal corresponding to the level difference calculated by the level difference calculation unit in each of the plurality of filters output by the filter output unit, and the determination unit includes the signal Among the level differences corresponding to the signals processed by the processing unit, the filter corresponding to the maximum level difference is determined as the inspection filter, and the signal for comparing the signal processed by the processing unit and a predetermined signal A comparison unit is provided, and the determination unit is configured to determine the quality of the device to be inspected based on the result of comparison by the signal comparison unit.
  • the level difference calculation unit calculates a level difference of a signal related to sound or vibration generated by the non-defective product and the plurality of defective products for each of the plurality of defective products
  • the signal processing unit is configured to process a signal corresponding to the level difference calculated for each of the plurality of defective products by the plurality of filters output from the filter output unit, and the signal processing unit For each of the plurality of defective products, and for each of the plurality of filters output by the filter output unit, the calculation unit is provided for each of the plurality of defective products.
  • a selection unit that selects a minimum value of the level difference among the calculated maximum values of the level difference, and the determination unit selects a maximum value among the minimum values of the level difference selected by the selection unit for each of the plurality of filters. No Characterized in that Le difference are so as to determine a filter for checking the corresponding filter.
  • the selection unit ends selection of the minimum value of the level difference, and the determination unit The filter corresponding to the minimum value of the level difference exceeding the value is determined as the inspection filter.
  • the filter determination method determines a test filter from a plurality of filters whose characteristics can be changed over time in order to inspect a device to be inspected that generates sound or vibration whose frequency changes over time.
  • a method for determining a filter wherein a level difference of a signal relating to sound or vibration generated by a non-defective product and one or a plurality of defective products prepared in advance is calculated for each of the plurality of defective products, and the calculated level difference is calculated.
  • a corresponding signal is processed by each of a plurality of filters, a maximum level difference corresponding to the processed signal is calculated for each of the one or more defective products, and each of the one or more defective products is calculated for each of the plurality of filters.
  • the minimum value of the level difference is selected from the maximum values of the level difference calculated in step S1, and the maximum level difference is selected from the minimum values of the level differences selected for each of the plurality of filters.
  • the program according to the present application is a program that causes a computer to execute a process of determining one filter from a plurality of filters whose characteristics can be changed over time based on a signal related to sound or vibration whose frequency changes over time. Calculating a level difference between one signal related to sound or vibration and one or a plurality of other signals for each of the one or a plurality of other signals, and changing a filter capable of changing the characteristics over time, Each of the plurality of changed filters processes a signal corresponding to the calculated level difference, calculates a maximum value of the level difference corresponding to the processed signal for each of the one or more other signals, and each of the plurality of filters.
  • the filter determination device is an inspection device for inspecting a device to be inspected that generates a sound or vibration whose frequency changes over time, from a plurality of filters whose characteristics can be changed over time.
  • a filter output unit that outputs a plurality of filters, and a level difference between signals relating to sound or vibration generated by a non-defective product and one or more defective products prepared in advance.
  • a level difference calculation unit that calculates for each of a plurality of defective products, a signal processing unit that processes a signal corresponding to the level difference calculated by the level difference calculation unit, and a plurality of filters that are output by the filter output unit;
  • a calculation unit that calculates a maximum level difference corresponding to the signal processed by the signal processing unit for each of the one or a plurality of defective products, and a plurality of frames output by the filter output unit; For each filter, the calculation unit selects a minimum value of the level difference among the maximum values of the level difference calculated for each one or a plurality of defective products, and the selection unit selects each of the plurality of filters.
  • a determination unit is provided that determines a filter corresponding to the maximum level difference among the minimum values of the level differences, as an inspection filter.
  • the optimal inspection filter is determined from a plurality of filters for each specification of the apparatus that generates sound or vibration.
  • the device can be inspected.
  • a device that generates a sound or vibration whose frequency changes with the passage of time is applied to the change of the sound or vibration. Can be inspected.
  • FIG. 10 is an explanatory diagram illustrating an outline of operation of an inspection system according to a second embodiment.
  • FIG. It is explanatory drawing which shows notionally the process which extracts defect detection capability value D (C) from detection value Over (ii, C) of n defective products. It is explanatory drawing which shows the example of a pattern. It is explanatory drawing which shows the example of a pattern. It is explanatory drawing which shows the example of a pattern. It is explanatory drawing which shows the image of the process which calculates
  • the inspection apparatus according to the present embodiment will be described with reference to the drawings.
  • the inspection apparatus according to the present embodiment inspects an apparatus that generates sound or vibration, and performs pass / fail determination of the apparatus.
  • the inspection apparatus according to the present embodiment will be described using a speaker as an example of an apparatus to be inspected.
  • FIG. 1 is a block diagram illustrating a hardware configuration example of an inspection system including an inspection apparatus according to the present application.
  • the inspection system includes a filter determination device 20, a measurement unit 30, a microphone 40, and an anechoic box 50.
  • An inspection device for inspecting the speaker 60 to be inspected includes a filter determination device 20 and a measurement unit 30.
  • the inspection apparatus uses the filter determination device 20 as a main part.
  • the filter determination device 20 uses a personal computer.
  • the filter determination device 20 outputs data corresponding to a sine wave sweep signal changing from a low frequency to a high frequency to the measurement unit 30.
  • the filter determination device 20 inputs a signal corresponding to the sound generated by the speaker 60 from the measurement unit 30.
  • the filter determination device 20 executes various processes based on the signal input from the measurement unit 30.
  • the filter determination device 20 includes a CPU 21, ROM 22, RAM 23, hard disk 24, disk drive 25, display unit 26, operation unit 27, communication unit 28, and I / O unit 29.
  • the CPU 21 is a processor that executes an instruction set described in the program 2P.
  • the ROM 22 stores a BIOS and the like executed by the CPU 21 when the filter determination device 20 is activated.
  • the RAM 23 is a main storage device.
  • the hard disk 24 is an auxiliary storage device.
  • the hard disk 24 stores a sound signal from the speaker 60 amplified by the measurement unit 30.
  • the disk drive 25 reads information from the optical disk 2d of the external storage medium and records information on the optical disk 2d.
  • the display unit 26 is a display device that displays an image.
  • the operation unit 27 is an input device that accepts various inputs from the user.
  • the communication unit 28 is an interface with the LAN 100.
  • the I / O unit 29 is an interface with the measurement unit 30.
  • the CPU 21 may read the program 2P from a portable storage medium such as the optical disk 2d via the disk drive 25.
  • the CPU 21 may read the program 2P from another information processing apparatus or storage device via the communication unit 28.
  • a semiconductor memory 2m such as a flash memory storing the program 2P may be mounted in the filter determination device 20.
  • the measurement unit 30 inputs a sweep signal to the speaker 60.
  • the measurement unit 30 includes a D / A converter 31, a power amplifier 32, a measurement amplifier 33, and an A / D converter 34.
  • the D / A converter 31 converts the digital sweep signal input from the filter determination device 20 into an analog signal, and outputs the converted analog signal to the power amplifier 32.
  • the power amplifier 32 amplifies the analog signal converted from the D / A converter 31.
  • the power amplifier 32 outputs a sweep signal to the speaker 60.
  • the anechoic box 50 has high sound insulation and sound absorption characteristics.
  • the upper wall of the anechoic box 50 facing the speaker 60 is provided with an opening for taking in sound generated by the speaker 60.
  • the microphone 40 is installed inside the anechoic box 50 so as to face the speaker 60.
  • the measurement amplifier 33 amplifies the level of the analog signal output from the microphone 40.
  • the A / D converter 34 converts the signal amplified by the measurement amplifier 33 into a digital signal.
  • the A / D converter 34 outputs the converted digital signal to the filter determination device 20.
  • FIG. 2 is an explanatory diagram showing an outline of the operation of the inspection system.
  • An auditory test is performed on a plurality of speakers 60 in advance, and one good product and a plurality of defective products are prepared as samples to be processed for determining a filter for inspection.
  • One prepared speaker 60 is attached to the anechoic box 50.
  • the filter determination apparatus 20 inputs data corresponding to a sweep signal of a sine wave that changes at a constant speed from a low frequency (for example, 10 Hz) to a high frequency (for example, 20 kHz) to the D / A converter 31 of the measurement unit 30.
  • the sweep time is, for example, 1.0 second.
  • the speaker 60 generates a sound corresponding to the sweep signal corresponding to the data input from the filter determination device 20 via the D / A converter 31 and the power amplifier 32.
  • the microphone 40 receives the sound generated by the sample speaker 60.
  • the microphone 40 outputs a signal to the filter determination device 20 via the measurement amplifier 33 and the A / D converter 34.
  • the CPU 21 inputs a signal from the measurement unit 30.
  • the CPU 21 stores the input signal in the hard disk 24 in association with the non-defective / defective product identification information related to the sample speaker 60.
  • the inspection system repeats the same storage process for all the prepared speakers 60 for the sample.
  • the CPU 21 reads the signal stored in the hard disk 24 together with the non-defective / defective product identification information related to the sample speaker 60 corresponding to the signal.
  • the CPU 21 calculates a three-dimensional spectrogram of time, frequency, and level from the read signal, and integrates the level at the time from the maximum frequency to the frequency for each frequency at each time.
  • FIG. 8 shows an example of a spectrogram
  • FIG. 9 shows a cross-section of a spectrogram at a certain time, that is, an example of a frequency distribution of levels
  • FIG. 10 shows an accumulation of levels from the maximum frequency to that frequency at that time. .
  • the integrated graph at all times is as shown in FIG. Details of FIGS. 8 to 11 will be described in the second embodiment.
  • the CPU 21 stores the calculated integrated graph in the hard disk 24 along with the good / bad classification.
  • the CPU 21 repeats this integrated graph calculation process for all the signals of the sample speaker 60 stored in the hard disk 24.
  • a defective speaker 60 generates an abnormal sound mainly at a high frequency with respect to an input signal. Abnormal noise is often lower in level than the fundamental, second, third, etc. harmonics of the input signal, so it is necessary to remove lower harmonics.
  • a high-pass filter that passes through is used. A level obtained by passing a signal of a speaker 60 for a sample through a high-pass filter having a frequency at a certain time as a cutoff frequency becomes a level at the frequency of the integrated graph of the speaker 60 at the time.
  • the CPU 21 determines the filter for processing the sound signal so that the level difference between the non-defective product and the defective product becomes large in order to clearly distinguish the speaker 60 to be inspected from the good product and the defective product. .
  • the CPU 21 dynamically changes the cutoff frequency of the high-pass filter so that the cutoff frequency of the high-pass filter follows the signal from the speaker 60 according to the sweep signal. That is, the filter used by the CPU 21 can be said to be a dynamic filter 1f.
  • the dynamic filter 1f is represented by a time change of the cutoff frequency.
  • the horizontal axis of the graph showing the dynamic filter 1f is time
  • the vertical axis is frequency. That is, the dynamic filter 1f can be expressed by a two-dimensional curve having time and frequency as elements.
  • the level graph obtained as a result of processing the signal of the speaker 60 for a sample by the dynamic filter 1f is a cut surface obtained by cutting the integrated graph of the speaker 60 by the curve of the dynamic filter 1f.
  • the CPU 21 uses the dynamic filter 1f to process components common to both good / defective samples.
  • CPU21 determines the filter for test
  • the CPU 21 inspects the speaker 60 to be inspected using the determined dynamic filter 1f. Since the dynamic filter 1f is a digital filter created by the CPU 21, the filter characteristics of the dynamic filter 1f are determined by filter parameters.
  • a filter determination process for inspecting the speaker 60 to be inspected will be described. If any characteristic change filter is allowed, it is sufficient to create a filter with a curve that passes through all the points of the maximum value of the difference graph of each defective product. It is difficult to realize such a filter.
  • the shape of the filter needs to be limited to a certain pattern that can be used in an actual inspection line.
  • a predetermined pattern of filters is stored in a plurality of hard disks 24.
  • the cut-off frequency continuously changes from a low frequency to a high frequency as time elapses.
  • the pattern includes a plurality of patterns that change following the frequency change of the fundamental wave or a plurality of harmonics from a low order to a high order.
  • the CPU 21 may generate a plurality of constant pattern filters in real time, and output each of the generated constant pattern filters for the following processing.
  • CPU21 selects one filter from the prepared filters.
  • CPU21 acquires the cut surface cut
  • the cut surface becomes a level graph obtained by filtering non-defective products with the filter.
  • FIG. 2 schematically shows level graphs 1g and 2g of the non-defective product and the defective product after processing, with the horizontal axis representing time and the vertical axis representing level.
  • the CPU 21 may store a non-defective level graph 1g that has been subjected to the filter processing in the hard disk 24.
  • the CPU 21 calculates the maximum value of the difference graph between the processed non-defective product level graph and the level graph of one defective product as a detected value of the defective product by the filter.
  • the CPU 21 stores the calculated detection value in the hard disk 24 in association with the identification information of the defective product.
  • the CPU 21 repeats the process of calculating the detection value for all defective samples.
  • FIG. 2 shows a list 200 of detection values stored in the hard disk 24. The larger the detected value, the farther the defective product is from the non-defective product.
  • the filter for discriminating whether the speaker 60 to be inspected is a non-defective product or a defective product needs to be a filter in which all defective products prepared as samples are determined as defective products. For that purpose, it is preferable from the viewpoint of safety to adopt a level difference corresponding to a defective product having the sound output characteristic closest to the non-defective product.
  • the CPU 21 extracts the minimum value as the defect detection capability value of the filter from the detection values calculated for all the defective products in the list 200.
  • the CPU 21 stores the extracted defect detection capability value in the hard disk 24 as a list 300 in association with the selected filter.
  • the CPU 21 repeatedly stores the defect detection capability value in the list 300 for each filter.
  • the CPU 21 selects the filter corresponding to the maximum value from the plurality of defect detection capability values stored as the list 300 in the hard disk 24 for inspection. Determine the filter.
  • the CPU 21 changes the filter and finds a filter having a defect detection capability value greater than or equal to a defect detection capability value that the user thinks is sufficient for separating a non-defective product and a defective product, for example, May be determined as a filter for inspection. As a result, the processing time for determining the inspection filter can be shortened.
  • FIG. 3 is a flowchart illustrating an example of an outline procedure of the filter determination process.
  • the CPU 21 selects an unprocessed filter (step S101).
  • the CPU 21 selects an unprocessed defective product integration graph (step S102).
  • the CPU 21 processes the integrated graph relating to the non-defective product and the defective product with the selected filter (step S103).
  • the CPU 21 calculates the maximum value of the difference graph between the processed non-defective product level graph and the defective product level graph as the detected value of the defective product by the filter (step S104).
  • the CPU 21 stores the calculated detection value in association with the defective product identification information in the list 200 in the hard disk 24 (step S105).
  • CPU21 determines whether the process was complete
  • CPU21 returns a process to step S102, when it determines with the process not being complete
  • the CPU 21 determines that the processing has been completed for the integration graph of all defective products (step S106: YES)
  • the CPU 21 determines the minimum value from the plurality of detection values stored in the list 200 in the hard disk 24 as the defect detection capability value of the filter. (Step S107).
  • the CPU 21 stores the extracted defect detection capability value in the list 300 in the hard disk 24 in association with the selected filter (step S108).
  • step S109 determines whether the process was completed about all the filters. If the CPU 21 determines that the processing has not been completed for all the filters (step S109: NO), the CPU 21 returns the processing to step S101. If the CPU 21 determines that the process has been completed for all filters (step S109: YES), the CPU 21 extracts the maximum value from the defect detection capability values stored in the list 300 in the hard disk 24 (step S110). The CPU 21 determines a filter corresponding to the extracted maximum value as an inspection filter (step S111), and ends the process.
  • the inspector attaches the speaker 60 to be inspected to the upper wall of the anechoic box 50.
  • the inspector operates the filter determination device 20 to generate a sound corresponding to the sweep signal on the speaker 60 to be inspected.
  • the microphone 40 receives sound from the speaker 60 to be inspected, and outputs the received sound signal to the filter determination device 20 via the measurement unit 30.
  • the filter determination device 20 inspects the speaker 60 to be inspected based on the input signal and the inspection filter.
  • the filter determination apparatus 20 displays the result of inspecting the inspection target speaker 60 on the display unit 26, for example.
  • the inspector removes the speaker 60 to be inspected from the upper wall of the anechoic box 50.
  • the inspector sorts the speaker 60 to be inspected into a non-defective product / defective product according to the result displayed on the display unit 26. Thereafter, the inspector repeats the above work. Of course, the work performed by the inspector may be automated.
  • FIG. 4 is a flowchart illustrating an example of the procedure of the inspection process.
  • the CPU 21 outputs data corresponding to the sweep signal to the measurement unit 30 (step S1).
  • the CPU 21 inputs a sound signal generated by the speaker 60 to be inspected from the measurement unit 30 from the measurement unit 30 (step S2).
  • the CPU 21 filters the input signal with the inspection filter determined in the filter determination process (step S3).
  • the CPU 21 compares the filtered signal with inspection reference data prepared in advance (step S4).
  • the CPU 21 determines whether or not the speaker 60 to be inspected is a non-defective product based on the comparison result (step S5).
  • CPU21 outputs a determination result to the display part 26 (step S6).
  • CPU21 determines whether the instruction
  • CPU21 returns a process to step S1, when it determines with the instruction
  • the case where the CPU 21 has not received an instruction to repeat the inspection process from the operation unit 27 is a case where the CPU 21 has received an instruction to stop the inspection process.
  • the filter determination device 20 inputs a sweep signal that changes from a low frequency to a high frequency to the speaker 60.
  • the filter determination device 20 may input a reverse sweep signal that changes from a high frequency to a low frequency to the speaker 60.
  • the filter determination device 20 inputs a sweep signal to the speaker 60.
  • the filter determination device 20 may input a step sine wave whose frequency is switched and a sine wave having a fixed frequency to the speaker 60. In such a case, the time change of the frequency may not be monotonously increasing or monotonically decreasing.
  • the filter determination device 20 may input a synthesized wave having a plurality of frequencies to the speaker 60.
  • the apparatus to be inspected is a speaker 60 that generates sound based on a signal input from the outside.
  • the device to be inspected may be a device that itself generates sound or vibration whose frequency changes over time.
  • the devices to be inspected are rotating members such as motors, gears and bearings, elevator cars, compressors, and the like.
  • the pattern prepared in advance for determining the filter for inspection may include a pattern in which the aspect of change in which the cutoff frequency increases changes with time.
  • FIG. 5 is an explanatory diagram showing an example of a filter in which the mode of change in which the cutoff frequency becomes higher changes with time. The time on the horizontal axis is a logarithm of milliseconds.
  • the cut-off frequency first changes following the frequency change of the higher-order harmonics, and maintains a constant frequency from the middle. Thereafter, the cut-off frequency transitions to the lower-order harmonic side, but deviates from the frequency change of the harmonic.
  • a plurality of filters having patterns as shown in FIG. 5 are prepared.
  • Patterns prepared in advance for determining a filter for inspection include not only the cutoff frequency that increases monotonously with time but also a portion that has a portion that changes to have a maximum value and a minimum value. But you can.
  • FIG. 6 is an explanatory diagram showing an example of a filter having a portion where the cutoff frequency becomes lower as time elapses. The time on the horizontal axis is a logarithm of milliseconds.
  • the cut-off frequency first changes following the frequency change of the higher order harmonics, becomes lower than the harmonic frequency from the middle, and changes so as to have a maximum value and a minimum value. . Therefore, the filter of FIG. 6 has a portion where the cutoff frequency becomes lower as time elapses.
  • a plurality of filters having patterns as shown in FIG. 6 are prepared.
  • the dynamic filter 1f employs a high-pass filter.
  • the dynamic filter 1f may be a band pass filter, a comb filter, a low pass filter, or the like.
  • the dynamic filter 1f may be a notch filter that cuts a specific frequency, a parametric equalizer that increases the level of a specific band, or the like.
  • the output extracted by the dynamic filter 1f may be an envelope output, a crest fact output, or the like.
  • the integration method for creating the integration graph varies depending on the type of filter.
  • Embodiment 2 In the second embodiment, after calculating the level difference between the signals related to the non-defective product and the defective speaker 60, the signal corresponding to the level difference is filtered by a plurality of dynamic filters 1f to determine the inspection filter. It relates to form.
  • the same reference numerals are assigned to the same components as those in the first embodiment, and detailed description thereof is omitted.
  • FIG. 7 is an explanatory diagram showing an outline of the operation of the inspection system according to the second embodiment.
  • a log sweep signal of a sine wave is input to a plurality of non-defective products and a plurality of defective products related to the sample speaker 60.
  • a sine wave log sweep signal is a sine wave signal whose frequency changes exponentially and continuously over time.
  • a signal related to the sound generated by the sample speaker 60 is input to the filter determination device 20.
  • the CPU 21 stores the input signal in the hard disk 24.
  • the CPU 21 reads the signal stored in the hard disk 24 together with the non-defective / defective product identification information related to the sample speaker 60 corresponding to the signal.
  • the CPU 21 calculates a three-dimensional spectrogram of time, frequency, and level from the read signal, and integrates the level at the time from the maximum frequency to the frequency for each frequency at each time.
  • FIG. 7 shows integrated graphs 1m and 2m calculated for non-defective products and defective products, with the horizontal axis representing time, the vertical axis representing frequency, and the height axis representing level.
  • the CPU 21 stores the calculated integrated graphs 1m and 2m in the hard disk 24 together with the good / bad classification.
  • the CPU 21 repeats this integrated graph calculation process for all the signals of the sample speaker 60 stored in the hard disk 24.
  • FIG. 8 is an example of a spectrogram 1 s calculated for a non-defective speaker 60.
  • the horizontal axis indicates time t
  • the vertical axis indicates the frequency f
  • the vertical axis with respect to the paper surface indicates the level. Note that the level in FIG. 8 is expressed by the difference in shading.
  • the signal waveform input from the measurement unit 30 to the filter determination device 20 is shown at the bottom of FIG.
  • FIG. 9 is an explanatory diagram showing a frequency distribution example of the level at a certain point in the spectrogram.
  • the horizontal axis in FIG. 9 is the level
  • the vertical axis is the frequency f.
  • the level at the point (t, f) at each time point / frequency distributed discretely on the spectrogram 1s is expressed by the following equation (1).
  • Level P (i) (t, f) is data of a two-dimensional array.
  • the CPU 21 calculates an integrated value of a level from the maximum frequency f max to each of a plurality of frequency points at each time point, and creates an integrated graph.
  • the integrated value of the level integrated up to a certain frequency point corresponds to the level of the output signal of the high-pass filter whose cutoff frequency is that frequency.
  • the integrated value g (i) (t, f) is expressed by the following equation (2).
  • FIG. 10 is an explanatory diagram showing a frequency distribution example of the integrated value of the level at a certain time point.
  • the horizontal axis in FIG. 10 is the level, and the vertical axis is the frequency.
  • the dotted line indicates the level, and the solid line indicates the integrated value of the level. That is, the black circle on the solid line in FIG. 10 indicates the integrated value of the level from the highest frequency f max to a certain frequency point at a certain time.
  • the integrated value g (i) (t, f) is two-dimensional array data.
  • FIG. 11 is an example of a three-dimensional graph of integrated values g (i) (t, f), that is, an integrated graph.
  • the vertical axis in FIG. 11 is the integrated value, and the axis in the plane direction is time t and frequency f, respectively.
  • the maximum value is extracted for each time point / frequency point (t, f) from the integrated value g (i) (t, f) of the level calculated for each of the non-defective speakers 60.
  • the extracted maximum value G (t, f) is the maximum value of the level integrated value for all the non-defective products, and is the level integrated value 1 m representing the non-defective product.
  • a filter that discriminates the speaker 60 to be inspected into a non-defective product and a defective product needs to be a filter that determines that all the non-defective products prepared as samples are determined to be good products. For that purpose, it is preferable from the viewpoint of safety to adopt the maximum value of the level corresponding to the non-defective product having the sound output characteristic closest to the defective product.
  • the maximum value G (t, f) is expressed by the following equation (3).
  • the integrated value 2m of the level is calculated in the same manner as the non-defective product.
  • the index for distinguishing a plurality of defective products is ii
  • the level at the point (t, f) at each time point / frequency on the spectrogram 2s is expressed by the following equation (4).
  • the ratio over (ii, t, f) of the level integrated value g (ii) (t, f) acquired from a plurality of defective products to the maximum value G (t, f) of the level integrated value acquired from a plurality of non-defective products f) is calculated for each of a plurality of defective products.
  • over (ii, t, f) is the decibel difference.
  • the ratio over (ii, t, f) corresponds to the level difference so that the logarithmic subtraction corresponds to the true division.
  • Over (ii, t, f) is expressed by the following equation (6).
  • FIG. 12 is an explanatory diagram conceptually showing the calculation of the ratio over (ii, t, f).
  • a graph of G (t, f) related to non-defective products is shown in the upper left of FIG. In the upper right of FIG. 12, g (ii) (t, f) related to the ii-th defective product is shown. Calculating the level difference between the non-defective product and the defective product corresponds to calculating the ratio over (ii, t, f).
  • the graph of over (ii, t, f) regarding the ii-th defective product is shown in the lower part of FIG.
  • the CPU21 calculates over (ii, t, f) corresponding to the difference 1d between the level integrated value 1m of the non-defective product and the level integrated value 2m of the plurality of defective products for each of the plurality of defective products.
  • the CPU 21 sequentially cuts the difference graph 1d of the integrated values calculated for each of the plurality of defective products using a plurality of curves of the dynamic filters 1f prepared in advance.
  • the cut surface is a graph of a level difference between the maximum value of the level when a plurality of non-defective products are processed by the dynamic filter 1f and the level of processing the defective products.
  • FIG. 13 is an explanatory diagram showing a curve showing a change over time of a cutoff frequency of a filter for discriminating a defective product from a non-defective product group on the graph of over (ii, t, f) relating to the ii-th defective product in FIG. It is.
  • the cutoff frequency of the filter changes with time t.
  • the change of the cut-off frequency in FIG. 13 is indicated by a bold line on the (t, f) plane of the graph of over (ii, t, f) in the lower part of FIG.
  • a curve indicating a temporal change in the cutoff frequency of the filter for discriminating defective products from the good product group is referred to as a cutoff curve C.
  • the cut-off curve C is indicated by a bold line.
  • FIG. 14 is a cross-sectional view of FIG. 13 cut along a cut-off curve C.
  • the vertical axis in FIG. 14 is the difference between the maximum value G (t, f) of the level integrated value of the non-defective product and the level integrated value g (ii) (t, f) of the ii-th defective product.
  • the horizontal axis in FIG. 14 is the frequency. Since the level integrated value is an integrated value of the level from the highest frequency to the frequency point, FIG. 14 shows a level integrated value that is filtered between a non-defective product and a defective product by a high-pass filter that filters a high frequency band. The level difference concerning is shown.
  • the arrows in FIG. 14 indicate the maximum value of the level difference, that is, the detected value of the defective product (ii).
  • the detected value Over (ii, C) of the defective product (ii) by the cut-off curve C is expressed by the following equation (7).
  • FIG. 7 shows the level difference 11 after the filtering process, with the vertical axis representing the level and the horizontal axis representing the frequency.
  • the CPU 21 calculates a detection value Over (ii, C) based on the cut-off curve C for each of a plurality of defective products.
  • the CPU 21 writes each detected value calculated for each of the plurality of defective products as a list 200 in the hard disk 24 in association with the identification information of the defective products.
  • FIG. 15 is an explanatory diagram conceptually showing a process of extracting the defect detection capability value D (C) from the detected value Over (ii, C) of n defective products.
  • the minimum detection value that is, the defect detection capability value D (C)
  • the minimum detection value can be divided into an integrated graph 1m for all non-defective products and an integrated graph 2m for all defective products by a cut-off curve C corresponding to D (C). It is the amount that can be.
  • the defect detection capability value D (C) is expressed by the following equation (8).
  • the CPU 21 writes the extracted defect detection capability value D (C) as a list 300 in the hard disk 24 in association with the filter identification information.
  • the CPU 21 changes the dynamic filter 1f prepared in advance and repeats the above processing.
  • the defect detection capability value D (C) is written in the list 300 by the number of prepared dynamic filters 1f.
  • the CPU 21 detects the maximum defect detection from each defect detection capability value D (C) written in the list 300.
  • the filter corresponding to the capability value D (C) is determined as the inspection filter.
  • the CPU 21 inspects the speaker 60 to be inspected with the determined inspection filter.
  • the CPU 21 inputs a signal related to the speaker 60 to be inspected.
  • the CPU 21 filters the input signal with the determined filter for inspection.
  • the CPU 21 compares the signal relating to the filtered speaker to be inspected with the inspection reference data prepared in advance.
  • the CPU 21 determines pass / fail of the speaker 60 to be inspected based on the comparison result.
  • the ratio over (ii, t, f) is used as reference data for inspection standard data for determining whether the level of the filtered speaker 60 to be inspected corresponds to a non-defective product or a defective product. Can do.
  • reference data representing non-defective products is obtained from level integrated values of spectrograms 1s in a plurality of non-defective speakers 60.
  • the spectrogram 1s of one good speaker 60 may be used as reference data.
  • the level integrated value at the point (t, f) at each time point / frequency related to one good product is used as reference data.
  • a method for obtaining reference data representing a non-defective product from level integrated values of spectrograms 1s in a plurality of non-defective speakers 60 may be applied to the first embodiment. In such a case, the level at the point (t, f) at each time point / frequency is obtained for each of a plurality of non-defective products, and the maximum level for each point (t, f) is used as reference data.
  • Embodiment 3 relates to a mode in which the filter used for the inspection of the speaker 60 is determined while suppressing the amount of calculation until the maximum defect detection capability value D (C) is obtained.
  • the filter used for the inspection of the speaker 60 is determined while suppressing the amount of calculation until the maximum defect detection capability value D (C) is obtained.
  • the defect detection capability value D (C) is calculated for many patterns prepared in advance, the amount of calculation becomes enormous, and the determination of the filter is a reality. It ’s not right. Therefore, an inspection filter is determined from three types of patterns that are highly practical. Note that determining the filter corresponds to determining the shape of the cut-off curve C.
  • the same reference numerals are assigned to the same components as those in the first and second embodiments, and detailed description thereof is omitted.
  • the dynamic filter 1f for determining the filter for inspection is a high-pass filter.
  • the cut-off frequency of the dynamic filter 1f changes following the frequency change of the fundamental wave or the harmonic wave of the signal based on the sound from the speaker 60 to which the sweep signal is input.
  • One of the parameters of the dynamic filter 1f is a magnification m of the harmonic frequency with respect to the fundamental frequency.
  • FIGS. 16A, 16B, and 16C are explanatory diagrams illustrating pattern examples.
  • the horizontal axis represents time
  • the vertical axis represents frequency.
  • the patterns in FIGS. 16A, 16B, and 16C are referred to as a pattern, b pattern, and c pattern for convenience.
  • the shape of the b pattern is a shape in which the cut-off frequency increases following high-order harmonics, and the cut-off frequency increases following low-order harmonics.
  • the shape of the b pattern is a shape in which the cut-off frequency is constant in an intermediate frequency band in which the higher order harmonics transition to the lower order harmonics.
  • the non-defective speaker 60 has an inspection specification that allows a plurality of harmonics in the vicinity of a frequency of 360 kHz to a time of 360 ms to 480 ms, the shape of the filter that blocks the signal common to the non-defective product and the defective product is b The pattern shape is reasonable.
  • the cut-off frequency follows the harmonic frequency of the order m0 while the fundamental frequency changes from the start frequency to the frequency f0.
  • the cut-off frequency maintains a constant frequency (f0 ⁇ m0) while the frequency of the fundamental wave changes from the frequency f0 to the frequency f1 (> f0) in the intermediate frequency band.
  • the C pattern is a modification of b pattern.
  • the shape of the c pattern is a shape in which the cutoff frequency is not constant in the frequency constant section of the b pattern (section in which the fundamental frequency changes from f0 to f1), but gradually changes.
  • the signal level difference between the non-defective product and the defective product is larger in the c pattern than in the b pattern. Sometimes. Therefore, the shape of the c pattern is considered.
  • the cut-off frequency monotonously increases in the interval in which the fundamental frequency changes from f0 to f1.
  • the shape of the c pattern is a straight line in a section where the frequency of the fundamental wave changes from f0 to f1.
  • f0 500 Hz
  • m0 7th order
  • f1 3 kHz
  • m1 2 order
  • the shape of the cut-off curve C is set to 3 of the a, b and c patterns here. It was made into a kind.
  • the dynamic filter 1f that is changed when the filter coefficient of the dynamic filter 1f is changed in order to reduce the amount of calculation until the maximum value of the defect detection capability value D (C) is calculated.
  • the order m of the harmonics followed by the cut-off frequency is limited to a certain range in advance.
  • the upper limit value and the lower limit value of the limit range of the order m in the harmonic be m max and mmin , respectively.
  • m max and m min are empirical values corresponding to the cut-off curve C that can discriminate between good and defective products.
  • a level difference at which one filter can distinguish the spectrogram 1s from the good product and the spectrogram 2s from the defective product is determined in advance.
  • the defect detection capability values D (C) corresponding to the upper limit value and the lower limit value of the level difference range are defined as D max and D min , respectively.
  • D max and D min are empirical values of level differences for discriminating between good products and defective products. For example, it has been found that when the level difference is larger than 6 to 10 dB, it is possible to discriminate between a good product and a defective product. Therefore, for example, 15 decibels is set as Dmax, and for example, 6 decibels is set as Dmin .
  • the filter determination process As a large flow of the filter determination process, the CPU 21 changes the shape of the filter from the low frequency side to the high frequency side, and determines a filter whose defect detection capability value D (C) exceeds Dmax as a filter for inspection. To do. At that time, the CPU 21 starts the filter determination process from the a pattern.
  • CPU21 performs a filter determination process about b pattern, when D (C) calculated about a pattern does not exceed Dmax .
  • CPU21 performs a filter determination process about c pattern, when D (C) calculated about b pattern does not exceed Dmax .
  • CPU21 receives the D max, D min, m max and m min.
  • CPU21 substitutes Dmin to the variable D0.
  • the variable D0 is a work variable that holds the maximum D (C) at the time of processing each time D (C) calculated for each of the plurality of defective products and the plurality of filters is updated to a larger value.
  • the CPU 21 terminates the filter search at that time. Thereby, shortening of processing time can be aimed at.
  • FIG. 17 is an explanatory diagram showing an image of processing for obtaining D (C) from the a pattern.
  • the CPU 21 extracts the minimum value D (C) from Over (ii, C) of all defective products.
  • the CPU 21 similarly extracts the minimum value D (C) for the pattern in which the order m is incremented by 1.
  • the CPU 21 substitutes the extracted minimum value D (C) into D0.
  • the CPU 21 repeats the same processing until the order m becomes m max .
  • the CPU 21 ends the search if the extracted minimum value D (C) exceeds D max during the process. At that time, the CPU 21 has determined a filter for inspecting the speaker 60 from the a-pattern filter.
  • the CPU 21 searches for the inspection filter for the b pattern. In that case, CPU21 performs the process which restrict
  • FIG. 18 is an explanatory diagram showing the limitation of the calculation range related to the filter.
  • the horizontal axis represents time, and the vertical axis represents frequency.
  • the frequency time change curve of the fundamental wave related to the sound generated by the speaker 60 is displayed by a broken line for reference.
  • the distribution of over (ii, t, f) regarding n defective products is indicated by white and gray areas.
  • the white area is an area where the value of over (ii, t, f) regarding each defective product is larger than D0 obtained by the filter determination in the a pattern.
  • the gray area is an area where the value of over (ii, t, f) for each defective product is equal to or less than D0 obtained by filter determination in the a pattern.
  • the white area corresponds to an area where the level of abnormal sound generated by the speaker 60 is high.
  • the lower part of FIG. 18 shows a region indicating a calculation range related to a filter on the distribution of over (ii, t, f) related to n defective products.
  • the outside of the area indicated by the broken line is an area where the signal level difference between the non-defective product and the defective product is small, and it is not necessary to change the shape of the pattern in order to determine the filter. Therefore, when changing the shape of the pattern, the CPU 21 excludes the outside of the broken line portion from the calculation process for determining the filter.
  • the CPU 21 determines the range of the order m (m L , m H ) that specifies the range in which the calculation is restricted based on the distribution of over (ii, t, f) regarding the defective product, and the range of the cutoff frequency. (Fmax, fmin) is determined. Thereby, the CPU 21 can improve the processing speed.
  • the CPU21 performs the filter determination for a test
  • the shape of the b pattern is determined by each shape of the harmonic part of the order m0, the part of the constant frequency section, and the harmonic part of the order m1 (see FIG. 16B).
  • the CPU 21 sequentially changes these three portions to dynamically change the pattern shape.
  • the harmonic part of order m0 corresponds to the part where the frequency of the fundamental wave changes from the lowest frequency to f0.
  • the portion of the constant frequency section corresponds to a portion where the frequency of the fundamental wave changes from f0 to f1.
  • the harmonic part of the order m1 corresponds to a part where the frequency of the fundamental wave changes from f1 to the highest frequency.
  • FIG. 19 is an explanatory diagram showing an image of the process of changing the shape of the b pattern.
  • the shape of the harmonic part of the order m0 is determined by the order of the harmonic.
  • the shape of the portion of the constant frequency section is determined by the positions of both ends of the line segment indicating the constant frequency.
  • One end position of the line segment is determined by the frequency f0 and the harmonic order m0.
  • the other end position of the line segment is determined by the frequency f1 and the harmonic order m1.
  • the shape of the harmonic part of the order m1 is determined by the order of the harmonic. Therefore, the shape of the b pattern is determined by the order m0, the frequency f0, the frequency f1, and the order m1.
  • the shape of the b pattern is determined by the order m0, the frequency f0, and the order m1.
  • the range in which the order m0 is changed is from mmin + 1 to mmax .
  • Range of changing the order m1 is from m min to m0-1.
  • the range in which f0 is changed is from the lowest frequency generated to fmax obtained as a calculation limit range based on the level difference regarding defective products.
  • the minimum frequency here is a frequency determined in consideration of the characteristics of the speaker 60, and is, for example, 20 Hz.
  • the order of changing each part of the b pattern is the order of the harmonic part of the order m1, the part of the constant frequency section, and the harmonic part of the order m0.
  • CPU21 is orders m0, by fixing the frequency f0, to vary the order m1 to m0-1 from m min.
  • the CPU 21 increases the frequency f0 by a constant value ⁇ f.
  • CPU21 is orders m0, by fixing the frequency f0 which is higher by a predetermined value Delta] f, changing the order m1 to m0-1 from m min.
  • CPU21 is as high as more constant value ⁇ f frequency f0, changing the order m1 to m0-1 from m min.
  • the CPU 21 increments the order m0 by 1 and sets it to mmin + 2. Then, CPU 21 may order m0, by fixing the frequency f0, to vary the order m1 to m0-1 from m min. As described above, the CPU 21 changes the shape of the b pattern by changing the order m1, the frequency f0, and the order m0 within the range sequentially determined in this order.
  • FIG. 20 is a flowchart illustrating an example of a processing procedure for sequentially changing each part of the b pattern.
  • CPU21 is the m min +1 in order m0, the lowest frequency in the frequency f0, substituting m min in order m1 (step S201).
  • CPU 21 is a m min +2 next degree m0, the lowest frequency in the frequency f0, substituting m min in order m1 (step S202).
  • the CPU 21 increments the order m1 (step S203).
  • the CPU 21 adds ⁇ f to f0 (step S206).
  • the CPU 21 increments the order m0 (step S209).
  • the CPU 21 limits the range in which the orders m0 and m1 are changed to the order limit range (m L , m H ) related to the calculation obtained based on the level difference regarding the defective product. Further, the CPU 21 limits the range in which f0 ⁇ m0 changes to the cutoff frequency limit range (fmax, fmin) related to the calculation obtained based on the level difference regarding the defective product.
  • FIG. 21 is an explanatory diagram showing three portions C0, C1, and C01 in the b pattern.
  • C0 is a portion corresponding to a section in which the frequency of the fundamental wave changes from the lowest frequency to f0 (the harmonic portion of the order m0).
  • C1 is a portion (harmonic portion of the order m1) corresponding to a section where the fundamental frequency changes from f1 to the maximum frequency.
  • C01 is a portion corresponding to a section in which the frequency of the fundamental wave changes from f0 to f1 (part of a constant frequency section).
  • the CPU 21 at the stage of the filter determining process in a pattern, when calculating the difference in the pattern corresponding to the harmonic of order m i, C0, C1 respectively corresponding Over0 (ii, f, mi) , Over1 (ii , F, mi) is calculated (see FIG. 12).
  • the CPU 21 stores the calculated Over0 (ii, f, mi) and Over1 (ii, f, mi) in the RAM 23 in advance.
  • FIG. 22A is an explanatory diagram showing a difference value on a pattern.
  • FIG. 22B is an explanatory diagram illustrating a process of calculating Over0 (ii, f, mi) corresponding to C0.
  • FIG. 22C is an explanatory diagram illustrating a process of calculating Over1 (ii, f, mi) corresponding to C1.
  • the time elapses from the left side to the right side.
  • the CPU 21 extracts the maximum value Over0 (ii, f, mi) corresponding to C0 by comparing the difference values in order from the start to the end of the frequency change of the input signal at the stage of the filter determination process in the a pattern. Keep it.
  • the CPU 21 extracts the maximum value Over1 (ii, f, mi) corresponding to C1 by comparing the difference values in order from the end of the frequency change of the input signal to the start at the stage of the filter determination process of the a pattern. Keep it.
  • the CPU 21 obtains Over (ii, Cx) for each of C0, C1, and C01 for each defective product, and sets the obtained maximum value of Over (ii, Cx) to Over (ii, C). .
  • Cx is any one of C0, C1, and C01.
  • FIG. 23 is an explanatory diagram showing an image in which the maximum value is set to Over (ii, C) from Over (ii, Cx) corresponding to C0, C1, and C01.
  • CPU21 substitutes Over0 (ii, f, mi) memorize
  • CPU21 substitutes Over1 (ii, f, mi) memorize
  • the CPU 21 obtains Over (ii, C01) from the calculated data of over (ii, t, F).
  • the CPU 21 sets the maximum value to Over (ii, C) among Over (ii, C0), Over (ii, C1), and Over (ii, C01).
  • the CPU 21 obtains the maximum value Over (ii, C) for all defective products, extracts the minimum value from the obtained maximum value Over (ii, C), and sets it to D (C).
  • the CPU 21 assigns the extracted minimum value D (C) to D0 and updates D0.
  • the CPU 21 changes the pattern shape, for example, by incrementing m1 by 1 for the harmonic portion of the order m1, and repeats the above processing.
  • the CPU 21 ends the filter determination process when the extracted minimum value D (C) exceeds D max during the process.
  • CPU21 performs a filter determination process about c pattern, when all the minimum value D (C) extracted about b pattern does not exceed Dmax . At that time, the CPU 21 executes the filter determination process for the c pattern after executing the restriction of the calculation range shown in FIG. 18 based on D0 which is the maximum value of the defect detection capability value calculated by the processes so far. .
  • FIG. 24 is an explanatory diagram showing an image of the process of changing the shape of the c pattern.
  • the shape of the c pattern is determined by each shape of the harmonic part of order m0, the part where the frequency gradually increases, and the harmonic part of order m1 (see FIG. 16C).
  • the harmonic part of order m0, the part where the frequency gradually increases, and the harmonic part of order m1 correspond to C0, C01 and C1 in FIG. 21, respectively. Therefore, C0, C01, and C1 are attached to the time axis of FIG. 24 for reference.
  • the CPU 21 sequentially changes the order m0, the frequency f0, and the order m1 to change the shape of the c pattern.
  • the CPU 21 since the frequency f1 is not fixed to one, the CPU 21 sequentially changes the frequency f1 to change the shape of the c pattern.
  • the range in which the frequency f1 is changed is from the lowest frequency generated to fmax.
  • the minimum frequency here is a frequency determined in consideration of the characteristics of the speaker 60, and is, for example, 20 Hz.
  • the CPU 21 changes the shape of the c pattern by changing the frequency f1, the order m1, the frequency f0, and the order m0 within a range sequentially determined in this order.
  • the change of the frequency f1 here corresponds to changing the position of the other end point of the portion where the frequency gradually increases (corresponding to C01) to the high frequency side along the harmonic of the order m1.
  • the CPU 21 changes the position of the other end point so that f1 ⁇ m1 changes from a value obtained by adding a constant value to f0 ⁇ m0 to fmax.
  • the CPU 21 obtains the maximum value Over (ii, C) for all defective products, extracts the minimum value from the obtained maximum value Over (ii, C), and sets it to D (C).
  • the CPU 21 assigns the extracted minimum value D (C) to D0 and updates D0.
  • the CPU 21 changes the shape of the c pattern and repeats the above processing. However, the CPU 21 ends the filter determination process when the extracted minimum value D (C) exceeds D max during the process.
  • the CPU 21 determines that the maximum value D0 of the defect detection capability value exceeds D min in the processing so far.
  • the filter corresponding to the maximum value D0 is determined as the inspection filter.
  • CPU21 receives the D max, D min, m max and m min (step S301).
  • CPU21 substitutes Dmin to the variable D0 (step S302).
  • CPU21 substitutes 1 to the variable P (step S303).
  • the variable P is a flag variable for determining the pattern.
  • the variable P is assigned 1 for the a pattern, 2 for the b pattern, and 3 for the c pattern.
  • step S307: YES YES
  • step S308 calculates D (C), when it determines with having calculated Over (ii, C) about all inferior goods (step S307: YES) (step S308).
  • the CPU 21 determines whether or not D (C)> D0 (step S309). If the CPU 21 determines that D (C)> D0 (step S309: YES), it substitutes D (C) for D0 (step S310).
  • the pattern corresponding to D (C) assigned to D0 in step S310 is a provisional pattern corresponding to the filter candidate determined for inspection. If the CPU 21 determines that D (C)> D0 is not satisfied (step S309: NO), the process proceeds to step S313.
  • the CPU 21 determines whether or not D (C) ⁇ Dmax (step S311). If the CPU 21 determines that D (C) ⁇ Dmax (step S311: YES), the CPU 21 outputs the filter coefficient of the filter corresponding to D (C) (step S312), and ends the process. When determining that D (C) ⁇ Dmax is not satisfied (step S311: NO), the CPU 21 determines whether or not processing has been executed for all filters included in the filter pattern being processed (step S313). If the CPU 21 determines that the process has not been executed for all the filters (step S313: NO), the CPU 21 changes the pattern shape according to the value of P (step S314), and returns the process to step S306.
  • the CPU 21 changes the filter to be determined from the a pattern to the b pattern and further from the b pattern to the c pattern when D (C) calculated for the pattern does not exceed D max . .
  • D (C) calculated for the a pattern does not exceed D max
  • the CPU 21 does not shift to the b pattern, determines the filter corresponding to D0 at that time as the inspection filter, and performs the filter determination process. May be terminated.
  • D (C) calculated for the b pattern does not exceed D max
  • the CPU 21 does not shift to the c pattern, determines the filter corresponding to D0 at that time as the inspection filter, and determines the filter. Processing may be terminated. Thereby, shortening of processing time can be aimed at.
  • the CPU 21 limits the calculation range related to the filter.
  • the CPU 21 may change the target of the filter determination process from the a pattern to the b pattern and from the b pattern to the c pattern without limiting the calculation range regarding the filter.
  • inspecting the said apparatus can be determined for every specification of the apparatus which generate
  • the CPU 21 can determine the filter that maximizes the signal level difference between the non-defective product and the defective product of the speaker 60 prepared in advance by repeating the process of dynamically changing the filter.
  • the filter determination device 20 it is possible to reduce the amount of calculation until the maximum defect detection capability value D (C) is obtained.
  • the filter determination device 20 executes filter determination within the accepted ranges of D max and D min . Therefore, the filter determination device 20 can reduce the number of D (C) to be calculated. Further, when the calculated D (C) exceeds D max , the filter determination device 20 ends the determination of the filter at that time, and selects a filter corresponding to D (C) exceeding D max as a test filter. And decide. Thereby, the filter determination apparatus 20 can aim at shortening of processing time.
  • the filter determination device 20 executes filter determination processing in three patterns of a pattern, b pattern, and c pattern.
  • the filter determination device 20 changes the shape of the pattern within the limited filter pattern range and within the accepted order ranges of m max and mmin . Such limiting conditions also contribute to a reduction in the amount of calculation.
  • the filter determination device 20 limits the calculation range based on D0 in determining the filters for the b pattern and the c pattern. The limitation of the calculation range can also reduce the calculation amount.

Landscapes

  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Otolaryngology (AREA)
  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Acoustics & Sound (AREA)
  • Signal Processing (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

 音又は振動を発生する装置の仕様毎に、最適な検査用のフィルタを複数のフィルタから決定して当該装置を検査することができる検査方法、検査装置、検査方法に使用するフィルタを決定することができるフィルタ決定方法、プログラム及びフィルタ決定装置を提供する。更に、特性を時間経過につれて変化させ得るフィルタを用いて、周波数が時間経過につれて変化する音又は振動を発生する装置を検査することができる検査方法及び検査装置を提供する。 検査対象の装置が発生した、周波数が時間経過につれて変化する音又は振動に係る信号を、特性を時間経過につれて変化させ得るフィルタを用いて処理して、前記装置を検査する。

Description

検査装置
 本願は、音又は振動を発生する装置を検査する検査方法、検査装置、検査方法に使用するフィルタを決定するフィルタ決定方法、プログラム及びフィルタ決定装置に関する。
 装置が発生する音又は振動から当該装置の異常を検出する技術がある。例えば、特許文献1には、スピーカが発生するびりつき音を検出して、当該スピーカの良否を判定する検出装置が記載されている。特許文献1に係る検出装置は、スピーカが発生した音に基づく信号から、予め最適化されたフィルタを用いてびりつき音の信号成分を抽出することにより、びりつき音を検出する。
特開平6-315197号公報
 しかしながら、特許文献1には、「予め最適化されたフィルタ」に関して具体的な説明が一切記載されていない。
 異常な装置が発生する音又は振動は、装置の仕様毎に異なる。そのため、様々な装置の良否を正確に判定するためには、仕様毎に最適化されたフィルタが必要である。しかし、装置の仕様毎にフィルタを予め最適化することは実用上困難である。
 本願は、かかる事情に鑑みてなされたものである。本願の目的は、音又は振動を発生する装置の仕様毎に、最適な検査用のフィルタを複数のフィルタから決定して当該装置を検査することができる検査方法、検査装置、検査方法に使用するフィルタを決定することができるフィルタ決定方法、プログラム及びフィルタ決定装置を提供することにある。
 更に、本願の目的は、特性を時間経過につれて変化させ得るフィルタを用いて、周波数が時間経過につれて変化する音又は振動を発生する装置を検査することができる検査方法及び検査装置を提供することにある。
 本願に係る検査方法は、検査対象の装置が発生した、周波数が時間経過につれて変化する音又は振動に係る信号を、特性を時間経過につれて変化させ得るフィルタを用いて処理して、前記装置を検査する。
 本願に係る検査方法は、検査対象の装置を検査するに先立ち、予め用意された装置の一又は複数の良品及び一又は複数の不良品夫々が発生した音又は振動に係る信号を複数のフィルタ各々で処理した複数の信号に基づいて、検査用のフィルタを決定しておき、決定した検査用のフィルタを用いて、検査対象の装置が発生した音又は振動に係る信号を処理し、処理した信号に基づいて、検査対象の装置の良否を判定することを特徴とする。
 本願に係る検査方法は、前記良品及び不良品夫々が発生した音又は振動に係る信号のレベル差を算出し、算出したレベル差に対応する信号を複数のフィルタ夫々で処理し、複数のフィルタ夫々で処理した信号に対応するレベル差のうち、最大のレベル差に対応するフィルタを検査用のフィルタと決定することを特徴とする。
 本願に係る検査方法は、前記良品及び複数の不良品夫々が発生した音又は振動に係る信号のレベル差を該複数の不良品毎に算出し、算出したレベル差に対応する信号を複数のフィルタ夫々で処理し、処理した信号に対応するレベル差の最大値を前記複数の不良品毎に算出し、複数のフィルタ夫々につき、前記複数の不良品毎に算出したレベル差の各最大値のうち、レベル差の最小値を選択し、複数のフィルタ夫々につき選択したレベル差の各最小値のうち、最大のレベル差に対応するフィルタを検査用のフィルタと決定することを特徴とする。
 本願に係る検査方法は、予め用意された複数の良品に係る信号から夫々算出された時間、周波数及びレベルの情報に基づく信号から、予め定められた間隔で画定された時間及び周波数毎に最大のレベルを各々抽出し、抽出した最大の各レベルを前記良品が発生した音又は振動に係る信号のレベルに設定し、設定した前記良品に係る信号のレベルに基づいて、レベル差の最大値を前記複数の不良品毎に算出することを特徴とする。
 本願に係る検査方法は、選択したレベル差の最小値が所定の値を超えた場合、レベル差の最小値の選択を終了し、前記所定の値を超えたレベル差の最小値に対応するフィルタを検査用のフィルタと決定することを特徴とする。
 本願に係る検査装置は、検査対象の装置が発生した、周波数が時間経過につれて変化する音又は振動に係る信号を、特性を時間経過につれて変化させ得るフィルタを用いて処理して、前記装置を検査する。
 本願に係る検査装置は、複数のフィルタを出力するフィルタ出力部と、予め用意された装置の一又は複数の良品及び一又は複数の不良品夫々が発生した音又は振動に係る信号に基づいて、前記フィルタ出力部が出力する複数のフィルタから検査用のフィルタを決定する決定部と、該決定部が決定した検査用のフィルタを用いて、検査対象の装置が発生した音又は振動に係る信号を処理する処理部と、該処理部が処理した信号に基づいて、検査対象の装置の良否を判定する判定部とを備えることを特徴とする。
 本願に係る検査装置は、前記良品に係る信号を記憶する良品信号記憶部と、該良品信号記憶部が記憶する良品に係る信号及び予め用意された前記不良品に係る信号のレベル差を算出するレベル差算出部と、前記フィルタ出力部が出力する複数のフィルタ夫々で、前記レベル差算出部が算出したレベル差に対応する信号を処理する信号処理部とを備え、前記決定部は、前記信号処理部が処理した信号に対応するレベル差のうち、最大のレベル差に対応するフィルタを検査用のフィルタと決定するようにしてあり、前記処理部が処理した信号及び所定の信号を比較する信号比較部を備え、前記判定部は前記信号比較部が比較した結果に基づいて、検査対象の装置の良否を判定するようにしてあることを特徴とする。
 本願に係る検査装置は、前記レベル差算出部は、前記良品及び複数の不良品夫々が発生した音又は振動に係る信号のレベル差を該複数の不良品毎に算出するようにしてあり、前記信号処理部は前記フィルタ出力部が出力する複数のフィルタ夫々で、前記レベル差算出部が前記複数の不良品毎に算出したレベル差に対応する信号を処理するようにしてあり、前記信号処理部が処理した信号に対応するレベル差の最大値を前記複数の不良品毎に算出する算出部と、前記フィルタ出力部が出力する複数のフィルタ夫々につき、前記算出部が前記複数の不良品毎に算出したレベル差の各最大値のうち、レベル差の最小値を選択する選択部とを備え、前記決定部は前記選択部が複数のフィルタ夫々につき選択したレベル差の各最小値のうち、最大のレベル差に対応するフィルタを検査用のフィルタと決定するようにしてあることを特徴とする。
 本願に係る検査装置は、前記選択部は、選択したレベル差の最小値が所定の値を超えた場合、レベル差の最小値の選択を終了するようにしてあり、前記決定部は前記所定の値を超えたレベル差の最小値に対応するフィルタを検査用のフィルタと決定するようにしてあることを特徴とする。
 本願に係るフィルタ決定方法は、周波数が時間経過につれて変化する音又は振動を発生する検査対象の装置を検査するために、特性を時間経過につれて変化させ得る複数のフィルタから検査用のフィルタを決定するフィルタ決定方法であって、予め用意された装置の良品及び一又は複数の不良品夫々が発生した音又は振動に係る信号のレベル差を該複数の不良品毎に算出し、算出したレベル差に対応する信号を複数のフィルタ夫々で処理し、処理した信号に対応するレベル差の最大値を前記一又は複数の不良品毎に算出し、複数のフィルタ夫々につき、前記一又は複数の不良品毎に算出したレベル差の各最大値のうち、レベル差の最小値を選択し、複数のフィルタ夫々につき選択したレベル差の各最小値のうち、最大のレベル差に対応するフィルタを検査用のフィルタと決定することを特徴とする。
 本願に係るプログラムは、周波数が時間経過につれて変化する音又は振動に係る信号に基づいて、特性を時間経過につれて変化させ得る複数のフィルタから一のフィルタを決定する処理をコンピュータに実行させるプログラムであって、音又は振動に係る一の信号と、一又は複数の他の信号とのレベル差を該一又は複数の他の信号毎に算出し、特性を時間経過につれて変化させ得るフィルタを変更し、変更した複数のフィルタ夫々で、算出したレベル差に対応する信号を処理し、処理した信号に対応するレベル差の最大値を前記一又は複数の他の信号毎に算出し、複数のフィルタ夫々につき、前記一又は複数の他の信号毎に算出したレベル差の各最大値のうち、レベル差の最小値を選択し、複数のフィルタ夫々につき選択したレベル差の各最小値のうち、最大のレベル差に対応するフィルタを一のフィルタと決定する処理をコンピュータに実行させることを特徴とする。
 本願に係るフィルタ決定装置は、周波数が時間経過につれて変化する音又は振動を発生する検査対象の装置を検査する検査装置のために、特性を時間経過につれて変化させ得る複数のフィルタから検査用のフィルタを決定するフィルタ決定装置であって、複数のフィルタを出力するフィルタ出力部と、予め用意された装置の良品及び一又は複数の不良品夫々が発生した音又は振動に係る信号のレベル差を該複数の不良品毎に算出するレベル差算出部と、前記フィルタ出力部が出力する複数のフィルタ夫々で、前記レベル差算出部が算出したレベル差に対応する信号を処理する信号処理部と、該信号処理部が処理した信号に対応するレベル差の最大値を前記一又は複数の不良品毎に算出する算出部と、前記フィルタ出力部が出力する複数のフィルタ夫々につき、前記算出部が前記一又は複数の不良品毎に算出したレベル差の各最大値のうち、レベル差の最小値を選択する選択部と、該選択部が複数のフィルタ夫々につき選択したレベル差の各最小値のうち、最大のレベル差に対応するフィルタを検査用のフィルタと決定する決定部とを備えることを特徴とする。
 本願発明に係る検査方法、検査装置、フィルタ決定方法、プログラム及びフィルタ決定装置によれば、音又は振動を発生する装置の仕様毎に、最適な検査用のフィルタを複数のフィルタから決定して当該装置を検査することができる。更に、本願発明に係る検査方法及び検査装置によれば、特性を時間経過につれて変化させ得るフィルタを用いるので、周波数が時間経過につれて変化する音又は振動を発生する装置を音又は振動の変化に適用させて検査することができる。
本願に係る検査装置を備える検査システムのハードウェア構成例を示すブロック図である。 検査システムの動作の概要を示す説明図である。 フィルタ決定処理の概要手順の一例を示すフローチャートである。 検査処理の手順の一例を示すフローチャートである。 カットオフ周波数が高くなる変化の態様が時間経過につれて変化するフィルタの一例を示す説明図である。 時間経過につれてカットオフ周波数が低くなる部分を有するフィルタの一例を示す説明図である。 実施の形態2に係る検査システムの動作の概要を示す説明図である。 良品のスピーカについて算出されたスペクトログラムの一例である。 ある時点におけるレベルの周波数分布例を示す説明図である。 ある時点におけるレベルの積算値の周波数分布例を示す説明図である。 積算値g(i)(t,f)の3次元グラフ例である。 比over(ii,t,f)の計算を概念的に示す説明図である。 図12におけるii番目の不良品に係るover(ii,t,f)のグラフ上で、不良品を良品群から判別するフィルタのカットオフ周波数の時間変化を示すカーブを示す説明図である。 図13をカットオフカーブCに沿って切断した断面図である。 n個の不良品の検出値Over(ii,C)から、不良検出能力値D(C)を抽出する処理を概念的に示す説明図である。 パターン例を示す説明図である。 パターン例を示す説明図である。 パターン例を示す説明図である。 aパターンからD(C)を求める処理のイメージを示す説明図である。 フィルタに関する演算範囲の制限を示す説明図である。 bパターンの形状を変更する処理のイメージを示す説明図である。 bパターンの各部分を順次変更する処理の手順の一例を示すフローチャートである。 bパターンにおける3つの部分C0、C1、C01を示す説明図である。 パターン上の差分値を示す説明図である。 C0に対応するOver0(ii,f,mi)を算出する処理を示す説明図である。 C1に対応するOver1(ii,f,mi)を算出する処理を示す説明図である。 C0、C1及びC01に対応するOver(ii,Cx)から、最大値をOver(ii,C)に設定するイメージを示す説明図である。 cパターンの形状を変更する処理のイメージを示す説明図である。 フィルタ決定処理の手順の一例を示すフローチャートである。 フィルタ決定処理の手順の一例を示すフローチャートである。
 本実施の形態に係る検査装置を、図面を参照して説明する。本実施の形態に係る検査装置は、音又は振動を発生する装置を検査して、当該装置の良否判定を実行する。以下、スピーカを検査対象の装置の例に挙げて、本実施の形態に係る検査装置を説明する。
 実施の形態1
 図1は、本願に係る検査装置を備える検査システムのハードウェア構成例を示すブロック図である。検査システムは、フィルタ決定装置20、測定部30、マイクロホン40及び無響箱50を含む。検査対象のスピーカ60を検査する検査装置は、フィルタ決定装置20及び測定部30から構成されている。検査装置は、フィルタ決定装置20を要部とする。
 フィルタ決定装置20は、パーソナルコンピュータを用いてなる。フィルタ決定装置20は、測定部30に低周波数から高周波数に変化する正弦波のスイープ信号に対応するデータを出力する。フィルタ決定装置20は、測定部30からスピーカ60が発生した音に対応する信号を入力する。フィルタ決定装置20は、測定部30から入力した信号に基づいて、各種処理を実行する。
 フィルタ決定装置20は、CPU21、ROM22、RAM23、ハードディスク24、ディスクドライブ25、表示部26、操作部27、通信部28及びI/O部29を含む。
 CPU21は、プログラム2Pに記述された命令セットを実行するプロセッサである。ROM22は、フィルタ決定装置20の起動時にCPU21が実行するBIOS等を記憶している。
 RAM23は、主記憶装置である。ハードディスク24は、補助記憶装置である。ハードディスク24は、測定部30が増幅したスピーカ60からの音の信号を記憶する。
 ディスクドライブ25は、外部記憶媒体の光ディスク2dから情報を読み出し、光ディスク2dに情報を記録する。表示部26は、画像を表示する表示装置である。操作部27は、ユーザから各種の入力を受け付ける入力デバイスである。通信部28は、LAN100とのインタフェースである。I/O部29は、測定部30とのインタフェースである。
 なお、CPU21は、ディスクドライブ25を介して、プログラム2Pを光ディスク2dなどの可搬記憶媒体から読み込んでもよい。CPU21は、通信部28を介して、プログラム2Pを他の情報処理装置又は記憶装置から読み込んでもよい。さらに、プログラム2Pを記憶したフラッシュメモリ等の半導体メモリ2mが、フィルタ決定装置20内に実装されていてもよい。
 測定部30は、スピーカ60にスイープ信号を入力する。
 測定部30は、D/A変換器31、パワーアンプ32、計測増幅器33及びA/D変換器34を含む。D/A変換器31は、フィルタ決定装置20から入力したデジタルのスイープ信号をアナログ信号に変換し、変換したアナログ信号をパワーアンプ32に出力する。パワーアンプ32は、D/A変換器31から変換されたアナログ信号を増幅する。パワーアンプ32は、スイープ信号をスピーカ60に出力する。
 無響箱50は、高い遮音性と、吸音特性を持つ。スピーカ60と対向する無響箱50の上壁には、スピーカ60が発生する音を内部に取り入れる開口が設けられている。マイクロホン40は、無響箱50の内部にスピーカ60に正対するように設置されている。
 計測増幅器33は、マイクロホン40から出力されるアナログの信号のレベルを増幅させる。A/D変換器34は、計測増幅器33が増幅した信号をデジタル信号に変換する。A/D変換器34は、変換したデジタル信号をフィルタ決定装置20へ出力する。
 図2は、検査システムの動作の概要を示す説明図である。予め、複数のスピーカ60について聴感テストを実施し、1つの良品と、複数の不良品とを、検査用のフィルタを決定するための処理対象のサンプルとして用意する。用意した1つのスピーカ60を無響箱50に取り付ける。フィルタ決定装置20は、測定部30のD/A変換器31へ、低い周波数(例えば10Hz)から高い周波数(例えば20kHz)に一定速度で変化する正弦波のスイープ信号に対応するデータを入力する。スイープ時間は、例えば1.0秒である。スピーカ60は、D/A変換器31及びパワーアンプ32を介して、フィルタ決定装置20から入力されたデータに対応するスイープ信号に対応する音を発生する。
 マイクロホン40は、サンプル用のスピーカ60が発生した音を受音する。マイクロホン40は、計測増幅器33及びA/D変換器34を介して、信号をフィルタ決定装置20に出力する。CPU21は、測定部30から信号を入力する。CPU21は、入力した信号をサンプル用のスピーカ60に係る良品/不良品の識別情報と対応付けてハードディスク24に記憶する。検査システムは、同様の記憶処理を、用意した全てのサンプル用のスピーカ60について繰り返す。
 CPU21は、ハードディスク24に記憶した信号を、当該信号に対応するサンプル用のスピーカ60に係る良品/不良品の識別情報と共に読み出す。CPU21は、読み出した信号から、時間、周波数、レベルの3次元からなるスペクトログラムを算出し、各時間における周波数毎に最大周波数から該周波数までの該時間におけるレベルの積算をする。図8はスペクトログラムの例、図9は、ある時点におけるスペクトログラムの断面、つまりレベルの周波数分布例を示し、図10はその時点における最大周波数から該周波数までのレベルの積算をしたものを示している。すべての時間での積算グラフは図11のようになる。図8~図11の詳細は、実施の形態2で説明する。CPU21は、算出した積算グラフを良/不良の別と共にハードディスク24に記憶する。CPU21は、この積算グラフ算出処理を、ハードディスク24に記憶したサンプル用のスピーカ60の信号全てについて繰り返す。
 一般的に不良品のスピーカ60は、入力信号に対して主に高域で異常音を発生する。異常音は、入力信号の基本波、2次、3次…の高調波に比べて、レベルが小さい場合が多いく低次の高調波を取り除く事が必要になるため、原則として高次高調波を通過させるハイパスフィルタを用いる。あるサンプル用のスピーカ60の信号をある時間にある周波数をカットオフ周波数に持つハイパスフィルタを通過させたレベルは、該スピーカ60の前記積算グラフの該時間の該周波数でのレベルになる。
 本発明では、CPU21は、検査対象のスピーカ60を良品と不良品と明確に判別するため、良品と不良品との間でレベル差が大きくなるように、音の信号を処理するフィルタを決定する。
 CPU21は、スイープ信号に応じたスピーカ60からの信号にハイパスフィルタのカットオフ周波数を追従させるために、ハイパスフィルタのカットオフ周波数を動的に変化させる。つまり、CPU21が使用するフィルタは動的フィルタ1fといえる。動的フィルタ1fは、カットオフ周波数の時間変化で表される。図2において、動的フィルタ1fを示すグラフの横軸は時間、縦軸は周波数である。すなわち、動的フィルタ1fは、時間及び周波数を要素とする2次元上の曲線で表現できる。また、あるサンプル用のスピーカ60の信号を該動的フィルタ1fで処理した結果のレベルグラフは、該スピーカ60の前記積算グラフを該動的フィルタ1fの曲線で切断した切断面となる。
 CPU21は、動的フィルタ1fを用いて、サンプル用の良品/不良品の双方に共通する成分を処理する。CPU21は、サンプル用の良品及び不良品に係る処理後のレベル差に基づいて、検査対象のスピーカ60を検査するためのフィルタを決定する。CPU21は、決定した動的フィルタ1fを用いて、検査対象のスピーカ60を検査する。動的フィルタ1fはCPU21で作成されるデジタルフィルタであるので、動的フィルタ1fのフィルタ特性はフィルタパラメータによって決定される。
 検査対象のスピーカ60を検査するためのフィルタの決定処理について説明する。
 どのような特性変化のフィルタも許容されるのであれば、夫々の不良品の差分グラフの最大値の点の全てを通る曲線でフィルタを作ればよいが、実際には様々な制約があるのでそのようなフィルタを実現することは困難となる。フィルタの形状は、実際の検査ラインで利用可能な一定のパターンに限定して考える必要がある。一定のパターンのフィルタを複数ハードディスク24に記憶しておく。本実施例でのパターンは、時間経過につれてカットオフ周波数が低周波数から高周波数へ連続的に変化する。また、パターンには、基本波又は低次数から高次数までの複数の高調波の周波数変化に夫々追従して変化する複数のパターンがある。
 なお、CPU21は、一定のパターンのフィルタをリアルタイムで複数生成し、生成した一定のパターンのフィルタ夫々を、以下の処理のために出力してもよい。
 CPU21は、用意されたフィルタから1つのフィルタを選択する。CPU21は、選択したフィルタでサンプル用の良品の積算グラフに対してフィルタ曲線で切断した切断面を取得する。該切断面は、該フィルタで良品をフィルタ処理したレベルグラフになる。
 CPU21は、選択したフィルタでサンプル用の複数の不良品の積算グラフをフィルタ曲線で切断した切断面を取得する。該切断面は、該フィルタで該不良品を処理したレベルグラフになる。図2に、処理後の良品及び不良品のレベルグラフ1g、2gを、横軸を時間、縦軸をレベルとして、模式的に示す。
 なお、CPU21は、フィルタ処理をした良品のレベルグラフ1gを、ハードディスク24に記憶しておいてもよい。
 CPU21は、処理した良品のレベルグラフと、1つの不良品のレベルグラフとの差分グラフの最大値を該不良品の該フィルタによる検出値として算出する。CPU21は、算出した検出値を不良品の識別情報と対応付けて、ハードディスク24に記憶する。CPU21は、全てのサンプル用の不良品について、検出値を算出する処理を繰り返す。図2には、ハードディスク24に記憶された検出値のリスト200が示されている。検出値が大きいほど、不良品は良品から懸け離れていることになる。
 複数の不良品の間には、音出力特性にバラツキがある。検査対象のスピーカ60を良品と不良品とを判別するフィルタは、サンプルとして用意した不良品の全てが不良品に判定されるフィルタである必要がある。そのためには、最も良品寄りの音出力特性を有する不良品に対応するレベル差を採用することが安全側に立ち、好ましい。
 そこで、CPU21は、リスト200にある、全ての不良品夫々について算出した検出値から、最小値を該フィルタの不良検出能力値として抽出する。CPU21は、抽出した不良検出能力値を、選択したフィルタと対応付けてハードディスク24にリスト300として記憶する。CPU21は、各フィルタについて、リスト300への不良検出能力値の記憶を繰り返す。
 CPU21は、予め用意したフィルタの全てを用いて不良検出能力値を計算し尽くした場合、ハードディスク24にリスト300として記憶した複数の不良検出能力値から、最大の値に対応するフィルタを、検査用のフィルタに決定する。あるいは、CPU21は、フィルタを変更して上記の処理途中で例えば利用者が良品と不良品の分離には充分と考える不良検出能力値以上の不良検出能力値を持つフィルタが見つかった場合、そのフィルタを、検査用のフィルタに決定してもよい。これにより、検査用のフィルタを決定する処理時間の短縮を図ることができる。
 図3は、フィルタ決定処理の概要手順の一例を示すフローチャートである。図3の手順では、複数の不良品に係る積算グラフの算出は予め終了しているものとしている。
 CPU21は、未処理のフィルタを選択する(ステップS101)。CPU21は、未処理の不良品の積算グラフを選択する(ステップS102)。CPU21は、良品及び不良品に係る積算グラフを、選択したフィルタで処理する(ステップS103)。CPU21は、処理した良品のレベルグラフと、不良品のレベルグラフとの差分グラフの最大値を該不良品の該フィルタによる検出値として算出する(ステップS104)。
 CPU21は、算出した検出値を不良品の識別情報と対応付けて、ハードディスク24内のリスト200に記憶する(ステップS105)。CPU21は、全ての不良品の積算グラフについて処理が終了したか否かを判定する(ステップS106)。CPU21は、全ての不良品の積算グラフについて処理が終了していないと判定した場合(ステップS106:NO)、ステップS102へ処理を戻す。CPU21は、全ての不良品の積算グラフについて処理が終了したと判定した場合(ステップS106:YES)、ハードディスク24内のリスト200に記憶した複数の検出値から最小値を該フィルタの不良検出能力値として抽出する(ステップS107)。CPU21は、抽出した不良検出能力値を、選択したフィルタと対応付けて、ハードディスク24内のリスト300に記憶する(ステップS108)。
 CPU21は、全てのフィルタについて処理が終了したか否かを判定する(ステップS109)。CPU21は、全てのフィルタについて処理が終了していないと判定した場合(ステップS109:NO)、ステップS101に処理を戻す。CPU21は、全てのフィルタについて処理が終了したと判定した場合(ステップS109:YES)、ハードディスク24内のリスト300に記憶した不良検出能力値から最大値を抽出する(ステップS110)。CPU21は、抽出した最大値に対応するフィルタを検査用のフィルタに決定し(ステップS111)、処理を終了する。
 次に、決定したフィルタを用いて、検査対象のスピーカ60を検査する処理について説明する。
 検査員は、検査対象のスピーカ60を、無響箱50の上壁に取り付ける。検査員は、フィルタ決定装置20を操作して、検査対象のスピーカ60にスイープ信号に対応する音を発生させる。マイクロホン40は、検査対象のスピーカ60からの音を受音し、受音した音の信号を、測定部30を介して、フィルタ決定装置20へ出力する。フィルタ決定装置20は、入力した信号及び検査用のフィルタに基づいて、検査対象のスピーカ60を検査する。フィルタ決定装置20は、検査対象のスピーカ60を検査した結果を、例えば表示部26に表示する。検査員は、検査対象のスピーカ60を無響箱50の上壁から取り外す。検査員は、表示部26に表示された結果に応じて、検査対象のスピーカ60を良品/不良品に分別する。以下、検査員は、上記の作業を繰り返す。
 なお、検査員が行なう作業は、自動化してもよいことは勿論である。
 図4は、検査処理の手順の一例を示すフローチャートである。
 CPU21は、測定部30へスイープ信号に対応するデータを出力する(ステップS1)。CPU21は、測定部30から検査対象のスピーカ60が発生した音の信号を測定部30から入力する(ステップS2)。CPU21は、入力した信号を、フィルタ決定処理で決定した検査用のフィルタでフィルタ処理する(ステップS3)。CPU21は、フィルタ処理した信号を、予め用意した検査基準データと比較する(ステップS4)。CPU21は、比較した結果に基づいて、検査対象のスピーカ60が良品か否かを判定する(ステップS5)。
 CPU21は、判定結果を表示部26に出力する(ステップS6)。CPU21は、上記の検査処理を繰り返すか否かの指示を操作部27から受け付けたか否かを判定する(ステップS7)。CPU21は、上記の検査処理を繰り返す指示を操作部27から受け付けたと判定した場合(ステップS7:YES)、ステップS1へ処理を戻す。CPU21は、上記の検査処理を繰り返す指示を操作部27から受け付けなかった場合(ステップS7:NO)、処理を終了する。なお、CPU21が検査処理を繰り返す指示を操作部27から受け付けなかった場合とは、CPU21が検査処理を停止する指示を受け付けた場合である。
 本実施の形態では、フィルタ決定装置20は、スピーカ60に低い周波数から高い周波数に変化するスイープ信号を入力している。しかし、フィルタ決定装置20は、スピーカ60に高い周波数から低い周波数に変化する逆スイープ信号を入力してもよい。
 本実施の形態では、フィルタ決定装置20は、スピーカ60にスイープ信号を入力している。しかし、フィルタ決定装置20は、固定周波数の正弦波であり、周波数が切り替わるステップ正弦波をスピーカ60に入力してもよい。かかる場合、周波数の時間変化は、単調増加又は単調減少でなくてもよい。あるいは、フィルタ決定装置20は、複数の周波数を有する合成波をスピーカ60に入力してもよい。
 本実施の形態に係る検査対象の装置は、外部から入力した信号に基づいて音を発生するスピーカ60であった。しかし、検査対象の装置は、時間経過につれて周波数が変化する音又は振動を自ら発生する装置でもよい。当該検査対象の装置は、モータ、ギヤ・ベアリング等の回転部材、エレベータの乗車かご、コンプレッサ等である。
 検査用のフィルタを決定するために、予め用意しておくパターンは、カットオフ周波数が高くなる変化の態様が時間経過につれて変化するものを含んでもよい。
 図5は、カットオフ周波数が高くなる変化の態様が時間経過につれて変化するフィルタの一例を示す説明図である。なお、横軸の時間は、ミリ秒の対数である。カットオフ周波数は、最初、高次数の高調波の周波数変化に追従して変化し、途中から約一定周波数を維持する。その後、カットオフ周波数は、低次数の高調波側へ遷移するが、高調波の周波数変化からはずれている。図5に示すようなパターンのフィルタを複数用意しておく。
 検査用のフィルタを決定するために、予め用意しておくパターンは、カットオフ周波数が時間経過につれて単調に高くなるだけでなく、極大値及び極小値を有するように変化する部分を有するものを含んでもよい。
 図6は、時間経過につれてカットオフ周波数が低くなる部分を有するフィルタの一例を示す説明図である。なお、横軸の時間は、ミリ秒の対数である。図6のフィルタの場合、カットオフ周波数は、最初、高次数の高調波の周波数変化に追従して変化し、途中から高調波の周波数より低くなり、極大値及び極小値を有するように変化する。そのため、図6フィルタは、カットオフ周波数が時間経過につれて低くなる部分を有している。図6に示すようなパターンのフィルタを複数用意しておく。
 本実施の形態では、動的フィルタ1fはハイパスフィルタを採用している。しかし、動的フィルタ1fは、バンドパスフィルタ、櫛形フィルタ、ローパスフィルタ等でもよい。動的フィルタ1fは、特定の周波数をカットするノッチフィルタ、特定の帯域のレベルを上げるパラメトリックイコライザ等でもよい。動的フィルタ1fにより抽出される出力は、エンベロープ出力、クレストファクト出力等でもよい。フィルタの種類に応じて、積算グラフを作成する積算方法は変化する。
 実施の形態2
 実施の形態2は、良品及び不良品のスピーカ60に係る信号のレベル差を算出してから、レベル差に対応する信号を複数の動的フィルタ1fでフィルタ処理し、検査用のフィルタを決定する形態に関する。
 なお、実施の形態2において、実施の形態1と同様である構成要素には同一の参照番号を付してその詳細な説明を省略する。
 図7は、実施の形態2に係る検査システムの動作の概要を示す説明図である。
 サンプル用のスピーカ60に係る複数の良品及び複数の不良品へ、正弦波のログスイープ信号が入力される。正弦波のログスイープ信号は、周波数が指数関数的かつ連続的に時間変化する正弦波信号である。サンプル用のスピーカ60が発生した音に係る信号がフィルタ決定装置20に入力される。CPU21は、入力した信号をハードディスク24に記憶する。
 CPU21は、ハードディスク24に記憶した信号を、当該信号に対応するサンプル用のスピーカ60に係る良品/不良品の識別情報と共に読み出す。CPU21は、読み出した信号から、時間、周波数、レベルの3次元からなるスペクトログラムを算出し、各時間における周波数毎に最大周波数から該周波数までの該時間におけるレベルの積算をする。図7に、良品及び不良品について夫々算出された積算グラフ1m、2mが、横軸を時間、縦軸を周波数、高さ方向の軸をレベルとして、示されている。CPU21は、算出した積算グラフ1m、2mを良/不良の別と共にハードディスク24に記憶する。CPU21は、この積算グラフ算出処理を、ハードディスク24に記憶したサンプル用のスピーカ60の信号全てについて繰り返す。
 実施の形態1では、聴感テストにより選別した1つの良品のスピーカ60における積算グラフ1mを良品の参照データとした。しかし、実施の形態2では、複数の良品のスピーカにおける積算グラフ1mから、良品を代表する参照データを算出する。
 図8は、良品のスピーカ60について算出されたスペクトログラム1sの一例である。図8のスペクトログラム1sにおいて、横軸は時間t、縦軸は周波数f、紙面に対する垂直軸はレベルを夫々示している。なお、図8におけるレベルの大小は濃淡の違いで表現されている。図8の下部には、参考にフィルタ決定装置20へ測定部30から入力された信号波形が示されている。
 図9は、スペクトログラムのある時点におけるレベルの周波数分布例を示す説明図である。図9の横軸はレベル、縦軸は周波数fである。
 複数の良品を区別するためのインデックスをiとした場合、スペクトログラム1s上で離散的に分布する各時点/周波数での点(t,f)におけるレベルは、次の(1)式で表される。
Figure JPOXMLDOC01-appb-M000001
 レベルP(i)(t,f)は、2次元配列のデータである。
 CPU21は、各良品のスペクトログラム1sにおいて、夫々の時点で最高周波数f maxから複数の周波数の点夫々までのレベルの積算値を夫々算出して積算グラフを作成する。ある周波数の点まで積算されたレベルの積算値は、その周波数をカットオフ周波数とするハイパスフィルタの出力信号のレベルに対応する。
 積算値g(i)(t,f)は、次の(2)式で表される。
Figure JPOXMLDOC01-appb-M000002
 図10は、ある時点におけるレベルの積算値の周波数分布例を示す説明図である。図10の横軸はレベル、縦軸は周波数である。図10において、点線はレベル、実線はレベルの積算値を示している。すなわち、図10における実線上の黒丸は、ある時点における最高周波数f maxからある周波数の点までのレベルの積算値を示している。積算値g(i)(t,f)は、2次元配列のデータである。
 図11は、積算値g(i)(t,f)の3次元グラフ、つまり積算グラフの例である。図11の縦軸は積算値、平面方向の軸は夫々時間t及び周波数fである。
 複数の良品のスピーカ60毎に算出したレベルの積算値g(i)(t,f)夫々から、各時点/周波数点(t,f)毎に最大値を夫々抽出する。抽出した最大値G(t,f)は、複数の良品全てについてのレベル積算値の最大値であり、良品を代表するレベル積算値1mである。
 検査対象のスピーカ60を良品と不良品とに判別するフィルタは、サンプルとして用意した良品の全てが良品に判定されるフィルタである必要がある。そのためには、最も不良品寄りの音出力特性を有する良品に対応するレベルの最大値を採用することが安全側に立ち、好ましい。
 最大値G(t,f)は、次の(3)式で表される。
Figure JPOXMLDOC01-appb-M000003
 複数の不良品のスペクトログラム2sから、良品と同様にレベルの積算値2mを算出する。複数の不良品を区別するためのインデックスをiiとした場合、スペクトログラム2s上の各時点/周波数での点(t,f)におけるレベルは、次の(4)式で表される。
Figure JPOXMLDOC01-appb-M000004
 各不良品のスペクトログラム2sから、ある時点における各周波数までのレベルの積算値を夫々算出する。積算値g(ii)(t,f)は、次の(5)式で表される。
Figure JPOXMLDOC01-appb-M000005
 次に、複数の良品から取得したレベル積算値の最大値G(t,f)に対する複数の不良品から取得したレベル積算値g(ii)(t,f)との比over(ii,t,f)を、複数の不良品毎に算出する。over(ii,t,f)は、デシベル差分である。対数の減算が真数の除算に対応するように、比over(ii,t,f)はレベル差に対応する。
 over(ii,t,f)は、次の(6)式で表される。
Figure JPOXMLDOC01-appb-M000006
 図12は、比over(ii,t,f)の計算を概念的に示す説明図である。図12の左上に、良品に係るG(t,f)のグラフを示す。図12の右上に、ii番目の不良品に係るg(ii)(t,f)を示す。良品と不良品とに係るレベルの差分を算出することは、比over(ii,t,f)を計算することに対応する。図12の下部に、ii番目の不良品に関するover(ii,t,f)のグラフを示す。
 CPU21は、良品のレベル積算値1mと、複数の不良品のレベル積算値2m夫々との差分1dに対応するover(ii,t,f)を、複数の不良品毎に算出する。図7に示すように、CPU21は、複数の不良品毎に算出した積算値の差分グラフ1dを、予め用意した複数の動的フィルタ1fの曲線で順次切断処理する。該切断面は、該動的フィルタ1fで複数の良品を処理したときのレベルの最大値と該不良品を処理したレベルとのレベル差のグラフになる。
 図13は、図12におけるii番目の不良品に係るover(ii,t,f)のグラフ上で、不良品を良品群から判別するフィルタのカットオフ周波数の時間変化を示すカーブを示す説明図である。フィルタのカットオフ周波数は、時間tの経過に伴って変化していく。図13におけるカットオフ周波数の変化は、図12の下部におけるover(ii,t,f)のグラフの(t,f)平面上に、太線で示されている。ここで、不良品を良品群から判別するフィルタのカットオフ周波数の時間変化を示すカーブをカットオフカーブCと呼ぶ。図13では、カットオフカーブCが太線で示されている。
 図14は、図13をカットオフカーブCに沿って切断した断面図である。図14の縦軸は、良品のレベル積算値の最大値G(t,f)と、ii番目の不良品のレベル積算値g(ii)(t,f)との差分である。図14の横軸は、周波数である。レベル積算値は最高周波数から周波数の点までのレベルの積算値であるため、図14は、高周波数帯域を濾波するハイパスフィルタにより、良品と不良品との間において、フィルタ処理されたレベル積算値に係るレベル差を示している。図14における矢印は、レベル差の最大値、つまり不良品(ii)の検出値を示している。カットオフカーブCによる不良品(ii)の検出値Over(ii,C)は、次の(7)式で表される。
Figure JPOXMLDOC01-appb-M000007
 図7に戻り、説明を続ける。図7に、フィルタ処理後のレベル差1lが、縦軸をレベル、横軸を周波数として示されている。CPU21は、カットオフカーブCによる検出値Over(ii,C)を複数の不良品毎に算出する。CPU21は、複数の不良品毎に算出した各検出値を、不良品の識別情報と対応付けてハードディスク24にリスト200として書き込む。
 CPU21は、リスト200に書き込んだ各検出値の最小値、つまり不良検出能力値を抽出する。
 図15は、n個の不良品の検出値Over(ii,C)から、不良検出能力値D(C)を抽出する処理を概念的に示す説明図である。
 最小の検出値、すなわち不良検出能力値D(C)は、D(C)に対応するカットオフカーブCで全ての良品の積算グラフ1mと、全ての不良品の積算グラフ2mとを切り分けることができる量である。不良検出能力値D(C)は、次の(8)式で表される。
Figure JPOXMLDOC01-appb-M000008
 再び図7に戻り、説明を続ける。CPU21は、抽出した不良検出能力値D(C)を、フィルタの識別情報と対応付けて、ハードディスク24にリスト300として書き込む。CPU21は、予め用意した動的フィルタ1fを変更して、上記の処理を繰り返す。その結果、リスト300には、用意した動的フィルタ1fの数だけ不良検出能力値D(C)が書き込まれる。
 CPU21は、予め用意した全ての動的フィルタ1fについて不良検出能力値D(C)の抽出処理が終了した場合、リスト300に書き込まれた各不良検出能力値D(C)から、最大の不良検出能力値D(C)に対応するフィルタを検査用のフィルタと決定する。
 CPU21は、決定した検査用のフィルタで、検査対象のスピーカ60を検査する。
 CPU21は、検査対象のスピーカ60に係る信号を入力する。CPU21は、入力した信号を決定した検査用のフィルタでフィルタ処理する。CPU21は、フィルタ処理した検査対象のスピーカに係る信号と、予め用意した検査基準データとを比較する。CPU21は、比較結果に基づいて、検査対象のスピーカ60の良否を判定する。
 なお、比over(ii,t,f)は、フィルタ処理した検査対象のスピーカ60に係るレベルが良品と不良品とのいずれに対応するか判定するための検査基準データの参考データとして利用することができる。
 実施の形態2では、複数の良品のスピーカ60におけるスペクトログラム1sのレベル積算値から、良品を代表する参照データを求めた。しかし、実施の形態1と同様に、1つの良品のスピーカ60におけるスペクトログラム1sを参照データとしてもよい。かかる場合、1つの良品に係る各時点/周波数での点(t,f)におけるレベル積算値を参照データとする。
 また、実施の形態1に、複数の良品のスピーカ60におけるスペクトログラム1sのレベル積算値から、良品を代表する参照データを求める手法を適用してもよい。かかる場合、複数の良品毎に各時点/周波数での点(t,f)におけるレベルを求め、各点(t,f)毎に最大のレベルを参照データとする。
 実施の形態3
 実施の形態3は、最大の不良検出能力値D(C)を求めるまでの演算量を抑えて、スピーカ60の検査に用いるフィルタを決定する形態に関する。不良検出能力値D(C)が最大となるフィルタを決定するため、予め用意した多くのパターンについて不良検出能力値D(C)を算出した場合、演算量が膨大になり、フィルタの決定が現実的でなくなる。そこで、実用性が高い3種類のパターンから検査用のフィルタを決定する。なお、フィルタを決定することは、カットオフカーブCの形状を決定することに対応する。
 なお、実施の形態3において、実施の形態1、2と同様である構成要素には同一の参照番号を付してその詳細な説明を省略する。
 検査用のフィルタを決定するための動的フィルタ1fは、ハイパスフィルタである。動的フィルタ1fのカットオフ周波数は、スイープ信号を入力したスピーカ60からの音に基づく信号の基本波又は高調波の周波数変化に追従して変化するものである。動的フィルタ1fのパラメータの1つは、基本波の周波数に対する高調波の周波数の倍率mである。
 図16A、図16B及び図16Cは、パターン例を示す説明図である。図16A、図16B及び図16C夫々における横軸は時間、縦軸は周波数である。図16A、図16B及び図16Cのパターンを便宜上aパターン、bパターン及びcパターンと呼ぶ。
 aパターンの形状は、カットオフ周波数が次数m=一定整数値である高調波の周波数の時間変化と同一の形状である。
 bパターンの形状は、高次数の高調波に追随してカットオフ周波数が高くなり、その後より低次数の高調波に追随してカットオフ周波数が高くなる形状である。そして、bパターンの形状は、高次数の高調波からより低次数の高調波へ遷移する中間の周波数帯域で、カットオフ周波数が一定となる形状である。良品として判別されるスピーカ60であっても、例えば、図8に示す良品のスペクトログラム1sにおいて、時間360ms~480ms、周波数10kHzの付近に、複数の高調波で、その高周波数側における境界の周波数がほぼ一定に現れている。良品のスピーカ60であっても、時間360ms~480ms、周波数10kHzの付近における複数の高調波を許容する検査仕様がある場合、良品と不良品とに共通する信号を遮断するフィルタの形状は、bパターンの形状が合理的である。
 カットオフ周波数は、基本波の周波数が開始周波数から周波数f0まで変化する間、次数m0の高調波の周波数に追従する。カットオフ周波数は、中間の周波数帯域において、基本波の周波数が周波数f0から周波数f1(>f0)まで変化する間、一定周波数(f0×m0)を保つ。カットオフ周波数は、基本波の周波数が周波数f1から終了周波数までの間、次数m1(<m0)の高調波の周波数に追従する。なお、f0×m0=f1×m1である。例えば、f0=500Hz、m0=7次、f1=1750Hz、m1=2次である場合、500Hz×7=1750Hz×2である。
 cパターンは、bパターンの変形例である。cパターンの形状は、カットオフ周波数がbパターンの周波数一定区間(基本波の周波数がf0からf1へ変化する区間)において一定ではなく、徐々に変化する形状である。中間の周波数帯域において、多少のノイズを発生してもよいという検査仕様のスピーカ60の場合、bパターンよりもcパターンの方が、良品と不良品との間の信号のレベル差がより大きくなることがある。そこで、cパターンの形状を考慮する。
 cパターンの場合、基本波の周波数がf0からf1へ変化する区間において、カットオフ周波数は単調増加する。時間軸と周波数軸とを有する図において、時間及び周波数を対数で表した場合、基本波の周波数がf0からf1へ変化する区間では、cパターンの形状は直線である。なお、f0×m0<f1×m1である。例えば、例えば、f0=500Hz、m0=7次、f1=3kHz、m1=2次である場合、500Hz×7<3kHz×2である。
 カットオフカーブCの形状は、他にも多くのパターンが考えられる。しかし、処理が複雑になることと、スピーカ60の異常音を検査する上では上記の3パターンで多くの場合最適化できることから、ここではカットオフカーブCの形状をa、b及びcパターンの3種類にとした。
 スピーカ60は良品であっても、スピーカ60における2次歪、3次歪等の低次歪はレベルが高い場合が多い。そのため、スピーカ60の良否を検査するためのハイパスフィルタにおけるカットオフ周波数が追従する高調波の次数mは、あまり低いものは適さない。また、この次数mが高すぎても、ハイパスフィルタのカットオフ周波数が高くなりすぎるため、ハイパスフィルタは不良品のスピーカ60が発生する異常音の周波数の信号を遮断してしまう。そこで、実施の形態3では、不良検出能力値D(C)の最大値を算出するまでの計算量を低減させるために、動的フィルタ1fのフィルタ係数を変化させるにあたり、変化させる動的フィルタ1fのカットオフ周波数が追従する高調波の次数mを予め一定の範囲に制限する。ここで、当該高調波における次数mの制限範囲の上限値及び下限値を夫々m max及びm minとする。なお、m max及びm minは、良品と不良品とを判別することができるカットオフカーブCに対応する経験値である。
 計算量の低減のために、検査対象のスピーカ60における良品と不良品とのバラツキも考慮して、1つのフィルタが良品のスペクトログラム1sと、不良品のスペクトログラム2sとを区別できるレベル差を、予め一定の範囲に制限する。ここで、レベル差の範囲の上限値及び下限値に対応する不良検出能力値D(C)を夫々D max及びD minとする。
 D max及びD minは、良品と不良品とを判別するレベル差の経験値である。例えば、レベル差が6~10デシベルより大きい場合、良品と不良品とを判別することができることが判明している。そこで、D maxとして例えば15デシベルが設定され、D minとして例えば6デシベルが設定される。
 次に、フィルタ決定処理について説明する。
 フィルタ決定処理の大きな流れとして、CPU21は、フィルタの形状を低周波数側から高周波数側へ向かって変化させて、不良検出能力値D(C)がD maxを超えるフィルタを検査用のフィルタに決定する。その際、CPU21は、aパターンからフィルタ決定処理を開始する。CPU21は、aパターンについて算出したD(C)がD maxを超えなかった場合、bパターンについてフィルタ決定処理を実行する。CPU21は、bパターンについて算出したD(C)がD maxを超えなかった場合、cパターンについてフィルタ決定処理を実行する。
 CPU21は、D max、D min、m max及びm minを受け付ける。CPU21は、変数D0にD minを代入する。変数D0は、複数の不良品及び複数のフィルタについて夫々算出されたD(C)がより大きな値に更新される都度、その処理時点における最大のD(C)を保持する作業変数である。CPU21は、算出したD(C)が受け付けたD maxを超えた場合、その時点でフィルタの検索を打ち切る。これにより、処理時間の短縮を図ることができる。
 図17は、aパターンからD(C)を求める処理のイメージを示す説明図である。
 CPU21は、aパターンについて、次数mにm minを設定してD(C)を算出する。つまり、CPU21は、サンプル用の各不良品iiについてover(ii,t,f)平面上におけるm=m minのパターン上のデータからOver(ii,C)を求める。CPU21は、全ての不良品のOver(ii,C)の中から最小値D(C)を抽出する。
 CPU21は、次数mを1だけインクリメントしたパターンについて同様に最小値D(C)を抽出する。CPU21は、抽出した最小値D(C)がD0より大きい場合、抽出した最小値D(C)をD0に代入する。CPU21は、同様の処理を次数mがm maxになるまで繰り返す。ただし、CPU21は、処理の途中で、抽出した最小値D(C)がD maxを超えた場合、検索を終了する。その時点で、CPU21は、aパターンのフィルタからスピーカ60を検査するフィルタを決定したことになる。
 CPU21は、aパターンについて抽出した最小値D(C)がD maxを超えなかった場合、bパターンについて検査用のフィルタの検索を行う。その際、CPU21は、フィルタに関する演算範囲を制限する処理を実行する。
 図18は、フィルタに関する演算範囲の制限を示す説明図である。図18に示す各図は、いずれもが横軸は時間、縦軸は周波数である。図18の各図には、参考のためにスピーカ60が発生した音に係る基本波の周波数時間変化曲線が破線で表示されている。図18の上部の図には、n個の不良品に関するover(ii,t,f)の分布が白抜き及び灰色の領域で示されている。白抜きの領域は、各不良品に関するover(ii,t,f)の値がaパターンにおけるフィルタ決定で得られたD0より大きい領域である。灰色の領域は、各不良品に関するover(ii,t,f)の値がaパターンにおけるフィルタ決定で得られたD0以下の領域である。白抜きの領域は、スピーカ60が発生する異常音のレベルが大きい領域に対応する。
 図18の下部の図には、n個の不良品に関するover(ii,t,f)の分布を重畳したものに、フィルタに関する演算範囲を示す領域が示されている。図18において、破線で示されている領域の外側は、良品と不良品との信号のレベル差が小さく、フィルタを決定するために、パターンの形状を変化させる必要がない領域である。従って、CPU21は、パターンの形状を変化させる場合、破線部分の外側を、フィルタを決定する演算処理から除外する。そのために、CPU21は、不良品に関するover(ii,t,f)の分布に基づいて、演算が制限される範囲を特定する次数mの範囲(mL,mH)と、カットオフ周波数の範囲(fmax,fmin)を決定する。これにより、CPU21は、処理速度を向上させることができる。
 CPU21は、図18に示した演算範囲の制限を実行した後、bパターンのフィルタについて検査用のフィルタ決定を行う。bパターンの形状は、次数m0の高調波の部分、周波数一定区間の部分及び次数m1の高調波の部分の各形状で決定される(図16B参照)。CPU21は、これらの3つの部分を順次変更して、パターンの形状を動的に変化させる。次数m0の高調波の部分は、基本波の周波数が最低周波数からf0まで変化する部分に対応する。周波数一定区間の部分は、基本波の周波数がf0からf1まで変化する部分に対応する。次数m1の高調波の部分は、基本波の周波数がf1から最高周波数まで変化する部分に対応する。
 図19は、bパターンの形状を変更する処理のイメージを示す説明図である。次数m0の高調波の部分の形状は、当該高調波の次数で決定される。周波数一定区間の部分の形状は、一定周波数を示す線分の両端位置で決定される。当該線分の一端位置は、周波数f0と、高調波の次数m0とで決定される。当該線分の他端位置は、周波数f1と、高調波の次数m1とで決定される。次数m1の高調波の部分の形状は、当該高調波の次数で決定される。よって、bパターンの形状は、次数m0、周波数f0、周波数f1及び次数m1で決定される。bパターンの場合、f0×m0=f1×m1であることから、f1はf0×m0/m1で求められる。従って、bパターンの形状は、次数m0、周波数f0及び次数m1で決定される。
 次数m0を変化させる範囲は、m min+1からm maxまでである。次数m1を変化させる範囲は、m minからm0-1までである。f0を変化させる範囲は、発生させた最低周波数から、不良品に関するレベル差に基づいて演算の制限範囲として求めたfmaxまでである。ここでの最低周波数は、スピーカ60の特性を考慮して、予め決められた周波数であり、例えば20Hzである。
 bパターンの各部分を変更する考え方を説明する。bパターンの各部分を変更する順序は、次数m1の高調波の部分、周波数一定区間の部分、次数m0の高調波の部分の順である。CPU21は、次数m0、周波数f0を固定して、次数m1をm minからm0-1まで変化させる。次にCPU21は、周波数f0を一定値Δfだけ高くする。CPU21は、次数m0、一定値Δfだけ高くした周波数f0を固定して、次数m1をm minからm0-1まで変化させる。次にCPU21は、周波数f0を更に一定値Δfだけ高くして、次数m1をm minからm0-1まで変化させる。CPU21は周波数f0がfmaxに達した場合、次数m0を1だけインクリメントして、m min+2に設定する。そして、CPU21は、次数m0、周波数f0を固定して、次数m1をm minからm0-1まで変化させる。このように、CPU21は、次数m1、周波数f0、次数m0をこの順で順次決められた範囲の中で変化させることにより、bパターンの形状を変更する。
 図20は、bパターンの各部分を順次変更する処理の手順の一例を示すフローチャートである。
 CPU21は、次数m0にm min+1を、周波数f0に最低周波数を、次数m1にm minを代入する(ステップS201)。この場合、既にm1=m0-1なので、CPU21は、次に次数m0にm min+2を、周波数f0に最低周波数を、次数m1にm minを代入する(ステップS202)。CPU21は、次数m1をインクリメントする(ステップS203)。CPU21は、m1=m0-1か否かを判定する(ステップS204)。CPU21は、m1=m0-1でないと判定した場合(ステップS204:NO)、ステップS203へ処理を戻す。CPU21は、m1=m0-1であると判定した場合(ステップS204:YES)、次数m1にm minを代入する(ステップS205)。CPU21は、f0にΔfを加算する(ステップS206)。
 CPU21は、f0=fmaxか否かを判定する(ステップS207)。CPU21は、f0=fmaxでないと判定した場合(ステップS207:NO)、ステップS203へ処理を戻す。CPU21は、f0=fmaxであると判定した場合(ステップS207:YES)、周波数f0に最低周波数を、次数m1にm minを代入する(ステップS208)。CPU21は、次数m0をインクリメントする(ステップS209)。CPU21はm0=m maxか否かを判定する(ステップS210)。CPU21は、m0=m maxでないと判定した場合(ステップS210:NO)、ステップS203へ処理を戻す。CPU21はm0=m maxであると判定した場合(ステップS210:YES)、処理を終了する。
 なお、CPU21は、次数m0、m1を夫々変化させる範囲を、不良品に関するレベル差に基づいて求めた演算に係る次数の制限範囲(mL,mH)に制限する。また、CPU21は、f0×m0が変化する範囲を、不良品に関するレベル差に基づいて求めた演算に係るカットオフ周波数の制限範囲(fmax,fmin)に制限する。
 CPU21がbパターンについて最小値D(m)を求める処理を説明する。CPU21は、bパターンのフィルタ(m0,f0,m1,f1)において、パターンを3つの部分C0、C1、C01に分割する。
 図21は、bパターンにおける3つの部分C0、C1、C01を示す説明図である。C0は、基本波の周波数が最低周波数からf0まで変化する区間に対応する部分(次数m0の高調波の部分)である。C1は、基本波の周波数がf1から最高周波数まで変化する区間に対応する部分(次数m1の高調波の部分)である。C01は、基本波の周波数がf0からf1まで変化する区間に対応する部分(周波数一定区間の部分)である。
 CPU21は、aパターンにおけるフィルタ決定処理の段階で、次数miの高調波に対応するパターン上の差分を算出する際、C0、C1に夫々対応するOver0(ii,f,mi)、Over1(ii,f,mi)を計算する(図12参照)。CPU21は、計算したOver0(ii,f,mi)、Over1(ii,f,mi)を予めRAM23に記憶しておく。
 図22Aは、パターン上の差分値を示す説明図である。図22Bは、C0に対応するOver0(ii,f,mi)を算出する処理を示す説明図である。図22Cは、C1に対応するOver1(ii,f,mi)を算出する処理を示す説明図である。図22A、図22B及び図22Cにおいて、夫々時間は左側から右側へ経過している。CPU21は、aパターンにおけるフィルタ決定処理の段階で、C0に対応する最大値Over0(ii,f,mi)を、入力信号の周波数変化の開始から終了へ向かって順に差分値を比較して抽出しておく。CPU21は、aパターンのフィルタ決定処理の段階で、C1に対応する最大値Over1(ii,f,mi)を、入力信号の周波数変化の終了から開始へ向かって順に差分値を比較して抽出しておく。
 ここで、CPU21は、不良品毎に、C0、C1、C01の夫々について、Over(ii,Cx)を求め、求めたOver(ii,Cx)の最大値をOver(ii,C)に設定する。Cxは、C0、C1及びC01のいずれかである。
 図23は、C0、C1及びC01に対応するOver(ii,Cx)から、最大値をOver(ii,C)に設定するイメージを示す説明図である。CPU21は、C0の区間では、Over(ii,C0)に、RAM23に記憶したOver0(ii,f,mi)を代入する。CPU21は、C1の区間では、Over(ii,C1)に、RAM23に記憶したOver1(ii,f,mi)を代入する。CPU21は、C01の区間では、算出済みのover(ii,t,F)のデータからOver(ii,C01)を求める。CPU21は、Over(ii,C0)、Over(ii,C1)及びOver(ii,C01)のうち、最大値をOver(ii,C)に設定する。
 CPU21は、全ての不良品について最大値Over(ii,C)を求め、求めた最大値Over(ii,C)から、最小値を抽出してD(C)に設定する。CPU21は、抽出した最小値D(C)がD0より大きい場合、抽出した最小値D(C)をD0に代入して、D0を更新する。CPU21は、抽出した最小値D(C)がD0以下である場合、例えば次数m1の高調波の部分についてm1を1つインクリメントすることにより、パターンの形状を変化させて、上記の処理を繰り返す。ただし、CPU21は、処理の途中で、抽出した最小値D(C)がD maxを超えた場合、フィルタ決定処理を終了する。
 CPU21は、bパターンについて抽出した最小値D(C)の全てがD maxを超えなかった場合、cパターンについてフィルタ決定処理を実行する。その際、CPU21は、ここまでの処理により算出した不良検出能力値の最大値であるD0に基づいて、図18に示した演算範囲の制限を実行した後、cパターンについてフィルタ決定処理を実行する。
 図24は、cパターンの形状を変更する処理のイメージを示す説明図である。
 cパターンの形状は、次数m0の高調波の部分、周波数が徐々に高くなる部分及び次数m1の高調波の部分の各形状で決定される(図16C参照)。次数m0の高調波の部分、周波数が徐々に高くなる部分及び次数m1の高調波の部分は、夫々図21のC0、C01及びC1に対応する。そこで、図24の時間軸には、参考にC0、C01及びC1が添えられている。
 CPU21は、bパターンの場合と同様に、次数m0、周波数f0、次数m1を順次変化させて、cパターンの形状を変化させる。ただし、cパターンの場合、周波数f1は1つに定まらないため、CPU21は周波数f1も順次変化させて、cパターンの形状を変化させる。周波数f1を変化させる範囲は、発生させた最低周波数からfmaxまでである。ここでの最低周波数は、スピーカ60の特性を考慮して、予め決められた周波数であり、例えば20Hzである。CPU21は、周波数f1、次数m1、周波数f0、次数m0をこの順で順次決められた範囲の中で変化させることにより、cパターンの形状を変更する。
 ここでの周波数f1の変化とは、周波数が徐々に高くなる部分(C01に対応する)の他端点の位置を次数m1の高調波に沿って、高周波数側へ変化させることに対応する。具体的には、CPU21は、f1×m1がf0×m0に一定値を加算した値からfmaxまで変化するように、他端点の位置を変化させる。
 cパターンの場合におけるD(C)の算出について説明する。
 CPU21は、bパターンの場合と同様に、cパターンを3つの部分C0、C1、C01に分割する。CPU21は、Over(ii,Cx)を求め、その最大値をOver(ii、C)に設定する。
 CPU21は、全ての不良品について最大値Over(ii,C)を求め、求めた最大値Over(ii,C)から、最小値を抽出してD(C)に設定する。CPU21は、抽出した最小値D(C)がD0より大きい場合、抽出した最小値D(C)をD0に代入して、D0を更新する。CPU21は、抽出した最小値D(C)がD0以下である場合、cパターンの形状を変化させて、上記の処理を繰り返す。ただし、CPU21は、処理の途中で、抽出した最小値D(C)がD maxを超えた場合、フィルタ決定処理を終了する。CPU21は、aパターン、bパターン及びcパターンの全てについて、D maxを超える最小値D(C)が見つからなかった場合、ここまでの処理で不良検出能力値の最大値D0がD minを超えているとき、最大値D0に対応するフィルタを検査用のフィルタと決定する。
 図25及び図26は、フィルタ決定処理の手順の一例を示すフローチャートである。
 CPU21は、D max、D min、m max及びm minを受け付ける(ステップS301)。CPU21は、変数D0にD minを代入する(ステップS302)。CPU21は、変数Pに1を代入する(ステップS303)。変数Pは、パターンを判別するためのフラグ変数である。変数Pには、aパターンの場合には1、bパターンの場合には2、cパターンの場合には3が夫々代入される。
 CPU21は、P=1か否かを判定する(ステップS304)。CPU21は、P=1でないと判定した場合(ステップS304:NO)、D0に基づいて、フィルタに関する演算範囲を制限し(ステップS305)、ステップS306へ処理を進める。ステップS305において、CPU21は、演算範囲を特定する次数mの範囲(mL,mH)と、カットオフ周波数の範囲(fmax,fmin)とを決定する。CPU21は、P=1であると判定した場合(ステップS304:YES)、未処理の不良品(ii)についてOver(ii、C)を算出する(ステップS306)。CPU21は、全ての不良品についてOver(ii、C)を算出したか否かを判定する(ステップS307)。CPU21は、全ての不良品についてOver(ii、C)を算出していないと判定した場合(ステップS307:NO)、ステップS306へ処理を戻す。
 CPU21は、全ての不良品についてOver(ii、C)を算出したと判定した場合(ステップS307:YES)、D(C)を算出する(ステップS308)。CPU21は、D(C)>D0か否かを判定する(ステップS309)。CPU21は、D(C)>D0であると判定した場合(ステップS309:YES)、D0にD(C)を代入する(ステップS310)。ステップS310でD0に代入されたD(C)に対応するパターンは、検査用として決定されるフィルタの候補に対応する暫定パターンである。CPU21は、D(C)>D0でないと判定した場合(ステップS309:NO)、ステップS313へ処理を進める。
 CPU21は、D(C)≧D maxか否かを判定する(ステップS311)。CPU21は、D(C)≧D maxであると判定した場合(ステップS311:YES)、D(C)に対応するフィルタのフィルタ係数を出力し(ステップS312)、処理を終了する。CPU21は、D(C)≧D maxでないと判定した場合(ステップS311:NO)、処理中のフィルタパターンに含まれる全てのフィルタについて処理を実行したか否かを判定する(ステップS313)。CPU21は、全てのフィルタについて処理を実行していないと判定した場合(ステップS313:NO)、Pの値に応じて、パターンの形状を変更し(ステップS314)、ステップS306へ処理を戻す。
 CPU21は、処理中のフィルタパターンに含まれる全てのフィルタについて処理を実行したと判定した場合(ステップS313:YES)、P=1か否かを判定する(ステップS315)。CPU22は、P=1であると判定した場合(ステップS315:YES)、Pに2を代入し(ステップS316)、ステップS304へ処理を戻す。CPU21は、P=1でないと判定した場合(ステップS315:NO)、P=2か否かを判定する(ステップS317)。CPU21は、P=2であると判定した場合(ステップS317:YES)、Pに3を代入し(ステップS318)、ステップS304へ処理を戻す。
CPU21は、P=2でないと判定した場合(ステップS317:NO)、D0に暫定パターンに対応するD(C)が代入されているか否かを判定する(ステップS319)。CPU21は、D0に暫定パターンが代入されていると判定した場合(ステップS319:YES)、ステップS312に処理を進める。CPU21は、D0に暫定パターンに対応するD(C)が代入されていないと判定した場合(ステップS319:NO)、良品/不良品を分離可能なカットオフカーブCはない旨を出力し(ステップS320)、処理を終了する。
 本実施の形態では、CPU21は、パターンについて算出したD(C)がD maxを超えなかった場合、決定対象のフィルタをaパターンからbパターンへ、更にはbパターンからcパターンへ変更している。しかし、CPU21は、aパターンについて算出したD(C)がD maxを超えなかった場合、bパターンへ移行せず、その時点のD0に対応するフィルタを検査用のフィルタと決定し、フィルタ決定処理を終了してもよい。同様に、CPU21は、bパターンについて算出したD(C)がD maxを超えなかった場合、cパターンへ移行せず、その時点のD0に対応するフィルタを検査用のフィルタと決定し、フィルタ決定処理を終了してもよい。これにより、処理時間の短縮を図ることができる。
 本実施の形態では、CPU21は、フィルタに関する演算範囲を限定した。しかし、CPU21は、フィルタに関する演算範囲を限定せずに、aパターンからbパターンへ、bパターンからcパターンへフィルタ決定処理の対象を変更してもよい。
 本願に係る検査装置によれば、音又は振動を発生する装置の仕様毎に、当該装置を検査するためのフィルタを決定することができる。
 CPU21は、フィルタを動的に変更する処理を繰り返すことにより、予め用意したスピーカ60の良品及び不良品の間における信号のレベル差が最大となるフィルタを決定することができる。
 フィルタ決定装置20によれば、最大の不良検出能力値D(C)を求めるまでの演算量を抑えることができる。
 フィルタ決定装置20は、受け付けたD max及びD minの範囲内で、フィルタの決定を実行する。そのため、フィルタ決定装置20は、計算するD(C)の数を低減することができる。また、フィルタ決定装置20は、算出したD(C)がD maxを超えた場合、その時点でフィルタの決定を終了し、D maxを超えたD(C)に対応するフィルタを検査用のフィルタと決定する。これにより、フィルタ決定装置20は、処理時間の短縮を図ることができる。
 フィルタ決定装置20は、aパターン、bパターン及びcパターンの3パターン内でフィルタの決定処理を実行する。フィルタ決定装置20は、限定したフィルタパターンの範囲内、かつ受け付けたm max及びm minの次数範囲内でパターンの形状を変更する。このような限定条件も演算量の低減に資する。更には、フィルタ決定装置20は、bパターン及びcパターンにおけるフィルタの決定において、D0に基づいて演算範囲を制限する。この演算範囲の制限も演算量を低減させることができる。
 開示された実施の形態は、全ての点で例示であって制限的なものではないと考えられるべきである。本発明の範囲は上述の説明ではなくて特許請求の範囲によって示され、特許請求の範囲と均等の意味及び範囲内での全ての変更が含まれることが意図される。
 各実施の形態で記載されている技術的特徴(構成要件)はお互いに組合せ可能であり、組み合わせすることにより、新しい技術的特徴を形成することができる。
 20 フィルタ決定装置
 21 CPU
 23 RAM
 29 I/O部
 24 ハードディスク
 30 測定部
 40 マイクロホン
 50 無響箱
 60 スピーカ
 2P プログラム

Claims (14)

  1.  検査対象の装置が発生した、周波数が時間経過につれて変化する音又は振動に係る信号を、特性を時間経過につれて変化させ得るフィルタを用いて処理して、前記装置を検査する検査方法。
  2.  検査対象の装置を検査するに先立ち、予め用意された装置の一又は複数の良品及び一又は複数の不良品夫々が発生した音又は振動に係る信号を複数のフィルタ各々で処理した複数の信号に基づいて、検査用のフィルタを決定しておき、
     決定した検査用のフィルタを用いて、検査対象の装置が発生した音又は振動に係る信号を処理し、
     処理した信号に基づいて、検査対象の装置の良否を判定する
     ことを特徴とする請求項1に記載の検査方法。
  3.  前記良品及び不良品夫々が発生した音又は振動に係る信号のレベル差を算出し、
     算出したレベル差に対応する信号を複数のフィルタ夫々で処理し、
     複数のフィルタ夫々で処理した信号に対応するレベル差のうち、最大のレベル差に対応するフィルタを検査用のフィルタと決定する
     ことを特徴とする請求項2に記載の検査方法。
  4.  前記良品及び複数の不良品夫々が発生した音又は振動に係る信号のレベル差を該複数の不良品毎に算出し、
     算出したレベル差に対応する信号を複数のフィルタ夫々で処理し、
     処理した信号に対応するレベル差の最大値を前記複数の不良品毎に算出し、
     複数のフィルタ夫々につき、前記複数の不良品毎に算出したレベル差の各最大値のうち、レベル差の最小値を選択し、
     複数のフィルタ夫々につき選択したレベル差の各最小値のうち、最大のレベル差に対応するフィルタを検査用のフィルタと決定する
     ことを特徴とする請求項2又は請求項3に記載の検査方法。
  5.  予め用意された複数の良品に係る信号から夫々算出された時間、周波数及びレベルの情報に基づく信号から、予め定められた間隔で画定された時間及び周波数毎に最大のレベルを各々抽出し、
     抽出した最大の各レベルを前記良品が発生した音又は振動に係る信号のレベルに設定し、
     設定した前記良品に係る信号のレベルに基づいて、レベル差の最大値を前記複数の不良品毎に算出する
     ことを特徴とする請求項4に記載の検査方法。
  6.  選択したレベル差の最小値が所定の値を超えた場合、レベル差の最小値の選択を終了し、
     前記所定の値を超えたレベル差の最小値に対応するフィルタを検査用のフィルタと決定する
     ことを特徴とする請求項4又は請求項5に記載の検査方法。
  7.  検査対象の装置が発生した、周波数が時間経過につれて変化する音又は振動に係る信号を、特性を時間経過につれて変化させ得るフィルタを用いて処理して、前記装置を検査する検査装置。
  8.  複数のフィルタを出力するフィルタ出力部と、
     予め用意された装置の一又は複数の良品及び一又は複数の不良品夫々が発生した音又は振動に係る信号に基づいて、前記フィルタ出力部が出力する複数のフィルタから検査用のフィルタを決定する決定部と、
     該決定部が決定した検査用のフィルタを用いて、検査対象の装置が発生した音又は振動に係る信号を処理する処理部と、
     該処理部が処理した信号に基づいて、検査対象の装置の良否を判定する判定部と
     を備える
     ことを特徴とする請求項7に記載の検査装置。
  9.  前記良品に係る信号を記憶する良品信号記憶部と、
     該良品信号記憶部が記憶する良品に係る信号及び予め用意された前記不良品に係る信号のレベル差を算出するレベル差算出部と、
     前記フィルタ出力部が出力する複数のフィルタ夫々で、前記レベル差算出部が算出したレベル差に対応する信号を処理する信号処理部と
     を備え、
     前記決定部は、前記信号処理部が処理した信号に対応するレベル差のうち、最大のレベル差に対応するフィルタを検査用のフィルタと決定するようにしてあり、
     前記処理部が処理した信号及び所定の信号を比較する信号比較部を備え、
     前記判定部は前記信号比較部が比較した結果に基づいて、検査対象の装置の良否を判定するようにしてある
     ことを特徴とする請求項8に記載の検査装置。
  10.  前記レベル差算出部は、前記良品及び複数の不良品夫々が発生した音又は振動に係る信号のレベル差を該複数の不良品毎に算出するようにしてあり、
     前記信号処理部は前記フィルタ出力部が出力する複数のフィルタ夫々で、前記レベル差算出部が前記複数の不良品毎に算出したレベル差に対応する信号を処理するようにしてあり、
     前記信号処理部が処理した信号に対応するレベル差の最大値を前記複数の不良品毎に算出する算出部と、
     前記フィルタ出力部が出力する複数のフィルタ夫々につき、前記算出部が前記複数の不良品毎に算出したレベル差の各最大値のうち、レベル差の最小値を選択する選択部と
     を備え、
     前記決定部は前記選択部が複数のフィルタ夫々につき選択したレベル差の各最小値のうち、最大のレベル差に対応するフィルタを検査用のフィルタと決定するようにしてある
     ことを特徴とする請求項9に記載の検査装置。
  11.  前記選択部は、選択したレベル差の最小値が所定の値を超えた場合、レベル差の最小値の選択を終了するようにしてあり、
     前記決定部は前記所定の値を超えたレベル差の最小値に対応するフィルタを検査用のフィルタと決定するようにしてある
     ことを特徴とする請求項10に記載の検査装置。
  12.  周波数が時間経過につれて変化する音又は振動を発生する検査対象の装置を検査するために、特性を時間経過につれて変化させ得る複数のフィルタから検査用のフィルタを決定するフィルタ決定方法であって、
     予め用意された装置の良品及び一又は複数の不良品夫々が発生した音又は振動に係る信号のレベル差を該複数の不良品毎に算出し、
     算出したレベル差に対応する信号を複数のフィルタ夫々で処理し、
     処理した信号に対応するレベル差の最大値を前記一又は複数の不良品毎に算出し、
     複数のフィルタ夫々につき、前記一又は複数の不良品毎に算出したレベル差の各最大値のうち、レベル差の最小値を選択し、
     複数のフィルタ夫々につき選択したレベル差の各最小値のうち、最大のレベル差に対応するフィルタを検査用のフィルタと決定する
     ことを特徴とするフィルタ決定方法。
  13.  周波数が時間経過につれて変化する音又は振動に係る信号に基づいて、特性を時間経過につれて変化させ得る複数のフィルタから一のフィルタを決定する処理をコンピュータに実行させるプログラムであって、
     音又は振動に係る一の信号と、一又は複数の他の信号とのレベル差を該一又は複数の他の信号毎に算出し、
     特性を時間経過につれて変化させ得るフィルタを変更し、
     変更した複数のフィルタ夫々で、算出したレベル差に対応する信号を処理し、
     処理した信号に対応するレベル差の最大値を前記一又は複数の他の信号毎に算出し、
     複数のフィルタ夫々につき、前記一又は複数の他の信号毎に算出したレベル差の各最大値のうち、レベル差の最小値を選択し、
     複数のフィルタ夫々につき選択したレベル差の各最小値のうち、最大のレベル差に対応するフィルタを一のフィルタと決定する
     処理をコンピュータに実行させることを特徴とするプログラム。
  14.  周波数が時間経過につれて変化する音又は振動を発生する検査対象の装置を検査する検査装置のために、特性を時間経過につれて変化させ得る複数のフィルタから検査用のフィルタを決定するフィルタ決定装置であって、
     複数のフィルタを出力するフィルタ出力部と、
     予め用意された装置の良品及び一又は複数の不良品夫々が発生した音又は振動に係る信号のレベル差を該複数の不良品毎に算出するレベル差算出部と、
     前記フィルタ出力部が出力する複数のフィルタ夫々で、前記レベル差算出部が算出したレベル差に対応する信号を処理する信号処理部と、
     該信号処理部が処理した信号に対応するレベル差の最大値を前記一又は複数の不良品毎に算出する算出部と、
     前記フィルタ出力部が出力する複数のフィルタ夫々につき、前記算出部が前記一又は複数の不良品毎に算出したレベル差の各最大値のうち、レベル差の最小値を選択する選択部と、
     該選択部が複数のフィルタ夫々につき選択したレベル差の各最小値のうち、最大のレベル差に対応するフィルタを検査用のフィルタと決定する決定部と
     を備える
     ことを特徴とするフィルタ決定装置。
PCT/JP2014/061479 2014-04-23 2014-04-23 検査装置 WO2015162736A1 (ja)

Priority Applications (2)

Application Number Priority Date Filing Date Title
PCT/JP2014/061479 WO2015162736A1 (ja) 2014-04-23 2014-04-23 検査装置
JP2014550576A JP5896434B1 (ja) 2014-04-23 2014-04-23 検査装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/JP2014/061479 WO2015162736A1 (ja) 2014-04-23 2014-04-23 検査装置

Publications (1)

Publication Number Publication Date
WO2015162736A1 true WO2015162736A1 (ja) 2015-10-29

Family

ID=54331925

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2014/061479 WO2015162736A1 (ja) 2014-04-23 2014-04-23 検査装置

Country Status (2)

Country Link
JP (1) JP5896434B1 (ja)
WO (1) WO2015162736A1 (ja)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6419392B1 (ja) * 2017-12-22 2018-11-07 三菱電機株式会社 音響計測システム及びパラメータ生成装置

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR102682417B1 (ko) * 2021-08-09 2024-07-05 (주)비올 고주파 피부 미용기용 핸드 피스의 탈조 검사 시스템 및 방법

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS5973737A (ja) * 1982-10-20 1984-04-26 Matsushita Electric Ind Co Ltd スピ−カシステム検査装置
JPH036400U (ja) * 1989-06-08 1991-01-22
JPH0334700A (ja) * 1989-06-29 1991-02-14 Matsushita Electric Ind Co Ltd 音質評価装置
JPH042299A (ja) * 1990-04-19 1992-01-07 Matsushita Electric Ind Co Ltd 音質評価装置
JPH04167699A (ja) * 1990-10-26 1992-06-15 Hosiden Corp ぴりつき音検出装置
JPH06315197A (ja) * 1993-04-28 1994-11-08 Hosiden Corp びりつき音検出装置
JP2007306410A (ja) * 2006-05-12 2007-11-22 Fujitsu Ltd 音再生装置検査方法、音再生装置検査プログラムおよび音再生装置検査システム
JP2009017490A (ja) * 2007-07-09 2009-01-22 Mitsubishi Electric Corp 音響診断システムおよびその診断機器

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS5973737A (ja) * 1982-10-20 1984-04-26 Matsushita Electric Ind Co Ltd スピ−カシステム検査装置
JPH036400U (ja) * 1989-06-08 1991-01-22
JPH0334700A (ja) * 1989-06-29 1991-02-14 Matsushita Electric Ind Co Ltd 音質評価装置
JPH042299A (ja) * 1990-04-19 1992-01-07 Matsushita Electric Ind Co Ltd 音質評価装置
JPH04167699A (ja) * 1990-10-26 1992-06-15 Hosiden Corp ぴりつき音検出装置
JPH06315197A (ja) * 1993-04-28 1994-11-08 Hosiden Corp びりつき音検出装置
JP2007306410A (ja) * 2006-05-12 2007-11-22 Fujitsu Ltd 音再生装置検査方法、音再生装置検査プログラムおよび音再生装置検査システム
JP2009017490A (ja) * 2007-07-09 2009-01-22 Mitsubishi Electric Corp 音響診断システムおよびその診断機器

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6419392B1 (ja) * 2017-12-22 2018-11-07 三菱電機株式会社 音響計測システム及びパラメータ生成装置
KR20200075884A (ko) * 2017-12-22 2020-06-26 미쓰비시덴키 가부시키가이샤 음향 계측 시스템 및 파라미터 생성 장치
CN111480196A (zh) * 2017-12-22 2020-07-31 三菱电机株式会社 声音计测系统和参数生成装置
KR102174553B1 (ko) 2017-12-22 2020-11-05 미쓰비시덴키 가부시키가이샤 음향 계측 시스템 및 파라미터 생성 장치
CN111480196B (zh) * 2017-12-22 2023-08-25 三菱电机株式会社 声音计测系统和参数生成装置

Also Published As

Publication number Publication date
JP5896434B1 (ja) 2016-03-30
JPWO2015162736A1 (ja) 2017-04-13

Similar Documents

Publication Publication Date Title
US9025781B2 (en) Sound quality evaluation apparatus and method thereof
EP2730099B1 (en) Estimating nonlinear distortion and parameter tuning for boosting sound
JP2007093221A (ja) 信号の波形分析方法とそのプログラム及び信号の波形分析方法を使用した車両運転特性の解析方法とそのプログラム
US10021483B2 (en) Sound capture apparatus, control method therefor, and computer-readable storage medium
KR100831539B1 (ko) 검사 기준 결정 지원 장치
JP5896434B1 (ja) 検査装置
US20090304195A1 (en) Method and device for diagnosing the operating state of a sound system
KR20070078784A (ko) 이음(異音) 추출 방법, 이음 추출 장치 및 이음 검사 장치
JP6822294B2 (ja) 評価プログラム、情報処理装置、及び評価方法
KR101293472B1 (ko) 소음 진동 분석 장치 및 방법
JP2010096547A (ja) 異音判定装置および異音判定方法
CN112700399A (zh) 缺陷检测视觉化方法及其系统
JP4630183B2 (ja) 音声信号分析装置、音声信号分析方法及び音声信号分析プログラム
JP5539437B2 (ja) 純音検査装置及びその制御方法
CN115460527A (zh) 基于心理声学的耳机失真测量方法、存储介质、耳机
JP2016188791A (ja) 音源探査装置および音源探査方法
JP2006226706A (ja) 欠陥検出方法及びプログラム
JP3688455B2 (ja) モータの作動音検査方法及びその検査装置
JP2877901B2 (ja) 騒音判別方法
Uchiyama et al. Real-Time Noise Suppression Using Harmonic/Percussive Separation with Morphological Operations for Hammering Test
JP7026026B2 (ja) びりつき音源探査システムおよびくし型フィルタの設定方法
JP6845544B2 (ja) 異音評価装置及び異音評価方法
TWI595791B (zh) 音訊檢測方法
WO2017130417A1 (ja) 生体音解析装置及び生体音解析方法、並びにコンピュータプログラム及び記録媒体
JP4471780B2 (ja) 音声信号処理装置及びその方法

Legal Events

Date Code Title Description
ENP Entry into the national phase

Ref document number: 2014550576

Country of ref document: JP

Kind code of ref document: A

121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 14890068

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 14890068

Country of ref document: EP

Kind code of ref document: A1