CN109870231A - Automatic Metrology Verification System and Identification Algorithm of Vibration Measuring Instrument - Google Patents

Automatic Metrology Verification System and Identification Algorithm of Vibration Measuring Instrument Download PDF

Info

Publication number
CN109870231A
CN109870231A CN201910199211.5A CN201910199211A CN109870231A CN 109870231 A CN109870231 A CN 109870231A CN 201910199211 A CN201910199211 A CN 201910199211A CN 109870231 A CN109870231 A CN 109870231A
Authority
CN
China
Prior art keywords
frequency
interface
measurement result
image
standard
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201910199211.5A
Other languages
Chinese (zh)
Other versions
CN109870231B (en
Inventor
金晖
何洁
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hangzhou City University
Original Assignee
Zhejiang University City College ZUCC
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 Zhejiang University City College ZUCC filed Critical Zhejiang University City College ZUCC
Priority to CN201910199211.5A priority Critical patent/CN109870231B/en
Publication of CN109870231A publication Critical patent/CN109870231A/en
Application granted granted Critical
Publication of CN109870231B publication Critical patent/CN109870231B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

本发明涉及振动测量仪器自动化计量检定系统,包括标准振源发生系统、测量结果识别系统、计算机程控系统;标准声源发生系统包括:ZigBee收发单元、ARM系统、程控信号源、电机控制单元和振动台;测量结果识别系统包括:被测仪器显示屏、摄像头、图像电信号接口单元、ARM系统、检测电路、AC接口、DC接口和ZigBee收发单元;计算机程控系统包括:ZigBee收发单元、ARM系统和计算机上位机。本发明的有益效果是:结合标准振源发生系统、测量结果识别系统、计算机程控系统形成的振动测量仪器自动化计量检定系统,实现了振动测量仪器检定的自动化,大大节省了检定过程的人力输出和时间成本。

The invention relates to an automatic metrology verification system for vibration measuring instruments, including a standard vibration source generation system, a measurement result identification system, and a computer program control system; the standard sound source generation system includes: a ZigBee transceiver unit, an ARM system, a program-controlled signal source, a motor control unit and a vibration source. The measurement result identification system includes: the display screen of the instrument under test, the camera, the image electrical signal interface unit, the ARM system, the detection circuit, the AC interface, the DC interface and the ZigBee transceiver unit; the computer program control system includes: ZigBee transceiver unit, ARM system and computer host computer. The beneficial effects of the invention are: the automatic measurement and verification system of vibration measuring instruments formed by combining the standard vibration source generating system, the measurement result identification system and the computer program control system, realizes the automation of the verification of the vibration measuring instruments, and greatly saves the manpower output and the cost of the verification process. Time costs.

Description

Vibration measuring instrument automatic gauge verification system and recognizer
Technical field
The present invention relates to automatic gauge verification system more particularly to vibration measuring instrument automatic gauge verification systems.
Background technique
Vibration measuring instrument includes analysis of vibration signal and vibration level size measuring instrument, these instruments are much for ring The fields such as border protection, belong to national measurement calibration equipment, according to national verification standard, analyze for vibration measurement, such as vibrate The 1/3oct of sound level is analyzed, and national regulations need to carry out the multiple spot centre frequency of analysis instrument the test of bandwidth performance, is needed The amplitude and frequency of artificial constantly setting measuring signal, while the index after test is subjected to artificial contrast, workload is very Greatly.For vibration level (including vertical (W.B.z) the weighted vibration level of whole body, whole body level (W.B.x-y) weighted vibration level, and not weighted (Lin) measure) measurement of size, it is provided according to the detection of ISO8041 or GB/10071, needs to test a variety of vibration frequency weighteds Artificial calibrating and index comparison are carried out with the linearity and decibel sizes values of time weighted, workload is also very big.
Summary of the invention
The purpose of the present invention is to overcome the above shortcomings and to provide a kind of structurally reasonable, vibration measurements of high degree of automation Analysis instrument measurement verification system and recognizer.To achieve the above object, the invention adopts the following technical scheme:
Vibration measuring instrument automatic gauge verification system, including standard vibration source generating system, measurement result identifying system, Computer program-control system;
Standard source generating system include: ZigBee Transmit-Receive Unit, ARM system, programme-controlled signal source, motor control unit and Shake table;ARM system receives the standard electric signal for the requirement output that computer program-control system is sent by Zigbee Transmit-Receive Unit Signal waveform, frequency and amplitude or standard vibration level size, all the way pass through circuit system carry out program control signal output, another way The frequency of shake table is hit by control motor control unit output and amplitude realizes the output of vibration level, while passing through ARM system To the signal frequency and height analysis of vibration level signal, the requirement for whether reaching standard vibration level seen, if be not achieved, adjusts ARM system The control to motor control unit output frequency and associated motor of uniting exports, and realizes motor control unit guidance striker to vibration The shock dynamics and frequency of platform, until the output of standard vibration level;The self-calibration process of realization standard vibration level;Programme-controlled signal source passes through ARM system controls AD9850 and exports corresponding signal waveform, amplitude and frequency by direct digital frequency synthesis technology.When output is marked When the vibration level signal or electric signal of standard, the ready order of signal is sent by ZigBee Transmit-Receive Unit and identifies system to measurement result System;
Measurement result identifying system includes: tested instrument display screen, camera, electric image signal interface unit, ARM system System, detection circuit, AC interface, DC interface and ZigBee Transmit-Receive Unit;ARM system receives standard by ZigBee Transmit-Receive Unit The ready order that vibration level generating system is sent, ARM system realize camera to tested by image telecommunication interface unit The measurement result of instrument display screen carries out image recognition, is surveyed by AC interface, DC interface of the detection circuit to tested instrument Amount;ARM system gets and sends the measurements to computer program-control system by ZigBee Transmit-Receive Unit after measurement result;
Computer program-control system includes: ZigBee Transmit-Receive Unit, ARM system and computer host computer;ARM system passes through The measurement result that ZigBee Transmit-Receive Unit obtains measurement result identifying system passes through USB interface again and data is sent to sequence controlled computer Host computer procedure, measurement result is sent in the output entry of preset probation redport by host computer procedure, then, passes through USB Waveform, amplitude and the frequency that interface sends the signal source output requirement required in next vertification regulation are identified to measurement result is System, measurement result identifying system pass through ZigBee Transmit-Receive Unit again and send standard vibration source system for signal source output requirement System.
The recognizer of the measurement result identifying system of vibration measuring instrument automatic gauge verification system, including following step It is rapid:
S1. vertical (W.B.z) the weighted vibration level of vibration level signal whole body, whole body level (W.B.x- are carried out by AC interface first Y) digital assay of weighted vibration level, specific as follows:
The 1/3oct analysis for completing 1Hz to 80Hz first, carries out according to following steps:
1) a 8 rank Butterworth filter devices are designed with the design method of analog filter;
The low pass chi square function of butterworth filter is expressed as follows:
Interim N takes 8,8 rank low-pass filters is designed, to obtain system function H (S);
2) predistortion is carried out to analog frequency;
It is being that, in order to guarantee that each edge frequency point is the frequency pre-established, filter is cut using two-wire row converter technique Only frequencies omegaCWith stopband side frequency ω s, it is pre- abnormal before determining simulation low-pass filter system function must to carry out as the following formula frequency Become:
H (s)=Ha(s/Ωc) (5)
3) transmission function of digital filtering is found out using Bilinear transformation method;
It is summarized as follows with the process of Bilinear transformation method design digital filter:
Normalize the frequency transformation of analog low-pass prototype to digital band-pass filter:
4) inputs a sampled value again, calculates an output valve, seeks the root mean square of each output valve;
5) the frequency band root mean square acceleration analyzed by 1/3oct is multiplied by the weighted factor, then takes r.m.s. and obtain Corresponding spectrum overlapping acceleration:
In formula: ai--- the accelerated speed effective value (m/s that the i-th frequency range is surveyed in third-octave frequency spectrum2);
ki--- the corresponding spectrum overlapping factor of the i-th frequency range in third-octave frequency spectrum;
It is analyzed by the calculating of above formula and 1/3oct, so that it may measure corresponding spectrum overlapping acceleration;
Weighted vibration level can be obtained as the following formula by weighted acceleration:
In formula: aw--- weighted acceleration virtual value (m/s2);
a0--- reference acceleration, a0=10-6m/s2
It is calculated by the digital collection of (10) formula, so that it may corresponding weighted vibration level;
Corresponding character match module is formed by these measurement results;
S2 obtains the measurement result of corresponding not weighted through ARM system sampling analysis by DC interface, which forms The respective symbols matching module of not weighted in optical identification measurement result;
The measurement result display screen of tested instrument is realized ARM system by RS232 interface by digital camera by S3 It reads, carries out optical identification.
As preferred: the algorithm of the optical identification of the step S3 are as follows:
It identifies that the image of object obtains and gray processing is realized by camera combination electric image signal unit, know for character Not, select median filtering technology realize image smoothing, by threshold value set realize image binaryzation, from binaryzation data into Line character segmentation, the method for selection is template matching method, then carries out character knowledge by extracting after image thinning using framework characteristic Not;It is extracted using framework characteristic, reduces corresponding calculation amount, only in the character feature and above-mentioned S1, S2, S3 step extracted The feature of the measurement result Character mother plate of formation is matched;An image to be identified is given, by micronization processes Character skeleton image, or the bianry image for eliminating noise not refined, it is horizontal and vertical all to cross at regular intervals, so The number across character picture, that is, number of crossings are successively write down afterwards;Template image is gathered in advance, extracts each template Pass through number, store with a matrix type, as eigenmatrix, after the feature for having extracted each character, as database It preserves;After obtaining single character picture to be identified, traversing times are equally extracted, are loaded into template data and mould The traversing times matrix of plate, successively makes comparisons, and calculates related coefficient, maximum related coefficient is exactly recognition result;Long and width etc. It is divided into five parts, is unable to equal part, round numbers part calculates traversing times in each position of halving, intersects each time, is denoted as primary It passes through, shares ten characteristic values in this way;Calculate poor method for distinguishing:
X can represent binary character to be identified in formula, and P is the eigenmatrix of template, and obtained matrix of differences is D;D= d11 2+…+d25 2For the number that the two correlation is measured, that maximum corresponding templates of the number mean that recognition result.
The beneficial effects of the present invention are:
1) vibration measurement that combined standard vibration source generating system, measurement result identifying system, computer program-control system are formed Automation equipment measurement verification system realizes the automation of vibration measuring instrument calibrating, the people of verification process is greatly saved Power output and time cost.
2) vibration level self-calibration technique, programme-controlled signal source design and Radio Transmission Technology is combined to devise the generation of standard vibration source System provides the technological means of innovation for the realization of vibration measuring instrument automatic gauge verification system.
3) method for combining digital image understanding technology and hardware interface circuit to detect greatly improves measurement result reading The accuracy and resolving ability taken also provides the technology of innovation for the realization of vibration measuring instrument automatic gauge verification system Means.
4) combining wireless transmission technology, embedded system technology and computer programming devise computer program-control system System.The system is the calibrating measurement that vibration measuring instrument automatic gauge verification system realizes tested calibration regulatory requirements The acquisition of the output order and verification result of signal and the printout of probation redport.
Detailed description of the invention
Fig. 1 is the system block diagram of standard vibration source generating system;
Fig. 2 is measurement result identifying system;
Fig. 3 is computer program-control system block diagram;
Fig. 4 is vibration measuring instrument automatic gauge verification system functional block diagram;
Fig. 5 is ARM system design general frame figure;
Fig. 6 is ARM unit circuit diagram;
Fig. 7 is DSP unit circuit diagram;
Fig. 8 is HPI bus connection figure;
Fig. 9 is A/D Acquisition Circuit figure;
Figure 10 is that ARM extends out SDRAM;
Figure 11 is that ARM extends out FLASH;
Figure 12 is that DSP extends out FLASH and SDRAM;
Figure 13 is usb circuit figure of the present invention;
Figure 14 is RS232 interface circuit diagram;
Figure 15 is RJ45 interface circuit figure;
Figure 16 is power system circuit figure;
Figure 17 is programme-controlled signal source circuit diagram;
Figure 18 is motor control unit figure;
Figure 19 is image recognition principle block diagram.
Specific embodiment
The present invention is described further below with reference to embodiment.The explanation of following embodiments is merely used to help understand this Invention.It should be pointed out that for those skilled in the art, without departing from the principle of the present invention, also Can be with several improvements and modifications are made to the present invention, these improvement and modification also fall into the protection scope of the claims in the present invention It is interior.
The system includes three part compositions, including standard vibration source generating system, measurement result identifying system, computer journey Control system.Its system design principle block diagram is as shown in Figure 4.Implementation steps are as follows:
Computer program-control system requires to receive and dispatch by ZigBee according to each in the vertification regulation for being detected instrument first Waveform, frequency and the amplitude or standard vibration level size that unit sends corresponding measuring signal are to standard vibration source generating system, standard Vibration source generating system is believed standard vibration level according to requiring to generate corresponding standard vibration level signal or electric signal to tested instrument Number or the instruction that has sent of electric signal measurement result identifying system, measurement result identification system be sent to by ZigBee Transmit-Receive Unit After system receives the instruction, starting ARM system measures the identification of result, detects phase according to image recognition technology and circuit In conjunction with method, obtain correct measurement result, the result be then sent to computer program-control by ZigBee Transmit-Receive Unit System, computer program-control system obtain after related calibrating entry results through the host computer procedure write-in corresponding result of calibration certificate Column.Then start waveform, frequency and the amplitude of the measuring signal of next vertification regulation requirement to standard vibration source generating system.It is logical Cross constantly circulation, finally realize the calibrating project of all vertification regulation requirements, simultaneous computer stored program controlled according to these as a result, The calibration certificate for forming specification comes from printer output.The automatic gauge inspection of vibration measurement analysis instrument can thus be completed It is fixed.
ARM system in above system designs general frame using following system, as shown in Figure 5.The general frame was both ARM can be given full play to, the real-time control of system and good man-machine interface are realized and played DSP and image and audio are believed Number analysis processing capacity.
Wherein ARM is using S3C2410A as central control module.S3C2410A is the very high ARM of a cost performance core SOC.Dominant frequency with 203MHz has the processing speed requirement needed to audio signal digit detection, there is USB interface, It can be communicated with computer, there is multiple serial ports, the control to ZigBee Transmit-Receive Unit may be implemented.Physical circuit figure is as schemed Shown in 6:
Wherein DSP is designed using TMS320VC5502 chip, which is a high-performance, low-power consumption, fixed-point number Word signal processor fully meets audio signal and image signal process ability.Its physical circuit figure is as shown in Figure 7.
The connection of ARM unit and DSP unit uses HPI bus, and connection type is as shown in Figure 8.
The acquisition of vibration level signal is realized using high-speed a/d, select TMS320AIC23B, which uses advanced Σ mono- △ oversampling technique can provide the hits of 16bit, 20bit, 24bit and 32bit under the sample rate of 8kHz to 96kHz According to.The output signal-to-noise ratio of ADC and DAC is respectively up to 90dB and 100dB.The performance meets the data acquisition energy of audio signal enough Power.DSP to the control of A/D is realized by I2C interface.Its circuit is as seen in figure 9.
It is used for DSP and ARM and extends out FLASH and SDRAM, increase its program capacity and Nonvolatile data result Storage capacity.
SDRAM model K4S561632C-TC75, the capacity 4Banks*16Bit*4M that ARM is extended out, that is, 256Mbit converts into 32M byte.Memory of two pieces of K4S561632C-TC75 as ARM is used in this system, circuit theory is such as Shown in Figure 10.
S3C2410 extends out the K9F1208U0B that NAND FLASH is one piece of Samsung company, and capacity is 64M byte, electricity Road principle is as shown in figure 11.
TMS320VC5502 extends out FLASH and the circuit theory of SDRAM is as shown in figure 12.
ARM system realizes data exchange and includes tri- type of RJ45, RS232 and USB to the control of ZIGBEE wireless module Type interface, usb circuit principle are as shown in figure 13.
The RS232 serial ports controller used in this system is MAX3232, and circuit theory is as shown in figure 14.
The circuit of S3C2410A and RJ45 interface is comparatively more complicated, it is necessary to use 802.3 Ethernet of IEEE Controller is attached.We use CS8900A to this ethernet controller, this is also more commonly used one piece at present Ethernet control chip, circuit theory are as shown in figure 15.
There are four the power supplys used in ARM system, including 5V, 3.3V, 1.8V and 1.26V.Due to 3.3V and 1.26V electricity The order of pressure power supply requires, it is necessary to and 1.26V powers prior to 3.3V, therefore devises the circuit powered in order here, Use the power management chip of Liang Kuai TI company, TPS75733 (3.3V) and TPS76801 (1.26V).
Program control signal source circuit uses AD9850 for the circuit design of core, and physical circuit is as shown in figure 17:
Program control signal source unit ARM system inputs frequency control word to AD9850 and phase control words realize sine wave Frequency and phase controlling, while being controlled by AD7521 12 D/A amplitudes, 1/4096 amplitude precision control can be reached System, can realize the technical requirements of programme-controlled signal source in standard source generating system well.
The calculating of frequency control word is as follows:
If the frequency of output signal is CLKIN, the frequency control word of AD9850 is △ PHASE, then the relationship between three Are as follows:
△ PHASE=(fouT × 232)/CLKIN.
The calculating of phase control words is as follows:
There are 5 bits for phase controlling in AD9850, so, the precision of phase controlling is 360./ 25=11.25., with two System is expressed as 00001, according to actual needs, different phase control words is arranged, so that it may realize accurate phase controlling.Table 1 gives the corresponding relationship between phase shift and phase control words.
Corresponding relationship between 1 phase shift of table and phase control words
Phase shift (degree) 0 22.5 45 67.5
Phase control words 00000 00010 00100 0010
Phase shift (degree) 180 202.5 225 247.5
Phase control words 10000 10010 10100 10110
Phase shift (degree) 90 112.5 135 157.5
Phase control words 01000 01010 01100 01110
Phase shift (degree) 270 292.5 315 337.5
Phase control words 11000 11010 111000 11110
The frequency/phase control word one of AD9850 shares 40 bits, wherein 32 bits are frequency control word, 5 bits are phase Position control word, 1 bit are power supply suspend mode control, and last 2 bit is work method control, control 1 power supply suspend mode in, 2 work method controls are set as " 000 ".The every allocation table of AD9850 parallel mode loading frequency/phase control words such as table 2 It is shown.
Every allocation table of 2 frequency/phase control word of table
Control word Data[7] Data[6] Data[5] Data[4]
W0 Phase—b4(MSB) Phase—b3 Phase—b2 Phase—bl
Wl Freq—b31(MSB) Freq—b30 Freq—b29 Freq—b28
W2 Freq—b23 Freq—b22 Freq—b21 Freq—b20
W3 Freq—b15 Freq—b14 Freq—b13 Freq—b12
W4 Freq—b7 Freq—b6 Freq—b5 Freq—b4
Control word Data[3] Data[2] Data[1] Data[0]
WO Phase—bO(LsB) Power D0wn Control (=0) Control (=0)
Wl Freq—b27 Freq—b26 Freq—b25 Freq—b24
W2 Freq—b19 Freq—b18 Freq—b17 Freq—b16
W3 Freq—b11 Freq—b10 Freq—b9 Freq—b8
W4 Freq—b3 Freq—b2 Freq—bl Freq—bo(LsB)
In loaded in parallel mode, by 8 position datawires, continuous 5 times AD9850 is written in 40 bit datas by load every time, WCLK and FQUD is used to determine address and load data sequence, and Wx (X=0,1,2,3,4) is written in the rising edge of WCLK, so The rising edge of WCLK, DATA should be ready to and keep stable, and the rising edge of FQUD is by 40 bit data write frequencies/phase data Register, while address pointer is directed toward first register WO.When initialization, high level is set by FQUD, WCLK is set as Low level, parallel mode writing process are as follows: firstly, FQUD switchs to low level by high level, AT89C51 exports data WO; Later, control WCLK switchs to high level by low level, then switchs to low level by high level, writes control word WO, AT89C51 at this time According to the process of write-in WO, it is sequentially written in W1, W2, W3, W4;Finally, AT89C51 controls FQUD, it is made to switch to height by low level Level, completes the write-in of 40 bit datas, and address pointer is directed toward W0, performs standard for next write frequency/phase control words It is standby.
Motor control unit controls multiple LN298N units using ARM and realizes that the rotation of multiple motors controls multiple strikers pair The control of the shock dynamics and velocity interpolation vibration level amplitude and speed of shake table, single LN298N element circuit figure such as Figure 18 institute Show.
The frequency of impact and width that motor drives striker are controlled by the frequency that ARM system controls the PWN wave of L298N input Degree control, realizes the control of vibration level amplitude and speed.
ZigBee Transmit-Receive Unit uses CC2420 for the ZigBee Transmit-Receive Unit module of core chips.Using ZigBee networking Technology realizes the wireless communication of three points well.Its channel from way of search, the topological structure of StarNet and entirely shake hands (Fully handshaked) agreement ensure that the stability of data transmission.Two-way is realized using RS232 interface with ARM system Letter.
It can be realized well in vibration measuring instrument automatic gauge verification system by the design of said units circuit The function of three systems (standard vibration source generating system, measurement result identifying system, computer program-control system).For three systems Circuit design be specifically described it is as follows:
Standard vibration source generating system: RS232 interface is controlled by ARM system (Fig. 5) and is realized to ZigBee transceiver module Communication receives the specific requirement of the standard vibration level signal or standard electric signal of the requirement output that computer program-control system is sent, On the one hand pass through ARM system control AD9850 unit (Figure 17) realization standard electric signal corresponding signal amplitude, frequency, phase Standard electric signal interface is given in control output, on the other hand controls motor-drive circuit (Figure 18) by ARM system, realizes multiple hit For the amplitude and frequency of shake table vibratory impulse, for the accurate control of amplitude and frequency that motor control controls multiple strikers System, is connected on shake table by standard vibration sensor and is acquired, by the high speed of the signal input ARM system after acquisition A/DTMS320AIC23B carries out digital detection analysis, and whether analysis vibration level reaches the requirement of standard, if do not reached, pass through The control for changing ARM system control motor-drive circuit, realizes the Oscillation Amplitude and frequency of shake table, by constantly feeding back, most Vibration level signal output definitely is realized eventually.The circuit system has following innovative design: (1) accurate program controlled signal source realizes Optional frequency control, the accurate amplitude to 1/4096 in multiple waveforms, audio range control contour performance signal source and fill The design requirement set.(2) by vibration level self calibration feedback technique, the output of standard vibration level signal is realized.(3) pass through ZigBee Networking technology, the multiple spot for realizing standard acoustic signal generating system and computer program-control system and measurement result identifying system are logical Letter.(4) rigorous analysis of vibration level signal is realized by digital detection technology.Digital detection includes vertical (W.B.z) weighted of whole body The vibration level digitalization index of vibration level, the digitlization of whole body level (W.B.x-y) weighted vibration level and AC signal is analyzed.
Vertical (W.B.z) the weighted vibration level of whole body, the digitized process of whole body level (W.B.x-y) weighted vibration level are as follows:
The 1/3oct analysis for completing 1Hz to 80Hz first, carries out according to following steps:
(1) the fertile hereby bandpass filter of a 8 rank Barts is designed with the design method of analog filter;
(2) predistortion is carried out to analog frequency;
(3) transmission function of digital filtering is found out using Bilinear transformation method;
(4) inputs a sampled value, calculates an output valve, seeks the root mean square of each output valve.
The frequency band root mean square acceleration analyzed by 1/3oct is multiplied by the weighted factor, then takes r.m.s. and obtain phase The spectrum overlapping acceleration answered:
In formula: ai--- the accelerated speed effective value (m/s that the i-th frequency range is surveyed in third-octave frequency spectrum2);
ki--- the corresponding spectrum overlapping factor of the i-th frequency range in third-octave frequency spectrum.
It is analyzed by the calculating of above formula and 1/3oct, so that it may measure corresponding spectrum overlapping acceleration.
Weighted vibration level can be obtained as the following formula by weighted acceleration:
In formula: aw--- weighted acceleration virtual value (m/s2);
a0--- reference acceleration, a0=10-6m/s2
It is calculated by the digital collection of (2) formula, so that it may corresponding weighted vibration level.
Measurement result identifying system: controlling communication of the RS232 interface realization to ZigBee transceiver module by ARM system, After receiving the ready order of signal that standard acoustic signal generating system method is come.ARM system completes following work: (1) The AC interface that tested instrument provides is docked to high-speed a/d TMS320AIC23B acquisition interface and carries out digital detection analysis.(2) will The built-in A/D translation interface that the DC interface that tested instrument provides is connected to ARM2410 offer carries out DC analysis.(3) by tested instrument The measurement result display screen of device passes through I by digital camera2C bus realizes the reading of ARM system, passes through optical character identification Technology combination AC interface and DC interface analysis are as a result, measure the identification of result.
Algorithm is analyzed as follows in the identification process:
S1 passes through AC interface first and carries out vertical (W.B.z) the weighted vibration level of vibration level signal whole body, whole body level (W.B.x-y) The digital assay of weighted vibration level, specific as follows:
The 1/3oct analysis for completing 1Hz to 80Hz first, carries out according to following steps:
1) a 8 rank Butterworth filter devices are designed with the design method of analog filter;
The low pass chi square function of butterworth filter is expressed as follows:
Interim N takes 8,8 rank low-pass filters is designed, to obtain system function H (S).
2) predistortion is carried out to analog frequency;
It is being that, in order to guarantee that each edge frequency point is the frequency pre-established, filter is cut using two-wire row converter technique Only frequencies omegaCWith stopband side frequency ω s, it is pre- abnormal before determining simulation low-pass filter system function must to carry out as the following formula frequency Become:
H (s)=Ha(s/Ωc) (5)
3) transmission function of digital filtering is found out using Bilinear transformation method;
It is summarized as follows with the process of Bilinear transformation method design digital filter:
Normalize the frequency transformation of analog low-pass prototype to digital band-pass filter:
4) inputs a sampled value again, calculates an output valve, seeks the root mean square of each output valve.
5) the frequency band root mean square acceleration analyzed by 1/3oct is multiplied by the weighted factor, then takes r.m.s. and obtain Corresponding spectrum overlapping
Acceleration:
In formula: ai--- the accelerated speed effective value (m/s2) that the i-th frequency range is surveyed in third-octave frequency spectrum;
ki--- the corresponding spectrum overlapping factor of the i-th frequency range in third-octave frequency spectrum.
It is analyzed by the calculating of above formula and 1/3oct, so that it may measure corresponding spectrum overlapping acceleration.
Weighted vibration level can be obtained as the following formula by weighted acceleration:
In formula: aw--- weighted acceleration virtual value (m/s2);
a0--- reference acceleration, a0=10-6m/s2
It is calculated by the digital collection of (10) formula, so that it may corresponding weighted vibration level.
Corresponding character match module is formed by these measurement results.
S2 obtains the measurement result of corresponding not weighted through ARM system sampling analysis by DC interface, which forms The respective symbols matching module of not weighted in optical identification measurement result.
The measurement result display screen of tested instrument is realized ARM system by RS232 interface by digital camera by S3 It reads, carries out optical identification, the algorithm such as Figure 19.
It identifies that the image of object obtains and gray processing is realized by camera combination electric image signal unit, know for character Not, select median filtering technology realize image smoothing, by threshold value set realize image binaryzation, from binaryzation data into Line character segmentation, the method for selection is template matching method, then carries out character knowledge by extracting after image thinning using framework characteristic Not.It is extracted using framework characteristic, can be reduced corresponding calculation amount, only in the character feature and above-mentioned S1, S2 step extracted The feature of the measurement result Character mother plate of formation is matched.An image to be identified is given, can be by refinement The character skeleton image of processing, be also possible to not refine eliminates the bianry image of noise, horizontal and vertical all with fixation Interval scribing line, then successively write down the number across character picture, that is, number of crossings.Template image is gathered in advance, Extract each template passes through number, stores with a matrix type, as eigenmatrix, extracted the feature of each character with Afterwards, it is preserved as database.After obtaining single character picture to be identified, traversing times are equally extracted, are loaded into The traversing times matrix of template data and template, successively makes comparisons, and calculates related coefficient, maximum related coefficient is exactly to identify As a result.In the design is long and width, is divided into five parts, is unable to equal part, and round numbers part is calculated in each position of halving Traversing times intersect each time, are denoted as and once pass through, and share ten characteristic values in this way.Calculate poor method for distinguishing:
X can represent binary character to be identified in formula, and P is the eigenmatrix of template, and obtained matrix of differences is D.D= d11 2+…+d25 2For the number that the two correlation is measured, that maximum corresponding templates of the number mean that recognition result.
In above process, corresponding character mould is obtained in conjunction with the measurement data that AC interface, DC interface measure acquisition Block carries out module matching, can play the role that the fast area positioning (2) of (1) measurement result makes up image recognition processing The slow and not accurate enough disadvantage of speed.By the algorithm, it is able to carry out measurement result identification fast and accurately.
After identification 30 seconds, if there is image recognition region still no-fix, illustrate the AC interface or DC interface of instrument The interface signal of offer is wrong, and the wrong conclusion of measurement result is sent to computer program-control system by ZigBee module, It determined AC interface or whether DC interface be correct.If correct identification, correct recognition result is sent to by ZigBee module Computer program-control system.
There is the circuit system following innovative design to be embodied in by passing through image on the basis of AC, DC analysis combination Also tested instrument is able to verify that while the shortcomings that identification technology can quickly and accurately find measurement result, compensate for image recognition The correctness of AC, DC interface of offer.
Computer program-control system: controlling communication of the RS232 interface realization to ZigBee transceiver module by ARM system, real Now the reading of the measurement result of measurement result identifying system is tied measurement by the USB interface or RJ45 interface of ARM system Fruit is read into computer system, and measurement result is written in the detection corresponding data base entries of regulation by computer system It goes.Then the standard signal output for starting next detection regulation requires to be sent to the generation of standard sound by ZigBee group net unit System.Computer programming realizes above-mentioned function using VB language combination SQL database technology.The innovative major embodiment of the circuit Realize that computer program-control system sends out standard sound by ZigBee networking technology and computer programming in computer program-control The two-way communication of raw system and measurement result identifying system and the formation and output of calibration certificate.

