math - Octave / Matlab: 多項式の根をプロットする方法

Question

多項式の根をプロットしようとしていますが、取得できません。

まず、多項式を作成します

p5 = [1 0 0 0 0 -1] %x^5 - 1
r5 = roots(p5)
stem (p5)

関数を使用していstemますが、茎を取り除き、根の周りに円を描きたいと思います。

これは可能ですか、ステムは正しいコマンドですか?

前もって感謝します、

PS: これは宿題ではありませんが、非常に近いので、リクエストがあればタグ付けします。

Answer 1

x 軸に実数部、y 軸に虚数部をプロットする複雑な根がある場合は、PLOT関数を使用できます。

plot(r5,'o');

関数とルートを一緒にプロットしたい場合は、複雑なルートを無視する必要があります (以下のコメントで yuk が言及しているように)。

p5 = [1 0 0 0 0 -1];
r5 = roots(p5);
realRoots = r5(isreal(r5));  %# Gets just the real roots
x = -2:0.01:2;               %# x values for the plot
plot(x,polyval(p5,x));       %# Evaluate the polynomial and plot it
hold on;                     %# Add to the existing plot
plot(realRoots,zeros(size(realRoots)),'o');  %# Plot circles for the roots