伪Wigner-Ville分布，可以运行。

这是把数据进行伪Wigner-Ville分布分布。在前面你可以先用分解方法分解，在求伪Wigner-Ville分布分布，可运行。

function [tfrtf] = tfrpwv（xtNhtrace）;
%TFRPWV	Pseudo Wigner-Ville time-frequency distribution.
%	[TFRTF]=TFRPWV（XTNHTRACE） computes the Pseudo Wigner-Ville  
%	distribution of a discrete-time signal X or the
%	cross Pseudo Wigner-Ville representation between two signals. 
% 
%	X     : signal if auto-PWV or [X1X2] if cross-PWV.
%	T     : time instant（s）          （default : 1:length（X））.
%	N     : number of frequency bins （default : length（X））.
%	H     : frequency smoothing window in the time-domain
%		H（0） being forced to 1   （default : Hamming（N/4））. 
%	TRACE : if nonzero the progression of the algorithm is shown
%	                                 （default : 0）.
%	TFR   : time-frequency representation. When called without 
%	        output arguments TFRPWV runs TFRQVIEW.
%	F     : vector of normalized frequencies.
%
%	Example:
%	 sig=fmlin（1280.10.4）; tfrpwv（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

[xrowxcol] = size（x）;
if （nargin < 1）
 error（‘At least 1 parameter is required‘）;
elseif （nargin <= 2）
 N=xrow;
end;

hlength=floor（N/4）;
hlength=hlength+1-rem（hlength2）;

if （nargin == 1）
 t=1:xrow; h = tftb_window（hlength）; trace=0;
elseif （nargin == 2）|（nargin == 3）
 h = tftb_window（hlength）; trace=0;
elseif （nargin == 4）
 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 & nargin==5）
 fprintf（‘For a faster computation N should be a power of two\n‘）;
end; 

[hrowhcol]=size（h）; Lh=（hrow-1）/2; h=h/h（Lh+1）;
if （hcol~=1）|（rem（hrow2）==0）
 error（‘H must be a smoothing window with odd length‘）;
end;

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

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

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

     文件       3285  2005-10-13 23:14  伪Wigner-Ville分布\tfrpwv.m

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

     文件        709  2018-11-12 21:16  伪Wigner-Ville分布\Untitled.m

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

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

              1127642                    5