2013-09-06 5 views
0

웨이브 렛 함수의 이름을 cwtft 또는 icwtft으로 설정하면 'parm'이란 무엇을 의미합니까? wave = {wname,[7.6]}. 'morl' 기능을 cwtft 또는 icwtft 변환으로 사용하면 Fb 및 Fc도 바꿀 수 있습니까? 그리고 그렇지 않다면 어떻게 cwt transform으로 신호를 재구성 할 수 있습니까? 원인에 따라 fb 및 fc (cmorfb-fc)에 대해 선택적 값을 선택하게하십시오. Matlab은 역 웨이브 렛 변환을위한 직접 함수를 가지고 있지 않습니다.역 연속 웨이브 렛 변환 및 cwtft의 [Parm]

N = 1024; 
t = linspace(0,1,N); 
y = sin(2*pi*8*t).*(t<=0.5)+sin(2*pi*16*t).*(t>0.5); 
dt = 0.05;s0 = 2*dt;ds = 0.4875;NbSc = 20; 
wname = 'morl';sig = {y,dt};sca = {s0,ds,NbSc}; 
wave = {wname,[7.6]}; 
cwtsig = cwtft(sig,'scales',sca,'wavelet',wave); 
sigrec = icwtft(cwtsig,'signal',sig,'plot'); 

답변

0

누군가가 대답을 필요로한다고 생각했습니다. cwtfticwtft은 웨이브 렛 함수의 푸리에 변환을 사용하여 신호를 재구성합니다. wname의 'morl'은 분석적 기능입니다. 그래서 그것은 정확히 복잡한 도박이고 신호에 대한 위상과 크기 정보를 줄 것입니다. wave={‘morl’,[parm]}의 'parm'은 wo 또는 2*pi*fc입니다. 그래서 그것은 중심 주파수와 일치합니다. 'parm'의 기본값은 6이므로 fc=6/2*pi. 잔물결 웨이브 렛 함수는 psi(t,fc)=exp(j*2*pi*fc*t)*exp(-t^2/2이고 푸리에 변환은 psi^(k)=sqrt(2*pi)exp(-0.5(2*pi*k-ko)^2)입니다. ko= parm = 2*pi*fc. 그래서 parm을 변경하여 morlet wavelet을 설정할 수 있습니다. 모르겠어, 제 작법을 멋지게 만드는 법. 누군가 그것을 편집 plz.