Claims (3)

1.一种振动测量仪器自动化计量检定系统,其特征在于:包括标准振源发生系统、测量结果识别系统和计算机程控系统;1. a vibration measuring instrument automatic metrology verification system is characterized in that: comprise standard vibration source generation system, measurement result identification system and computer program control system; 标准声源发生系统包括:ZigBee收发单元、ARM系统、程控信号源、电机控制单元和振动台;ARM系统通过Zigbee收发单元接收计算机程控系统发来的要求输出的标准电信号的信号波形、频率和幅度或标准振级的大小,一路通过系统电路进行程控信号输出,另一路通过控制电机控制单元输出撞击振动台的频率和幅度实现振级的输出,同时通过ARM系统对振级信号的信号频率及幅度分析,看是否达到标准振级的要求,如果达不到,调节ARM系统对电机控制单元输出频率以及相关电机的控制输出,实现电机控制单元引导撞针对振动台的撞击力度及频率,直到标准振级的输出;实现标准振级的自校准过程;程控信号源通过ARM系控制AD9850通过直接数字频率合成技术输出相应的信号波形、幅度和频率;当输出标准的振级信号或电信号时,通过ZigBee收发单元发送信号准备就绪命令给测量结果识别系统;The standard sound source generation system includes: ZigBee transceiver unit, ARM system, program-controlled signal source, motor control unit and vibration table; ARM system receives the signal waveform, frequency and frequency of the standard electrical signal sent by the computer program control system through the Zigbee transceiver unit. Amplitude or standard vibration level, one channel outputs program-controlled signals through the system circuit, and the other channel outputs the vibration level by controlling the motor control unit to output the frequency and amplitude of the impact vibration table. Amplitude analysis to see if it meets the requirements of the standard vibration level. If not, adjust the output frequency of the ARM system to the motor control unit and the control output of the related motors, so as to realize the impact strength and frequency of the motor control unit to guide the impact on the vibration table until the standard The output of vibration level; the self-calibration process of standard vibration level is realized; the program-controlled signal source controls AD9850 through ARM system to output corresponding signal waveform, amplitude and frequency through direct digital frequency synthesis technology; when outputting standard vibration level signal or electrical signal, Send a signal ready command to the measurement result identification system through the ZigBee transceiver unit; 测量结果识别系统包括:被测仪器显示屏、摄像头、图像电信号接口单元、ARM系统、检测电路、AC接口、DC接口和ZigBee收发单元;ARM系统通过Zigbee收发单元接收到标准振级发生系统发来的准备就绪的命令,ARM系统通过图像电信接口单元实现摄像头对被测仪器显示屏的测量结果进行图像识别,通过检测电路对被测仪器的AC接口、DC接口进行测量;ARM系统获取到测量结果后通过ZigBee收发单元将测量结果发送到计算机程控系统;The measurement result identification system includes: the display screen of the instrument under test, the camera, the image electrical signal interface unit, the ARM system, the detection circuit, the AC interface, the DC interface and the ZigBee transceiver unit; the ARM system receives the standard vibration level generation system through the Zigbee transceiver unit. When the ready command comes, the ARM system realizes image recognition of the measurement results of the display screen of the instrument under test through the image telecommunication interface unit, and measures the AC interface and DC interface of the instrument under test through the detection circuit; the ARM system obtains the measurement results. After the result, the measurement result is sent to the computer program control system through the ZigBee transceiver unit; 计算机程控系统包括:ZigBee收发单元、ARM系统和计算机上位机;ARM系统通过ZigBee收发单元获取测量结果识别系统的测量结果再通过USB接口将数据送到程控计算机的上位机程序,上位机程序将测量结果送到预置的检定报告的输出条目中,然后,通过USB接口发送下一条检定规程中要求的输出信号源要求的波形、幅度和频率到测量结果识别系统,测量结果识别系统再通过ZigBee收发单元将输出信号源要求发送到标准振源发生系统。The computer program control system includes: ZigBee transceiver unit, ARM system and computer host computer; the ARM system obtains the measurement results through the ZigBee transceiver unit and identifies the measurement results of the system, and then sends the data to the host computer program of the program control computer through the USB interface, and the host computer program will measure the measurement results. The result is sent to the output entry of the preset verification report, and then, the waveform, amplitude and frequency required by the output signal source required in the next verification procedure are sent to the measurement result identification system through the USB interface, and the measurement result identification system is then sent and received through ZigBee The unit sends the output signal source request to the standard source generator system. 2.一种如权利要求1所述的振动测量仪器自动化计量检定系统的测量结果识别系统的识别算法,其特征在于,包括以下步骤:2. the identification algorithm of the measurement result identification system of the vibration measuring instrument automatic metrological verification system as claimed in claim 1, is characterized in that, comprises the following steps: S1.首先通过AC接口进行振级信号全身垂向(W.B.z)计权振级、全身水平(W.B.x-y)计权振级的数字化分析,具体如下:S1. First, perform the digital analysis of the vibration level signal whole body vertical (W.B.z) weighted vibration level and whole body horizontal (W.B.x-y) weighted vibration level through the AC interface, as follows: 首先完成1Hz到80Hz的1/3oct分析,根据以下几个步骤进行:First complete the 1/3 oct analysis from 1Hz to 80Hz, according to the following steps: 1)用模拟滤波器的设计方法设计一个8阶巴特沃思低通滤波器;1) Design an 8th-order Butterworth low-pass filter using the analog filter design method; 巴特沃兹滤波器的低通平方函数表示如下:The low-pass squared function of the Butterworth filter is expressed as follows: 期中N取8,设计8阶低通滤波器,从而得到系统函数H(S);In the period, N is taken as 8, and an 8th-order low-pass filter is designed to obtain the system function H(S); 2)对模拟频率进行预畸变;2) Pre-distort the analog frequency; 在利用双线行变换法是,为了保证各边界频率点为预先制定的频率,滤波器的截止频率ωC与阻带边频ωs,在确定模拟低通滤波器系统函数之前必须按下式进行频率预畸变:When using the bilinear line transformation method, in order to ensure that each boundary frequency point is a predetermined frequency, the cut-off frequency ω C of the filter and the stop-band side frequency ωs must be determined by the following formula before the analog low-pass filter system function. Frequency Predistortion: H(s)=Ha(s/Ωc) (5)H(s)=H a (s/Ω c ) (5) 3)利用双线性变换法求出数字滤波的传递函数;3) Use the bilinear transformation method to obtain the transfer function of digital filtering; 用双线性变换法设计数字滤波器的过程归纳如下:The process of designing a digital filter using the bilinear transformation method is summarized as follows: 归一化模拟低通原型到数字带通滤波器的频率变换:Frequency transformation of normalized analog low-pass prototype to digital band-pass filter: 4)再.输入一个采样值,计算出一个输出值,求各输出值的均方根;4) Then. Input a sample value, calculate an output value, and find the root mean square of each output value; 5)通过1/3oct分析得到的频带均方根加速度乘以计权因子,再取方均根和,获得相应的频率计权加速度:5) Multiply the frequency band root mean square acceleration obtained by 1/3oct analysis by the weighting factor, and then take the root mean square sum to obtain the corresponding frequency-weighted acceleration: 式中:ai——1/3倍频程频谱中第i频段实测的加速度有效值(m/s2);In the formula: a i ——the actual measured acceleration effective value of the i-th frequency band in the 1/3 octave frequency spectrum (m/s 2 ); ki——1/3倍频程频谱中第i频段相应的频率计权因数;k i ——the frequency weighting factor corresponding to the i-th frequency band in the 1/3 octave frequency spectrum; 通过上式的计算以及1/3oct分析,就可以测得相应的频率计权加速度;Through the calculation of the above formula and 1/3 oct analysis, the corresponding frequency-weighted acceleration can be measured; 由计权加速度可以按下式得到计权振级:From the weighted acceleration, the weighted vibration level can be obtained as follows: 式中:aw——计权加速度有效值(m/s2);where: a w — effective value of weighted acceleration (m/s 2 ); a0——参考加速度,a0=10-6m/s2a 0 ——reference acceleration, a 0 =10 -6 m/s 2 ; 通过(10)式的数字化采集计算,就可以相应的计权振级;Through the digital acquisition and calculation of formula (10), the corresponding weighted vibration level can be obtained; 通过这些测量结果形成相应的字符匹配模块;The corresponding character matching modules are formed through these measurement results; S2通过DC接口,经ARM系统采样分析,获得相应的不计权的测量结果,该结果形成光学识别测量结果中的不计权的相应字符匹配模块;S2 obtains the corresponding unweighted measurement result through the DC interface and the ARM system sampling and analysis, and the result forms the unweighted corresponding character matching module in the optical recognition measurement result; S3将被测仪器的测量结果显示屏通过数字摄像头通过RS232接口实现ARM系统的读取,进行光学识别。S3 realizes the reading of the ARM system through the digital camera through the RS232 interface on the display screen of the measurement results of the instrument under test, and performs optical identification. 3.根据权利要求2所述的测量结果识别系统的识别算法,其特征在于:所述步骤S3的光学识别的算法为:3. the identification algorithm of the measurement result identification system according to claim 2, is characterized in that: the algorithm of the optical identification of described step S3 is: 识别对象的图像获取通过摄像头结合图像电信号单元实现灰度化,针对字符识别,选用中值滤波技术实现图像平滑,通过阈值设定实现图像的二值化,从二值化数据中进行字符分割,选用的方法是模板匹配法,再通过图像细化后采用骨架特征提取进行字符识别;采用骨架特征提取,减少相应的计算量,只把提取出的字符特征与上述S1、S2步骤中形成的测量结果字符模板的特征进行匹配即可;给定一个待识别的图像,经过细化处理的字符骨架图像,或没有细化的去除了噪声的二值图像,横向和纵向都以固定的间隔划线,然后依次记下穿过字符图像的次数,也就是交叉次数;事先收集好模板图像,提取出每个模板的穿过次数,以矩阵的形式存放,作为特征矩阵,提取完每个字符的特征以后,作为数据库保存下来;获取了待识别的单个字符图像以后,同样提取出穿越次数,载入模板数据,和模板的穿越次数矩阵,依次作比较,计算相关系数,相关系数最大的就是识别结果;长和宽等分成五份,不能等分的,取整数部分,在每个等分位置计算穿越次数,每一次交叉,记为一次穿越,这样共有十个特征值;计算差别的方法:The image of the recognized object is obtained by combining the camera with the image electrical signal unit to achieve grayscale. For character recognition, median filtering technology is used to achieve image smoothing, and threshold setting is used to achieve image binarization, and character segmentation is performed from the binarized data. , the selected method is the template matching method, and then the skeleton feature extraction is used for character recognition after image refinement; the skeleton feature extraction is used to reduce the corresponding calculation amount, and only the extracted character features are compared with those formed in the above S1 and S2 steps. The features of the character template of the measurement result can be matched; given an image to be recognized, a character skeleton image that has been refined, or a binary image that has not been refined with noise removed, the horizontal and vertical lines are drawn at fixed intervals. Line, and then record the times of crossing the character image in turn, that is, the number of crossings; collect the template images in advance, extract the crossing times of each template, store them in the form of a matrix, and use them as a feature matrix. After the feature, it is saved as a database; after obtaining the single character image to be recognized, the number of crossings is also extracted, the template data is loaded, and the matrix of crossing times of the template is compared in turn, and the correlation coefficient is calculated. The largest correlation coefficient is recognition. Result; the length and width are divided into five equal parts, and if they cannot be divided equally, the integer part is taken, and the number of crossings is calculated at each equal position. 式中X可以代表待识别的二值字符,P为模板的特征矩阵,得到的差值矩阵为D;d=d11 2+…+d25 2为两者相关性衡量的数,该数最大的那个对应模板就表示识别结果。In the formula, X can represent the binary character to be recognized, P is the feature matrix of the template, and the obtained difference matrix is D; d=d 11 2 +...+d 25 2 is the number measured by the correlation between the two, and the number is the largest The corresponding template of is the recognition result.
CN201910199211.5A 2019-03-15 2019-03-15 Vibration measuring instrument automated metrology calibration system and identification algorithm Expired - Fee Related CN109870231B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910199211.5A CN109870231B (en) 2019-03-15 2019-03-15 Vibration measuring instrument automated metrology calibration system and identification algorithm

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910199211.5A CN109870231B (en) 2019-03-15 2019-03-15 Vibration measuring instrument automated metrology calibration system and identification algorithm

