求实验数据的spwvd分布，包括功率谱，二维以及三维spwvd分布，你指的拥有，完美运行

把你实验数据进行求spwvd分布，然后在进行你想要的运算，得到你想要的东西，完美运行，你指的拥有。

function [tfrtf] = tfrspwv（xtNghtrace）;
%TFRSPWV Smoothed Pseudo Wigner-Ville time-frequency distribution.
%	[TFRTF]=TFRSPWV（XTNGHTRACE） computes the Smoothed Pseudo 
%	Wigner-Ville distribution of a discrete-time signal X or the 
%	cross Smoothed Pseudo Wigner-Ville representation between two 
%	signals. 
% 
%	X     : signal if auto-SPWV or [X1X2] if cross-SPWV.
%	T     : time instant（s）          （default : 1:length（X））		.
%	N     : number of frequency bins （default : length（X））.
%	G     : time smoothing window G（0） being forced to 1. 
%	                                 （default : Hamming（N/10））. 
%	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 TFRSPWV runs TFRQVIEW.
%	F     : vector of normalized frequencies.
%
%	Example :
%	 sig=fmlin（1280.050.15）+fmlin（1280.30.4）;   
%	 g=tftb_window（15‘Kaiser‘）; h=tftb_window（63‘Kaiser‘）; 
%	 tfrspwv（sig1:12864gh1）;
% 
%	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 1 parameter required‘）;
end;
[xrowxcol] = size（x）;
if （xcol==0）|（xcol>2）
 error（‘X must have one or two columns‘）;
end

if （nargin <= 2）
 N=xrow;
elseif （N<0）
 error（‘N must be greater than zero‘）;
elseif （2^nextpow2（N）~=N）
 fprintf（‘For a faster computation N should be a power of two\n‘）;
end;

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

if （nargin == 1）
 t=1:xrow; g = tftb_window（glength）; h = tftb_window（hlength）; trace = 0;
elseif （nargin == 2）|（nargin == 3）
 g = tftb_window（glength）; h = tftb_window（hlength）; trace = 0;
elseif （nargin == 4）
 h = tftb_window（hlength）; trace = 0;
elseif （nargin == 5）
 trace = 0;
end;

[trowtcol] = size（t）;
if （trow~=1）
 error（‘T must only have one row‘）; 
end; 

[growgcol]=size（g）; Lg=（grow-1）/2; % g=g/sum（g）;
if （gcol~=1）|（rem（grow2）==0）
 error（‘G must be a smoothing window with odd l

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

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

     文件       4368  2005-10-13 23:14  平滑伪Wigner-Ville分布\平滑伪Wigner-Ville分布\tfrspwv.m

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

     文件        898  2018-11-28 15:58  平滑伪Wigner-Ville分布\平滑伪Wigner-Ville分布\Untitled.m

     目录          0  2018-11-12 21:12  平滑伪Wigner-Ville分布\平滑伪Wigner-Ville分布

     目录          0  2018-11-28 15:35  平滑伪Wigner-Ville分布

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

              1128914                    6