KR20070069637A - 직교필터를 이용한 주파수 추정 장치 및 그 방법 - Google Patents

직교필터를 이용한 주파수 추정 장치 및 그 방법 Download PDF

Info

Publication number
KR20070069637A
KR20070069637A KR1020050131977A KR20050131977A KR20070069637A KR 20070069637 A KR20070069637 A KR 20070069637A KR 1020050131977 A KR1020050131977 A KR 1020050131977A KR 20050131977 A KR20050131977 A KR 20050131977A KR 20070069637 A KR20070069637 A KR 20070069637A
Authority
KR
South Korea
Prior art keywords
filter
frequency
sine
cosine
digital signal
Prior art date
Application number
KR1020050131977A
Other languages
English (en)
Other versions
KR100758492B1 (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 KR1020050131977A priority Critical patent/KR100758492B1/ko
Publication of KR20070069637A publication Critical patent/KR20070069637A/ko
Application granted granted Critical
Publication of KR100758492B1 publication Critical patent/KR100758492B1/ko

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
    • G01R23/16—Spectrum analysis; Fourier analysis
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
    • G01R23/16—Spectrum analysis; Fourier analysis
    • G01R23/165—Spectrum analysis; Fourier analysis using filters

Landscapes

  • Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • General Physics & Mathematics (AREA)
  • Measuring Frequencies, Analyzing Spectra (AREA)

Abstract

본 발명은 직교필터를 이용한 주파수 추정 장치 및 그 방법에 관한 것이다.
본 발명에 따른 주파수 추정 장치는, 입력되는 전류 또는 전압값을 디지털신호로 변환하는 A/D 변환부와, A/D변환부에서 변환된 디지털신호로부터 코사인(cosine) 신호를 추출해 내는 코사인 필터 및 A/D변환부에서 변환된 디지털신호로부터 사인(sine) 신호를 추출해 내는 사인 필터로 구성된 직교필터와, 코사인 필터 및 사인 필터에서 출력되는 신호를 일정 시간 간격으로 샘플링하고, 샘플링한 신호를 이용하여 입력된 신호에 대한 주파수를 추정하는 주파수 추정부를 포함하여 이루어진 것을 특징으로 한다.
주파수 추정, 직교필터, 코사인, 사인, 샘플링

Description

직교필터를 이용한 주파수 추정 장치 및 그 방법{Apparatus and method for estimating frequency using orthogonal filter}
도 1a 및 1b는 종래 영점 교차법을 이용한 주파수 추정 방법을 설명하기 위해 도시한 도면,
도 2는 종래 영점 교차법을 이용하여 주파수 추정시 발생하는 문제점을 설명하기 위해 도시한 도면,
도 3은 본 발명의 바람직한 실시예에 따른 주파수 추정 장치에 대한 개략적인 블록도,
도 4는 도 3에 도시된 코사인 필터와 사인 필터의 시간영역에서의 특성을 도시한 도면,
도 5는 도 3에 도시된 코사인 필터와 사인 필터의 주파수 응답특성을 설명하기 위해 도시한 도면, 그리고,
도 6 및 도 7은 도 3에 도시된 코사인 필터와 사인 필터에서 출력되는 신호의 궤적을 복소평면상에 나타낸 도면이다.
* 도면의 주요 부분에 대한 부호의 설명 *
100 : 주파수 추정 장치 110 : A/D 변환부
120 : 직교필터 122 : 코사인 필터
124 : 사인 필터 130 : 주파수 추정부
본 발명은 주파수 추정 장치 및 그 방법에 관한 것으로, 보다 상세하게는, 직교필터를 이용하여 전력계통에 입력되는 신호에 대한 주파수를 정밀하게 추정할 수 있도록 한 주파수 추정 장치 및 그 방법에 관한 것이다.
일반적으로, 보호 계전기나 페이저 측정 장치(PMU : Phasor Measurement Unit) 등과 같은 전력기기는 전원장치로부터 입력되는 전압 또는 전류에 대한 주파수를 추정하기 위해 주파수 추정 장치를 구비하고 있다. 상기 보호 계전기에서 상기 주파수 추정 장치는 계통의 고장, 안정도 판단 및 고장 판단을 위해 이용되며, 페이저 측정 장치에서는 측정된 주파수 정보를 이용하여 계통의 상태 모니터링 및 안정적인 전력계통 운영에 필요한 기본적인 데이터로 사용한다.
종래 주파수 추정장치는 영점 교차법(Zero-crossing Method)을 입력되는 신호에 대한 주파수를 추정하였다. 종래 영점 교차법을 이용한 주파수 추정 방법으로는, 도 1a에 도시된 바와 같이, 입력되는 전압 또는 전류의 신호파형에서 영점을 교차하는 시점(T1, T2)을 검출하고, 상기 검출된 영점 교차 시점(T1, T2) 사이의 간격을 계산하여 주파수를 계산하는 방식과, 도 1b에 도시된 바와 같이 입력되는 신호파형을 일정한 주기로 샘플링(Tn-1, Tn)하고, 상기 샘플링한 두 점(Tn-1, Tn) 사이의 파형을 선형이라고 가정한 후, 직선의 방정식을 이용하여 영점 교차점(Tzc) 을 계산하는 방식이 있다.
그런데, 상기 영점 교차법을 이용한 주파수 추정 방식은 정확한 영점 교차점을 검출하기 어려워 계산된 주파수와 실제 주파수 사이에 큰 오차가 발생하는 문제점이 있었다. 상기와 같은 오차를 줄이기 위해서는 상기 샘플링된 두 점(Tn-1, Tn) 사이의 파형을 2차 함수 또는 3차 함수로 가정하는 2차 또는 3차 보간법을 이용하면 된다. 그러나, 상기 2차 또는 3차 보간법을 이용하는 경우 차수의 증가에 따라 연산량이 기하급수적으로 증가하게 된다. 이러한 연산량의 증가는 주파수 추정 장치에 큰 부담으로 작용한다.
또한, 상기 영점 교차법은 원 신호에 포함되어 입력되는 잡음에 대해 매우 취약한 구조를 갖는다. 즉, 노이즈가 포함된 신호가 인가되는 경우, 도 2에 도시된 바와 같이, 여러 개의 영점 교차점(Tzc1, Tzc2, Tzc3)이 발생될 수 있으므로 영점 교차점이 검출이 어려운 문제점이 있다.
따라서, 본 발명은 상기와 같은 문제점을 해결하기 위해 창안된 것으로, 본 발명의 목적은 직교필터를 이용하여 입력되는 신호를 직교성분을 갖는 두 개의 신호로 분리한 후 일정한 주기로 샘플링하고, 샘플링된 두 개의 신호를 이용하여 주파수를 정밀하게 추정할 수 있는 직교필터를 이용한 주파수 추정 장치 및 그 방법을 제공하는 데 있다.
상기와 같은 기술적 과제를 해결하기 위한, 본 발명에 따른 주파수 추정 장치는, 입력되는 전류 또는 전압값을 디지털신호로 변환하는 A/D 변환부; 상기 A/D변환부에서 변환된 디지털신호로부터 코사인(cosine) 신호를 추출해 내는 코사인 필터; 상기 A/D변환부에서 변환된 디지털신호로부터 사인(sine) 신호를 추출해 내는 사인 필터; 및 상기 코사인 필터 및 상기 사인 필터에서 출력되는 신호를 일정 시간 간격으로 샘플링하고, 상기 샘플링한 신호를 이용하여 상기 디지털신호에 대한 주파수를 추정하는 주파수 추정부를 포함하여 이루어진 것을 특징으로 한다.
상기 주파수 추정부는, 상기 코사인 필터 및 상기 사인 필터에서 출력되는 신호에서 각각 두 개의 샘플값을 취하고, 상기 두 개의 샘플값과 상기 두 개의 샘플값을 취한 시간 정보를 이용하여 상기 디지털신호에 대한 주파수를 추정하는 것을 특징으로 한다.
상기 디지털 신호에 대한 주파수는 상기 두 개의 샘플값의 위상차에 따라 다음의 수학식 1 또는 수학식 2로 계산되는 것을 특징으로 한다.
수학식 1
Figure 112005077332663-PAT00001
,
Figure 112005077332663-PAT00002
Figure 112005077332663-PAT00003
수학식 2
Figure 112005077332663-PAT00004
,
Figure 112005077332663-PAT00005
Figure 112005077332663-PAT00006
상기 수학식 1 및 2에서, N은 윈도우 크기, t1은 첫 번째 샘플링 시간, t2는 두 번째 샘플링 시간, k는 0이상의 정수값, Xre[n1]은 상기 t1에 샘플링된 상기 코사인 필터의 출력값, Xre[n2]는 t2에 샘플링된 상기 코사인 필터의 출력값, Xim[n1]은 t1에 샘플링된 사인 필터의 출력값, Xim[n2]는 t2에 샘플링된 상기 사인 필터의 출력값, T는 상기 두 개의 샘플값의 위상차.
한편, 상기와 같은 기술적 과제를 해결하기 위한, 본 발명에 따른 주파수 추정 방법은, 입력되는 전류 또는 전압값을 디지털신호로 변환하는 단계; 코사인(cosine) 필터 및 사인(sine) 필터를 이용하여 상기 디지털신호로부터 코사인 신호 및 사인 신호를 각각 추출해 내는 단계; 및 상기 코사인 필터 및 상기 사인 필터에서 출력되는 신호에서 각각 두 개의 샘플값을 취하고, 상기 두 개의 샘플값과 상기 두 개의 샘플값을 취한 시간 정보를 이용하여 상기 디지털신호에 대한 주파수를 추정하는 단계를 포함하여 이루어진 것을 특징으로 한다.
상기 디지털 신호에 대한 주파수는 상기 두 개의 샘플값의 위상차에 따라 다음의 수학식 1 또는 수학식 2로 계산되는 것을 특징으로 한다.
수학식 1
Figure 112005077332663-PAT00007
,
Figure 112005077332663-PAT00008
Figure 112005077332663-PAT00009
수학식 2
Figure 112005077332663-PAT00010
,
Figure 112005077332663-PAT00011
Figure 112005077332663-PAT00012
상기 수학식 1 및 2에서, N은 윈도우 크기, t1은 첫 번째 샘플링 시간, t2는 두 번째 샘플링 시간, k는 0이상의 정수값, Xre[n1]은 상기 t1에 샘플링된 상기 코사인 필터의 출력값, Xre[n2]는 t2에 샘플링된 상기 코사인 필터의 출력값, Xim[n1]은 t1에 샘플링된 사인 필터의 출력값, Xim[n2]는 t2에 샘플링된 상기 사인 필터의 출력값, T는 상기 두 개의 샘플값의 위상차.
이하에서는 첨부된 도면들을 참조하여 본 발명을 보다 상세하게 설명한다. 다만, 본 발명을 설명함에 있어, 관련된 공지 기능 혹은 구성에 대한 구체적인 설명이 본 발명의 요지를 불필요하게 흐릴 수 있다고 판단되는 경우 그에 대한 상세한 설명은 생략한다.
도 3은 본 발명의 바람직한 실시예에 따른 주파수 추정 장치에 대한 개략적인 블록도이다.
도 3에 도시된 바와 같이, 본 발명에 따른 주파수 추정 장치(100)는 A/D 변환부(110), 직교필터(120) 및 주파수 추정부(130)를 포함하여 구성된다.
상기 A/D 변환부(110)는 외부로부터 입력되는 아날로그 전압 또는 전류 신호를 디지털신호로 변환하여 출력한다.
상기 직교필터(120)는 상기 A/D 변환부(110)로부터 입력되는 디지털신호를 직교성분을 갖는 두 개의 신호로 분리하여 출력한다. 상기 직교필터(120)는 코사인 필터(cosine Filter)(122)와 사인 필터(Sine Filter)(124)로 구성되어 있다.
상기 코사인 필터(122)는 상기 A/D 변환부(110)에 의해 변환된 디지털신호로부터 코사인 신호를 추출하여 출력한다.
상기 사인 필터(124)는 상기 A/D 변환부(110)에 의해 변환된 디지털신호로부터 사인 신호를 추출하여 출력한다. 상기 코사인 필터(122)와 사인 필터(124)에서 출력되는 신호는 90도의 위상차를 갖는다.
상기 코사인 필터(122) 및 사인 필터(124)에서 출력되는 신호 Xre[n]과 Xim[n]은 아래 수학식 1과 같이, 유한 임펄스 응답(FIR : Finite Impulse Response) 특성을 갖는 필터이다.
Figure 112005077332663-PAT00013
Figure 112005077332663-PAT00014
여기서, x는 입력신호이고, N은 필터의 윈도우 크기이다. 도 4에서와 같이 N 의 크기를 가진 윈도우는 시간에 따라 이동하며 윈도우 내의 샘플값 x를 이용하여 상기 수학식 1 및 2를 계산한다.
상기 코사인 필터(122) 및 상기 사인 필터(124)에서 출력되는 신호 Xre[n]과 Xim[n]은 도 5에서와 같은 주파수 응답특성을 갖는다. 즉, 주파수가 기본 주파수인 60Hz일 때에는 상기 코사인 필터(122)와 사인 필터(124)의 주파수 응답특성의 크기가 1로 동일하나, 주파수에 따라 크기 응답특성은 달라진다. 예를 들면, 주파수가 기본 주파수보다 작을 때는 상기 사인 필터(124)의 크기 응답특성이 상기 코사인 필터(122)의 크기 응답특성보다 큰 반면, 주파수가 기본 주파수보다 클 때에는 상기 코사인 필터(122)의 크기 응답특성이 상기 사인 필터(124)의 주파수 응답특성 보다 크다.
이하에서는 상기 코사인 필터(122) 및 사인 필터(124)로 입력되는 신호, x[n]이 수학식 2와 같은 주파수 f이고, 위상각 φ인 코사인 신호라고 가정한다.
Figure 112005077332663-PAT00015
여기서, A는 입력되는 신호의 크기이고, N은 윈도우 크기이다.
상기 코사인 필터(122) 및 사인 필터(124)로 상기 수학식 3과 같은 신호가 입력될 경우, 코사인 필터(122)와 사인 필터(124)에서 출력되는 신호 Xre[n], Xim[n]은 다음 수학식 3과 같다.
Figure 112005077332663-PAT00016
Figure 112005077332663-PAT00017
여기서, Hc(f)는 코사인 필터(122)의 크기 응답특성이고, Hs(f)는 사인 필터(124)의 주파수 응답특성 크기이다.
상기 코사인 필터(122)와 사인 필터(124)에서 출력되는 신호 Xre[n], Xim[n]은 90도의 위상차를 가지므로, 아래의 수학식 4가 성립된다.
Figure 112005077332663-PAT00018
상기 수학식 4를 상기 수학식 3에 대입하여 정리하면 아래 수학식 5와 같은 코사인 함수 및 사인 함수로 나타낼 수 있다.
Figure 112005077332663-PAT00019
Figure 112005077332663-PAT00020
상기에서 코사인 함수와 사인 함수의 내부 변수를 Ω이라고 두고, 코사인 필터(122)와 사인 필터(124)에서 출력되는 신호 Xre[n]과 Xim[n]의 궤적을 복소평면 상에 도시하면 도 6과 같다. 도 6에 도시된 바와 같이, 궤적과 실수축(Re) 및 허수 축(Im)과 만나는 지점은 코사인 필터(122)의 크기 응답특성(Hc(f))과 사인 필터(124)의 크기 응답특성(Hs(f))에 따라 달라진다. 그리고, 상기 코사인 필터(122)의 크기 응답특성(Hc(f))과 사인 필터(124)의 크기 응답특성(Hs(f))은 주파수에 따라 달라지며, 그에 따라 복소평면 상의 궤적 또한 달라진다. 예를 들면, 코사인 필터(122)와 사인 필터(124)에서 출력되는 신호 Xre[n]과 Xim[n]은 주파수가 기본 주파수(60Hz)일 경우에는 복소평면 상에 원의 궤적을 그리며, 주파수가 기본 주파수(60Hz)가 아닐 경우에는 타원의 궤적을 그린다.
상기 주파수 추정부(130)는 상기 코사인 필터(122) 및 사인 필터(124)에서 출력되는 신호로부터 두 개의 샘플값(X1, X2)을 취한다. 상기 두 개의 샘플값(X1, X2)을 복소평면 상에 나타내면 도 7과 같다. 도 7에서 X1은 첫 번째 샘플링 시간(t1)에 샘플링된 값이고, X2는 두 번째 샘플링 시간(t2)에 샘플링된 값이다. 상기 주파수 추정부(130)는 상기 샘플링한 두 개의 샘플값(X1, X2)을 이용하여 코사인 필터(122) 및 사인 필터(124)에서 출력되는 신호의 궤적 및 주파수를 추정할 수 있다.
상기 두 샘플링 시간(t1, t2)에 샘플링된 두 개의 샘플값 X1과 X2는 각각 아래의 수학식 6과 7을 만족한다.
Figure 112005077332663-PAT00021
Figure 112005077332663-PAT00022
상기에서 Xre[n1]은 상기 첫 번째 샘플링 시간(t1)에 샘플링된 코사인 필터(122)의 출력값이고, Xim[n1]은 상기 첫 번째 샘플링 시간(t1)에 샘플링된 사인 필터(124)의 출력값이다.
Figure 112005077332663-PAT00023
Figure 112005077332663-PAT00024
상기에서 Xre[n2]는 상기 두 번째 샘플링 시간(t2)에 샘플링된 코사인 필터(122)의 출력값이고, Xim[n2]는 상기 두 번째 샘플링 시간(t2)에 샘플링된 사인 필터(124)의 출력값이다.
상식 수학식 6 및 7에서 미지수는 A, f, φ,φc으로, 연립방정식을 풀면 주파수(f)를 계산할 수 있다. 상기 수학식 6, 7을 정리하면 아래 수학식 8과 같다.
Figure 112005077332663-PAT00025
Figure 112005077332663-PAT00026
Figure 112005077332663-PAT00027
,
Figure 112005077332663-PAT00028
상기 수학식 8에서 샘플링 시간(t1, t2)에 따른 샘플값 Xre,와 Xim을 나누면 다음의 수학식 9와 같다.
Figure 112005077332663-PAT00029
Figure 112005077332663-PAT00030
상기 수학식 9에서 tanΦ를 소거하여 식을 정리하면 아래 수학식 10과 같다.
Figure 112005077332663-PAT00031
Figure 112005077332663-PAT00032
상기 수학식 10을 정리하면 입력 신호에 대한 주파수(f)를 구할 수 있다.
상기 주파수 추정부(130)는 코사인 필터(122) 및 사인 필터(124)로부터 취한 상기 두 개의 샘플값의 위상차(T)에 따라 아래의 수학식 11 또는 12를 이용하여 입력된 신호에 대한 주파수(f)를 계산하게 된다.
Figure 112005077332663-PAT00033
Figure 112005077332663-PAT00034
,
Figure 112005077332663-PAT00035
Figure 112005077332663-PAT00036
,
상기 수학식 11 및 12에서, N은 윈도우 크기, t1은 첫 번째 샘플링 시간, t2는 두 번째 샘플링 시간, k는 0이상의 정수값, T는 두 개의 샘플값의 위상차이다.
상기에서 cos-1α값은 일정한 주기성을 갖고 반복되므로, 상기 cos-1α값이 어느 범위에 속하는지를 알아내야 한다. 상기 cos-1α값이 어느 범위에 속하는지의 여부는 상기 두 개의 샘플값의 위상차(T)를 통해 알 수 있다. 그리고, 상기 두 개의 샘플값은 상기 주파수 추정부(130)에서 취해지는 값이므로, 상기 주파수 추정부(130)는 추정을 통해 상기 두 개의 샘플값의 위상차(T)를 알 수 있다. 상기 두 개의 샘플값의 위상차(T)가 결정되면, 상기 k값이 정해지면 상기 cos-1α값이 코사인 함수의 어느 범위에 해당되는지 알 수 있으므로, 주파수(f)를 계산할 수 있다.
한편, 상기에서 두 개의 샘플값을 취하는 간격(t2-t1)을 증가시키는 경우 주파수 추정시 발생하는 오차를 감소시킬 수 있다. 여기서, 상기 주파수 추정시 발생하는 오차는 주파수 추정을 위해 이용되는 요소들 계산에 따른 오차로, 상기 수학식 11 또는 12에서는 'cos-1 α'값 계산시 많은 오차가 발생한다. 상기 'cos-1 α'값 계산에 따른 오차는 상기 두 개의 샘플값을 취하는 간격(t2-t1)이 클수록 감소하게 된다. 이에 따라 주파수 추정 정밀도는 증가한다. 그러나, 상기 두 개의 샘플값을 취하는 간격(t2-t1)이 증가하게 되면 주파수 추정 속도는 감소하게 되므로, 주파수 추정 밀도와 추정 속도를 고려하여 적절한 샘플링 주기를 설정하여야 한다.
상기와 같은 과정을 통해 추정된 주파수는, 보호 계전기를 구성하는 주파수 계전요소(예컨대, Underfrequency, Overfrequency, df/dt 등)의 동작을 위한 기본 동작값으로 사용하거나, 페이저 측정 시스템에서 페이저 측정수단에 의해 측정된 페이저를 보정하는 데 사용할 수 있다.
한편, 이상에서 대표적인 실시예를 통하여 본 발명에 대하여 상세하게 설명하였으나, 본 발명이 속하는 기술분야에서 통상의 지식을 가진 자는 상술한 실시예에 대하여 본 발명의 범주에서 벗어나지 않는 한도 내에서 다양한 변형이 가능함을 이해할 것이다. 그러므로 본 발명의 권리범위는 설명된 실시예에 국한되어 정해져서는 안 되며, 후술하는 특허청구범위 뿐만 아니라 이 특허청구범위와 균등한 것들에 의해 정해져야 한다.
지금까지 설명한 바와 같이, 본 발명에 따르면, 직교필터를 이용하여 필터링한 신호를 이용하여 주파수를 추정함으로써, 신호에 포함된 잡음에 의해 발생하는 오차를 감소시킬 수 있다. 또한, 두 개의 샘플값을 취하는 간격을 변화시키는 경우 주파수 추정 정밀도와 추정 속도를 변화시킬 수 있다. 따라서, 설계조건에 따라 주파수 추정 특성을 가변시킬 수 있다.

Claims (5)

  1. 입력되는 전류 또는 전압값을 디지털신호로 변환하는 A/D 변환부;
    상기 A/D변환부에서 변환된 디지털신호로부터 코사인(cosine) 신호를 추출해 내는 코사인 필터;
    상기 A/D변환부에서 변환된 디지털신호로부터 사인(sine) 신호를 추출해 내는 사인 필터; 및
    상기 코사인 필터 및 상기 사인 필터에서 출력되는 신호를 일정 시간 간격으로 샘플링하고, 상기 샘플링한 신호를 이용하여 상기 디지털신호에 대한 주파수를 추정하는 주파수 추정부;를 포함하여 이루어진 주파수 추정장치.
  2. 제 1항에 있어서,
    상기 주파수 추정부는, 상기 코사인 필터 및 상기 사인 필터에서 출력되는 신호에서 각각 두 개의 샘플값을 취하고, 상기 두 개의 샘플값과 상기 두 개의 샘플값을 취한 시간 정보를 이용하여 상기 디지털신호에 대한 주파수를 추정하는 것을 특징으로 하는 주파수 추정 장치.
  3. 제 2항에 있어서,
    상기 디지털 신호에 대한 주파수는 상기 두 개의 샘플값의 위상차에 따라 다음의 수학식 1 또는 수학식 2로 계산되는 것을 특징으로 하는 주파수 추정 장치.
    수학식 1
    Figure 112005077332663-PAT00037
    ,
    Figure 112005077332663-PAT00038
    Figure 112005077332663-PAT00039
    수학식 2
    Figure 112005077332663-PAT00040
    ,
    Figure 112005077332663-PAT00041
    Figure 112005077332663-PAT00042
    상기 수학식 1 및 2에서, N은 윈도우 크기, t1은 첫 번째 샘플링 시간, t2는 두 번째 샘플링 시간, k는 0이상의 정수값, Xre[n1]은 상기 t1에 샘플링된 상기 코사인 필터의 출력값, Xre[n2]는 t2에 샘플링된 상기 코사인 필터의 출력값, Xim[n1]은 t1에 샘플링된 사인 필터의 출력값, Xim[n2]는 t2에 샘플링된 상기 사인 필터의 출력값, T는 상기 두 개의 샘플값의 위상차.
  4. 입력되는 전류 또는 전압값을 디지털신호로 변환하는 단계;
    코사인(cosine) 필터 및 사인(sine) 필터를 이용하여 상기 디지털신호로부터 코사인 신호 및 사인 신호를 각각 추출해 내는 단계; 및
    상기 코사인 필터 및 상기 사인 필터에서 출력되는 신호에서 각각 두 개의 샘플값을 취하고, 상기 두 개의 샘플값과 상기 두 개의 샘플값을 취한 시간 정보를 이용하여 상기 디지털신호에 대한 주파수를 추정하는 단계;를 포함하여 이루어진 주파수 추정 방법.
  5. 제 4항에 있어서,
    상기 디지털 신호에 대한 주파수는 상기 두 개의 샘플값의 위상차에 따라 다음의 수학식 1 또는 수학식 2로 계산되는 것을 특징으로 하는 주파수 추정 방법.
    수학식 1
    Figure 112005077332663-PAT00043
    ,
    Figure 112005077332663-PAT00044
    Figure 112005077332663-PAT00045
    수학식 2
    Figure 112005077332663-PAT00046
    ,
    Figure 112005077332663-PAT00047
    Figure 112005077332663-PAT00048
    상기 수학식 1 및 2에서, N은 윈도우 크기, t1은 첫 번째 샘플링 시간, t2는 두 번째 샘플링 시간, k는 0이상의 정수값, Xre[n1]은 상기 t1에 샘플링된 상기 코사인 필터의 출력값, Xre[n2]는 t2에 샘플링된 상기 코사인 필터의 출력값, Xim[n1]은 t1에 샘플링된 사인 필터의 출력값, Xim[n2]는 t2에 샘플링된 상기 사인 필터의 출력값, T는 상기 두 개의 샘플값의 위상차.
KR1020050131977A 2005-12-28 2005-12-28 직교필터를 이용한 주파수 추정 장치 및 그 방법 KR100758492B1 (ko)

Priority Applications (1)

Application Number Priority Date Filing Date Title
KR1020050131977A KR100758492B1 (ko) 2005-12-28 2005-12-28 직교필터를 이용한 주파수 추정 장치 및 그 방법

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020050131977A KR100758492B1 (ko) 2005-12-28 2005-12-28 직교필터를 이용한 주파수 추정 장치 및 그 방법

Publications (2)

Publication Number Publication Date
KR20070069637A true KR20070069637A (ko) 2007-07-03
KR100758492B1 KR100758492B1 (ko) 2007-09-12

Family

ID=38505227

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020050131977A KR100758492B1 (ko) 2005-12-28 2005-12-28 직교필터를 이용한 주파수 추정 장치 및 그 방법

Country Status (1)

Country Link
KR (1) KR100758492B1 (ko)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100823832B1 (ko) * 2006-12-29 2008-04-21 엘에스산전 주식회사 전력계통의 주파수 추정장치 및 그 방법
KR100839436B1 (ko) * 2006-10-25 2008-06-19 명지대학교 산학협력단 코사인 필터와 사인 필터의 이득차를 이용한 주파수 측정방법
CN102928665A (zh) * 2012-11-01 2013-02-13 南京国睿安泰信科技股份有限公司 一种中频数字化的频谱分析仪及其方法
KR101538738B1 (ko) * 2014-05-21 2015-07-27 와이피피 주식회사 삼단 이산 푸리에 변환을 이용한 주파수 측정 장치 및 방법
RU189666U1 (ru) * 2019-03-29 2019-05-30 Евгений Борисович Колесников Устройство для измерения частоты трехфазного синусоидального напряжения
CN114047683A (zh) * 2021-11-15 2022-02-15 星汉时空科技(长沙)有限公司 一种基于正交采样内插的时间间隔测量方法和装置

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100964114B1 (ko) 2008-04-25 2010-06-16 엘에스산전 주식회사 페이저 측정 장치
KR101681145B1 (ko) * 2015-08-27 2016-12-01 명지대학교 산학협력단 전력 계통의 비동기 신호에 대한 기본 주파수 추정 장치 및 방법

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR860000054A (ko) * 1984-06-26 1986-01-25 모리이 하지메 무연소배장치(無煙燒焙漿置)
US4983906A (en) * 1989-08-17 1991-01-08 Hewlett-Packard Company Frequency estimation system
KR100377248B1 (ko) 2001-04-11 2003-03-26 최한규 디지털 자기 상관 분광기
KR100403870B1 (ko) 2001-10-09 2003-11-03 명지대학교 사고 전류 또는 전압 신호의 기본파 성분 검출방법

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100839436B1 (ko) * 2006-10-25 2008-06-19 명지대학교 산학협력단 코사인 필터와 사인 필터의 이득차를 이용한 주파수 측정방법
KR100823832B1 (ko) * 2006-12-29 2008-04-21 엘에스산전 주식회사 전력계통의 주파수 추정장치 및 그 방법
CN102928665A (zh) * 2012-11-01 2013-02-13 南京国睿安泰信科技股份有限公司 一种中频数字化的频谱分析仪及其方法
KR101538738B1 (ko) * 2014-05-21 2015-07-27 와이피피 주식회사 삼단 이산 푸리에 변환을 이용한 주파수 측정 장치 및 방법
RU189666U1 (ru) * 2019-03-29 2019-05-30 Евгений Борисович Колесников Устройство для измерения частоты трехфазного синусоидального напряжения
CN114047683A (zh) * 2021-11-15 2022-02-15 星汉时空科技(长沙)有限公司 一种基于正交采样内插的时间间隔测量方法和装置
CN114047683B (zh) * 2021-11-15 2022-05-24 星汉时空科技(长沙)有限公司 一种基于正交采样内插的时间间隔测量方法和装置

Also Published As

Publication number Publication date
KR100758492B1 (ko) 2007-09-12

Similar Documents

Publication Publication Date Title
KR100758492B1 (ko) 직교필터를 이용한 주파수 추정 장치 및 그 방법
KR100823832B1 (ko) 전력계통의 주파수 추정장치 및 그 방법
Phadke et al. Synchronized phasor and frequency measurement under transient conditions
Sidhu et al. An iterative technique for fast and accurate measurement of power system frequency
US9658259B2 (en) System frequency measurement method, synchrophasor measurement method and device thereof
Kusljevic A simple recursive algorithm for frequency estimation
US8145443B2 (en) Fourier transform-based phasor estimation method and apparatus capable of eliminating influence of exponentially decaying DC offsets
CN103869162B (zh) 一种基于时域准同步的动态信号相量测量方法
Denys et al. Measurement of voltage phase for the French future defence plan against losses of synchronism
CN105044456A (zh) 一种基于正交子带的电网瞬时频率测量与跟踪方法
KR100186888B1 (ko) 전자유량계
Yamada High-accuracy estimations of frequency, amplitude, and phase with a modified DFT for asynchronous sampling
US20210389405A1 (en) Method, apparatus and device for detecting abnormity of energy metering chip, and medium
Macii et al. Design criteria of digital filters for synchrophasor estimation
CN104155521A (zh) 相位差的确定方法和装置
Kumar et al. A novel sub-cycle-based method for estimation of decaying DC component and fundamental phasor
Thilakarathne et al. Static performance comparison of prominent synchrophasor algorithms
ElRefaie et al. A novel technique to eliminate the effect of decaying DC component on DFT based phasor estimation
Kasztenny A new method for fast frequency measurement for protection applications
Belega et al. Amplitude estimation by a multipoint interpolated DFT approach
Petrović Frequency and parameter estimation of multi-sinusoidal signal
Ferrero et al. Employment of interpolated DFT-based PMU algorithms in three-phase systems
SE513954C2 (sv) Förfarande och system för bearbetning av signaler från en givare driven med en växelströmsexcitationssignal
CN111505375A (zh) 一种频率测量方法及装置
Kušljević A simultaneous estimation of frequency, magnitude, and active and reactive power by using decoupled modules

Legal Events

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

Payment date: 20120524

Year of fee payment: 6

FPAY Annual fee payment

Payment date: 20130515

Year of fee payment: 7

FPAY Annual fee payment

Payment date: 20141217

Year of fee payment: 8

FPAY Annual fee payment

Payment date: 20150907

Year of fee payment: 9

FPAY Annual fee payment

Payment date: 20161004

Year of fee payment: 10

FPAY Annual fee payment

Payment date: 20170906

Year of fee payment: 11

FPAY Annual fee payment

Payment date: 20180818

Year of fee payment: 12