| 查看: 1135 | 回复: 1 | ||
wangshitong银虫 (小有名气)
|
[求助]
hankel逆变换
|
| 哪位用过Hankel逆变换，我从网上搜到一些matlab的代码，但还是搞不定。希望交流一下，谢谢啦！ |
回复此楼» 猜你喜欢
323求调剂
已经有6人回复
-
一志愿北京化工大学 070300 学硕 336分 求调剂
已经有4人回复
-
352求调剂
已经有3人回复
-
一志愿东华大学化学070300，求调剂
已经有8人回复
-
277材料科学与工程080500求调剂
已经有7人回复
-
317求调剂
已经有18人回复
-
293求调剂
已经有5人回复
-
280分求调剂 一志愿085802
已经有7人回复
-
0854电子信息求调剂
已经有3人回复
-
263求调剂
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 关于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
