版塊導(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í)熱門版塊排行榜

返回列表

xiaolobster

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 140
帖子: 10
在線: 37.3小時(shí)
蟲號(hào): 964557

[交流] 用matlab求bessel函數(shù)的根

問題看上去很簡(jiǎn)單，方程如下所示：
v*sqrt（1-b）*J1（v*sqrt(1-b)）/J0(v*sqrt(1-b))=-v*sqrt(b)*K1(v*sqrt(b))/K0(v*sqrt(b))
其中b為未知數(shù)。v為常數(shù)，J1我第一類第1階貝塞爾函數(shù)，J0為第一類第0階貝塞爾函數(shù)；K1為第二類第1階貝塞爾函數(shù)，K2為第二類第0階貝塞爾函數(shù)。

問題：通過解上述方程，求解b的值。
我在matlab中嘗試了fzero，和solve，但是都沒能給出正確的結(jié)果。
求高手指點(diǎn)！

回復(fù)此樓

» 猜你喜歡

268求調(diào)劑已經(jīng)有3人回復(fù)
302求調(diào)劑已經(jīng)有5人回復(fù)
311求調(diào)劑已經(jīng)有10人回復(fù)
被我言中：新模板不強(qiáng)調(diào)格式了，假專家開始管格式了已經(jīng)有4人回復(fù)
303求調(diào)劑已經(jīng)有5人回復(fù)
考研求調(diào)劑已經(jīng)有4人回復(fù)
085601專碩，總分342求調(diào)劑，地區(qū)不限已經(jīng)有4人回復(fù)
085601材料工程專碩求調(diào)劑已經(jīng)有5人回復(fù)
一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑已經(jīng)有6人回復(fù)
能源材料化學(xué)課題組招收碩士研究生8-10名已經(jīng)有6人回復(fù)

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

matlab如何求解含有bessel方程中的參數(shù) 已經(jīng)有5人回復(fù)
求助：Matlab中如何繪制三維隱函數(shù)圖像已經(jīng)有20人回復(fù)
誰給解釋一下擴(kuò)散系數(shù)的單位問題？？已經(jīng)有5人回復(fù)
數(shù)值求解N-S方程一般用什么編程��？？？已經(jīng)有9人回復(fù)
Matlab畫隱函數(shù)，ezplot 已經(jīng)有3人回復(fù)
求助matlab的發(fā)fgoalattain函數(shù) 已經(jīng)有4人回復(fù)
哪些書里面有講解關(guān)于含有積分方程的未知數(shù)求解問題�。糠浅８兄x！已經(jīng)有9人回復(fù)
函數(shù)有第一類間斷點(diǎn)，則不存在原函數(shù) 已經(jīng)有3人回復(fù)
求matlab畫一個(gè)函數(shù)的圖已經(jīng)有6人回復(fù)
如何用Matlab畫z關(guān)于x和y的函數(shù)圖？謝謝大家！已經(jīng)有8人回復(fù)
求助：能不能幫忙寫一個(gè)用matlab求一個(gè)函數(shù)的傅里葉變換的程序已經(jīng)有16人回復(fù)
如何根據(jù)已知函數(shù)畫函數(shù)曲線？謝謝已經(jīng)有4人回復(fù)
【求助】請(qǐng)問怎么求函數(shù)方程的函數(shù)表達(dá)式？已經(jīng)有11人回復(fù)
【求助】關(guān)于matlab中的函數(shù)代換已經(jīng)有10人回復(fù)

» 搶金幣啦！回帖就可以得到:

查看全部散金貼

1樓 2013-11-28 10:33:28

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

bruceleem

銅蟲 (著名寫手)

應(yīng)助: 3 (幼兒園)
金幣: 3710
帖子: 1033
在線: 284.7小時(shí)
蟲號(hào): 946615

★ ★ ★ ★ ★ ★
xiaolobster: 金幣+5, ★有幫助 2013-11-29 09:03:39
xzhdty: 金幣+1, 謝謝參與 2013-11-30 10:51:42

function bessel_problem
%-----------------------------------------
function f=besselfunc(b,v)
f=v.*sqrt(1-b).*besselj(1,v.*sqrt(1-b))./besselj(0,v.*sqrt(1-b))-(-v.*sqrt(b).*bessely(1,(v.*sqrt(b)))./bessely(0,(v.*sqrt(b))));
end
%-----------------------------------------
v=10;
b0=0.1;
[b,yb,exitflag]=fsolve(@(b)besselfunc(b,v),b0);
b
yb
exitflag
end

將代碼保存為m文件，在命令窗口中輸入bessel_problem即可求解。給v賦不同的值得到的結(jié)果不大一樣，這里給v賦值為10。

贊一下

回復(fù)此樓

2樓2013-11-28 19:12:56

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

xiaolobster

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 140
帖子: 10
在線: 37.3小時(shí)
蟲號(hào): 964557

引用回帖:

