版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

返回列表

當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

jjgy_2002

木蟲 (著名寫手)

應(yīng)助: 16 (小學(xué)生)
金幣: 1285.3
散金: 1557
紅花: 4
帖子: 1949
在線: 984.6小時(shí)
蟲號(hào): 646930
注冊(cè): 2008-11-05
專業(yè): 信號(hào)理論與信號(hào)處理

[求助] [求助] 幫看此第一類貝塞爾函數(shù)求解代碼根據(jù)哪篇文獻(xiàn)的計(jì)算方法編寫？

[求助] 幫看此第一類貝塞爾函數(shù)求解代碼根據(jù)哪篇文獻(xiàn)的計(jì)算方法編寫？
代碼如下：
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

回復(fù)此樓

» 猜你喜歡

290求調(diào)劑已經(jīng)有7人回復(fù)
一志愿中南化學(xué)（0703）總分337求調(diào)劑已經(jīng)有8人回復(fù)
北科281學(xué)碩材料求調(diào)劑已經(jīng)有5人回復(fù)
求調(diào)劑一志愿南京航空航天大學(xué)289分已經(jīng)有3人回復(fù)
A區(qū)線材料學(xué)調(diào)劑已經(jīng)有5人回復(fù)
材料學(xué)碩297已過(guò)四六級(jí)求調(diào)劑推薦已經(jīng)有11人回復(fù)
316求調(diào)劑已經(jīng)有5人回復(fù)
一志愿西南交通專碩材料355 本科雙非求調(diào)劑已經(jīng)有5人回復(fù)
一志愿武漢理工材料工程專碩調(diào)劑已經(jīng)有8人回復(fù)
329求調(diào)劑已經(jīng)有8人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

關(guān)于貝塞爾函數(shù)的問(wèn)題已經(jīng)有25人回復(fù)
一個(gè)高一函數(shù)題求解已經(jīng)有7人回復(fù)
不同階第一類貝塞爾函數(shù)的正交性已經(jīng)有11人回復(fù)
【求助】貝塞爾函數(shù)求解已經(jīng)有4人回復(fù)

1樓 2011-06-25 14:16:19

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

jjgy_2002

木蟲 (著名寫手)

應(yīng)助: 16 (小學(xué)生)
金幣: 1285.3
散金: 1557
紅花: 4
帖子: 1949
在線: 984.6小時(shí)
蟲號(hào): 646930
注冊(cè): 2008-11-05
專業(yè): 信號(hào)理論與信號(hào)處理

怎沒(méi)人響應(yīng)呢?

回復(fù)此樓

3樓2011-06-27 09:58:13

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 9 個(gè)回答

vivian5251

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 35.3
帖子: 29
在線: 2.7小時(shí)
蟲號(hào): 1335737
注冊(cè): 2011-07-02
性別: MM
專業(yè): 信息安全

【答案】應(yīng)助回帖

想幫助,但不會(huì),希望有人能幫到你

贊一下

回復(fù)此樓

大家好

4樓2011-07-03 01:10:06

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

jjgy_2002

木蟲 (著名寫手)

應(yīng)助: 16 (小學(xué)生)
金幣: 1285.3
散金: 1557
紅花: 4
帖子: 1949
在線: 984.6小時(shí)
蟲號(hào): 646930
注冊(cè): 2008-11-05
專業(yè): 信號(hào)理論與信號(hào)處理

引用回帖:

5樓: Originally posted by tong_xiao at 2013-11-12 21:14:14
金幣有點(diǎn)少哈

呵呵，能幫忙嗎？我多給金幣就是了

贊一下

回復(fù)此樓

6樓2013-11-13 13:02:23

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

tong_xiao

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 11
帖子: 3
在線: 3.2小時(shí)
蟲號(hào): 2603090
注冊(cè): 2013-08-18
專業(yè): 應(yīng)用數(shù)學(xué)方法

引用回帖:

6樓: Originally posted by jjgy_2002 at 2013-11-13 13:02:23
呵呵，能幫忙嗎？我多給金幣就是了

...

我在一本書上看到過(guò)Vb編程的，可以嗎？

贊一下

回復(fù)此樓

7樓2013-11-19 21:03:13

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 9 個(gè)回答

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] A區(qū)線材料學(xué)調(diào)劑 +5	周周無(wú)極 2026-03-20	5/250	2026-03-20 21:33 by laoshidan
[考研] 材料學(xué)碩297已過(guò)四六級(jí)求調(diào)劑推薦 +11	adaie 2026-03-19	11/550	2026-03-20 21:30 by laoshidan
[考研] 求調(diào)劑 +5	Mqqqqqq 2026-03-19	5/250	2026-03-20 20:46 by zhukairuo
[考研] 279分求調(diào)劑一志愿211 +10	chaojifeixia 2026-03-19	11/550	2026-03-20 20:32 by JourneyLucky
[考研] 一志愿西安交通大學(xué) 學(xué)碩 354求調(diào)劑 +3	我想要讀研究生 2026-03-20	3/150	2026-03-20 20:13 by JourneyLucky
[基金申請(qǐng)] 學(xué)校已經(jīng)提交到NSFC，還能修改嗎？ 40+4	babangida 2026-03-19	8/400	2026-03-20 15:58 by babero
[考研] 【考研調(diào)劑】化學(xué)專業(yè) 281分，一志愿四川大學(xué)，誠(chéng)心求調(diào)劑 +6	吃吃吃才有意義 2026-03-19	6/300	2026-03-20 10:47 by 盡舜堯1
[論文投稿] 申請(qǐng)回稿延期一個(gè)月，編輯同意了。但系統(tǒng)上的時(shí)間沒(méi)變，給編輯又寫郵件了，沒(méi)回復(fù) 10+3	wangf9518 2026-03-17	4/200	2026-03-19 23:55 by babero
[考研] 復(fù)試調(diào)劑 +4	z1z2z3879 2026-03-14	6/300	2026-03-19 17:18 by fei626-918
[考研] 0703化學(xué)調(diào)劑 +4	18889395102 2026-03-18	4/200	2026-03-19 16:13 by 30660438
[考研] 化學(xué)求調(diào)劑 +3	臨澤境llllll 2026-03-17	4/200	2026-03-19 13:59 by houyaoxu
[考研] 0703化學(xué)調(diào)劑 +5	pupcoco 2026-03-17	8/400	2026-03-19 13:58 by houyaoxu
[考研] 330求調(diào)劑 +3	小材化本科 2026-03-18	3/150	2026-03-18 21:55 by 無(wú)懈可擊111
[考研] 311求調(diào)劑 +4	冬十三 2026-03-18	4/200	2026-03-18 21:47 by 盡舜堯1
[考研] 354求調(diào)劑 +4	Tyoumou 2026-03-18	7/350	2026-03-18 21:45 by Tyoumou
[考研] 085601專碩，總分342求調(diào)劑，地區(qū)不限 +5	share_joy 2026-03-16	5/250	2026-03-18 14:48 by haxia
[考研] 收復(fù)試調(diào)劑生 +4	雨后秋荷 2026-03-18	4/200	2026-03-18 14:16 by elevennnne
[考研] 280求調(diào)劑 +6	咕嚕曉曉 2026-03-18	7/350	2026-03-18 11:25 by 無(wú)際的草原
[考研] 275求調(diào)劑 +4	太陽(yáng)花天天開(kāi)心 2026-03-16	4/200	2026-03-17 10:53 by 功夫瘋狂
[考研] 318求調(diào)劑 +3	Yanyali 2026-03-15	3/150	2026-03-16 16:41 by houyaoxu