| 查看: 582 | 回复: 0 | |||
jjgy_2002木虫 (著名写手)
|
[求助]
[求助] 第一类贝塞尔函数求解代码的计算复杂度如何计算？
|
|
[求助] 第一类贝塞尔函数求解代码的计算复杂度如何计算？ 代码如下： function Jx = bessel(n,x) format long; %n = 1 仅计算一阶贝塞尔函数 c = [7.2362614232e10;-7.895059235e9;2.423968531e8; -2.972611439e6;1.570448260e4;-3.016036606e1]; d = [1.44725228443e11;2.300535178e9;1.858330474e7; 9.944743994e4;3.769991397e2;1.0]; p = [1.0;0.183105e-2;-0.3516396496e-4; 0.2457520174e-5;-0.240337019e-6]; q = [0.4687499995e-1;-0.2002690873e-3; 0.8449199096e-5;-0.88228987e-6; 0.105787412e-6]; if abs(x) < 8 J1_u = c(1); J1_d = d(1); for i=2:6 J1_u = J1_u + c(i)*power(x,2*(i-1)); J1_d = J1_d + d(i)*power(x,2*(i-1)); end J1 = x*J1_u/J1_d; else z = 8/abs(x); fi = abs(x) - 3*pi/4; P1 = p(1); Q1 = q(1); for i=2:5 P1 = P1 + p(i)*power(z,2*(i-1)); Q1 = Q1 + q(i)*power(z,2*(i-1)); end J1 = (P1*cos(fi)-z*Q1*sin(fi))*sqrt(2/pi/abs(x)); if x<0 J1 = -J1; end end Jx = J1; |
回复此楼» 猜你喜欢
之前让一硕士生水了7个发明专利，现在这7个获批发明专利的维护费可从哪儿支出哈？
已经有9人回复
-
心脉受损
已经有6人回复
-
博士申请都是内定的吗？
已经有8人回复
-
博士读完未来一定会好吗
已经有32人回复
-
读博
已经有5人回复
-
投稿精细化工
已经有4人回复
-
高职单位投计算机相关的北核或SCI四区期刊推荐，求支招！
已经有4人回复
-
导师想让我从独立一作变成了共一第一
已经有9人回复
-
Springer期刊投稿求助
已经有4人回复
找到一些相关的精华帖子，希望有用哦~
请教关于本机的Dmol计算问题——望高手指点
已经有6人回复- CASTEP在计算是出现问题啦，求大侠帮忙啊。。。
已经有13人回复
- 用CHEM3D软件PM3理论水平计算分子轨道却无法显示轨道能量
已经有5人回复
- 学机械的小硕从事通信设备制造行业的结构工程师，个人发展前景和空间如何？
已经有6人回复
- 为什么含有稀土元素的计算就这么难呢
已经有16人回复
- 关键词：钻孔、配钻、工装，告诉求解啊！！！！！！
已经有9人回复
- 天然气燃烧的能量利用效率应该怎么计算？
已经有4人回复
- 纯理论计算投稿的问题
已经有21人回复
- Matlab求解偏微分方程组
已经有22人回复
- 求助DFT计算中CuCl2在γ-Al2O3 表面的吸附结构？
已经有5人回复
- [求助] 帮看此第一类贝塞尔函数求解代码根据哪篇文献的计算方法编写？
已经有8人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人