besselh.m

可用于对信号的频谱进行分析,希望站长能整理好,以便大家互相学习

% BESSELH   戞3庬偺Bessel娭悢(Hankel娭悢)
% 
% K = 1傑偨偼2偵懳偟偰丄H = BESSELH(NU,K,Z) 偼丄暋慺攝楍 Z 偺奺梫慺偵
% 懳偟偰丄Hankel 娭悢 H1_nu(Z) 傑偨偼 H2_nu(Z) 傪寁嶼偟傑偡丅
%
% H = BESSELH(NU,Z) 偼丄K = 1 傪巊偄傑偡丅
% H = BESSELH(NU,1,Z,1) 偼丄exp(-i*z) 偱 H1_nu(z) 傪僗働乕儕儞僌偟傑偡丅
% H = BESSELH(NU,2,Z,1) 偼丄exp(+i*z) 偱H2_nu(z) 傪僗働乕儕儞僌偟傑偡丅
%
% NU 偲 Z 偑摨偠戝偒偝偺攝楍偺応崌偼丄寢壥偺攝楍傕摨偠戝偒偝偵側傝傑偡丅
% 偳偪傜偐偺擖椡偑僗僇儔偺応崌偼丄傕偆堦曽偺擖椡偺僒僀僘偲摨偠戝偒偝偵
% 奼挘偝傟傑偡丅1偮偺擖椡偑峴儀僋僩儖偱丄懠偺擖椡偑楍儀僋僩儖偺応崌丄
% 寢壥偼丄娭悢抣偐傜側傞2師尦偺僥乕僽儖偵側傝傑偡丅
%
% [H,IERR] = BESSELH(NU,K,Z) 偼丄僄儔乕僼儔僌偺攝楍傕弌椡偟傑偡丅
%     ierr = 1   堷悢偑娫堘偭偰偄傑偡丅
%     ierr = 2   僆乕僶僼儘乕偱 Inf傪弌椡偟傑偡丅
%     ierr = 3   堷悢傪尭傜偡偙偲偵傛傝丄惛搙偑掅壓偟傑偡丅
%     ierr = 4   z傑偨偼NU偑戝偒夁偓傞偙偲偵傛傝丄惛搙偑掅壓偟傑偡丅
%     ierr = 5   廂懇偟傑偣傫丅NaN傪弌椡偟傑偡丅
%
% Hankel娭悢偲Bessel娭悢偺娭學偼丄偮偓偺傛偆偵側傝傑偡丅
%
%     besselh(nu,1,z) = besselj(nu,z) + i*bessely(nu,z)
%     besselh(nu,2,z) = besselj(nu,z) - i*bessely(nu,z)
%
% 椺戣:
% 偙偺椺偼丄Abramowitz 偲 Stegun挊"Handbook of Mathematical Functions"
% 偺359儁乕僕偵婰偝傟偰偄傞丄Hankel娭悢H1_0(z)偺儌僕儏儔僗偲埵憡偺僐儞
% 僞乕僾儘僢僩傪嶌惉偟傑偡丅
%
%     [X,Y] = meshgrid(-4:0.025:2,-1.5:0.025:1.5);
%     H = besselh(0,1,X+i*Y);
%     contour(X,Y,abs(H),0:0.2:3.2)丄hold on
%     contour(X,Y,(180/pi)*angle(H),-180:10:180); hold off
%
% 偙偺M-僼傽僀儖偼丄D. E. Amos偵傛傞Fortran儔僀僽儔儕偺MEX僀儞僞僼僃乕僗
% 傪巊梡偟傑偡丅
%
% 嶲峫  BESSELJ, BESSELY, BESSELI, BESSELK.


%   Copyright 1984-2004 The MathWorks, Inc.