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 PDFInfo
- 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
Links
Images
Classifications
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01J—ELECTRIC DISCHARGE TUBES OR DISCHARGE LAMPS
- H01J49/00—Particle spectrometers or separator tubes
- H01J49/0027—Methods for using particle spectrometers
- H01J49/0036—Step by step routines describing the handling of the data generated during a measurement
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01J—ELECTRIC DISCHARGE TUBES OR DISCHARGE LAMPS
- H01J49/00—Particle spectrometers or separator tubes
- H01J49/26—Mass spectrometers or separator tubes
- H01J49/34—Dynamic spectrometers
- H01J49/36—Radio frequency spectrometers, e.g. Bennett-type spectrometers, Redhead-type spectrometers
- H01J49/38—Omegatrons ; using ion cyclotron resonance
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N27/00—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
- G01N27/62—Investigating 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/622—Ion mobility spectrometry
- G01N27/623—Ion mobility spectrometry combined with mass spectrometry
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01J—ELECTRIC DISCHARGE TUBES OR DISCHARGE LAMPS
- H01J49/00—Particle spectrometers or separator tubes
- H01J49/02—Details
- H01J49/022—Circuit 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
Description
본 발명은 이온 싸이클로트론 공명 질량 분석기 신호처리방법에 관한 것이며, 더욱 상세하게는 노이즈 저감 알고리즘을 사용하여 노이즈 성분을 감쇄하도록 하는 신호 대 잡음 성분 개선을 위한 이온 싸이클로트론 공명 질량분석기의 신호처리방법에 관한 것이다.
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
여기된 신호는 이온트랩의 또 다른 전극을 통하여 도 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
종래기술에 따른 이온 싸이클로트론 공명 질량 분석기의 신호처리부는 도 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
이와 같이, 상기 컴퓨터(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
이와 같이 이루어진 본 발명의 실시 예에 따른 신호 대 잡음 개선을 위한 이온 싸이클로트론 공명 질량 분석기의 신호처리방법에 대하여 상세히 설명하면 다음과 같다.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,
다음 수학식 2는 n번의 미분에 대한 식이다.
결국, 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).
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
여기서, 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
도 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)를 통하여 얻어진 아날로그 신호를 디지털신호로 변환하는 데이터취득과정(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
상기 미분윈도우처리과정(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.
상기 미분윈도우처리과정(S30)은 다음 수학식 1과 같이 n번의 미분을 수행하는 것을 특징으로 하는 신호대 잡음 개선을 위한 이온 싸이클로트론 공명 질량 분석기 신호처리방법.
[수학식 1]
여기서, 는 시간영역에서의 신호 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]
here, Is the Fourier transform of the signal f (t) in the time domain.
상기 노이즈저감연산과정(S50)은 전체 원래 신호 데이터에 대한 평균값(nmean)을 산출하는 단계;
상기 전체 원래 신호 데이터에 대한 각 노이즈값(ni)를 다음 수학식 2에 의해 산출하는 단계;
[수학식 2]
여기서, ni은 현재의 노이즈값, ni-1는 직전의 노이즈 값, Mi는 현재의 처리하는 데이터 값, n0는 첫 노이즈 값이고, 상기 α값에 대한 결정은 n0를 계산할 때는 M1값이 nmean보다 크면 0 보다 크고 1과 같거나 작은 수 이며, 상기 산출된 노이즈 값들은 각 원래 신호 데이터 값에 대하여 다음 수학식 3에 의하여 노이즈 성분의 증감에 따라 가감 연산을 수행하는 단계;로 이루어진 것을 특징으로 하는 신호대 잡음 개선을 위한 이온 싸이클로트론 공명 질량 분석기 신호처리방법.
[수학식 3]
여기서, 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) "
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]
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 .
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)
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)
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 |
-
2010
- 2010-12-17 KR KR1020100129680A patent/KR101162178B1/en active IP Right Grant
-
2011
- 2011-10-31 WO PCT/KR2011/008230 patent/WO2012081825A1/en active Application Filing
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 |