| 查看: 162 | 回复: 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 ] |
回复此楼» 猜你喜欢
重庆交大26年硕士生招生拟调剂通知已出！欢迎加入光子学微结构与器件课题组。
已经有0人回复
-
UJN物理学专业调剂
已经有12人回复
-
物理学I论文润色/翻译怎么收费?
已经有152人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有0人回复
-
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组招收2026硕士研究生1名
已经有0人回复
-
山西大同大学物理学硕士研究生招收调剂生
已经有0人回复
-
[调剂信息]211智能人工感知方向国家青年特聘专家课题组招收调剂研究生
已经有0人回复
-
［调剂信息］211智能人工感知方向国家青年特聘专家课题组招收调剂研究生
已经有0人回复
-
[调剂信息]211智能人工感知方向国家青年特聘专家课题组招收调剂研究生
已经有0人回复
-
一志愿物理所凝聚态物理352求调剂
已经有0人回复
lzjuxmch
铁虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 299.3
- 帖子: 112
- 在线: 44.3小时
- 虫号: 490252
- 注册: 2008-01-03
- 专业: 凝聚态物性 II ：电子结构
2楼2009-04-16 16:58:38
