HT_decomposition
所属分类:matlab编程
开发工具:matlab
文件大小:132KB
下载次数:37
上传日期:2007-05-26 23:50:53
上 传 者:
nwei
说明: hibert变换程序,力学领域应用hht必用的程序。
(hibert transform procedures, mechanical application area hht used in the procedure.)
文件列表:
decomp.p (8629, 2006-03-09)
demo (0, 2006-03-10)
demo\ammod.mat (11888, 2005-12-07)
demo\duf.mat (38956, 2005-12-07)
demo\hf.mat (2763, 2005-12-07)
demo\run.fig (11293, 2005-12-15)
demo\run.p (35754, 2006-03-09)
demo\square2.mat (17086, 2005-12-07)
demo\two.mat (16206, 2005-12-07)
hf.mat (2763, 2005-12-07)
pl.m (1768, 2005-12-14)
tilefigs.m (2303, 2004-03-02)
file "README.TXT"
The package includes the main Matlab program DECOMP for non-stationary
and nonlinear vibration signal decomposition together with the program
PL.M for plotting the obtained decomposed results. It also includes some
simulation examples and the demo program RUN (Matlab 7.x).
The Hilbert Vibration Decomposition method is based on the Hilbert transform
signal processing [*]. Run the program RUN and choose a vibration signal.
Press the button 'RUN' to start the decomposition procedure.
[*] M. Feldman. TIME-VARYING VIBRATION DECOMPOSITION AND ANALYSIS
BASED ON THE HILBERT TRANSFORM, Journal of Sound and Vibration.
2006, Vol 295/3-5 pp 518-530.
*****************************************************
Function DECOMP (Hilbert Vibration Decomposition)
The function DECOMP is a signal processing M-File (Matlab 7.x).
It extracts intrinsic quasi-harmonic components Y(t) from the initial
signal composition y(t).
[Y,A,F,dev]=decomp(y,2);
pl(Y,A,F,Fs)
Inputs:
y – the initial vibration signal in time domain, 2 – the
specified number of components
Outputs: Y – the matrix of row component (vectors), A – the matrix of
the component envelopes, F – the matrix of the relative component
instantaneous frequencies, dev – the vector of each component standard
deviation.
Conversion from the relative frequency F to the absolute frequency Fa
[Hz] is made with the formula: Fa=Fs*F where Fs – is the sampling
frequency [Hz].
Function PL (Plot)
The function PL will plot the results of the decomposition in four
subplots: The component waveforms; The sum of the component waveforms;
The component instantaneous frequencies and the envelopes; The 3D plot
of the instantaneous frequency and the envelope of each component.
pl(Y,A,F,Fs);
Inputs: Y – the matrix of row component (vectors), A – the matrix of the
component envelopes, F – the matrix of the relative component
instantaneous frequencies, Fs – the sampling frequency [Hz],
EXAMPLE:
[Y,A,F,dev]=decomp(y,2); pl(Y,F,A,1000)
******************************************************
LIMITATIONS: The sampling frequency Fs has to be in the range
Fs=(20-80)*f0. The minimum of the initial points in time domain is 721.
2006 - Michael Feldman, PhD
Technion, Haifa, Israel
http://hitech.technion.ac.il/feldman
bug-reports email to MFeldman@technion.ac.il
近期下载者:
相关文件:
收藏者: