| 查看: 596 | 回复: 4 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
怎样加快这个代码的运行？
|
||
|
Hi， 我要在Mathematica里做Hankel transform（它就是像傅立叶变换，但是以Bessel函数为base），具体的转换公式可以参考wiki，http://en.wikipedia.org/wiki/Hankel_transform。但是我在Mathematica里直接用数值积分的办法来进行转换，好像很慢。 我上传我的代码。其中 a : 我要转换的离散数值，可以是任何函数。第一个index是坐标r，第二个index是具体数值。 然后我做interpolation得到他的插值的函数。 然后用NIntegrate来算Hankel transform。 像傅立叶转换一样，Hankel transform也有逆转换。倒数第二行就是逆转换的定义，然后最后一行给出的数值应该是落于{0.9986238244, 0.0002427139}, {1.0101920314, 0.0002393707}之间，这个是从a里读出来的数据。 但是这个运行的时间好像太长，我运行了20分钟，还没有结果。不知道大家有没有更好的办法，我Mathematica不是很熟。谢谢。 |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : xiaomuchong.nb
2014-10-13 20:40:46, 45.25 K
» 猜你喜欢
调剂推荐
已经有5人回复
-
中国科学院深圳先进技术研究院-光纤传感课题组招生-中国科学院大学、深圳理工大学联培
已经有4人回复
-
总分322求生物学/生化与分子/生物信息学相关调剂
已经有5人回复
-
0703化学338求调剂！
已经有3人回复
-
284求调剂
已经有7人回复
-
271求调剂
已经有6人回复
-
312求调剂
已经有9人回复
-
352求调剂
已经有4人回复
-
085602化学工程求调剂。
已经有3人回复
-
085600材料与化工306
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 运行这段代码到出图片,你的电脑用了多久?
已经有33人回复
- fortran 源代码在powersation4.0运行问题
已经有5人回复
- for循环语句的加速优化求助！
已经有8人回复
- 关于手写字符的卷积神经网络代码运行错误
已经有5人回复
- 求助，怎样去掉双重for循环加快程序运行
已经有10人回复
- 用MPICH2运行并行程序，代码中的printf语句只有在所有进程都结束后才会打印出来么？
已经有5人回复
- 在linux下运行ADF出错
已经有18人回复
- linux下如何让一个程序运行n次？
已经有5人回复
- 用matlab运行代码出不来图像 ！！
已经有3人回复
- linux 终端如何运行.csh格式的程序代码
已经有4人回复
- 源代码如何生成可执行文件？
已经有5人回复
- 审稿人要求开源代码或提供可执行程序，如何拒绝？
已经有5人回复
- 【求助】如何将源代码在linux下运行？
已经有8人回复
- forcal 代码运行后出错
已经有4人回复
- 【求助】一个matlab窗口占25%CPU，如何提高？
已经有11人回复
- 【求助】此段C语言代码运行不了
已经有9人回复
- 【讨论】如何过写代码运行心理上这一关
已经有4人回复
walk1997
金虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 4676.2
- 红花: 22
- 帖子: 1066
- 在线: 798.1小时
- 虫号: 416039
- 注册: 2007-06-29
- 性别: GG
- 专业: 粒子物理学和场论
3楼2014-10-14 12:44:09
walk1997
金虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 4676.2
- 红花: 22
- 帖子: 1066
- 在线: 798.1小时
- 虫号: 416039
- 注册: 2007-06-29
- 性别: GG
- 专业: 粒子物理学和场论
|
几个问题 1. 数学，这个变换的存在对这个函数是有要求的吧， 要求这个函数的绝对值|f(r)|*r在0到无穷间的积分区域收敛吧. 看你的函数，原始差不多只是定义在小量到20左右，并且在这附近都是常数，一方面你插值以后外推到无穷远处不太合理，另外一方面要是在大x都取成常数的话，也不符合做这个变换的条件。我觉得你这个时候不应该用到变换吧，而应该用类似的无穷级数展开,比如.. f[w]=\sum_{1} BesselJ[0,\sqrt{mu_m \rho}] 2. 另外，如果直接截断第一个积分到20，给出的结果感觉差不多 －------------- r1 = 1; NIntegrate[ ifun[r] *r* BesselJ[0, k*r]*BesselJ[0, k*r1]*k, {k, 0, Infinity}, {r, 0, 20}, PrecisionGoal -> 2, MaxRecursion -> 40] // AbsoluteTiming －----------- {1092.899510,0.000242286} |
2楼2014-10-14 12:34:56
4楼2014-10-14 22:05:29
大家的诚
银虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 655.8
- 散金: 3
- 帖子: 87
- 在线: 141.8小时
- 虫号: 2514004
- 注册: 2013-06-19
- 性别: GG
- 专业: 半导体电子器件
5楼2015-12-10 09:09:13
6