Publications (2)

Publication Number Publication Date
CN109870231A true CN109870231A (en) 2019-06-11
CN109870231B CN109870231B (en) 2023-09-26

Family

ID=66920520

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910199211.5A Expired - Fee Related CN109870231B (en) 2019-03-15 2019-03-15 Vibration measuring instrument automated metrology calibration system and identification algorithm

Country Status (1)

Country Link
CN (1) CN109870231B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112888013A (en) * 2020-12-08 2021-06-01 福建星网天合智能科技有限公司 System and method for testing consistency of zigbee radio frequency parameters

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102288278A (en) * 2011-07-23 2011-12-21 浙江大学 Automatic calibrating device of digital display type vibration measurer and calibrating method thereof
CN102506994A (en) * 2011-11-21 2012-06-20 嘉兴中科声学科技有限公司 Digitized acoustic detection system
JP2012208043A (en) * 2011-03-30 2012-10-25 Railway Technical Research Institute Method and device for identifying vibration characteristic of railroad structure
CN102829941A (en) * 2012-08-16 2012-12-19 杭州亿恒科技有限公司 Automated method for verifying vibration system
CN105334460A (en) * 2015-11-27 2016-02-17 浙江大学城市学院 Machine running state online monitoring analysis system based on noise and vibration analysis
CN105606207A (en) * 2016-01-29 2016-05-25 湖南省计量检测研究院 Calibration apparatus and calibration method applied to vibration sensor
CN106404160A (en) * 2016-06-30 2017-02-15 国网重庆市电力公司电力科学研究院 Power transmission line monitoring device metering system and method in field environment
CN106645949A (en) * 2016-09-26 2017-05-10 武汉大学 Heterodyne sweep-frequency type spectrum analyzer based on low frequency detection

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012208043A (en) * 2011-03-30 2012-10-25 Railway Technical Research Institute Method and device for identifying vibration characteristic of railroad structure
CN102288278A (en) * 2011-07-23 2011-12-21 浙江大学 Automatic calibrating device of digital display type vibration measurer and calibrating method thereof
CN102506994A (en) * 2011-11-21 2012-06-20 嘉兴中科声学科技有限公司 Digitized acoustic detection system
CN102829941A (en) * 2012-08-16 2012-12-19 杭州亿恒科技有限公司 Automated method for verifying vibration system
CN105334460A (en) * 2015-11-27 2016-02-17 浙江大学城市学院 Machine running state online monitoring analysis system based on noise and vibration analysis
CN105606207A (en) * 2016-01-29 2016-05-25 湖南省计量检测研究院 Calibration apparatus and calibration method applied to vibration sensor
CN106404160A (en) * 2016-06-30 2017-02-15 国网重庆市电力公司电力科学研究院 Power transmission line monitoring device metering system and method in field environment
CN106645949A (en) * 2016-09-26 2017-05-10 武汉大学 Heterodyne sweep-frequency type spectrum analyzer based on low frequency detection

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
王志鹏;何闻;于梅;沈润杰;: "超低频振动校准自动控制系统的研究", 机电工程, no. 06, pages 693 - 697 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112888013A (en) * 2020-12-08 2021-06-01 福建星网天合智能科技有限公司 System and method for testing consistency of zigbee radio frequency parameters
CN112888013B (en) * 2020-12-08 2022-06-28 福建星网天合智能科技有限公司 System and method for testing radio frequency parameter consistency of zigbee

