| 查看: 1152 | 回复: 1 | ||
wangshitong银虫 (小有名气)
|
[求助]
hankel逆变换
|
| 哪位用过Hankel逆变换，我从网上搜到一些matlab的代码，但还是搞不定。希望交流一下，谢谢啦！ |
回复此楼» 猜你喜欢
有谁可曾问过你过的还好吗？
已经有22人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有5人回复
-
E0414, 我的本子有没有希望？
已经有7人回复
-
一篇论文同时出现在两个期刊，一模一样，这算不算学术不端，请各位老师斧正。
已经有12人回复
-
希望面上有个好结果
已经有7人回复
-
三区计算机方向期刊推荐
已经有5人回复
-
sci论文二审求助
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 关于matlab进行傅里叶变换和逆变换的程序，求大神指导！
已经有9人回复
- 零相位滤波器、傅里叶变换及逆变换滤波器
已经有4人回复
- 二阶偏微分方程Hankel变换后的结果
已经有3人回复
- 傅里叶逆变换 频域到时域
已经有3人回复
wangshitong
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1687.4
- 散金: 20
- 帖子: 214
- 在线: 219.1小时
- 虫号: 439970
- 注册: 2007-10-26
- 专业: 生物、医学光子学
|
说的再具体一下。 下面是我搜到的一个程序，这里面的参数应该怎么设置呢？有用过这个程序的大侠吗？请求指导，谢谢！ %[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