应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 数学 » 计算数学 » [求助] 帮看此第一类贝塞尔函数求解代码根据哪篇文献的计算方法编写？
北京石油化工学院2026年研究生招生接收调剂公告
91/1返回列表
查看: 1811  |  回复: 8
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2011-06-25 14:16:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jjgy_2002

木虫 (著名写手)
2楼2011-06-25 16:52:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jjgy_2002

木虫 (著名写手)
怎没人响应呢?
回复此楼
3楼2011-06-27 09:58:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

vivian5251

新虫 (初入文坛)

【答案】应助回帖

想帮助,但不会,希望有人能帮到你
赞一下 回复此楼
大家好
4楼2011-07-03 01:10:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tong_xiao

新虫 (初入文坛)
金币有点少哈
回复此楼
5楼2013-11-12 21:14:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jjgy_2002

木虫 (著名写手)
引用回帖:
5楼: Originally posted by tong_xiao at 2013-11-12 21:14:14
金币有点少哈

呵呵，能帮忙吗？我多给金币就是了
赞一下 回复此楼
6楼2013-11-13 13:02:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tong_xiao

新虫 (初入文坛)
引用回帖:
6楼: Originally posted by jjgy_2002 at 2013-11-13 13:02:23
呵呵，能帮忙吗？我多给金币就是了...

我在一本书上看到过Vb编程的，可以吗？
赞一下 回复此楼
7楼2013-11-19 21:03:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tong_xiao

新虫 (初入文坛)
就是一个visual basic 常用算法集里面有一个此种算法，不过是vb，跟你的matlab大同小异哦
赞一下 回复此楼
8楼2013-11-19 21:06:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jjgy_2002

木虫 (著名写手)
引用回帖:
7楼: Originally posted by tong_xiao at 2013-11-19 21:03:13
我在一本书上看到过Vb编程的，可以吗？...

可以啊
回复此楼
9楼2013-11-20 15:51:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 jjgy_2002 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 288求调剂 一志愿哈工大 材料与化工 +5 洛神哥哥 2026-03-31 5/250 2026-03-31 11:23 by jp9609
[考研] 一志愿085600中科院宁波所276分求调剂 +8 材料学257求调剂 2026-03-28 9/450 2026-03-31 10:27 by cal0306
[考研] 考研生物与医药调剂 +3 铁憨憨123425 2026-03-31 3/150 2026-03-31 09:42 by longlotian
[考研] 福建理工大学材料学院先进合金团队招收考研调剂学生 +3 大华金商都 2026-03-30 4/200 2026-03-31 01:04 by 方英俊602
[考研] 288资源与环境专硕求调剂，不限专业，有学上就行 +11 lllllos 2026-03-30 11/550 2026-03-30 21:19 by 研究僧导导
[考研] 327求调剂 +5 小卡不卡. 2026-03-29 5/250 2026-03-30 19:30 by Wang200018
[考研] 324求调剂 +9 hanamiko 2026-03-26 11/550 2026-03-30 14:27 by JourneyLucky
[基金申请] 面上5B能上会吗？ +7 redcom 2026-03-29 7/350 2026-03-30 12:44 by 爱搞研究的小孩
[考研] 一志愿中南大学化学0703总分337求调剂 +6 niko- 2026-03-27 6/300 2026-03-30 10:25 by herarysara
[考研] 295求调剂 +5 wei-5 2026-03-26 5/250 2026-03-30 08:34 by 探123
[考研] 考研调剂 +7 小蜡新笔 2026-03-29 7/350 2026-03-29 19:00 by 学员8dgXkO
[考研] 299求调剂 +10 15188958825 2026-03-25 10/500 2026-03-29 17:51 by 王亮_大连医科大
[考研] 330分求调剂 +5 qzenlc 2026-03-29 5/250 2026-03-29 07:37 by 无际的草原
[考研] 药学105500求调剂 +3 Ssun。。 2026-03-28 3/150 2026-03-28 11:24 by lxf170613
[考研] 安徽大学专硕生物与医药专业(086000)324分，英语已过四六级，六级521，求调剂 +4 美味可乐鸡翅 2026-03-26 4/200 2026-03-27 15:27 by 星空星月
[考研] 调剂求收留 +7 果然有我 2026-03-26 7/350 2026-03-27 00:26 by wxiongid
[考研] 环境专硕324分求调剂推荐 +5 轩小宁—— 2026-03-26 5/250 2026-03-26 12:05 by i_cooler
[考研] 化学调剂一志愿上海交通大学336分-本科上海211 +4 小鱼爱有机 2026-03-25 4/200 2026-03-26 10:19 by aa331100
[考研] 机械学硕总分317求调剂！！！！ +4 Acaciad 2026-03-25 4/200 2026-03-25 19:59 by hanserlol
[考研] 0854人工智能方向招收调剂 +4 章小鱼567 2026-03-24 4/200 2026-03-25 13:29 by 2177681040
信息提示
关闭
请填处理意见
关闭