| 查看: 1121 | 回复: 1 | ||
wangshitong银虫 (小有名气)
|
[求助]
hankel逆变换
|
| 哪位用过Hankel逆变换，我从网上搜到一些matlab的代码，但还是搞不定。希望交流一下，谢谢啦！ |
回复此楼» 猜你喜欢
不自信的我
已经有11人回复
-
北核录用
已经有3人回复
-
要不要辞职读博？
已经有6人回复
-
实验室接单子
已经有3人回复
-
磺酰氟产物，毕不了业了！
已经有8人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有10人回复
-
26申博（荧光探针方向，有机合成）
已经有4人回复
-
论文终于录用啦！满足毕业条件了
已经有26人回复
-
2026年机械制造与材料应用国际会议 （ICMMMA 2026）
已经有4人回复
-
Cas 72-43-5需要30g，定制合成，能接单的留言
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 关于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