 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1804 | 回復(fù): 8 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | ||
jjgy_2002木蟲 (著名寫手)
|
[求助]
[求助] 幫看此第一類貝塞爾函數(shù)求解代碼根據(jù)哪篇文獻(xiàn)的計算方法編寫？
|
|
|
[求助] 幫看此第一類貝塞爾函數(shù)求解代碼根據(jù)哪篇文獻(xiàn)的計算方法編寫？ 代碼如下： 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ù)此樓» 猜你喜歡
327求調(diào)劑
已經(jīng)有4人回復(fù)
-
085600材料與化工調(diào)劑
已經(jīng)有3人回復(fù)
-
309求調(diào)劑
已經(jīng)有8人回復(fù)
-
一志愿南航 335分 | 0856 | GPA 4.07 | 有科研經(jīng)歷
已經(jīng)有8人回復(fù)
-
求調(diào)劑
已經(jīng)有5人回復(fù)
-
求化學(xué)調(diào)劑
已經(jīng)有5人回復(fù)
-
085600，材料與化工321分求調(diào)劑
已經(jīng)有10人回復(fù)
-
0703 化學(xué) 求調(diào)劑，一志愿山東大學(xué) 342 分
已經(jīng)有5人回復(fù)
-
一志愿南開大學(xué)0710生物學(xué)359求調(diào)劑
已經(jīng)有3人回復(fù)
-
085600，專業(yè)課化工原理，320分求調(diào)劑
已經(jīng)有4人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 關(guān)于貝塞爾函數(shù)的問題
已經(jīng)有25人回復(fù)
- 一個高一函數(shù)題求解
已經(jīng)有7人回復(fù)
- 不同階第一類貝塞爾函數(shù)的正交性
已經(jīng)有11人回復(fù)
- 【求助】貝塞爾函數(shù)求解
已經(jīng)有4人回復(fù)
tong_xiao
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 11
- 帖子: 3
- 在線: 3.2小時
- 蟲號: 2603090
- 注冊: 2013-08-18
- 專業(yè): 應(yīng)用數(shù)學(xué)方法
8樓2013-11-19 21:06:47
jjgy_2002
木蟲 (著名寫手)
- 應(yīng)助: 16 (小學(xué)生)
- 金幣: 1285.3
- 散金: 1557
- 紅花: 4
- 帖子: 1954
- 在線: 986.4小時
- 蟲號: 646930
- 注冊: 2008-11-05
- 專業(yè): 信號理論與信號處理
3樓2011-06-27 09:58:13
vivian5251
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 35.3
- 帖子: 29
- 在線: 2.7小時
- 蟲號: 1335737
- 注冊: 2011-07-02
- 性別: MM
- 專業(yè): 信息安全
4樓2011-07-03 01:10:06
jjgy_2002
木蟲 (著名寫手)
- 應(yīng)助: 16 (小學(xué)生)
- 金幣: 1285.3
- 散金: 1557
- 紅花: 4
- 帖子: 1954
- 在線: 986.4小時
- 蟲號: 646930
- 注冊: 2008-11-05
- 專業(yè): 信號理論與信號處理
6樓2013-11-13 13:02:23
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085601材料工程找調(diào)劑 +14 | oatmealR 2026-03-29 | 15/750 |
|
|---|---|---|---|---|
|
[考研] 一志愿哈爾濱工業(yè)大學(xué)材料與化工方向336分 +10 | 辰沐5211314 2026-03-26 | 10/500 |
|
|
|
[基金申請] 面上5B能上會嗎？ +4 | redcom 2026-03-29 | 4/200 |
|
|
|
[考研] 070305高分子化學(xué)與物理 304分求調(diào)劑 +12 | c297914 2026-03-28 | 12/600 |
|
|
|
[考研] 329求調(diào)劑 +7 | 星野? 2026-03-26 | 7/350 |
|
|
|
[考研] 321求調(diào)劑 +7 | 璞玉~~ 2026-03-25 | 8/400 |
|
|
|
[考研] 322求調(diào)劑 +7 | 宋明欣 2026-03-27 | 7/350 |
|
|
|
[考研] 071000生物學(xué)求調(diào)劑，初試成績343 +7 | 小小甜面團(tuán) 2026-03-25 | 7/350 |
|
|
|
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +8 | 邱gl 2026-03-27 | 8/400 |
|
|
|
[考研] 266分，求材料冶金能源化工等調(diào)劑 +7 | 哇呼哼呼哼 2026-03-27 | 9/450 |
|
|
|
[考研] 286求調(diào)劑 +4 | 丟掉懶惰 2026-03-27 | 7/350 |
|
|
|
[考研] 張芳銘-中國農(nóng)業(yè)大學(xué)-環(huán)境工程專碩-298 +4 | 手機(jī)用戶 2026-03-26 | 4/200 |
|
|
|
[考研] 308求調(diào)劑 +7 | 墨墨漠 2026-03-25 | 7/350 |
|
|
|
[考研] 化學(xué)308分求調(diào)劑 +8 | 你好明天你好 2026-03-23 | 9/450 |
|
|
|
[考研] 325求調(diào)劑 +3 | Aoyijiang 2026-03-23 | 3/150 |
|
|
|
[考研] 一志愿河工大 081700 276求調(diào)劑 +4 | 地球繞著太陽轉(zhuǎn) 2026-03-23 | 4/200 |
|
|
|
[考研] 環(huán)境專碩324分求調(diào)劑推薦 +5 | 軒小寧—— 2026-03-26 | 5/250 |
|
|
|
[考研] 化學(xué)調(diào)劑一志愿上海交通大學(xué)336分-本科上海211 +4 | 小魚愛有機(jī) 2026-03-25 | 4/200 |
|
|
|
[考研] 0854人工智能方向招收調(diào)劑 +4 | 章小魚567 2026-03-24 | 4/200 |
|
|
|
[論文投稿] 急發(fā)核心期刊論文 +3 | 賢達(dá)問津 2026-03-23 | 5/250 |
|