Wavelets : 주파수에서 비늘을 결정하는 방법?

Question

다음 코드 샘플에서는 인쇄 된 출력 행이 적어도 대략 서로 일치해야합니다. 그러나 그들은 그렇지 않습니다.

% input data 
y = rand(10000,1); 
dt = 1/500; % 500Hz is the sampling rate 


% wavelet transform 
f = 1:50;     % frequencies at which to perform WT 
s = scal2frq(f,'morl',dt); % corresponding scales estimates 
wt = cwtft(y, 'wavelet', 'morl', 'scales', s); 

[ f; wt.frequencies/dt ] 

내가 매트랩 변환 연속 웨이블릿을 사용하려고, 나는 내 입력 신호에 다른 주파수에서 변환을 수행 할,하지만 해당 을 추정하는 신뢰할 수있는 방법을 찾기 위해 고군분투 가늠자. 이 코드 샘플에서

, 나는 그 규모를 추정 scal2freq을 사용하고,하지만 내가 cwtft에 의해 반환 된 주파수 추정치가 될 것이라고 기대 (가 거듭 제곱 인 것을 알고, scal2frq(scal2frq(x,'morl',dt),'morl',dt) == x 수치의 정밀도까지의 모든 x, dt에 대한 그 즉) 비슷한. 보시다시피 그들은 그렇지 않습니다.

누구나 cwtft이 이러한 빈도를 예측하는 방법을 알고 있으며, 필요한 저울을 예측하기 위해 역방향 추정을 사용할 수 있는지 여부를 아는 사람이 있습니까?

Answer 1

Morlet 웨이브 렛의 빈도와 스케일 간의 일치는 scale = 6/(2*pi * frequency*dt)입니다 (여기서 frequency*dt은 사이클/단위입니다).