| 7 | 1/1 | 返回列表 |
| 查看: 1911 | 回復: 6 | |||
[交流]
用matlab求bessel函數(shù)的根
|
|
問題看上去很簡單,方程如下所示: 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,但是都沒能給出正確的結果。 求高手指點! |
» 搶金幣啦!回帖就可以得到:
+1/280
+1/97
+1/81
+1/78
+1/71
+1/35
+1/28
+1/17
+1/14
+1/13
+1/7
+1/7
+1/7
+1/7
+1/5
+1/5
+1/5
+1/2
+1/1
+1/1
|
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賦不同的值得到的結果不大一樣,這里給v賦值為10。 |
| 7 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 311求調(diào)劑 +9 | 冬十三 2026-03-15 | 10/500 |
|
|---|---|---|---|---|
|
[考研] 材料,紡織,生物(0856、0710),化學招生啦 +3 | Eember. 2026-03-17 | 7/350 |
|
|
[考研] 299求調(diào)劑 +4 | △小透明* 2026-03-17 | 4/200 |
|
|
[考研] 085600材料與化工求調(diào)劑 +5 | 緒幸與子 2026-03-17 | 5/250 |
|
|
[考研] 材料與化工專碩調(diào)劑 +5 | heming3743 2026-03-16 | 5/250 |
|
|
[碩博家園] 深圳大學碩士招生(2026秋,傳感器方向,僅錄取第一志愿) +4 | xujiaoszu 2026-03-11 | 9/450 |
|
|
[考研] 274求調(diào)劑 +5 | 時間點 2026-03-13 | 5/250 |
|
|
[考研] 304求調(diào)劑 +5 | 素年祭語 2026-03-15 | 5/250 |
|
|
[考研] 中科大材料專碩319求調(diào)劑 +3 | 孟鑫材料 2026-03-13 | 3/150 |
|
|
[考研] 復試調(diào)劑 +3 | 呼呼?~+123456 2026-03-14 | 3/150 |
|
|
[考研] 333求調(diào)劑 +3 | 球球古力 2026-03-11 | 3/150 |
|
|
[考研] (081700)化學工程與技術-298分求調(diào)劑 +12 | 11啦啦啦 2026-03-11 | 35/1750 |
|
|
[考研] 329求調(diào)劑 +3 | miaodesi 2026-03-12 | 4/200 |
|
|
[考研] 材料與化工085600調(diào)劑求老師收留 +9 | jiaanl 2026-03-11 | 9/450 |
|
|
[考研] 0703化學求調(diào)劑 +7 | 綠豆芹菜湯 2026-03-12 | 7/350 |
|
|
[考研] 材料專碩350 求調(diào)劑 +4 | 王金科 2026-03-12 | 4/200 |
|
|
[考研] 土木第一志愿276求調(diào)劑,科研和技能十分豐富,求新興方向的導師收留 +3 | 土木小天才 2026-03-12 | 3/150 |
|
|
[考研] 328化工專碩求調(diào)劑 +4 | 。,。,。,。i 2026-03-12 | 4/200 |
|
|
[考研] 308求調(diào)劑 +3 | 是Lupa啊 2026-03-12 | 3/150 |
|
|
[考研] 270求調(diào)劑 085600材料與化工專碩 +3 | YXCT 2026-03-11 | 3/150 |
|