KR20120068181A - A signal process method for improving s/n of ion cyclotron resonance mass spectrometer - Google Patents

A signal process method for improving s/n of ion cyclotron resonance mass spectrometer Download PDF

Info

Publication number
KR20120068181A
KR20120068181A KR1020100129680A KR20100129680A KR20120068181A KR 20120068181 A KR20120068181 A KR 20120068181A KR 1020100129680 A KR1020100129680 A KR 1020100129680A KR 20100129680 A KR20100129680 A KR 20100129680A KR 20120068181 A KR20120068181 A KR 20120068181A
Authority
KR
South Korea
Prior art keywords
signal
noise
value
data
cyclotron resonance
Prior art date
Application number
KR1020100129680A
Other languages
Korean (ko)
Other versions
KR101162178B1 (en
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 KR1020100129680A priority Critical patent/KR101162178B1/en
Priority to PCT/KR2011/008230 priority patent/WO2012081825A1/en
Publication of KR20120068181A publication Critical patent/KR20120068181A/en
Application granted granted Critical
Publication of KR101162178B1 publication Critical patent/KR101162178B1/en

Links

Images

Classifications

    • HELECTRICITY
    • H01ELECTRIC ELEMENTS
    • H01JELECTRIC DISCHARGE TUBES OR DISCHARGE LAMPS
    • H01J49/00Particle spectrometers or separator tubes
    • H01J49/0027Methods for using particle spectrometers
    • H01J49/0036Step by step routines describing the handling of the data generated during a measurement
    • HELECTRICITY
    • H01ELECTRIC ELEMENTS
    • H01JELECTRIC DISCHARGE TUBES OR DISCHARGE LAMPS
    • H01J49/00Particle spectrometers or separator tubes
    • H01J49/26Mass spectrometers or separator tubes
    • H01J49/34Dynamic spectrometers
    • H01J49/36Radio frequency spectrometers, e.g. Bennett-type spectrometers, Redhead-type spectrometers
    • H01J49/38Omegatrons ; using ion cyclotron resonance
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N27/00Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
    • G01N27/62Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating the ionisation of gases, e.g. aerosols; by investigating electric discharges, e.g. emission of cathode
    • G01N27/622Ion mobility spectrometry
    • G01N27/623Ion mobility spectrometry combined with mass spectrometry
    • HELECTRICITY
    • H01ELECTRIC ELEMENTS
    • H01JELECTRIC DISCHARGE TUBES OR DISCHARGE LAMPS
    • H01J49/00Particle spectrometers or separator tubes
    • H01J49/02Details
    • H01J49/022Circuit arrangements, e.g. for generating deviation currents or voltages ; Components associated with high voltage supply

Landscapes

  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Physics & Mathematics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Molecular Biology (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Electrochemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Other Investigation Or Analysis Of Materials By Electrical Means (AREA)

Abstract

PURPOSE: A method for processing an ion cyclotron resonance mass spectrograph is provided to reduce noise by performing a noise reduction operation process in regard to a noise component. CONSTITUTION: An analog signal obtained by a digitizer is transformed to a digital signal(S10). A DC voltage included in a signal component is eliminated(S20). Loss in regard to the signal component is reduced through a differentiation windowing signal processing process(S30). FFT calculation in regard to data of a time area is operated(S40). A noise component which is increased in the process that the FFT calculation is operated is eliminated(S50). A frequency area signal of the data in which noise is reduced is changed into a mass-to-electron ratio area(S60).

Description

신호 대 잡음 개선을 위한 이온 싸이클로트론 공명 질량 분석기 신호처리방법{A Signal Process Method for improving S/N of Ion Cyclotron Resonance Mass Spectrometer} A Signal Process Method for improving S / N of Ion Cyclotron Resonance Mass Spectrometer

본 발명은 이온 싸이클로트론 공명 질량 분석기 신호처리방법에 관한 것이며, 더욱 상세하게는 노이즈 저감 알고리즘을 사용하여 노이즈 성분을 감쇄하도록 하는 신호 대 잡음 성분 개선을 위한 이온 싸이클로트론 공명 질량분석기의 신호처리방법에 관한 것이다.
The present invention relates to an ion cyclotron resonance mass spectrometer signal processing method, and more particularly, to a signal processing method of an ion cyclotron resonance mass spectrometer for signal-to-noise component improvement to reduce noise components using a noise reduction algorithm. .

일반적인 이온 싸이클로트론 공명 질량 분석기의 제어장치에 대하여 첨부된 도1을 참고하여 상세히 설명하면 다음과 같다.A control device of a general ion cyclotron resonance mass spectrometer will be described in detail with reference to FIG. 1.

도 1은 일반적인 이온 싸이클로트론 공명 질량 분석기 제어장치의 전체 구성도이다.1 is an overall configuration diagram of a general ion cyclotron resonance mass spectrometer control device.

일반적인 이온 싸이클로트론 공명 질량 분석기는 도 1에 도시된 바와 같이, 주입된 시료를 이온화시켜 방출하는 시료주입 이온화부(1)와, 상기 시료주입 이온화부(1)로부터 방출된 이온을 전송하는 제1이온 전송부(2), 상기 제1이온 전송부(2)를 통해 전송된 이온을 특정 목적에 따라 선택하거나 분리하여 방출하는 이온선택(분리)부(3), 상기 이온선택(분리)부(3)에 의해 선택되거나 분리된 이온을 충돌 가스와 충돌시켜 더 작은 크기로 나누어 방출하는 이온 충돌부(4), 상기 이온 충돌부(4)에 의해 나누어진 이온을 전송하는 제2이온 전송부(5), 상기 제2이온 전송부(5)를 통해 전송된 이온을 이온 트랩에 수집한 후 특정 목적에 부합하는 이온의 질량을 나타내는 전기적 신호를 검출하는 이온트랩(6), 이온트랩(6)에서 신호 검출은 컴퓨터(10)의 제어 프로그램에 의하여 임의파형을 발생시키는 임의파형발생부(AWG)(8)와 발생된 임의 파형을 증폭시키는 고주파증폭기(RF Amp)(7), 증폭된 파형을 이온트랩(6)에 인가하여 이온을 여기시킨다. As shown in FIG. 1, a typical ion cyclotron resonance mass spectrometer includes a sample injection ionizer 1 for ionizing and releasing an injected sample, and a first ion for transmitting ions emitted from the sample injection ionizer 1. An ion selector (separation) unit 3 for selecting or separating and transmitting ions transmitted through the transmission unit 2 and the first ion transport unit 2 according to a specific purpose, and the ion selector (separation) unit 3 Ion bombardment part 4 which collides the ion selected or separated by impingement with the collision gas and emits it into smaller sizes, and a second ion transmitter 5 which transmits the ions divided by the ion bombardment part 4. In the ion trap 6 and the ion trap 6, which collects ions transmitted through the second ion transmitter 5 to an ion trap and detects an electrical signal representing a mass of ions meeting a specific purpose. Signal detection is carried out by the control program of the computer 10. An arbitrary waveform generator (AWG) 8 for generating an arbitrary waveform, an RF amplifier 7 for amplifying the generated arbitrary waveform, and an amplified waveform to the ion trap 6 to excite ions. .

여기된 신호는 이온트랩의 또 다른 전극을 통하여 도 1에서 나타나는 전치증폭기(Pre Amp)(11)를 통과하여 검출하기 적당한 신호크기로 증폭하여 디지타이저(A/D)(10)를 통과하여 디지털 신호가 되어 컴퓨터에서 신호처리가 이루어진다. 각 단계별로 이루어지는 제어는 컴퓨터에서 운영되는 제어프로그램에 의하여 제어 신호가 제어부(12)를 통하여 각 단계의 제어 장치에 신호가 전달된다.The excited signal is passed through the digitizer (A / D) 10 through another electrode of the ion trap and amplified to a signal size suitable for detection through the preamplifier 11 shown in FIG. Signal processing is performed on the computer. In each step of control, a control signal is transmitted to the control device of each step through the control unit 12 by a control program operated by a computer.

종래기술에 따른 이온 싸이클로트론 공명 질량 분석기의 신호처리부는 도 2에 도시된 바완 같다.The signal processor of the ion cyclotron resonance mass spectrometer according to the prior art is as shown in FIG.

즉, 도 2는 이온 싸이클로트론 공명 질량 분석기의 신호처리 과정 흐름도로서, 종래의 이온 싸이클로트론 공명 질량 분석기는 디지타이저에서 데이터를 얻는 데이터 취득과정(S1), 얻은 데이터의 신호 대 잡음 비 향상을 위하여 윈도우잉 기법을 사용하는 윈도우잉 처리과정(S2), 시간 영역 신호를 주파수 영역신호로 처리하는 FFT(Fast Fourier Transformer) 처리과정(S3), 주파수 영역 신호를 m/z 영역으로 변환하는 질량영역 표시과정(S4)으로 구성되어, 신호처리 알고리즘을 수행한다. In other words, Figure 2 is a flow diagram of the signal processing process of the ion cyclotron resonance mass spectrometer, the conventional ion cyclotron resonance mass spectrometer is a data acquisition process (S1) for obtaining data from the digitizer, windowing technique for improving the signal-to-noise ratio of the obtained data Windowing processing using S2, Fast Fourier Transformer (F3) processing time domain signal into frequency domain signal (S3), Mass domain display process converting frequency domain signal into m / z domain (S4) ) To perform a signal processing algorithm.

여기서, 신호처리 과정인 윈도우잉 처리과정(S2)의 신호처리를 수행할 때 잡음을 줄일 수는 있지만 신호의 손실이 발생한다. Here, although the noise may be reduced when performing the signal processing of the windowing process S2, which is a signal processing process, signal loss occurs.

따라서, 종래기술에 따른 이온 싸이클로트론 공명 질량 분석기는 윈도우잉 자체가 제한된 영역만을 선택하는 과정이므로 신호 손실이 불가피한 문제점이 있었다.
Therefore, the ion cyclotron resonance mass spectrometer according to the prior art has a problem in that signal loss is inevitable because the windowing process selects only a limited area.

본 발명은 상기한 종래의 문제점을 해결하기 위한 것으로서, 이온 싸이클로트론 공명 질량 분석기의 윈도우잉 과정에서 생기는 신호손실을 미분 윈도우로 처리하여 신호손실을 줄이고 미분 윈도우잉의 단점인 노이즈 증가를 노이즈 저감 알고리즘을 사용하여 노이즈 성분을 감쇄하도록 하는 신호대 잡음 개선을 위한 이온 싸이클로트론 공명 질량 분석기 신호처리방법을 제공하는데 그 목적이 있다.The present invention is to solve the above-mentioned problems, the signal loss generated during the windowing process of the ion cyclotron resonance mass spectrometer by processing the differential window to reduce the signal loss and increase the noise which is a disadvantage of the differential windowing noise reduction algorithm An object of the present invention is to provide an ion cyclotron resonance mass spectrometer signal processing method for improving signal-to-noise that is used to attenuate noise components.

본 발명의 목적을 달성하기 위한 신호대 잡음 개선을 위한 이온 싸이클로트론 공명 질량 분석기 신호처리과정은 디지타이저에서 얻어진 데이터의 신호 대 잡음비 향상 후, 시간 영역 신호를 주파수 영역신호로 FFT(Fast Fourier Transformer) 처리하고, 주파수 영역신호를 질량대 전하비(m/z)영역으로 변환하는 이온 싸이클로트론 공명 질량분석기의 신호처리방법에 있어서, 상기 디지타이저를 통하여 얻어진 아날로그 신호를 디지털신호로 변환하는 데이터취득과정(S10); 상기 데이터취득과정(S10)으로부터 얻어진 신호 성분에 포함된 DC전압을 제거해 주는 DC전압처리과정(S20); 상기 DC전압처리과정(S20)을 통해 얻어진 신호를 미분 윈도우잉 신호처리를 통해 신호성분에 대한 손실을 줄이는 미분윈도우처리과정(S30); 상기 미분윈도우처리과정(S30)에서 얻어진 시간 영역의 데이터를 FFT(Fast Fourier Transformer) 연산을 수행하는 FFT처리과정(S40); 상기 FFT처리과정(S40)을 통해 연산되는 과정에서 증가된 노이즈성분을 제거하는 노이즈저감연산과정(S50); 및 상기 노이즈저감연산과정(S50)으로부터 노이즈 감소된 데이터의 주파수 영역 신호를 질량 대 전하비(m/z) 영역으로 변환하는 질량영역표시과정(S60);를 포함하여 이루어진 것을 특징으로 한다.
In order to achieve the object of the present invention, the ion cyclotron resonance mass spectrometer signal processing process for improving the signal-to-noise is performed after the signal-to-noise ratio of the data obtained from the digitizer is improved, and the time-domain signal is processed to the frequency-domain signal by FFT (Fast Fourier Transformer). A signal processing method of an ion cyclotron resonance mass spectrometer for converting a frequency domain signal into a mass-to-charge ratio (m / z) region, the method comprising: a data acquisition process of converting an analog signal obtained through the digitizer into a digital signal (S10); A DC voltage processing step (S20) of removing the DC voltage included in the signal component obtained from the data acquisition step (S10); A differential window processing step (S30) of reducing the loss of signal components through differential windowing signal processing of the signal obtained through the DC voltage processing step (S20); An FFT process (S40) for performing a fast fourier transformer (FFT) operation on the data in the time domain obtained in the differential window processing process (S30); A noise reduction operation (S50) for removing an increased noise component in the operation calculated through the FFT process (S40); And a mass region display process S60 for converting the frequency domain signal of the noise reduced data from the noise reduction operation S50 to a mass-to-charge ratio (m / z) region.

본 발명에 따른 신호 대 잡음 개선을 위한 이온 싸이클로트론 공명 질량 분석기의 신호처리방법은 신호 성분의 손실을 줄이기 위한 미분 차수 윈도우잉과 늘어난 노이즈 성분에 대한 노이즈 저감 연산과정을 수행하므로, 신호성분의 향상 및 노이즈 저감 효과가 있다.
The signal processing method of the ion cyclotron resonance mass spectrometer for improving the signal-to-noise according to the present invention performs differential order windowing to reduce the loss of signal components and a noise reduction operation for the increased noise components. It has a noise reduction effect.

도 1은 일반적인 이온 싸이클로트론 공명 질량 분석기 제어장치의 전체 구성도이고,
도 2은 종래기술에 따른 이온 싸이클로트론 공명 질량 분석기의 신호처리 과정 흐름도이고,
도 3는 본 발명의 실시예에 따른 신호 대 잡음 개선을 위한 이온 싸이클로트론 공명 질량 분석기의 신호처리과정의 흐름도이고,
도 4은 본 발명의 실시예에 따른 미분 윈도우과정 및 노이즈저감연산과정을 통해 신호 대 잡음비 산출 및 분해능(resolving power)에 대한 최적점을 찾기 위한 비교 도표이고,
도 5는 미분 윈도우잉 연산에 대한 이해를 돕기 위한 도면이다.
1 is an overall configuration diagram of a general ion cyclotron resonance mass spectrometer control device,
Figure 2 is a flow diagram of the signal processing process of the ion cyclotron resonance mass spectrometer according to the prior art,
3 is a flowchart of a signal processing procedure of an ion cyclotron resonance mass spectrometer for signal-to-noise improvement according to an embodiment of the present invention;
4 is a comparison chart for finding an optimal point for signal-to-noise ratio calculation and resolving power through a differential window process and a noise reduction operation process according to an embodiment of the present invention.
5 is a diagram to help understand the differential windowing operation.

본 발명의 실시예에 따른 신호 대 잡음 개선을 위한 이온 싸이클로트론 공명 질량 분석기의 신호처리방법을 첨부된 도 3 내지 도 5를 참조하여 상세히 설명하면 다음과 같다.The signal processing method of the ion cyclotron resonance mass spectrometer for signal-to-noise improvement according to an embodiment of the present invention will be described in detail with reference to FIGS. 3 to 5 as follows.

본 발명을 구현하기 위한 온 싸이클로트론 공명 질량 분석기 제어장치는 도 1에 도시된 바와 같이 배경기술에서 설명하였으며, 본 발명은 동일부호를 부여하고 그 상세한 설명은 생략한다. An on-cyclotron resonance mass spectrometer control device for implementing the present invention has been described in the background art as shown in FIG. 1, and the present invention is denoted by the same reference numerals and the detailed description thereof is omitted.

즉, 본 발명은 상기 디지타이저(10)를 통하여 얻어진 디지털신호를 수신한 컴퓨터(9)에서 신호처리 알고리즘에 의하여 발생된 제어 신호가 제어부(12)를 통하여 각 단계의 제어 장치에 신호가 전달하게 된다.That is, according to the present invention, the control signal generated by the signal processing algorithm in the computer 9 receiving the digital signal obtained through the digitizer 10 is transmitted to the control device of each stage through the control unit 12. .

이와 같이, 상기 컴퓨터(9)에서 신호 대 잡음 개선을 위한 이온 싸이클로트론 공명 질량 분석기의 신호처리 과정에 대하여 다음 도 3을 참조하여 상세히 설명하면 다음과 같다.As described above, the signal processing of the ion cyclotron resonance mass spectrometer for signal-to-noise improvement in the computer 9 will be described in detail with reference to FIG. 3.

도 3은 본 발명의 실시예에 따른 신호 대 잡음 개선을 위한 이온 싸이클로트론 공명 질량 분석기의 신호처리과정의 흐름도로서, 디지타이저(10)를 통하여 얻어진 아날로그 신호를 디지털신호로 변환하는 데이터취득과정(S10)과, 상기 데이터취득과정(S10)으로부터 얻어진 신호 성분에 포함된 DC전압을 제거해 주는 DC전압처리과정(S20)과, 상기 DC전압처리과정(S20)을 통해 얻어진 신호를 미분 윈도우잉 신호처리를 통해 신호성분에 대한 손실을 줄이는 미분윈도우처리과정(S30)과, 상기 미분윈도우처리과정(S30)에서 얻어진 시간 영역의 데이터를 FFT(Fast Fourier Transformer) 연산을 수행하는 FFT처리과정(S40)과, 상기 FFT처리과정(S40)을 통해 연산되는 과정에서 증가된 노이즈성분을 제거하는 노이즈저감연산과정(S50)과, 상기 노이즈저감연산과정(S50)으로부터 노이즈 감소된 데이터의 주파수 영역 신호를 질량 대 전하비(m/z) 영역으로 변환하는 질량영역표시과정(S60)으로 이루어진다.3 is a flowchart illustrating a signal processing process of an ion cyclotron resonance mass spectrometer for improving signal-to-noise according to an embodiment of the present invention, and a data acquisition process of converting an analog signal obtained through the digitizer 10 into a digital signal (S10). And a DC voltage process (S20) for removing a DC voltage included in the signal component obtained from the data acquisition process (S10), and a signal obtained through the DC voltage process (S20) through differential windowing signal processing. Differential window processing to reduce the loss of signal components (S30), FFT processing (S40) to perform the FFT (Fast Fourier Transformer) operation on the data in the time domain obtained in the differential window processing (S30), and Noise reduction operation (S50) for removing the increased noise component in the process calculated through the FFT process (S40), and noise reduction from the noise reduction operation (S50) A frequency domain signal of the data comprises a mass area display process (S60) to be converted to a mass to charge ratio (m / z) domain.

이와 같이 이루어진 본 발명의 실시 예에 따른 신호 대 잡음 개선을 위한 이온 싸이클로트론 공명 질량 분석기의 신호처리방법에 대하여 상세히 설명하면 다음과 같다.When the signal processing method of the ion cyclotron resonance mass spectrometer for improving the signal-to-noise according to the embodiment of the present invention made as described above will be described in detail.

상기 미분윈도우 처리과정(30)에서 중요한 것은 몇 번의 미분을 행하는 가에 있다. What is important in the differential window processing process 30 is how many different derivatives are performed.

시간 영역에서의 신호 f(t)를 푸리에 변환식으로 표현하면 다음 수학식 1과 같으며,When the signal f (t) in the time domain is expressed by a Fourier transform equation,

Figure pat00001
Figure pat00001

다음 수학식 2는 n번의 미분에 대한 식이다.Equation 2 is an expression for n derivatives.

Figure pat00002
Figure pat00002

결국, n번 미분을 행한다는 것은 시간영역에서 t를 n번 곱하는 것과 같다. After all, performing n derivatives is equivalent to multiplying t times in the time domain.

여기서, 문제는 몇 번 곱하는 것이 신호 성분과 노이즈 그 밖에 신호특성(질량분석기에서는 Resolving power 변수)에 최적화된 값을 찾아야 한다.Here, the problem is to multiply a few times to find a value that is optimized for signal components, noise and other signal characteristics (resolving power variables in mass spectrometers).

상기 노이즈 저감 연산과정(S50)은 다음 수학식 3을 통해 노이즈값을 산출하게 된다.The noise reduction calculation process (S50) is to calculate the noise value through the following equation (3).

Figure pat00003
Figure pat00003

ni : 현재의 노이즈 성분, ni-1 : 직전의 노이즈 성분,n i : current noise component, n i-1 : previous noise component,

nmean : 모든 데이터의 평균값, Mi : 현재의 처리하는 데이터 값,n mean : average of all data, M i : current data value

n0 : 첫 노이즈 성분 계산n 0 : First noise component calculation

먼저, 전체 데이터에 대한 평균값(nmean)을 산출한 다음, 각 데이터에 대한 노이즈 레벨을 계산한다. First, an average value (n mean ) for all data is calculated, and then a noise level for each data is calculated.

처음 데이터에 대한 처음 노이즈 성분(n0)을 구한다. 위의 식에 따르거나 아니면 적절한 값을 준다. α값에 대한 결정은 n0를 계산할 때, M1값이 nmean보다 크면 0보다 크고 1과 같거나 작은 수로 결정한다. 바람직하게는 1로 결정하며, 반대인 경우는 0보다 작고 -1과 같거나 큰 수, 바람직하게는 -1로 결정한다. 이와 같은 방법으로 전체 각 데이터에 대한 ni를 모두 구한다.Find the first noise component (n 0 ) for the first data. Follow the above formula or give an appropriate value. When determination of the value of α calculate the n 0, and determines the number equal to or smaller than zero is greater 1 M 1 n value is greater than the mean. Preferably it is determined to be 1, and vice versa to a number less than 0 and equal to or greater than -1, preferably -1. In this way, we find all n i for each data.

즉, 신호상태에 따라 다르나, 증가하면 +1, 감소하면 -1로 연산하나, 정밀도가 필요한 경우 알고리즘에 의하여 증가하면 0보다 크고 1과 같거나 작은 수, 감소하면 0보다 작거나 -1보다 크거나 같은 수로 결정한다. That is, it depends on the signal condition, but if it is increased, it is +1 and if it is decreased, it is -1.However, if the precision is needed, it is bigger than 0 and equal to or less than 1 if it is increased by the algorithm. Or the same number.

상기 수학식 3에서 산출된 노이즈 값들은 다음 수학식 4를 통해 감산을 행하면 산출된 새로운 데이터값 Si가 구해진다. When the noise values calculated in Equation 3 are subtracted through the following Equation 4, the calculated new data value S i is obtained.

Figure pat00004
Figure pat00004

여기서, Mi 는 현재의 처리하는 데이터 값이고, Si는 계산된 현재의 데이터 값이고, ni 는 현재의 노이즈 성분 계산이다.Where M i is the current data value to be processed, S i is the calculated current data value, and n i is the current noise component calculation.

단, 여기서 Mi 값이 nmean 보다 3~5 배 큰 경우, 바람직하게는 5로 결정한다. Si는 Mi를 적용한다.However, when the M i value is 3 to 5 times larger than n mean , it is preferably determined as 5. S i applies to M i .

즉, 상기 수학식 4에서와 같은 방식의 노이즈 감산은 신호에 포함된 노이즈 값이 증가되면 노이즈 레벨도 올라가고, 노이즈값이 감소하면 노이즈 레벨도 감소하여 노이즈 성분의 증감을 따라가며 가감 연산을 하는 특징이 있다.
That is, the noise subtraction in the same manner as in Equation 4 increases the noise level when the noise value included in the signal increases, and decreases the noise level when the noise value decreases to follow the increase and decrease of the noise component. There is this.

도 4는 본 발명의 실시예에 따른 미분윈도우 처리과정 및 노이즈저감과정을 통해 신호 대 잡음비 산출 및 분해능(resolving power)에 대한 비교 도표로서, 상기 DC전압처리과정(S20)을 통해 신호 성분에 포함된 DC전압을 제거한 원 데이터와 상기 미분윈도우처리과정(S30)을 통해 1차~7차 미분 값들에 대한 신호 대 잡음비(S/N)를 산출하고, 상기 노이즈저감연산과정(S50)을 통한 신호 대 잡음비(Noise reduction S/N)를 산출하고, 질량신호를 특징짓는 분해능(Resolving power)에 대한 비교를 통하여 최적의 미분 차수를 결정지을 수 있다. 4 is a comparison chart of signal-to-noise ratio calculation and resolving power through differential window processing and noise reduction according to an embodiment of the present invention, which is included in a signal component through the DC voltage processing (S20). The signal-to-noise ratio (S / N) for the first to seventh order derivatives is calculated through the raw data from which the DC voltage is removed and the differential window processing process (S30), and the signal through the noise reduction operation process (S50). Noise reduction S / N can be calculated and the optimum differential order can be determined by comparing the resolving power that characterizes the mass signal.

여기서, 도 4에서와 같이 1차 미분을 행하였을 때가 최적의 신호 향상을 가져옴을 알 수 있다. 즉, 최적 미분 연산점은 1차 미분 연산점이 된다.Here, as shown in FIG. 4, it can be seen that the optimal signal improvement is obtained when the first derivative is performed. In other words, the optimal differential calculation point becomes the first derivative calculation point.

도 5a, 5b는 미분 윈도우잉에 따른 신호손실 저감을 위한 설명도로서, 미분 윈도우잉을 행할 때 도 5b에서처럼 곱하는 함수의 시작점이 x=0일 때, y=1에서 시작하면 원신호를 살리는 결과가 나온다. 하지만 도 5a에서 처럼 x=0일 때, y=0인 원점에서 시작하여도 충분히 신호 손실은 있지만 결과값을 얻는데는 동일한 결과가 나옴을 알 수 있다.
5A and 5B are explanatory diagrams for reducing signal loss due to differential windowing. When differential windowing is performed, starting at y = 1 when the starting point of the multiplication function as in FIG. 5B is x = 0, the original signal is saved. Comes out. However, when x = 0, as shown in FIG. 5a, even when starting from the origin of y = 0, there is a sufficient signal loss, but the same result is obtained.

본 발명은 상기한 실시 예에 한정되지 않고, 이하의 특허청구 범위에서 청구하는 본 발명의 요지를 벗어남이 없이 본 발명이 속하는 분야에서 통상의 지식을 가진 자라면 누구든지 다양한 변경 실시가 가능한 범위까지 본 발명의 기술적 정신이 있다고 할 것이다.
The present invention is not limited to the above-described embodiments, and any person having ordinary skill in the art to which the present invention pertains can make various changes without departing from the gist of the present invention as claimed in the following claims. It will be said that there is a technical spirit of the present invention.

9 : 컴퓨터 10 : 디지타이저
12 : 제어부
9: computer 10: digitizer
12: control unit

Claims (4)

디지타이저(10)에서 얻어진 데이터를 수신한 컴퓨터(9)에 의해 신호 대 잡음비 향상 후, 시간 영역 신호를 주파수 영역신호로 FFT(Fast Fourier Transformer) 처리하고, 주파수 영역신호를 질량대 전하비(m/z)영역으로 변환하는 이온 싸이클로트론 공명 질량분석기의 신호처리방법에 있어서,
상기 디지타이저(10)를 통하여 얻어진 아날로그 신호를 디지털신호로 변환하는 데이터취득과정(S10);
상기 컴퓨터(10)에 의하여 데이터취득과정(S10)으로부터 얻어진 신호 성분에 포함된 DC전압을 제거해 주는 DC전압처리과정(S20);
상기 DC전압처리과정(S20)을 통해 얻어진 호를 미분 윈도우잉 신호처리를 통해 신호성분에 대한 손실을 줄이는 미분윈도우처리과정(S30);
상기 미분윈도우처리과정(S30)에서 얻어진 시간 영역의 데이터를 FFT(Fast Fourier Transformer) 연산을 수행하는 FFT처리과정(S40);
상기 FFT처리과정(S40)을 통해 연산되는 과정에서 증가된 노이즈성분을 제거하는 노이즈저감연산과정(S50); 및
상기 노이즈저감연산과정(S50)으로부터 노이즈 감소된 데이터의 주파수 영역 신호를 질량 대 전하비(m/z) 영역으로 변환하는 질량영역표시과정(S60);를 포함하여 이루어진 것을 특징으로 하는 신호대 잡음 개선을 위한 이온 싸이클로트론 공명 질량 분석기 신호처리방법.
After the signal-to-noise ratio is improved by the computer 9 receiving the data obtained from the digitizer 10, the time-domain signal is processed into a frequency domain signal by FFT (Fast Fourier Transformer), and the frequency-domain signal is mass-to-charge ratio (m / In the signal processing method of the ion cyclotron resonance mass spectrometer to convert to z) region,
A data acquisition process of converting an analog signal obtained through the digitizer 10 into a digital signal (S10);
DC voltage processing (S20) for removing the DC voltage included in the signal component obtained from the data acquisition process (S10) by the computer (10);
Differential windowing process to reduce the loss of the signal component of the signal obtained from the DC voltage processes (S20) through differential signal windowing process (S30);
An FFT process (S40) for performing a fast fourier transformer (FFT) operation on the data in the time domain obtained in the differential window processing process (S30);
A noise reduction operation (S50) for removing an increased noise component in the operation calculated through the FFT process (S40); And
And a mass domain display process (S60) for converting the frequency domain signal of the noise reduced data from the noise reduction operation process (S50) into a mass to charge ratio (m / z) region. Ion cyclotron resonance mass spectrometer signal processing method for
제 1 항에 있어서,
상기 미분윈도우처리과정(S30) 및 노이즈저감연산과정(S50)은 신호성분에 따른 노이즈 값의 증감을 통해 새로운 신호성분을 연산하여 신호 대 잡음비를 향상시킨 데이터와 n차의 미분연산에 대한 최적 미분 연산점을 산출하여 노이즈성분을 제거하는 것을 특징으로 하는 신호대 잡음 개선을 위한 이온 싸이클로트론 공명 질량 분석기 신호처리방법.
The method of claim 1,
The differential window processing process (S30) and the noise reduction operation process (S50) calculate the new signal component by increasing or decreasing the noise value according to the signal component to improve the signal-to-noise ratio and the optimal derivative for the nth-order differential operation. Ion cyclotron resonance mass spectrometer signal processing method for improving the signal-to-noise, characterized in that to remove the noise components by calculating the operating point.
제 1 항에 있어서,
상기 미분윈도우처리과정(S30)은 다음 수학식 1과 같이 n번의 미분을 수행하는 것을 특징으로 하는 신호대 잡음 개선을 위한 이온 싸이클로트론 공명 질량 분석기 신호처리방법.
[수학식 1]
Figure pat00005

여기서,
Figure pat00006
는 시간영역에서의 신호 f(t)를 푸리에 변환식으로 표현한 것임.
The method of claim 1,
The differential window processing process (S30) is an ion cyclotron resonance mass spectrometer signal processing method for improving the signal-to-noise, characterized in that performing the differential n times as shown in the following equation (1).
[Equation 1]
Figure pat00005

here,
Figure pat00006
Is the Fourier transform of the signal f (t) in the time domain.
제 1 항에 있어서,
상기 노이즈저감연산과정(S50)은 전체 원래 신호 데이터에 대한 평균값(nmean)을 산출하는 단계;
상기 전체 원래 신호 데이터에 대한 각 노이즈값(ni)를 다음 수학식 2에 의해 산출하는 단계;
[수학식 2]
Figure pat00007

여기서, ni은 현재의 노이즈값, ni-1는 직전의 노이즈 값, Mi는 현재의 처리하는 데이터 값, n0는 첫 노이즈 값이고, 상기 α값에 대한 결정은 n0를 계산할 때는 M1값이 nmean보다 크면 0 보다 크고 1과 같거나 작은 수 이며, 상기 산출된 노이즈 값들은 각 원래 신호 데이터 값에 대하여 다음 수학식 3에 의하여 노이즈 성분의 증감에 따라 가감 연산을 수행하는 단계;로 이루어진 것을 특징으로 하는 신호대 잡음 개선을 위한 이온 싸이클로트론 공명 질량 분석기 신호처리방법.
[수학식 3]
Figure pat00008

여기서, Mi 는 현재의 처리하는 데이터 값이고, Si는 계산된 현재의 데이터 값이고, ni 는 현재의 노이즈 값이고, 여기서 Mi 값이 nmean 보다 3~5 배 큰 경우, Si는 Mi를 적용함.
The method of claim 1,
The noise reduction operation S50 may include calculating a mean value n of all original signal data;
Calculating each noise value (n i ) for the entire original signal data by the following equation (2);
&Quot; (2) "
Figure pat00007

Where, n i is the current noise value, n i-1 is the noise value of the immediately preceding, M i is the data value for the current processing, n 0 is the initial noise value, determined for the α value when calculating a n 0 If the M 1 value is greater than n mean, the number is greater than 0 and equal to or less than 1, and the calculated noise values are subjected to the arithmetic operation according to the increase and decrease of the noise component according to the following equation 3 for each original signal data value. Ion cyclotron resonance mass spectrometer signal processing method for improving the signal-to-noise, characterized in that consisting of.
[Equation 3]
Figure pat00008

Where M i is the current data value, S i is the calculated current data value, n i is the current noise value, and where M i is 3 to 5 times larger than n mean , S i Applies M i .
KR1020100129680A 2010-12-17 2010-12-17 A Signal Process Method for improving S/N of Ion Cyclotron Resonance Mass Spectrometer KR101162178B1 (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
KR1020100129680A KR101162178B1 (en) 2010-12-17 2010-12-17 A Signal Process Method for improving S/N of Ion Cyclotron Resonance Mass Spectrometer
PCT/KR2011/008230 WO2012081825A1 (en) 2010-12-17 2011-10-31 Method for processing the signal of an ion cyclotron resonance mass spectrometer for signal-to-noise improvement

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020100129680A KR101162178B1 (en) 2010-12-17 2010-12-17 A Signal Process Method for improving S/N of Ion Cyclotron Resonance Mass Spectrometer

Publications (2)

Publication Number Publication Date
KR20120068181A true KR20120068181A (en) 2012-06-27
KR101162178B1 KR101162178B1 (en) 2012-07-04

Family

ID=46244890

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020100129680A KR101162178B1 (en) 2010-12-17 2010-12-17 A Signal Process Method for improving S/N of Ion Cyclotron Resonance Mass Spectrometer

Country Status (2)

Country Link
KR (1) KR101162178B1 (en)
WO (1) WO2012081825A1 (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111949919B (en) * 2020-06-14 2022-04-01 杭州谱育科技发展有限公司 Method for reducing noise in AES detection

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3741563B2 (en) * 1999-04-15 2006-02-01 日本電子株式会社 Data collection system for mass spectrometer
KR100947869B1 (en) * 2007-12-31 2010-03-18 한국기초과학지원연구원 Method for controlling signal generation of Fourier Transform Ion Cyclotron Resonance Mass Spectrometer
KR100962550B1 (en) * 2007-12-31 2010-06-14 한국기초과학지원연구원 Apparatus for acquiring ion signal of Fourier Transform Ion Cyclotron Resonance Mass Spectrometer and method for the same

Also Published As

Publication number Publication date
KR101162178B1 (en) 2012-07-04
WO2012081825A1 (en) 2012-06-21

Similar Documents

Publication Publication Date Title
US8289029B2 (en) Application of wideband sampling for arc detection with a probabilistic model for quantitatively measuring arc events
JP6090479B2 (en) Mass spectrometer
US8637809B2 (en) Tandem quadrupole mass spectrometer
JPH10508386A (en) Methods for reducing interference in elemental mass spectrometers
JP2006236795A (en) Mass spectroscope and mass spectrometry
JPWO2008129929A1 (en) Gas analyzer
US8477956B2 (en) Howling suppression device, howling suppression method, program, and integrated circuit
JP2008020459A (en) System and method of removing noise from spectral data
US8264237B2 (en) Application of wideband sampling for arc detection with a probabilistic model for quantitatively measuring arc events
CN112037816A (en) Voice signal frequency domain frequency correction, howling detection and suppression method and device
JP2018004298A (en) Mass spectrometer
KR101162178B1 (en) A Signal Process Method for improving S/N of Ion Cyclotron Resonance Mass Spectrometer
US9196467B2 (en) Mass spectrum noise cancellation by alternating inverted synchronous RF
EP2674961A3 (en) Electron microscope and method of operating the same
US20190325003A1 (en) Noise reduction apparatus and noise suppressing method
JP6232012B2 (en) Detection device, reception device, and detection method
JP6236755B2 (en) Passive sonar device, transient signal processing method and signal processing program thereof
US11769517B2 (en) Signal processing apparatus, signal processing method, and signal processing program
JP7047698B2 (en) Analysis equipment
US20210174820A1 (en) Signal processing apparatus, voice speech communication terminal, signal processing method, and signal processing program
JPH1164431A (en) Partial discharge measurement method
US20140044266A1 (en) Real-time noise reduction apparatus for radio monitoring system
US20120267524A1 (en) Apparatus and method for controlling a pipeline-type ion cyclotron resonance mass spectrometer
JP5990422B2 (en) Peak suppression circuit
US10416130B2 (en) Mass spectrometry data processing apparatus, mass spectrometry system, and method for processing mass spectrometry data

Legal Events

Date Code Title Description
A201 Request for examination
E902 Notification of reason for refusal
E701 Decision to grant or registration of patent right
E701 Decision to grant or registration of patent right
GRNT Written decision to grant
FPAY Annual fee payment

Payment date: 20150731

Year of fee payment: 4

FPAY Annual fee payment

Payment date: 20160613

Year of fee payment: 5

FPAY Annual fee payment

Payment date: 20170609

Year of fee payment: 6

FPAY Annual fee payment

Payment date: 20180627

Year of fee payment: 7