| 查看: 1129 | 回复: 1 | |||
wangshitong银虫 (小有名气)
|
[求助]
hankel逆变换
|
| 哪位用过Hankel逆变换，我从网上搜到一些matlab的代码，但还是搞不定。希望交流一下，谢谢啦！ |
回复此楼» 猜你喜欢
青椒八年已不青，大家都被折磨成啥样了？
已经有7人回复
-
为什么nbs上溴 没有产物点出现呢
已经有10人回复
-
救命帖
已经有11人回复
-
招博士
已经有5人回复
-
青年基金C终止
已经有3人回复
-
26申博求博导推荐-遥感图像处理方向
已经有4人回复
-
限项规定
已经有7人回复
-
西南交通大学国家级人才团队2026年博士研究生招生（考核制）—机械、材料、力学方向
已经有3人回复
-
英文综述是否需要润色及查重
已经有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