wavelets

DWT 필터 뱅크 시간 영역 웨이블릿

구문

psi = wavelets(fb)

[psi,t] = wavelets(fb)

설명

psi = wavelets(fb)는 DWT(이산 웨이블릿 변환) 필터 뱅크 fb의 웨이블릿 통과대역 필터에 대응하는 시간 영역/기준의 웨이블릿을 반환합니다.

[psi,t] = wavelets(fb)는 샘플링 시점 t를 반환합니다.

예제

DWT 필터 뱅크 웨이블릿

Daubechies db2 웨이블릿과 샘플링 주파수 1kHz를 사용하여 신호 길이가 1000개 샘플인 7 레벨 DWT 필터 뱅크를 만듭니다.

wv = "db4";
len = 1000;
lev = 7;
Fs = 1e3;
fb = dwtfilterbank('Wavelet',wv,'SignalLength',len, ...
    'Level',lev,'SamplingFrequency',Fs);

웨이블릿 대역통과 필터에 대응하는 시간 영역/기준의 웨이블릿을 플로팅합니다.

[psi,t] = wavelets(fb);
plot(t,psi')
grid on
title('Time-domain Wavelets')

가장 세밀한 스케일 시간 영역 웨이블릿과 대응하는 웨이블릿 대역통과 필터의 단측 크기 주파수 응답을 플로팅합니다.

sc = 1;
[psidft,f] = freqz(fb);
tiledlayout(2,1)
nexttile
plot(t,psi(sc,:))
grid on
xlabel('Time (sec)')
ylabel('Magnitude')
title(['Level ',num2str(sc),' Time-Domain Wavelet'])
nexttile
plot(f(len/2:end),abs(psidft(sc,len/2:end)))
grid on
xlabel('Hz')
ylabel('Magnitude')
title('Magnitude Frequency Response')

입력 인수

fb — 이산 웨이블릿 변환 필터 뱅크

dwtfilterbank 객체

DWT(이산 웨이블릿 변환) 필터 뱅크로, dwtfilterbank 객체로 지정됩니다.

출력 인수

psi — 시간 기준 웨이블릿

실수 값 행렬

웨이블릿 통과대역 필터에 대응하는 시간 기준 웨이블릿으로, L×N 행렬로 반환됩니다. 여기서 L은 필터 뱅크 Level이고 N은 SignalLength입니다. psi에서 웨이블릿은 가장 세밀한 스케일 분해능에서 가장 성긴 스케일 분해능의 순서로 정렬됩니다.

t — 샘플링 시점

실수 값 벡터

샘플링 시점으로, 길이가 N인 실수 값 벡터 t로 반환됩니다. 여기서 N은 필터 뱅크 SignalLength입니다. 샘플링 시점은 구간 $$\left[-\mathrm{½}NDT,\mathrm{½}NDT\right)$$ 내에 있습니다. 여기서 DT는 필터 뱅크 샘플링 주기(필터 뱅크 샘플링 주파수의 역수)입니다.

버전 내역

R2018a에 개발됨

참고 항목

dwtfilterbank | scalingfunctions | freqz