資源簡(jiǎn)介
頻率切片小波變換程序(FTWT),matlab的m文件。從一種新的角度出發(fā),通過自由選擇頻率切片函數(shù)、引進(jìn)新尺度參數(shù),在頻率域?qū)崿F(xiàn)小波變換,該變換能夠很好地刻畫信號(hào)各成分之間的相對(duì)能量關(guān)系。此外,頻率切片小波變換的時(shí)頻窗中心頻率就是觀測(cè)頻率,而無需進(jìn)行尺度換算。并在適當(dāng)?shù)胤阶龊昧藗渥ⅰ?img width=510 src=http://www.oy-car.com/pic/58172.jpg />
代碼片段和文件信息
function?[a]?=?GetFSWT(sFsfpkapaTn)
%%?Tn:?Time?Resample?Point?Number
%%?fp?:?Observed?frequency?range
N=length(s);??%%%?Original?Point?Number
Y?=?fft(sN);?%%?FFT?
Y(1)=0;?%%?cut?the?constant?part
nl=length(fp);?%%?Observation?frequency?length
a=zeros(Tnnl);
s1=zeros(Tn1);
TNN=fix(Tn/2);
k0=0;
for??p=1:nl;
???s1=s1*0;??
???k0=0;
??for?k=1:Tn;
?????s1(k)=0;
?????if?fp(p)==0?continue;end;
?????m=fp(p)+k-1-TNN;
????if?m>=0?&&??abs(k-1-TNN)<.5*TNN?&&?m ?????ct=kapa*(k-1-TNN)/fp(p);
??????s1(k)=Y(m+1)/(1+ct*ct);
?????%s1(k)=Y(m+1)*exp(-ct*ct/2);??%%%?Here?you?can?use?your?Frequency?slice?function.
?????end;
??end
???s1(1)=0;
???s1=conj(s1);
???Y1=fft(s1Tn);
???for?k=1:Tn;?if?mod(k-12)==1?Y1(k)=-Y1(k);?end;??end??
???a(:p)=Y1;
end
?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????文件??????56745??2010-04-30?15:13??FSWT\FSWT-Example.jpg
?????文件????????784??2009-08-31?09:41??FSWT\GetFSWT.asv
?????文件????????792??2010-04-30?12:55??FSWT\GetFSWT.m
?????文件????????451??2010-04-30?15:11??FSWT\GetInvFSWT.asv
?????文件????????557??2010-04-30?14:44??FSWT\GetInvFSWT.m
?????文件??????33390??2010-04-30?14:56??FSWT\TestFSWT.asv
?????文件??????33413??2010-04-30?14:57??FSWT\TestFSWT.m
?????文件???????4335??2013-02-16?00:57??FSWT\wtspecdecom.m
?????目錄??????????0??2012-08-15?13:40??FSWT
-----------?---------??----------?-----??----
???????????????130467????????????????????9
評(píng)論
共有 條評(píng)論