이 번역 페이지는 최신 내용을 담고 있지 않습니다. 최신 내용을 영문으로 보려면 여기를 클릭하십시오.
hilbert
R2026a힐베르트 변환을 사용한 이산시간 해석적 신호
설명
예제
시퀀스를 정의하고 hilbert를 사용하여 이 시퀀스의 해석적 신호를 계산합니다.
xr = [1 2 3 4]'; x = hilbert(xr)
x = 4×1 complex
1.0000 + 1.0000i
2.0000 - 1.0000i
3.0000 - 1.0000i
4.0000 + 1.0000i
x의 허수부는 xr의 힐베르트 변환이고, 실수부는 xr 그 자체입니다.
imx = imag(x)
imx = 4×1
1
-1
-1
1
rex = real(x)
rex = 4×1
1
2
3
4
x의 이산 푸리에 변환(DFT)의 마지막 절반은 0입니다. (이 예제에서 변환의 마지막 절반은 마지막 요소입니다.) fft(x)의 DC 요소와 나이퀴스트 요소는 순수 실수입니다.
dft = fft(x)
dft = 4×1 complex
10.0000 + 0.0000i
-4.0000 + 4.0000i
-2.0000 + 0.0000i
0.0000 + 0.0000i
hilbert 함수는 유한한 데이터 블록에 대해 정확한 해석적 신호를 구합니다. 또한 유한 임펄스 응답(FIR) 힐베르트 변환기 필터를 사용하여 허수부에 대한 근삿값을 계산하는 방법으로 해석적 신호를 생성할 수도 있습니다.
주파수가 203Hz, 721Hz, 1001Hz인 3개의 정현파로 구성된 시퀀스를 생성합니다. 이 시퀀스는 약 1초 동안 10kHz로 샘플링됩니다. hilbert 함수를 사용하여 해석적 신호를 계산합니다. 0.01초와 0.03초 사이의 해석적 신호를 플로팅합니다.
fs = 1e4; t = 0:1/fs:1; x = 2.5 + cos(2*pi*203*t) + sin(2*pi*721*t) + cos(2*pi*1001*t); y = hilbert(x); plot(t,real(y),t,imag(y)) xlim([0.01 0.03]) legend("real","imaginary") title("Hilbert Function") xlabel("Time (s)")
원래 시퀀스와 해석적 신호의 파워 스펙트럼 밀도에 대한 Welch 추정값을 계산합니다. 시퀀스를 길이가 256이며 해밍 윈도우가 적용된 겹치지 않는 조각으로 나눕니다. 음수 주파수에서는 해석적 신호가 전력을 가지지 않는 것을 확인합니다.
pwelch([x;y].',256,0,[],fs,"centered") legend("Original","Hilbert")
designfilt 함수를 사용하여 60차 힐베르트 변환기 FIR 필터를 설계합니다. 천이 폭을 400Hz로 지정합니다. 필터의 주파수 응답을 시각화합니다.
fo = 60; d = designfilt("hilbertfir",FilterOrder=fo, ... TransitionWidth=400,SampleRate=fs); freqz(d,1024,fs)
정현파 시퀀스에 필터를 적용하여 해석적 신호의 허수부에 대한 근삿값을 구합니다.
hb = filter(d,x);
필터의 군지연 grd는 필터 차수의 절반과 동일합니다. 이 지연을 보정합니다. 허수부의 처음 grd개 샘플과 실수부와 시간 벡터의 마지막 grd개 샘플을 제거합니다. 0.01초와 0.03초 사이의 결과를 플로팅합니다.
grd = fo/2; y2 = x(1:end-grd) + 1j*hb(grd+1:end); t2 = t(1:end-grd); plot(t2,real(y2),t2,imag(y2)) xlim([0.01 0.03]) legend("real","imaginary") title("FIR Filter") xlabel("Time (s)")
해석적 신호 근삿값의 파워 스펙트럼 밀도(PSD)를 추정하고 이를 hilbert 결과와 비교합니다.
pwelch([y;[y2 zeros(1,grd)]].',256,0,[],fs,"centered") legend("Hilbert","FIR Filter")
입력 인수
출력 인수
세부 정보
알고리즘
시퀀스 xr의 해석적 신호는 단측 푸리에 변환을 가집니다. 즉, 음수 주파수에 대해 이 변환은 0이 됩니다. 해석적 신호의 근삿값을 구하기 위해 hilbert는 입력 시퀀스의 FFT를 계산하고, 음수 주파수에 대응되는 이러한 FFT 계수를 0으로 바꾸고, 그 결과의 역 FFT를 계산합니다.
hilbert는 다음 4단계의 알고리즘을 사용합니다.
입력 시퀀스의 FFT를 계산하고 벡터
x에 결과를 저장합니다.요소
h(i)가 다음과 같은 값을 가지는 벡터h를 생성합니다.i= 1,(n/2)+1의 경우 1i= 2, 3, … ,(n/2)의 경우 2i=(n/2)+2, … ,n의 경우 0
x와h의 요소별 곱을 계산합니다.3단계에서 구한 시퀀스의 역 FFT를 계산하고 결과의 처음
n개 요소를 반환합니다.
이 알고리즘은 [2]에서 처음으로 소개되었습니다. 이 기법은 입력 신호 x가 유한한 데이터 블록인 것으로 가정합니다. 이러한 가정을 통해 함수가 x에서 중복된 스펙트럼 요소를 정확하게 제거할 수 있습니다. FIR 필터 방식을 기반으로 한 방법은 해석적 신호의 근삿값만 구할 수 있지만, 데이터에 대해 지속적으로 동작하는 강점이 있습니다. FIR 필터를 사용하여 계산된 힐베르트 변환에 대한 또 다른 예제는 Single-Sideband Amplitude Modulation 항목을 참조하십시오.
참고 문헌
[1] Claerbout, Jon F. Fundamentals of Geophysical Data Processing with Applications to Petroleum Prospecting. Oxford, UK: Blackwell, 1985.
[2] Marple, S. L. “Computing the Discrete-Time Analytic Signal via FFT.” IEEE® Transactions on Signal Processing. Vol. 47, 1999, pp. 2600–2603.
[3] Oppenheim, Alan V., Ronald W. Schafer, and John R. Buck. Discrete-Time Signal Processing. 2nd Ed. Upper Saddle River, NJ: Prentice Hall, 1999.