2樓: Originally posted by bruceleem at 2013-11-28 19:12:56
function bessel_problem
%-----------------------------------------
function f=besselfunc(b,v)
f=v.*sqrt(1-b).*besselj(1,v.*sqrt(1-b))./besselj(0,v.*sqrt(1-b))-(-v.*sqrt(b).*bessely(1,(v.*sqrt(b))) ...

非常感謝。主要是我的版本有點(diǎn)低，升級(jí)版本之后，用solve能解決。不過仍然十分感謝您的熱情回復(fù)。5個(gè)金幣要給，必須的。

贊一下

回復(fù)此樓

3樓2013-11-29 09:04:49

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

bruceleem

銅蟲 (著名寫手)

應(yīng)助: 3 (幼兒園)
金幣: 3710
帖子: 1033
在線: 284.7小時(shí)
蟲號(hào): 946615

祝樓主好運(yùn)！

回復(fù)此樓

4樓2013-11-29 09:23:25

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

真心妹夫

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 0.5
帖子: 1
在線:
蟲號(hào): 2974012

★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖

請(qǐng)問樓主最后用solve怎么解決的？

贊一下

回復(fù)此樓

5樓2014-02-18 04:21:59

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

九垓

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 1
帖子: 2
在線:
蟲號(hào): 3878447

★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖

樓主問題解決了嗎/

贊一下

回復(fù)此樓

6樓2015-05-19 14:17:14

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

九垓

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 1
帖子: 2
在線:
蟲號(hào): 3878447

★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖

引用回帖:

3樓: Originally posted by xiaolobster at 2013-11-29 09:04:49
非常感謝。主要是我的版本有點(diǎn)低，升級(jí)版本之后，用solve能解決。不過仍然十分感謝您的熱情回復(fù)。5個(gè)金幣要給，必須的。...

對(duì)復(fù)雜貝塞爾形式的化簡(jiǎn)，，樓主可好

小弟求助

贊一下

回復(fù)此樓

7樓2015-05-19 14:20:12

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

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 xiaolobster 的主題更新

返回列表

普通表情龍兔虎貓高級(jí)回復(fù) (可上傳附件)

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

[考研] 085601材料工程專碩求調(diào)劑 +5	慕寒mio 2026-03-16	5/250	2026-03-17 21:31 by hmn_wj
[考研] 344求調(diào)劑 +4	knight344 2026-03-16	4/200	2026-03-17 17:27 by ruiyingmiao
[考研] 332求調(diào)劑 +6	Zz版 2026-03-13	6/300	2026-03-17 17:03 by ruiyingmiao
[考研] 302求調(diào)劑 +4	小賈同學(xué)123 2026-03-15	8/400	2026-03-17 10:33 by 小賈同學(xué)123
[考研] 304求調(diào)劑 +5	素年祭語 2026-03-15	5/250	2026-03-16 17:00 by 我的船我的海
[考研] 318求調(diào)劑 +3	Yanyali 2026-03-15	3/150	2026-03-16 16:41 by houyaoxu
[考研] 304求調(diào)劑 +3	曼殊2266 2026-03-14	3/150	2026-03-16 16:39 by houyaoxu
[考研] 070300化學(xué)學(xué)碩求調(diào)劑 +6	太想進(jìn)步了0608 2026-03-16	6/300	2026-03-16 16:13 by kykm678
[考研] 283求調(diào)劑 +10	小樓。 2026-03-12	14/700	2026-03-16 16:08 by 13811244083
[考研] 0703化學(xué)調(diào)劑，求各位老師收留 +8	秋有木北 2026-03-14	8/400	2026-03-16 15:21 by 哦哦123
[考研] 070303 總分349求調(diào)劑 +3	LJY9966 2026-03-15	5/250	2026-03-16 14:24 by xwxstudy
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +3	自由煎餅果子 2026-03-16	3/150	2026-03-16 14:10 by 運(yùn)氣yunqi
[考研] 中科大材料專碩319求調(diào)劑 +3	孟鑫材料 2026-03-13	3/150	2026-03-14 18:10 by houyaoxu
[考研] 材料與化工求調(diào)劑一志愿 985 總分 295 +8	dream…… 2026-03-12	8/400	2026-03-13 22:17 by 星空星月
[考研] ［0860］321分求調(diào)劑，ab區(qū)皆可 +4	寶貴熱 2026-03-13	4/200	2026-03-13 22:01 by 星空星月
[考研] 一志愿西南交大，材料專碩317求調(diào)劑 +5	lx8568 2026-03-11	5/250	2026-03-13 21:43 by peike
[考研] 26調(diào)劑/材料科學(xué)與工程/總分295/求收留 +9	2026調(diào)劑俠 2026-03-12	9/450	2026-03-13 20:46 by 18595523086
[考研] 工科278分求調(diào)劑 +5	周慢熱啊 2026-03-12	7/350	2026-03-13 15:49 by JourneyLucky
[考研] 一志愿山大07化學(xué) 332分四六級(jí)已過本科山東雙非求調(diào)劑！ +3	不想理你 2026-03-12	3/150	2026-03-13 14:18 by JourneyLucky
[考研] 0856化學(xué)工程280分求調(diào)劑 +4	shenzxsn 2026-03-11	4/200	2026-03-13 11:55 by ymwdoctor

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

[交流] 用matlab求bessel函數(shù)的根

» 猜你喜歡

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

» 搶金幣啦！回帖就可以得到: 查看全部散金貼

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

» 搶金幣啦！回帖就可以得到:

查看全部散金貼