如何在给定区间内求多项式的根

Question

我想写一个函数，它的输入是一个多项式(p)和一个范围(a，b)，它给出了多项式在这个范围内的根的个数。但是我不知道如何在函数中设置范围。

x=poly(0,"x")
a = 0;
b = 5;
p=-6+11*x-6*x^2+x^3

function [n_roots] = numroots (p,a,b)
    rootn = roots(p)
//todo
    n_roots = length(rootn)
endfunction

下一步是在图中表示多项式和根。

Answer 1

也许你应该在计算完根后选择哪些根在这个范围内。请注意，如果要将根分类到一个范围内，则根应该是实数。

function [n_roots] = numroots(p,a,b)
    rootn = roots(p)
    rootn = rootn(find(real(rootn) > a & real(rootn) < b))
    n_roots = length(rootn)
endfunction