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