Wigner-Ville分布分布

把你的实验数据变成Wigner-Ville分布分布，可以运行，你指的拥有。

function [tfrtf] = tfrwv（xtNtrace）;
%TFRWV	Wigner-Ville time-frequency distribution.
%	[TFRTF]=TFRWV（XTNTRACE） computes the Wigner-Ville distribution
%	of a discrete-time signal X 
%	or the cross Wigner-Ville representation between two signals. 
% 
%	X     : signal if auto-WV or [X1X2] if cross-WV.
%	T     : time instant（s）          （default : 1:length（X））.
%	N     : number of frequency bins （default : length（X））.
%	TRACE : if nonzero the progression of the algorithm is shown
%	                                 （default : 0）.
%	TFR   : time-frequency representation. When called without 
%	        output arguments TFRWV runs TFRQVIEW.
%	F     : vector of normalized frequencies.
%
%	Example :
%	 sig=fmlin（1280.10.4）; tfrwv（sig）;
% 
%	See also all the time-frequency representations listed in
%	the file CONTENTS （TFR*）

%	F. Auger May-August 1994 July 1995.
%	Copyright （c） 1996 by CNRS （France）.
%
%  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（‘At least one parameter required‘）;
end;
[xrowxcol] = size（x）;

if （nargin == 1）
 t=1:xrow; N=xrow ; trace=0;
elseif （nargin == 2）
 N=xrow ; trace=0;
elseif （nargin == 3）
 trace = 0;
end;

if （N<0）
 error（‘N must be greater than zero‘）;
end;

[trowtcol] = size（t）;
if （xcol==0）|（xcol>2）
 error（‘X must have one or two columns‘）;
elseif （trow~=1）
 error（‘T must only have one row‘）; 
elseif （2^nextpow2（N）~=N）
 fprintf（‘For a faster computation N should be a power of two\n‘）;
end; 

tfr= zeros （Ntcol）;  
if trace disp（‘Wigner-Ville distribution‘）; end;
for icol=1:tcol
 ti= t（icol）; taumax=min（[ti-1xrow-tiround（N/2）-1]）;
 tau=-taumax:taumax; indices= rem（N+tauN）+1;
 tfr（indicesicol） = x（ti+tau1） .* conj（x（ti-tauxcol））;
 tau=round（N/2）; 
 if （ti<=xrow-tau）&（ti>=tau+1）
  tfr（tau+1icol） = 0.5 * （x（ti+tau1） * conj（x（ti-tauxcol））  + ...
                           x（ti-tau1） * conj（x（ti+tauxcol））） ;
 end;
 if trace disprog（icoltcol10）; end;
end; 
tfr= fft（tfr）; 
if （xcol==1） tfr=real（tfr）; end ;

if （nargout==0）
 tfrqview（tfrxt‘tfrwv‘）;
elseif （nargout==3）
 f=（0.5*（0:N-1）/N）‘;
end;

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

     文件    1118292  2018-11-12 20:45  Wigner-Ville分布\1224_2_2磨损击穿C-OUT.xlsx

     文件       2784  2005-10-13 23:14  Wigner-Ville分布\tfrwv.m

     文件       5356  2005-10-13 23:14  Wigner-Ville分布\tftb_window.m

     文件        708  2018-11-12 21:19  Wigner-Ville分布\Untitled.m

     目录          0  2018-11-12 21:19  Wigner-Ville分布

----------- ---------  ---------- -----  ----

              1127140                    5