Matlab 时频分析（TFTB） 工具箱 源代码.m文件

此工具箱为 .m文件的Matlab源代码，直接拷贝toolbox根目录，添加路径即可使用，具体EMD,HHT下载安装步骤在我个人博客有介绍，祝好。 www.cnblogs.com/BrowserSnake

function x=altes（Nfminfmaxalpha） ;
%ALTES	Altes signal in time domain.
%	X=ALTES（NFMINFMAXALPHA） generates the Altes signal in
%	the time domain.
%
%	N     : number of points in time   
%	FMIN  : lower frequency bound （value of the hyperbolic
%	        instantaneous frequency law at the sample N） 
%	        in normalized frequency             （default : .05） 
%	FMAX  : upper frequency bound （value of the hyperbolic
%	        instantaneous frequency law at the first sample） 
%	        in normalized frequency             （default : 0.5）
%	ALPHA : attenuation factor of the envelope （default : 300）
%	X     : time row vector containing the Altes signal samples.
%
%	Example: 
%	 x=altes（1280.10.45）; plot（x）;
%
%	See also KLAUDER ANASING ANAPULSE ANASTEP DOPPLER.

%	P. Goncalves - September 1995.
%	Copyright （c） 1995 Rice University
%
%  This program is free software; you can redistribute it and/or modify
%  it under the terms of the GNU General Public License as published by
%  the Free Software Foundation; either version 2 of the License or
%  （at your option） any later version.
%
%  This program is distributed in the hope that it will be useful
%  but WITHOUT ANY WARRANTY; without even the implied warranty of
%  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
%  GNU General Public License for more details.
%
%  You should have received a copy of the GNU General Public License
%  along with this program; if not write to the Free Software
%  Foundation Inc. 51 Franklin St Fifth Floor Boston MA  02110-1301  USA

if （nargin == 0）
 error （ ‘The number of parameters must be at least 1.‘ ）;
elseif （nargin == 1）
 fmin=0.05; fmax=0.5; alpha=300;
elseif （nargin == 2）
 fmax=0.5; alpha=300;
elseif （nargin == 3）
 alpha=300 ;
end;

if （N <= 0）
 error （‘The signal length N must be strictly positive‘ ）;
elseif （fmin > 0.5） | （fmin < 0）
 error （ ‘FMIN must be in ]0  0.5]‘ ） ;
elseif （fmax > 0.5） | （fmax < 0）
 error （ ‘FMAX must be between 0 and 0.5‘ ） ;
elseif （alpha <= 1）
 error （ ‘ALPHA must be > 1‘ ） ;
else
 g = exp（（log（fmax/fmin））^2/（8*log（alpha））） ;
 nu0 = sqrt（fmin*fmax） ;
 beta=sqrt（2*log（g）*log（alpha））;
 t0 = N/（exp（beta）-exp（-beta）） ;
 t1 = t0*exp（-beta）; t2 = t0*exp（beta） ;
 b = -t0*nu0*g*log（g） ;
 t = linspace（t1t2N+1） ; t = t（1:N） ;
 x = （exp（-（log（t./t0）.^2）/（2*log（g））））.*cos（2*pi*b*log（t./t0）/log（g）） ;
 x = x.‘/norm（x） ;
end

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

     文件       2409  2005-10-13 16:14  TFTB\altes.m

     文件       2884  2005-10-13 16:14  TFTB\ambifunb.m

     文件       5455  2005-10-13 16:14  TFTB\ambifuwb.m

     文件       1620  2005-10-13 16:14  TFTB\amexpo1s.m

     文件       1628  2005-10-13 16:14  TFTB\amexpo2s.m

     文件       1654  2005-10-13 16:14  TFTB\amgauss.m

     文件       1601  2005-10-13 16:14  TFTB\amrect.m

     文件       1634  2005-10-13 16:14  TFTB\amtriang.m

     文件       1971  2005-10-13 16:14  TFTB\anaask.m

     文件       2060  2005-10-13 16:14  TFTB\anabpsk.m

     文件       2100  2005-10-13 16:14  TFTB\anafsk.m

     文件       1397  2005-10-13 16:14  TFTB\anapulse.m

     文件       2217  2005-10-13 16:14  TFTB\anaqpsk.m

     文件       1878  2005-10-13 16:14  TFTB\anasing.m

     文件       1390  2005-10-13 16:14  TFTB\anastep.m

     文件       5976  2005-10-13 16:14  TFTB\atoms.m

     文件       8126  2005-10-13 16:14  TFTB\Contents.m

     文件       3805  2005-10-13 16:14  TFTB\contwtgn.m

     文件       2688  2005-10-13 16:14  TFTB\contwtgnmir.m

     文件       2242  2005-10-13 16:14  TFTB\correlmx.m

     文件       5573  2008-02-15 15:38  TFTB\CVS\Entries

     文件         12  2008-02-15 15:38  TFTB\CVS\Repository

     文件         57  2008-02-15 15:38  TFTB\CVS\Root

     文件       1478  2005-10-13 16:14  TFTB\d2statio.m

     文件       1392  2005-10-13 16:14  TFTB\disprog.m

     文件       1134  2005-10-13 16:14  TFTB\divider.m

     文件       2430  2005-10-13 16:14  TFTB\dopnoise.m

     文件       2680  2005-10-13 16:14  TFTB\doppler.m

     文件       1315  2005-10-13 16:14  TFTB\dwindow.m

     文件       1878  2005-10-13 16:14  TFTB\fmconst.m

............此处省略111个文件信息