| 查看: 1592 | 回复: 8 | ||
jjgy_2002木虫 (著名写手)
|
[求助]
[求助] 帮看此第一类贝塞尔函数求解代码根据哪篇文献的计算方法编写？
|
|
[求助] 帮看此第一类贝塞尔函数求解代码根据哪篇文献的计算方法编写？ 代码如下： function Jx = bessel(n,x) format long; if n == 0 a = [5.7568490574e10;-1.3362590354e10;6.516196407e8; -1.121442418e7;7.739233017e4;-1.849052456e2]; b = [5.7568490411e10;1.029532985e9;9.494680718e6; 5.927264853e4;2.678532712e2;1.0]; r =[1.0;-0.1098628627e-2;0.2734510407e-4; -0.2073370639e-5;0.2093887211e-6]; s =[-0.1562499995e-1;0.1430488765e-3;-0.6911147651e-5; 0.7621095161e-6;-0.934945152e-7]; if abs(x) < 8 J0_u = a(1); J0_d = b(1); for i=2:6 J0_u = J0_u + a(i)*power(x,2*(i-1)); J0_d = J0_d + b(i)*power(x,2*(i-1)); end J0 = J0_u/J0_d; else z = 8/abs(x); sita = abs(x) - pi/4; R0 = r(1); S0 = s(1); for i=2:5 R0 = R0 + r(i)*power(z,2*(i-1)); S0 = S0 + s(i)*power(z,2*(i-1)); end J0 = (R0*cos(sita)-z*S0*sin(sita))*sqrt(2/pi/abs(x)); end Jx = J0; else if 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; else if abs(x) > n y = 2/x; bess1 = bessel(0,x); bess2 = bessel(1,x); for j=1:n-1 Jx = j*y*bess2 - bess1; bess1 = bess2; bess2 = Jx; end else M = 2*floor(((n + floor(sqrt(40*n))))/2); JS = zeros(M+2 ,1); JS(M+2) = 0; JS(M+1) = 1; y = 2/x; bSum = 0; for k=M:-1:1 JS(k) = JS(k+1)*y*k - JS(k+2); end BK = JS(1); for k=1:floor(M/2) BK = BK + 2*JS(2*k+1); end Jx = JS(n+1)/BK; end end end |
回复此楼» 猜你喜欢
读博
已经有5人回复
-
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗
已经有13人回复
-
博士申请都是内定的吗？
已经有6人回复
-
之前让一硕士生水了7个发明专利，现在这7个获批发明专利的维护费可从哪儿支出哈？
已经有5人回复
-
博士读完未来一定会好吗
已经有29人回复
-
投稿精细化工
已经有4人回复
-
高职单位投计算机相关的北核或SCI四区期刊推荐，求支招！
已经有4人回复
-
导师想让我从独立一作变成了共一第一
已经有9人回复
-
心脉受损
已经有5人回复
-
Springer期刊投稿求助
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
关于贝塞尔函数的问题
已经有25人回复- 一个高一函数题求解
已经有7人回复
- 不同阶第一类贝塞尔函数的正交性
已经有11人回复
- 【求助】贝塞尔函数求解
已经有4人回复
 |
2楼2011-06-25 16:52:32
jjgy_2002
木虫 (著名写手)
- 应助: 16 (小学生)
- 金币: 1285.3
- 散金: 1556
- 红花: 4
- 帖子: 1947
- 在线: 984.2小时
- 虫号: 646930
- 注册: 2008-11-05
- 专业: 信号理论与信号处理
3楼2011-06-27 09:58:13
4楼2011-07-03 01:10:06
5楼2013-11-12 21:14:14
jjgy_2002
木虫 (著名写手)
- 应助: 16 (小学生)
- 金币: 1285.3
- 散金: 1556
- 红花: 4
- 帖子: 1947
- 在线: 984.2小时
- 虫号: 646930
- 注册: 2008-11-05
- 专业: 信号理论与信号处理
6楼2013-11-13 13:02:23
7楼2013-11-19 21:03:13
8楼2013-11-19 21:06:47
jjgy_2002
木虫 (著名写手)
- 应助: 16 (小学生)
- 金币: 1285.3
- 散金: 1556
- 红花: 4
- 帖子: 1947
- 在线: 984.2小时
- 虫号: 646930
- 注册: 2008-11-05
- 专业: 信号理论与信号处理
9楼2013-11-20 15:51:13