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')

Figure contains an axes object. The axes object with title Time-domain Wavelets contains 7 objects of type line.

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

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')

Figure contains 2 axes objects. Axes object 1 with title Level 1 Time-Domain Wavelet, xlabel Time (sec), ylabel Magnitude contains an object of type line. Axes object 2 with title Magnitude Frequency Response, xlabel Hz, ylabel Magnitude contains an object of type line.

입력 인수

모두 축소

`fb` — 이산 웨이블릿 변환 필터 뱅크
`dwtfilterbank` 객체

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

출력 인수

모두 축소

`psi` — 시간 기준 웨이블릿
실수 값 행렬

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

`t` — 샘플링 시점
실수 값 벡터

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