如何更新复根图

2024-09-27 07:33:40 发布

您现在位置:Python中文网/ 问答频道 /正文

我想连续绘制几个雅可比多项式的(复数)根。这些多项式的参数和阶数随尺寸因子n的变化而变化。到目前为止,我已经有了计算根并在复平面上绘制它们的代码,但我希望每个绘图都覆盖上一个绘图,以便查看是否会形成根的模式。你知道吗

我写的代码现在生成了这些图,但它们都是独立的。你知道吗

import sympy
from sympy import *
import numpy
x = sympy.symbols('x')

def show_roots(n,a,b,c):
    for k in range (1,n+1):
        p = sympy.jacobi(a*k,-(a+b)*k,(b+c)*k,x)
        coeff = sympy.poly(p).coeffs()
        roots = numpy.roots(coeff)
        X = [x.real for x in roots]
        Y = [x.imag for x in roots]
        plt.scatter(X,Y, color='red')
        plt.show()
        plt.pause(0.5)

Tags: 代码inimportnumpy绘图for参数show

热门问题