synchdem matlab,希尔伯特振动分解(HVD)的matlab程序

CODE:

function [Y,A,om_r,dev]=hvd(x,n,fp); % % x - initial signal, n - number of decomposed components % Y - decomposed components, A - component envelopes , % F_r - component relative angular frequency % F=Fs*om_r/2/pi - Absolute frequecy [Hz], Fs -sampling frequency, % dev=std(Y_i)/std(Y_1)) - relative standard deviation of the decomposed component % % Example:  [Y,A,om_r,dev]=hvd(x,2,0.02); % % LIMITATIONS: %    The sampling frequency Fs has to be in the range Fs=(20-80)*f0. %    The minimum of points in time domain is 230*3+1 = 691 % % © 2011 Michael Feldman % For use with the book "HILBERT TRANSFORM APPLICATION % IN MECHANICAL VIBRATION", John Wiley & Sons, 2011 %  if n>7; disp('Max number of components not greater than 7'); end if n<=0; disp('Number of components less than 1');Y=[];A=[];F_r=[];dev=[];return;end x=x(:); s(1)=std(x); if s(1)==0,Y=[];A=[];F_r=[];dev=[];disp('Zero signal');return,end; for k=1:n;     [At,Ft,phit]=inst(x,1);     omf=2*pi*lpf(Ft,fp);% Angular Frequency lowpass filtering (Smoothing)     [yi,Ai,phi]=synchdem(x,omf,fp);     Y(:,k)=yi; A(:,k)=Ai;     om_r(:,k)=omf;            % Angular Frequency, [Radians]         x=x-yi;         s(k)=std(x)/s(1);         if k == 7, dev=[1 diff(s)]; return,  end end dev=s;  % Relative standard deviation of the components  return

hvd相关程序集合congr.m – Congruent envelope and instantaneous frequency coph.m – Congruent phase estimation diffir.m – Remez differentiator with filtering procedure forcevib.m – Identification of modal parameters of vibration SDOF system under forced excitation freevib.m – Identification of modal parameters of SDOF system free vibration hilbfft.m – Hilbert transform (based on FFT) and filtering procedure hilbfir.m – Remez Hilbert transformer and filtering procedure hilbturner.m – Hilbert transform via Turner filter and filtering procedure hvd.m – Hilbert Vibration Decomposition ilpf.m – An ideal lowpass filter and filtering inst.m – Instantaneous characteristics (envelope and frequency) of the signal integ.m – Signal integration procedure for vibration lpf.m – Remez lowpass filter and filtering procedure phaseh.m – Instantaneous phase between two signals pl.m – Hilbert spectrum presentation of decomposed components plfreq.m – Hilbert spectrum presentation of decomposed components with frequency arranging plfor.m – Presentation of the Hilbert transform forced vibration identification plfree.m – Presentation of the Hilbert transform free vibration identification print1.m – Print Matlab figure in Portable Network Graphic 24-bit truecolor image synchdem.m – Synchronous demodulation of a component with the known frequency from composition tilefigs.m – Automatically tiles Matlab figures (by Charles Plum) duffod.mat – Duffing equation forced excitation and vibration solution duffrd.mat – Duffing equation free vibration solution hilbTurner.mat – the coefficients for the Hilbert quadrature filters (by Clay S. Turner)
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值