| 查看: 165 | 回复: 1 | |||
| 当前主题已经存档。 | |||
| 【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 lzjuxmch 的 10 个金币 | |||
lzjuxmch铁虫 (小有名气)
|
[交流]
【求助】hankel变换翻译成C语言
|
||
|
本人急需hankel变换的C程序，目前只能找到matlab中hankel变换的原代码，但是本人不懂matlab，请各位帮个忙能否帮把下列代码翻译成C或者C++语言？提供下列程序的思路也可。万分感谢！ %[H,I]=ht(h,r,k,I) %----------------- % %Hankel transform of order 0. % %Input: % h Signal h(r) % r Radial positions [m] {0:numel(h)-1} % k Spatial frequencies [rad/m] {pi/numel(h)*(0:numel(h)-1)} % I Integration kernel ? {default} % %Output: % H Spectrum H(k) % I Integration kernel % % % ? If the integration kernel is missing, it is % recomputed from the Bessel functions (slow). % % Marcel Leutenegger ?June 2006 % function [H,I]=ht(h,r,k,I) if sum(size(h) > 1) > 1 error('Signal must be a vector.'); end if nargin < 2 | isempty(r) r=0:numel(h)-1; else [r,w]=sort(r(  .'); h=h(w); end if nargin < 3 | isempty(k) k=pi/numel(h)*(0:numel(h)-1); end if nargin < 4 | isempty(I) r=[(r(2:end) + r(1:end-1))/2 r(end)]; I=2*pi./k( *r.*besselj(1,k(*r); I(k == 0, =pi*r.*r; I=I - [zeros(numel(k),1) I(:,1:end-1)]; elseif exist('w','var') I=I(:,w); end H=reshape(I*h( ,size(k));[ Last edited by mingdong on 2009-4-16 at 17:53 ] |
回复此楼» 猜你喜欢
【修正版】单环用磁约束低速电子实验方案（简化版）
已经有0人回复
-
桂林理工大学物理学专业招收调剂，还有三个名额！！！
已经有22人回复
-
物理学I论文润色/翻译怎么收费?
已经有74人回复
-
考博自荐
已经有4人回复
-
山东大学第二批博士研究生招生
已经有0人回复
-
中国科学院东莞材料科学与技术研究所-2026年博士招生-吴昊研究员-磁学与自旋电子学
已经有0人回复
-
《电磁学》教材推荐
已经有1人回复
-
【急招】合肥工大核聚变材料计算方向2026级工程博士生
已经有4人回复
-
大豆异黄酮分离
已经有0人回复
-
湖南大学材料学院急招2026年博士生，临时增加一名博士联培指标
已经有10人回复
-
天津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生1-2名
已经有0人回复
lzjuxmch
铁虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 299.3
- 帖子: 112
- 在线: 44.3小时
- 虫号: 490252
- 注册: 2008-01-03
- 专业: 凝聚态物性 II ：电子结构
2楼2009-04-16 16:58:38
30