24小时热门版块排行榜    

查看: 1192  |  回复: 1

wangshitong

银虫 (小有名气)

[求助] hankel逆变换

哪位用过Hankel逆变换,我从网上搜到一些matlab的代码,但还是搞不定。希望交流一下,谢谢啦!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangshitong

银虫 (小有名气)

说的再具体一下。
下面是我搜到的一个程序,这里面的参数应该怎么设置呢?有用过这个程序的大侠吗?请求指导,谢谢!

%[h,I]=ihat(H,k,r,n;I)
%---------------------
%
%Inverse Hankel transform of order n.
%
%Input:
% H      Spectrum K(k)
% k      Spatial frequencies [rad/m]   {pi/numel(H)*(0:numel(H)-1)}
% r      Radial positions [m]          {0:numel(H)-1}
% n      Transform order               {0}
%   or
% I      Integration kernel �         {default}
%
%Output:
% h      Signal h(r)
% I      Integration kernel
%
%
% �  If the integration kernel is missing, it is
%     recomputed from the Bessel functions (slow).
%

%     Marcel Leutenegger �June 2006
%
function [h,I]=ihat(H,k,r,n)
if sum(size(H) > 1) > 1
   error('Spectrum must be a vector.');
end
if nargin < 2 | isempty(k)
   k=pi/numel(H)*(0:numel(H)-1).';
else
   [k,w]=sort(k();
   H=H(w);
end
if nargin < 3 | isempty(r)
   r=0:numel(H)-1;
end
if nargin < 4 | isempty(n)
   n=0;
end
if numel(n) > 1
   if exist('w','var')
      I=n(w,;
   else
      I=n;
   end
else
   I=besselj(n,k*r(.');
end
h=reshape(frdr(H,k).'*I/(2*pi)^2,size(r));
2楼2013-10-10 11:44:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wangshitong 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 23jxep3nCNZb 2026-09-14 6/300 2026-09-18 11:30 by HjTo14L7GlnE
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 23jxep3nCNZb 2026-09-14 5/250 2026-09-18 11:28 by HjTo14L7GlnE
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 s3fFTmArrBt6 2026-09-14 6/300 2026-09-18 10:28 by yS4CkeWDhD6V
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 s3fFTmArrBt6 2026-09-14 7/350 2026-09-18 10:18 by yS4CkeWDhD6V
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +8 3n8v2C8RimXI 2026-09-13 9/450 2026-09-18 09:16 by yS4CkeWDhD6V
[找工作] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 vZfe6xYu34yj 2026-09-14 4/200 2026-09-18 02:28 by cITaGg3p5edr
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 Aj1rhIDL5ixY 2026-09-14 3/150 2026-09-18 01:40 by cITaGg3p5edr
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +5 s3fFTmArrBt6 2026-09-13 7/350 2026-09-17 23:28 by cuPZTDXS3VOd
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +8 LwdutQ8HoqWP 2026-09-13 8/400 2026-09-17 22:53 by cuPZTDXS3VOd
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 Aj1rhIDL5ixY 2026-09-14 4/200 2026-09-17 14:27 by ObVzqQQh1rrE
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 6F5UbRU2I5hL 2026-09-14 5/250 2026-09-17 13:17 by ObVzqQQh1rrE
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 vZfe6xYu34yj 2026-09-14 4/200 2026-09-17 05:50 by k06BKNblrGRH
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 vZfe6xYu34yj 2026-09-14 3/150 2026-09-17 05:38 by k06BKNblrGRH
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-17 03:16 by Tql5LQhh5rLK
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-16 17:00 by i3KBRD19ERkG
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 QUjhNVAcOSff 2026-09-13 6/300 2026-09-16 14:27 by i3KBRD19ERkG
[考博] 上海工程技术大学激光智能制造课题组招收博士研究生 +8 两三岁ss 2026-09-14 8/400 2026-09-16 09:46 by 水士口
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +7 3n8v2C8RimXI 2026-09-13 7/350 2026-09-16 02:35 by DgNGHc3h5tPl
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 08:26 by C79jjtjAKjEn
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 07:50 by C79jjtjAKjEn
信息提示
请填处理意见