Wigner-Ville分布分布
把你的实验数据变成Wigner-Ville分布分布,可以运行,你指的拥有。
代码片段和文件信息
属性 大小 日期 时间 名称
----------- --------- ---------- ----- ----
文件 1118292 2018-11-12 20:45 Wigner-Ville分布1224_2_2磨损击穿C-OUT.xlsx
文件 2784 2005-10-13 23:14 Wigner-Ville分布 frwv.m
文件 5356 2005-10-13 23:14 Wigner-Ville分布 ftb_window.m
文件 708 2018-11-12 21:19 Wigner-Ville分布Untitled.m
目录 0 2018-11-12 21:19 Wigner-Ville分布
----------- --------- ---------- ----- ----
1127140 5
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
‘);
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分布 frwv.m
文件 5356 2005-10-13 23:14 Wigner-Ville分布 ftb_window.m
文件 708 2018-11-12 21:19 Wigner-Ville分布Unti
目录 0 2018-11-12 21:19 Wigner-Ville分布
----------- --------- ---------- ----- ----
1127140 5
版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌抄袭侵权/违法违规的内容, 请发送邮件举报,一经查实,本站将立刻删除。
评论列表(条)