KR100852195B1 - 적응 필터 적용 여부 판단 시스템 및 그 방법 - Google Patents

적응 필터 적용 여부 판단 시스템 및 그 방법 Download PDF

Info

Publication number
KR100852195B1
KR100852195B1 KR1020060132214A KR20060132214A KR100852195B1 KR 100852195 B1 KR100852195 B1 KR 100852195B1 KR 1020060132214 A KR1020060132214 A KR 1020060132214A KR 20060132214 A KR20060132214 A KR 20060132214A KR 100852195 B1 KR100852195 B1 KR 100852195B1
Authority
KR
South Korea
Prior art keywords
adaptive filter
signal
adaptive
apply
determining whether
Prior art date
Application number
KR1020060132214A
Other languages
English (en)
Other versions
KR20080058098A (ko
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 KR1020060132214A priority Critical patent/KR100852195B1/ko
Priority to US11/819,477 priority patent/US8005882B2/en
Publication of KR20080058098A publication Critical patent/KR20080058098A/ko
Application granted granted Critical
Publication of KR100852195B1 publication Critical patent/KR100852195B1/ko

Links

Images

Classifications

    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03HIMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
    • H03H17/00Networks using digital techniques
    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03HIMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
    • H03H21/00Adaptive networks
    • H03H21/0012Digital adaptive filters

Landscapes

  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
  • Noise Elimination (AREA)
  • Filters That Use Time-Delay Elements (AREA)

Abstract

본 발명은 검출 신호를 센싱하는 신호 센싱부와, 상기 센싱된 신호를 적응적으로 필터링하는 적응 필터부 및 상기 센싱된 신호 및 상기 필터링된 신호를 분석하고, 상기 분석된 결과에 기초하여 상기 센싱된 신호에 대한 상기 적응 필터부의 적용 여부를 판단하는 적응 필터 적용 여부 판단부를 포함하는 적응 필터 적용 판단 시스템을 제공한다.
적응 필터, 적용, 판단, PPG, NW, ACC

Description

적응 필터 적용 여부 판단 시스템 및 그 방법{system for determined application of and Method thereof}
도 1은 본 발명의 일실시예에 따른 적응 필터 적용 여부 판단 시스템의 구성도를 나타내는 도면이다.
도 2는 적응 필터부에서 PPG 신호에 포함되어 있는 움직임 노이즈가 제거된 결과를 나타내는 도면이다.
도 3은 PPG 신호와 움직임 노이즈 신호의 패턴이 동일한 구간에서의 적응 필터부의 출력 결과를 나타내는 도면이다.
도 4는 본 발명에 따른 적응 필터 적용 여부 판단부의 동작에 대한 적응 필터 적용 여부 판별 알고리즘을 나타내는 도면이다.
도 5는 측정하고자 하는 신호와 노이즈 신호의 주기가 일치하는 경우 R(n)값의 변화를 나타내는 도면이다.
도 6은 R(n)에 대한 변화 그래프 및 NW에 대한 변화 그래프의 일례를 나타내는 도면이다.
도 7은 ACC 비율에 대한 변화 그래프 및 NW에 대한 변화 그래프의 일례를 나타내는 도면이다.
도 8은 도 4에 도시된 적응 필터 적용 여부 판단 방법의 알고리즘을 수행한 결과의 일례를 나타내는 도면이다.
<도면의 주요 부분에 대한 부호의 설명>
100: 적응 필터 적용 여부 판단 시스템
111~11n: 제1 내지 제N 신호 센싱부
121, 122: 제1, 제2 전처리부
130: 적응 필터 적용 여부 판단부
140: 적응 필터부
150: 신호 처리부
160: 신호 출력부
본 발명은 적응 필터 적용 여부를 판단하는 시스템 및 그 방법에 관한 것으로, 더욱 상세하게는 측정하고자 하는 신호에 영향을 주는 복수 개의 노이즈 신호에 대한 영향을 파악함으로써 노이즈를 제거하는 적용 필터에 대한 적용 여부를 판단하는 시스템 및 그 방법에 관한 것이다.
일반적으로 생체 신호는 저주파 신호로서 전극의 움직임에 의해 노이즈가 다량 발생되므로 이를 제거하기 위해 다양한 전극 및 신호 처리 방법이 적용되고 있다. 특히 가속도 신호를 이용하여 전극의 움직임 노이즈를 제거하는 방법이 최근 많은 연구 기관에서 연구되고 있으며, Principal Company Analysis(PCA) 및 적응 필터(Adaptive filter)가 대표적으로 사용되고 있다. 그러나, PCA는 측정 신호가 2채널 이상이어야 한다는 제약 조건이 있기 때문에 입력 센서가 2개 이상 필요하고 데이터 처리 시에도 많은 양의 데이터가 필요하므로 실시간으로 처리하기에 어려움이 많은 문제점이 있다. 한편 적응 필터는 얻고자 하는 생체 신호와 제거하고자 하는 노이즈 패턴이 일정한 경우 출력값이 제로에 가깝게 나오는 문제점이 있다.
본 발명은 상기와 같은 종래기술의 문제점을 해결하기 위해 안출된 것으로서, 측정하고자 하는 신호에 영향을 주는 복수 개의 노이즈 신호의 영향을 파악하여 노이즈를 제거하는 적응 필터를 적용할지 여부를 판단하는 시스템 및 그 방법을 제공하는 것을 목적으로 한다.
본 발명의 다른 목적은 실시간 처리가 가능하고 입력 채널의 장점을 가지는 적응 필터를 사용할 때 생체 신호와 제거하고자 하는 노이즈 패턴이 동일한 경우 발생되는 문제점을 원천적으로 해결할 수 있는 적응 필터 적용 여부 판단 시스템 및 그 방법을 제공하는 것이다.
상기의 목적을 이루고 종래기술의 문제점을 해결하기 위하여, 본 발명은 검출 신호를 센싱하는 신호 센싱부와, 상기 센싱된 신호를 적응적으로 필터링하는 적응 필터부 및 상기 센싱된 신호 및 상기 필터링된 신호를 분석하고, 상기 분석된 결과에 기초하여 상기 센싱된 신호에 대한 상기 적응 필터부의 적용 여부를 판단하는 적응 필터 적용 여부 판단부를 포함하는 적응 필터 적용 판단 시스템을 제공한 다.
본 발명의 다른 일측에 따르는 적응 필터 적용 여부 판단 방법은 각종 센서를 통해 측정하고자 하는 신호를 검출하는 단계와, 적응 필터를 통해 상기 검출된 신호에 대해 적응 필터링 동작을 수행하는 단계와, 상기 검출된 신호 및 상기 적응 필터링된 신호를 분석하고, 상기 분석된 결과에 기초하여 상기 검출된 신호에 대한 적용 필터링 적용 여부를 판단하는 단계를 포함한다.
이하에서는 첨부된 도면들을 참조하여 본 발명의 실시예에 따른 적응 필터 적용 여부 판단 시스템 및 그 방법을 상세히 설명한다. 아울러 본 발명의 상세한 설명에서는 운동중의 맥박 측정 시스템에서 PPG 신호를 측정하는 경우를 일례로 들어 설명하기로 한다.
도 1은 본 발명의 일실시예에 따른 적응 필터에 대한 적용 여부를 판단하는 시스템의 구성도를 나타내는 도면이다.
도 1을 참조하면, 적응 필터 적용 여부 판단 시스템(100)은 제1 내지 N신호 센싱부(111~11n), 제1 및 제2 전처리부(121, 122), 적응 필터 적용 여부 판단부(130), 적응 필터부(140), 신호 처리부(150) 및 신호 출력부(160)를 포함한다.
제1 신호 센싱부(111)는 측정하고자 하는 제1 대상에서 발생되는 제1 신호를 센싱한다. 제1 신호 센싱부(110)는 일례로 운동중인 사용자의 생체 중 세동맥의 변화를 광을 이용하여 PPG(Photoplethysmography) 신호를 센싱(sensing)한다. 즉, 제1 신호 센싱부(111)는 상기 사용자가 안정 상태에서 상기 사용자의 손가락에 PPG 신호 측정 센서를 부착시켜 상기 PPG 신호를 센싱할 수 있다. 상기 사용자가 운동하는 동안에 상기 PPG 신호 센서가 고정되어 있지 않을 경우 상기 PPG 센서의 움직임으로 광의 경로가 변경되어 제대로 된 신호를 검출할 수 없다. 특히 상기 사용자의 귓볼과 같이 PPG 신호의 크기가 미약한 경우 움직임에 의한 노이즈가 상대적으로 원 신호에 크게 영향을 미치게 되는 문제점이 있다. 따라서, 이러한 노이즈를 제거하기 위한 다양한 시도가 이루어지고 있으며, 그중에서 대표적인 것으로 적응 필터부(140)에 의한 노이즈 제거이다.
제2 내지 제N 신호 센싱부(112~11n)도 제1 신호 센싱부(111)와 마찬가지로 각기 측정하고자 하는 대상에서 발생되는 제2 신호 내지 제N 신호를 센싱한다.
제1 전처리부(121)는 제1 신호 센싱부(111)에 의해 센싱된 제1 신호를 증폭시키는 전처리 동작 수행한다. 제2 전처리부(122)는 제2 내지 제N 신호 센싱부(112~11n)에 의해 센싱된 제2 신호 내지 제N 신호를 증폭시키는 전처리 동작을 수행한다.
도 2는 적응 필터를 이용한 PPG 신호에 포함되어 있는 움직임 노이즈 제거 결과의 일례를 나타내는 도면이다.
도 2를 참조하면, 참조부호(210)는 적응 필터부(140)에 의한 적응 필터링이 적용되기 전의 PPG 신호로서 노이즈의 영향으로 PPG 신호의 간격(A1, A2)이 서로 상이한 상태이다. 참조부호(220)는 적응 필터부(140)에 의한 적응 필터링이 적용된 PPG 신호로서 노이즈가 제거됨에 따라 참조부호(210)과 비교해서 보면 PPG 신호의 간격(T1, T2)이 동일한 상태이다.
참조부호(230)는 가속도 신호로서 가속도 신호 주기(A1', A2')가 적응 필터 부(140)에 의한 적응 필터링이 적용되기 전의 PPG 신호에 영향을 받은 상태를 나타낸다. 즉, 적응 필터부(140)에 의한 적응 필터링 동작이 적용되기 전의 PPG 신호의 간격(A1, A2)은 가속도 신호에 직접적으로 영향을 받아 원래 측정하고자 하는 PPG 신호에 가속도 신호 주기(A1', A2')가 그대로 실리게 된다.
참조부호(240)는 심전도 신호로서 적응 필터부(140)에 의해 적응 필터링이 수행됨에 따라 PPG 신호에 적응 필터링 동작이 수행되기 전에 실려 있던 신호의 간격(A1, A2)의 주기가 제거되고, 적응 필터링이 적용된 신호의 주기(T1, T2)와 같은 주기(T1', T2')를 나타낸다. 즉, 적응 필터부(140)에 의한 적응 필터링 동작이 적용되기 전에 실려 있던 신호의 간격(A1, A2)의 주기가 제거되고, 심전도의 R-Peak의 주기(T1', T2')가 적응 필터링 동작이 적용된 PPG 신호의 간격(T1, T2)와 동일하게 나타난다.
이와 같이, 적응 필터부(140)는 제1 신호 센싱부(111)에 의해 센싱된 PPG 신호에 대해 적응 필터링 동작을 수행함에 따라 상기 PPG 신호에 포함된 움직임 노이즈를 제거할 수 있다.
하지만, 적응 필터부(140)는 측정하고자 하는 입력 신호와 제거하고자 하는 노이즈 신호의 패턴이 일정할 경우 원하는 출력값을 얻을 수 없는 문제점이 있다.
도 3은 PPG 신호와 움직임 노이즈 신호의 패턴이 동일한 구간에서의 적응 필터 출력 결과의 일례를 나타내는 도면이다.
도 3을 참조하면, 참조부호(310)은 적응 필터부(140)에 의한 적응 필터링 동작이 수행되기 전의 PPG 신호와 가속도 신호를 각기 나타내며, 참조부호(320)는 적응 필터부(140)에 의한 적응 필터링 동작이 수행된 후의 PPG 신호와 가속도 신호를 각기 나타낸다.
이와 같이, 측정하고자 하는 PPG 신호와 움직임 노이즈 신호가 포함된 가속도 신호가 동일한 구간인 경우는 적응 필터부(140)에 의한 적응 필터링 동작이 수행된 결과로써 상기 적응 필터링 동작이 적용되었음에도 불구하고, 상기 PPG 신호를 복원하지 못하고 오히려 적응 필터링 동작이 적용되기 전보다 상기 PPG 신호 신호 패턴이 깨져서 S/N(Signal/Noise)비가 낮아진다.
따라서, 본 발명에 따른 적응필터 적용 여부 판단부(130)의 동작이 필요한 것이다.
도 4는 본 발명에 따른 적응 필터 적용 여부 판단부의 동작에 대한 적응 필터 적용 여부 판별 알고리즘을 나타내는 도면이다.
도 4를 참조하면, 적응 필터부(140)는 멀티 채널 적응 필터로 구성되며, 원 데이터(D(n))과 각 적응 필터의 기준(reference) 입력값(X1(n)~Xn(n))을 입력받는다. 그리고, 적응 필터부(140)는 각 적응 필터의 기준 입력값(X1(n)~Xn(n))에 대한 계수(coefficient)(W1~Wn)(401~40n)가 반영된 값을 합산하고(411), 상기 원 데이터(D(n))에서 상기 합산된 값을 차감(412)하여 출력값(Y(n))으로 출력한다.
적응 필터 적용 여부 판단부(130)는 적응 필터부(140)에서 출력된 상기 출력값(Y(n)), 적응 필터부(140)에 입력된 상기 원 데이터(D(n)) 및 적응 필터부(140)의 기준 입력값(X1(n)~Xn(n))에 대한 계수(W1~Wn)값을 입력받는다.
단계(431)에서 적응 필터 적용 여부 판단부(130)는 상기 원 데이터(D(n))에 대한 절대값(Abs[D(n)])을 구한다. 단계(432)에서 적응 필터 적용 여부 판단부(130)는 상기 출력값(Y(n))에 대한 절대값(Abs[Y(n)])을 구한다.
단계(441)에서 적응 필터 적용 여부 판단부(130)는 상기 원 데이터의 절대값(Abs[D(n)])에 대해 저주파 컷오프(cutoff)를 갖는 저주파 패스 필터(LPF: Low Pass Filter)를 통과시키거나 또는 이동 평균(Moving average)값(D(n)')을 구한다. 단계(442)에서 적응 필터 적용 여부 판단부(130)는 상기 출력값(Y(n))에 대한 절대값(Abs[Y(n)])에 대해 저주파 컷오프(cutoff)를 갖는 저주파 패스 필터(LPF)를 통과시키거나 또는 이동 평균(Moving average)값(Y(n)')을 구한다.
단계(450)에서 적응 필터 적용 여부 판단부(130)는 단계(441)에서 구한 값(D(n)')과 단계(442)에서 구한 값(Y(n)')에 대한 출력비(R(n)=Y(n)'/D(n)')를 구한다. 적응 필터 적용 여부 판단부(130)는 상기 출력비(R(n))를 이용하여 적응 필터부(140)에 대한 동작 적용 여부를 판단한다.
단계(460)에서 적응 필터 적용 여부 판단부(130)는 하기 수학식 1을 이용하여 NW를 구하거나 또는 하기 수학식 2를 이용하여 ACC 비율(ratio)를 구한다.
Figure 112006095189967-pat00001
ACC 비율: X1 비율, X2 비율,
여기서, X1 비율: ave(abs(X1-W))/ave(abs(X1)),
X2 비율: ave(abs(X2-W))/ave(abs(X2)), ...
Xn 비율: ave(abs(Xn-W))/ave(abs(Xn))
abs : 절대값(absolute), ave : 이동 평균값(Moving average)
단계(470)에서 적응 필터 적용 여부 판단부(130)는 상기 출력비(R(n))와 제1 임계값(threshold 1)을 비교하여 상기 출력비(R(n))가 상기 제1 임계값(threshold 1)보다 작은지 여부를 판단한다.
즉, 단계(470)에서 적응 필터 적용 여부 판단부(130)는 상기 출력비(R(n))와 제1 임계값(threshold 1)을 비교한 결과로 상기 출력비(R(n))가 제1 임계값 이상이면, 상기 적응 필터부에 의한 적응 필터링 동작을 적용시키는 것으로 판단할 수 있다.
한편, 상기 출력비(R(n))가 상기 제1 임계값(threshold 1)보다 작은 경우, 단계(475)에서 적응 필터 적용 여부 판단부(130)는 상기 NW 또는 상기 ACC 비율이 제2 임계값보다 큰지 여부를 판단한다.
즉, 적응 필터 적용 여부 판단부(130)는 상기 NW 또는 상기 ACC 비율이 상기 제2 임계값보다 크지 않으면, 적응 필터부(140)에 의한 적응 필터링 동작을 적용시키는 것으로 판단한다. 이와 달리, 적응 필터 적용 여부 판단부(130)는 NW 또는 상기 ACC 비율이 상기 제2 임계값보다 크면, 적응 필터부(140)에 의한 적응 필터링 동작을 적용시키지 않는 것으로 판단한다.
적응 필터부(140)에 의해 적응 필터링 동작을 적용시키는 것으로 판단하면, 단계(480)에서 적응 필터 적용 여부 판단부(130)는 적응 필터부(140)에 의해 적응 필터링 동작이 적용된 출력 신호(Y(n))가 출력되도록 한다.
적응 필터부(140)에 의해 적응 필터링 동작을 적용시키지 않는 것으로 판단하면, 단계(490)에서 적응 필터 적용 여부 판단부(130)는 적응 필터부(140)에 의해 적응 필터링 동작이 적용되지 않은 상기 원 신호(D(n))가 출력되도록 한다.
신호 처리부(150)는 적응 필터 적용 여부 판단부(130)로부터 출력된 상기 원 신호(D(n)) 또는 적응 필터부(140)로부터 출력된 상기 적응 필터링 동작이 적용된 출력 신호(Y(n))를 증폭시키는 신호 처리 동작을 수행한다.
신호 출력부(160)는 상기 증폭된 원 신호(D(n)) 또는 상기 증폭된 출력 신호(Y(n))를 출력한다.
도 5는 측정하고자 하는 신호와 노이즈 신호의 주기가 일치하는 경우 R(n)에 대한 변화를 나타내는 도면이다.
도 5를 참조하면, 참조부호(510)는 적응 필터부(140)에 의한 적응 필터링 동작이 적용되기 전 신호를 나타내고, 참조부호(520)는 적응 필터부(140)에 의한 적응 필터링 동작이 적용된 후 신호를 나타내고, 참조부호(530)는 상기 출력값(R(n))의 변화를 나타낸다. 상기 출력값(R(n))이 측정하고자 하는 신호와 노이즈 신호의 주기가 일치하지 않는 경우 1에 근접한 값을 가지게 되나 일치하는 경우 그 값이 현저하게 작아지게 된다.
그러므로, 본 발명에 따른 적응 필터 적용 여부 판단 시스템(100)은 상기 출력값(R(n))이 제1 임계값 이하로 내려갈 경우 적응 필터부(140)에 의한 적응 필터링 동작을 적용하지 않고, 상기 제1 임계값 이상인 경우 상기 적응 필터링 동작 을 적용함으로써 종래 적응 필터가 가지는 문제점을 극복할 수 있다.
한편, 측정하고자 하는 신호와 노이즈 신호의 패턴이 일정하지 않더라도 측정하고자 하는 신호보다 큰 노이즈 신호가 포함된 경우, 일례로 피검자가 뛰는 경우에 적응 필터 적용 시 R(n)값이 작아지는 문제점이 있다. 따라서, 적응 필터 적용 여부 판단부(130)는 각각의 입력 기준 신호에 연산되어진 계수(W1,W2,...Wn)의 Norm(
Figure 112006095189967-pat00002
)값 즉, 계수 벡터(coefficient vector)의 크기(NW)와 ACC 비율을 구한 후 상기 NW 또는 상기 ACC 비율이 제2 임계치를 넘으면, 적응 필터링 동작을 적용하지 않고, 상기 NW 또는 상기 ACC 비율이 상기 제2 임계치를 넘지 않으면 적응 필터링 동작을 적용시킴으로써 상기 문제점을 보완할 수 있다.
도 6은 R(n)에 대한 변화 그래프 및 NW에 대한 변화 그래프의 일례를 나타내는 도면이다.
도 6을 참조하면, 참조부호(610)는 상기 출력비(R(n))에 대한 변화 그래프를 나타내고, 참조부호(620)는 상기 NW에 대한 변화 그래프를 나타낸다. 여기서, 빗금친 구간이 상기 측정하고자 하는 신호와 노이즈 신호의 패턴이 일치되는 구간이다.
도 7은 ACC 비율에 대한 변화 그래프 및 NW에 대한 변화 그래프의 일례를 나타내는 도면이다.
도 7을 참조하면, 참조부호(710)는 상기 ACC 비율에 대한 변화 그래프를 나타내고, 참조부호(720)는 상기 NW에 대한 변화 그래프를 나타낸다. 여기서, 굵은 점선 사이 구간이 상기 가속도 신호가 동기되는 구간이다.
도 8은 적응 필터 적용 여부 판단 방법에 의한 알고리즘을 수행한 결과의 일례를 나타내는 도면이다.
도 8을 참조하면, 참조부호(810)는 원(Raw) PPG 신호를 나타내고, 참조부호(820)는 상기 PPG 신호와 상기 가속도 신호가 동일한 주기일 경우, 적응 필터링 동작이 적용됨에 따라 상기 원 PPG 신호가 복원되지 않은 상태를 나타낸다. 참조부호(830)는 본 발명에 따른 적응 필터 적용 여부 판단 방법의 알고리즘이 실제로 적용됨에 따라 상기 PPG 신호와 상기 가속도 신호가 중첩될 경우, 적응 필터링 동작이 적용되지 않음으로써 상기 원 PPG 신호가 유지된 상태를 나타낸다. 참조부호(831)는 적응 필터부(140)의 온/오프 동작에 의해 발생되는 노이즈를 나타낸다.
이와 같이, 본 발명에 따른 적응 필터 적용 여부 판단 방법은 종래 일률적으로 적응 필터링 동작을 적용하지 않고 상기 PPG 신호와 상기 가속도 신호를 분석한 결과에 따라 적응 필터링 동작에 대한 적용 여부를 판단함으로써 상기 PPG 신호와 상기 가속도 신호가 동일한 주기인 경우 원 신호를 그대로 반영함으로써 종래보다 PPG 신호에 대한 S/N비를 높일 수 있다.
상술한 바와 같이 본 발명의 상세한 설명에서는 운동 중의 맥박 측정 시스템을 기준으로 설명하였으나 본 발명은 맥박 측정 시스템에서 측정되는 PPG 신호뿐만 아니라 저주파 신호를 측정하는 시스템에 모두 적용될 수 있다.
한편 본 발명에 따른 적응 필터 적용 여부 방법은 다양한 컴퓨터로 구현되는 동작을 수행하기 위한 프로그램 명령을 포함하는 컴퓨터 판독 가능 매체를 포함 한다. 상기 컴퓨터 판독 가능 매체는 프로그램 명령, 데이터 파일, 데이터 구조 등을 단독으로 또는 조합하여 포함할 수 있다. 상기 매체는 프로그램 명령은 본 발명을 위하여 특별히 설계되고 구성된 것들이거나 컴퓨터 소프트웨어 당업자에게 공지되어 사용 가능한 것일 수도 있다. 컴퓨터 판독 가능 기록 매체의 예에는 하드 디스크, 플로피 디스크 및 자기 테이프와 같은 자기 매체(magnetic media), CD-ROM, DVD와 같은 광기록 매체(optical media), 플롭티컬 디스크(floptical disk)와 같은 자기-광 매체(magneto-optical media), 및 롬(ROM), 램(RAM), 플래시 메모리 등과 같은 프로그램 명령을 저장하고 수행하도록 특별히 구성된 하드웨어 장치가 포함된다. 프로그램 명령의 예에는 컴파일러에 의해 만들어지는 것과 같은 기계어 코드뿐만 아니라 인터프리터 등을 사용해서 컴퓨터에 의해서 실행될 수 있는 고급 언어 코드를 포함한다.
이상과 같이 본 발명은 비록 한정된 실시예와 도면에 의해 설명되었으나, 본 발명은 상기의 실시예에 한정되는 것은 아니며, 이는 본 발명이 속하는 분야에서 통상의 지식을 가진 자라면 이러한 기재로부터 다양한 수정 및 변형이 가능하다. 따라서, 본 발명 사상은 아래에 기재된 특허청구범위에 의해서만 파악되어야 하고, 이의 균등 또는 등가적 변형 모두는 본 발명 사상의 범주에 속한다고 할 것이다.
본 발명에 따르면, 측정하고자 하는 신호에 영향을 주는 복수 개의 노이즈 신호의 영향을 파악하여 노이즈를 제거하는 적응 필터를 적용할지 여부를 판단하는 시스템 및 그 방법을을 제공할 수 있다.
또한 본 발명에 따르면, 실시간 처리가 가능하고 입력 채널의 장점을 가지는 적응 필터를 사용할 때 생체 신호와 제거하고자 하는 노이즈 패턴이 동일한 경우 발생되는 문제점을 원천적으로 해결할 수 있는 적응 필터 적용 여부 판단 시스템 및 그 방법을 제공할 수 있다.

Claims (13)

  1. 검출 신호를 센싱하는 신호 센싱부;
    상기 센싱된 신호를 적응적으로 필터링하는 적응 필터부; 및
    상기 센싱된 신호 및 상기 필터링된 신호를 분석하고, 상기 분석된 결과에 기초하여 상기 센싱된 신호에 대한 상기 적응 필터부의 적용 여부를 판단하는 적응 필터 적용 여부 판단부
    를 포함하는 적응 필터 적용 판단 시스템.
  2. 제1항에 있어서,
    상기 적응 필터 적용 여부 판단부는,
    상기 적응 필터부의 출력값(Y(n)), 상기 적응 필터부에 입력되는 원 데이터(D(n)) 및 상기 적응 필터부의 기준 입력(X1(n)~Xn(n))의 계수(W1~Wn)값을 입력 받고, 상기 출력값(Y(n))과 상기 원 데이터(D(n))의 절대값(Abs[Y(n)], Abs[D(n)])을 구한 후 각각 저주파 컷오프(cutoff)를 갖는 저주파 패스 필터(LPF)를 적용시키거나 이동 평균값을 구한 결과값(Y(n)', D(n)')의 출력비(R(n)=Y(n)'/D(n)')을 구한 후 상기 출력비(R(n)에 기초하여 상기 적응 필터부에 의한 적응 필터링 동작에 대한 적용 여부를 판단하는 적응 필터 적용 판단 시스템.
  3. 제2항에 있어서,
    상기 적응 필터 적용 여부 판단부는,
    상기 출력비(R(n))가 제1 임계값 이상이면, 상기 적응 필터부에 의한 적응 필터링 동작을 적용시키는 것으로 판단하는 적응 필터 적용 판단 시스템.
  4. 제3항에 있어서,
    상기 적응 필터 적용 여부 판단부는,
    상기 출력비(R(n))가 상기 제1 임계값보다 작은 경우, NW 또는 ACC 비율이 제2 임계값보다 크지 않으면, 상기 적응 필터부에 의한 적응 필터링 동작을 적용시키는 것으로 판단하는 적응 필터 적용 판단 시스템.
  5. 제4항에 있어서,
    상기 적응 필터 적용 여부 판단부는,
    상기 NW 또는 상기 ACC 비율이 상기 제2 임계값 큰 경우, 상기 원 데이터(D(n))를 출력하도록 하기 위해 상기 적응 필터부에 의한 적응 필터링 동작을 적용시키지 않는 것으로 판단하는 적용 필터 적용 판단 시스템.
  6. 제2항 내지 제5항 중 어느 한 항에 있어서,
    상기 적응 필터 적용 여부 판단부는,
    상기 적용 필터링 적용 여부 판단 결과에 기초하여 상기 적응 필터에 의해 적용 필터링 동작이 적용된 신호(Y(n))를 출력하거나 또는 상기 원 데이터(D(n))를 출력하는 적용 필터 적용 판단 시스템.
  7. 각종 센서를 통해 측정하고자 하는 신호를 검출하는 단계;
    적응 필터를 통해 상기 검출된 신호에 대해 적응 필터링 동작을 수행하는 단계; 및
    상기 검출된 신호 및 상기 적응 필터링된 신호를 분석하고, 상기 분석된 결과에 기초하여 상기 검출된 신호에 대한 적용 필터링 적용 여부를 판단하는 단계
    를 포함하는 적응 필터 적용 판단 방법.
  8. 제7항에 있어서,
    상기 검출된 신호에 대한 적용 필터링 적용 여부를 판단하는 단계는,
    상기 적응 필터의 출력값(Y(n)), 상기 적응 필터에 입력되는 원 데이터(D(n)) 및 상기 적응 필터의 레퍼런스 입력(X1(n)~Xn(n))의 계수(W1~Wn)값을 입력 받는 단계;
    상기 출력값(Y(n))의 절대값과 상기 원 데이터(D(n))의 절대값을 각각 구하는 단계;
    상기 절대값을 각각 저주파 컷오프(cutoff)를 갖는 저주파 패스 필터(LPF)를 통과시키는 단계;
    입력 신호 대비 상기 적응 필터의 출력비(R(n))을 구하는 단계; 및
    상기 출력비를 기초하여 상기 적응 필터에 의한 상기 검출된 신호에 대한 적 응 필터링 적용 여부를 판단하는 단계
    를 포함하는 적응 필터 적용 판단 방법.
  9. 제8항에 있어서,
    상기 검출된 신호에 대한 적용 필터링 적용 여부를 판단하는 단계는,
    상기 출력비(R(n))가 제1 임계값 이상이면, 상기 적응 필터에 의한 적응 필터링 동작을 적용시키는 것으로 판단하는 적응 필터 적용 판단 방법.
  10. 제9항에 있어서,
    상기 검출된 신호에 대한 적용 필터링 적용 여부를 판단하는 단계는,
    상기 출력비(R(n))가 상기 제1 임계값보다 작은 경우, NW 또는 ACC 비율이 제2 임계값보다 크지 않으면, 상기 적응 필터부에 의한 적응 필터링 동작을 적용시키는 것으로 판단하는 적응 필터 적용 판단 방법.
  11. 제10항에 있어서,
    상기 검출된 신호에 대한 적용 필터링 적용 여부를 판단하는 단계는,
    상기 NW 또는 상기 ACC 비율이 상기 제2 임계값보다 큰 경우, 상기 원 데이터(D(n))를 출력하도록 하기 위해 상기 적응 필터부에 의한 적응 필터링 동작을 적용시키지 않는 것으로 판단하는 적용 필터 적용 판단 방법.
  12. 제8항 내지 제11항 중 어느 한 항에 있어서,
    상기 적용 필터링 적용 여부 판단 결과에 기초하여 상기 적응 필터에 의해 적용 필터링 동작이 적용된 신호(Y(n))를 출력하거나 또는 상기 원 데이터(D(n))를 출력하는 단계를 더 포함하는 적용 필터 적용 판단 방법.
  13. 제7항 내지 제11항 중 어느 한 항의 방법을 컴퓨터에서 실행하기 위한 프로그램을 기록하는 컴퓨터 판독 가능한 기록매체.
KR1020060132214A 2006-12-21 2006-12-21 적응 필터 적용 여부 판단 시스템 및 그 방법 KR100852195B1 (ko)

Priority Applications (2)

Application Number Priority Date Filing Date Title
KR1020060132214A KR100852195B1 (ko) 2006-12-21 2006-12-21 적응 필터 적용 여부 판단 시스템 및 그 방법
US11/819,477 US8005882B2 (en) 2006-12-21 2007-06-27 System and method for determining application of adaptive filter

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020060132214A KR100852195B1 (ko) 2006-12-21 2006-12-21 적응 필터 적용 여부 판단 시스템 및 그 방법

Publications (2)

Publication Number Publication Date
KR20080058098A KR20080058098A (ko) 2008-06-25
KR100852195B1 true KR100852195B1 (ko) 2008-08-13

Family

ID=39542860

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020060132214A KR100852195B1 (ko) 2006-12-21 2006-12-21 적응 필터 적용 여부 판단 시스템 및 그 방법

Country Status (2)

Country Link
US (1) US8005882B2 (ko)
KR (1) KR100852195B1 (ko)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100995318B1 (ko) 2009-03-27 2010-11-19 서울과학기술대학교 산학협력단 신호 검출 장치 및 방법

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2009076523A1 (en) 2007-12-11 2009-06-18 Andrea Electronics Corporation Adaptive filtering in a sensor array system
US8150054B2 (en) * 2007-12-11 2012-04-03 Andrea Electronics Corporation Adaptive filter in a sensor array system
US9392360B2 (en) 2007-12-11 2016-07-12 Andrea Electronics Corporation Steerable sensor array system with video input
US20170164847A1 (en) * 2015-12-15 2017-06-15 Texas Instruments Incorporated Reducing Motion Induced Artifacts in Photoplethysmography (PPG) Signals
TWI819722B (zh) * 2022-07-26 2023-10-21 瑞昱半導體股份有限公司 用於光體積變化描記圖法裝置的數位濾波方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0316413A (ja) * 1989-06-14 1991-01-24 Matsushita Electric Ind Co Ltd 適応フィルタ装置
JPH07134061A (ja) * 1993-11-12 1995-05-23 Ishikawajima Harima Heavy Ind Co Ltd 適応フィルタを用いた信号計測方法
JPH10173489A (ja) 1996-12-12 1998-06-26 Kenwood Corp 適応フィルタ

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3522018B2 (ja) 1995-10-12 2004-04-26 富士通テン株式会社 騒音制御装置
FI973455A (fi) * 1997-08-22 1999-02-23 Nokia Mobile Phones Ltd Menetelmä ja järjestely melun vaimentamiseksi tilassa muodostamalla vastamelua
US5971930A (en) * 1997-10-17 1999-10-26 Siemens Medical Systems, Inc. Method and apparatus for removing artifact from physiological signals
US5908393A (en) * 1998-05-04 1999-06-01 Cambridge Heart, Inc. Reducing noise in a biological signal
US8275120B2 (en) * 2006-05-30 2012-09-25 Microsoft Corp. Adaptive acoustic echo cancellation
EP2197214A4 (en) * 2007-10-03 2010-12-29 Panasonic Corp METHOD AND DEVICE FOR REDUCING BLOCK DEFORMATION
US7860200B2 (en) * 2007-10-12 2010-12-28 Harris Corporation Communications system using adaptive filter that is selected based on output power

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0316413A (ja) * 1989-06-14 1991-01-24 Matsushita Electric Ind Co Ltd 適応フィルタ装置
JPH07134061A (ja) * 1993-11-12 1995-05-23 Ishikawajima Harima Heavy Ind Co Ltd 適応フィルタを用いた信号計測方法
JPH10173489A (ja) 1996-12-12 1998-06-26 Kenwood Corp 適応フィルタ

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100995318B1 (ko) 2009-03-27 2010-11-19 서울과학기술대학교 산학협력단 신호 검출 장치 및 방법

Also Published As

Publication number Publication date
US20080152161A1 (en) 2008-06-26
US8005882B2 (en) 2011-08-23
KR20080058098A (ko) 2008-06-25

Similar Documents

Publication Publication Date Title
Zhang et al. Heart sound classification based on scaled spectrogram and partial least squares regression
KR100852195B1 (ko) 적응 필터 적용 여부 판단 시스템 및 그 방법
WO2008036911A3 (en) System and method for acoustic detection of coronary artery disease
US9192311B2 (en) Pulse wave detection device and pulse wave detection method
DE60331842D1 (de) Verfahren zum nachweis von stress durch segmentierung und analyse eines herzschlagsignals
WO2001025802A3 (en) Method, apparatus and system for removing motion artifacts from measurements of bodily parameters
JPWO2011033784A1 (ja) 信号抽出装置
KR101301304B1 (ko) 이산 웨이블릿 변환을 이용한 운동 중 실시간 특정 점 검출 및 리듬 분석
KR101704491B1 (ko) 심박수 측정 장치 및 방법, 이를 수행하기 위한 기록매체
Jain et al. A novel method for suppression of motion artifacts from the seismocardiogram signal
JP2001095768A (ja) ペースメーカが定めたペースの心電図を自動的に検出しかつ解釈する方法および装置
US20200245875A1 (en) Heartbeat detection device, heartbeat detection method, and program
US9772368B2 (en) Detection of an abnormal signal in a compound sampled signal
CN112315440A (zh) 心率检测方法、可穿戴设备及可读存储介质
Ukil et al. Cardiac condition monitoring through photoplethysmogram signal denoising using wearables: can we detect coronary artery disease with higher performance efficacy?
WO2002096266A3 (en) Method and system for optically performing an assay to determine a medical condition
JP3687135B2 (ja) 音振動評価装置
JP2001070265A (ja) 脈波解析法
Yu et al. Heart beat detection method with estimation of regular intervals between ECG and blood pressure
US10993676B2 (en) Signal processing method and apparatus
Zaeni et al. Implementation of adaptive threshold for peak detection of photoplethysmography applied on microcontroller
JP2002224061A (ja) 電子血圧計
JP7224163B2 (ja) 心拍周期測定装置及び心拍周期測定プログラム
KR20170050033A (ko) 맥의 긴장도를 검출하는 장치 및 방법
JP7267002B2 (ja) 呼吸周期測定装置及び呼吸周期測定プログラム

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
GRNT Written decision to grant
FPAY Annual fee payment

Payment date: 20120716

Year of fee payment: 5

FPAY Annual fee payment

Payment date: 20130724

Year of fee payment: 6

FPAY Annual fee payment

Payment date: 20140721

Year of fee payment: 7

FPAY Annual fee payment

Payment date: 20150716

Year of fee payment: 8

FPAY Annual fee payment

Payment date: 20160718

Year of fee payment: 9

FPAY Annual fee payment

Payment date: 20170719

Year of fee payment: 10

FPAY Annual fee payment

Payment date: 20180717

Year of fee payment: 11

FPAY Annual fee payment

Payment date: 20190717

Year of fee payment: 12