求助:求用MATLAB实现hilbert变换过程的程序谁有共享一下哈谢谢

nihao00000 · 发表于 2006-6-4 22:24

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

如题 求用MATLAB实现hilbert变换过程的程序谁有共享一下哈谢谢

bainhome · 发表于 2006-6-5 01:49

function x = hilbert(xr,n) %HILBERT Discrete-time analytic signal via Hilbert transform. % X = HILBERT(Xr) computes the so-called discrete-time analytic signal % X = Xr + i*Xi such that Xi is the Hilbert transform of real vector Xr. % If the input Xr is complex, then only the real part is used: Xr=real(Xr). % If Xr is a matrix, then HILBERT operates along the columns of Xr. % % HILBERT(Xr,N) computes the N-point Hilbert transform. Xr is padded with % zeros if it has less than N points, and truncated if it has more. % % For a discrete-time analytic signal X, the last half of fft(X) is zero, % and the first (DC) and center (Nyquist) elements of fft(X) are purely real. % % Example: % Xr = [1 2 3 4]; % X = hilbert(Xr) % produces X=[1+1i 2-1i 3-1i 4+1i] such that Xi=imag(X)=[1 -1 -1 1] is the % Hilbert transform of Xr, and Xr=real(X)=[1 2 3 4]. Note that the last half % of fft(X)=[10 -4+4i -2 0] is zero (in this example, the last half is just % the last element). Also note that the DC and Nyquist elements of fft(X) % (10 and -2) are purely real. % % See also FFT, IFFT.
% Copyright 1988-2002 The MathWorks, Inc. % $Revision: 1.10 $ $Date: 2002/03/28 17:28:20 $
% References: % [1] Alan V. Oppenheim and Ronald W. Schafer, Discrete-Time % Signal Processing, 2nd ed., Prentice-Hall, Upper Saddle River, % New Jersey, 1998. % % [2] S. Lawrence Marple, Jr., Computing the discrete-time analytic % signal via FFT, IEEE Transactions on Signal Processing, Vol. 47, % No. 9, September 1999, pp.2600--2603.
if nargin<2, n=[]; end if ~isreal(xr) warning('HILBERT ignores imaginary part of input.') xr = real(xr); end % Work along the first nonsingleton dimension [xr,nshifts] = shiftdim(xr); if isempty(n) n = size(xr,1); end x = fft(xr,n,1); % n-point FFT over columns. h = zeros(n,~isempty(x)); % nx1 for nonempty. 0x0 for empty. if n>0 & 2*fix(n/2)==n % even and nonempty h([1 n/2+1]) = 1; h(2:n/2) = 2; elseif n>0 % odd and nonempty h(1) = 1; h(2:(n+1)/2) = 2; end x = ifft(x.*h(:,ones(1,size(x,2))));
% Convert back to the original shape. x = shiftdim(x,-nshifts);

nihao00000 · 发表于 2006-6-5 08:17

谢谢楼上的把要进行hilbert变换的信号代入xr就行吗?我怎么运行不出来

bainhome · 发表于 2006-6-5 12:45

实际这是MATLAB信号工具箱中的官方函数，可以查阅看看它的实例以及官方的在线FAQ，有详细的介绍 <a href="http://www.mathworks.com/cgi-bin/texis/webinator/search/?db=MSS&prox=page&rorder=750&rprox=750&rdfreq=500&rwfreq=500&rlead=250&sufs=0&order=r&is_summary_on=1&ResultCount=10&query=hilbert" target="_blank" >http://www.mathworks.com/cgi-bin/texis/webinator/search/?db=MSS&prox=page&rorder=750&rprox=750&rdfreq=500&rwfreq=500&rlead=250&sufs=0&order=r&is_summary_on=1&ResultCount=10&query=hilbert</A>

账号		自动登录	找回密码
密码			我要加入

求助:求用MATLAB实现hilbert变换过程的程序 谁有共享一下哈 谢谢

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

求助:求用MATLAB实现hilbert变换过程的程序谁有共享一下哈谢谢