All zeros for Bessel function
이전 댓글 표시


I dont want to find all zeros of this function by using 'bessel function' besselj(nu, Z) in matlab toolbox. I would like to create a matlab function to calculate all roots. I create a function on newton iteration method, but it calculated just one root. thanks for your helps in advance.
댓글 수: 10
David Goodmanson
2020년 2월 29일
Hi Zeynep,
J_1/2(x) equals zero if and only if the right hand side of that equation is zero. So you should consider when the r.h.s. is zero.
Zeynep Toprak
2020년 2월 29일
편집: Zeynep Toprak
2020년 2월 29일
Walter Roberson
2020년 2월 29일
You examine the formula and see that the first term is 0 only if x is infinite, which cannot occur for that range. You then examine the second term and see that it is zero when sin(x) is 0, which happens exactly every π/2 . No need to consult a graph or use Newton's method.
David Goodmanson
2020년 2월 29일
편집: David Goodmanson
2020년 2월 29일
Hi Walter, you meant to say every pi
Walter Roberson
2020년 2월 29일
You are right, should be every π
Zeynep Toprak
2020년 2월 29일
Walter Roberson
2020년 2월 29일
That is one way; there are other approaches.
Zeynep Toprak
2020년 2월 29일
편집: Zeynep Toprak
2020년 2월 29일
Walter Roberson
2020년 2월 29일
https://en.m.wikipedia.org/wiki/Root-finding_algorithm
There are also techniques that involve splitting up the interval into a number of subintervals and running a Newton type algorithm on the entire vector of starting points, and then at the end taking the unique values (taking into account round-off error). This approach only really works if you have information about the minium separation of the zeros.
Zeynep Toprak
2020년 3월 1일
답변 (0개)
카테고리
도움말 센터 및 File Exchange에서 Mathematics에 대해 자세히 알아보기
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!