式は I 0 (a*x)=b で、a、b は定数、I 0 (y) はゼロ次ベッセル関数です。
I 0 (x)=1/(2*pi)*integral(x*cos(t)) dt [0 ->2*pi から]
a と b が与えられたときの x の値を取得したい。結果が欲しいだけなので、Matlabで解く手順を全て実装する必要はありませんし、おおよその答えでも構いません。
式は I 0 (a*x)=b で、a、b は定数、I 0 (y) はゼロ次ベッセル関数です。
I 0 (x)=1/(2*pi)*integral(x*cos(t)) dt [0 ->2*pi から]
a と b が与えられたときの x の値を取得したい。結果が欲しいだけなので、Matlabで解く手順を全て実装する必要はありませんし、おおよその答えでも構いません。
使用できます(http://nl.mathworks.com/help/matlab/ref/besselj.html、http://nl.mathworks.com/help/optim/ug/fsolve.html)
x = fsolve(@(x)(besselj(0,a*x)-b),x0);
または、以下のコメントで指摘されているように、より単純なソルバーを使用します ( http://nl.mathworks.com/help/optim/ug/fzero.html )
x = fzero(@(x)(besselj(0,a*x)-b),x0);