Also Published As

Publication number Publication date
CN109870231B (en) 2023-09-26

Similar Documents

Publication Publication Date Title
CN104198829B (en) Electromagnetic radiation measurement device and measurement method based on ARM band temperature and humidity self-correction
CN105415191B (en) A kind of grinding status investigating method and device based on sound emission
CN105425128B (en) A kind of partial discharge of transformer ultrasound examination and accurate positioning device and method
CN103954903B (en) A kind of can real-time resolving multi-mode output circuit test system
CN106526332A (en) System and method for detecting capacitor of scanning line of capacitive touch screen
CN201083539Y (en) Portable digital sound and vibration analyzer
CN114676004B (en) A dynamic index measurement system and method of ADC chip based on FPGA
CN106771589A (en) A kind of implementation method of the real-time Waterfall plot of spectrum analysis
CN112612656A (en) Online debugging system for digital nuclear pulse processor
CN104407510A (en) Method and device for time service
CN106645942A (en) Low cost and high precision embedded type signal collection and analysis system and method
CN109870231A (en) Automatic Metrology Verification System and Identification Algorithm of Vibration Measuring Instrument
CN105955918B (en) An experimental universal instrument system
CN201047798Y (en) Virtual vibration table detecting apparatus
CN109101419A (en) A kind of test method, device, equipment and computer readable storage medium
CN109935239B (en) Automatic metering verification system and recognition algorithm for audio analysis instrument
CN105842566B (en) The system and its mapping method of automation mapping digital filter amplitude-versus-frequency curve
CN203164407U (en) Ultrahigh-frequency partial discharge state detection instrument calibrating device
CN118859310A (en) A method and device for predicting earthquake source parameters based on deep learning
CN106291753B (en) A kind of Special testing device of half aviation Transient Electromagnetic Receiver
CN209559329U (en) Automatic Metrology Verification System for Vibration Measuring Instruments
CN104568137B (en) The decision method of specific noise source contribution rate in noise testing
CN105891907A (en) Portable piezoelectric detector multi-parameter testing system and testing method
CN201199265Y (en) Full-automatic calibrating apparatus for voltage monitoring instrument
CN109061717A (en) A kind of one key of nuclear radiation Calibration Experiment room calibrates control system and its control method

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20230926

CF01 Termination of patent right due to non-payment of annual fee