频率切片小波变换MALTAB代码

频率切片小波变换程序（FTWT），matlab的m文件。从一种新的角度出发，通过自由选择频率切片函数、引进新尺度参数，在频率域实现小波变换，该变换能够很好地刻画信号各成分之间的相对能量关系。此外，频率切片小波变换的时频窗中心频率就是观测频率，而无需进行尺度换算。并在适当地方做好了备注。

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

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件      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