我画了一条椭圆曲线。我想沿着P,Q,R
画一条线(其中P
和Q
将独立于此问题确定)。P
的主要问题是sympy solve()
返回另一个方程,它需要返回一个值,以便用于绘制P
的x值。据我所知,solve()
应该返回一个值,所以我显然在这里做了一些错误的事情,只是我完全没有看到。作为参考,下面是P+Q=R
的外观:
我已经检查了文档和other material,这是我能让自己陷入麻烦的最远距离:
from mpl_toolkits.axes_grid.axislines import SubplotZero
from pylab import *
import numpy as np
import matplotlib.pyplot as plt
from matplotlib.path import Path
import matplotlib.patches as patches
from matplotlib import rc
import random
from sympy.solvers import solve
from sympy import *
def plotGraph():
fig = plt.figure(1)
#ax = SubplotZero(fig, 111)
#fig.add_subplot(ax)
#for direction in ["xzero", "yzero"]:
#ax.axis[direction].set_axisline_style("-|>")
#ax.axis[direction].set_visible(True)
#ax.axis([-10,10,-10,10])
a = -2; b = 1
y, x = np.ogrid[-10:10:100j, -10:10:100j]
xlist = x.ravel(); ylist = y.ravel()
elliptic_curve = pow(y, 2) - pow(x, 3) - x * a - b
plt.contour(xlist, ylist, elliptic_curve, [0])
#rand = random.uniform(-5,5)
randmid = random.randint(30,70)
#y = ylist[randmid]; x = xlist[randmid]
xsym, ysym = symbols('x ylist[randmid]')
x_result = solve(pow(ysym, 2) - pow(xsym, 3) - xsym * a - b, xsym) # 11/5/13 needs to return a value
plt.plot([-1.5,5], [-1,8], color = "c", linewidth=1) # plot([x1,x2,x3,...],[y1,y2,y3,...])
plt.plot([xlist[randmid],5], [ylist[randmid],8], color = "m", linewidth=1)
#rc('text', usetex=True)
text(-9,6,' size of xlist: %s \n size of ylist: %s \n x_coord: %s \n random_y: %s'
%(len(xlist),len(ylist),x_result,ylist[randmid]),
fontsize=10, color = 'blue',bbox=dict(facecolor='tan', alpha=0.5))
plt.annotate('$P+Q=R$', xy=(2, 1), xytext=(3, 1.5),arrowprops=dict(facecolor='black', shrink=0.05))
## verts = [(-5, -10),(5, 10)] # [(x,y)startpoint,(x,y)endpoint] #,(0, 0)]
## codes = [Path.MOVETO,Path.LINETO] # related to verts[] #,Path.STOP]
## path = Path(verts, codes)
## patch = patches.PathPatch(path, facecolor='none', lw=2)
## ax.add_patch(patch)
plt.grid(True)
plt.show()
def main():
plotGraph()
if __name__ == '__main__':
main()
最后,我想画一条线来显示P+Q=R
,所以如果有人还想增加一些关于如何编码得到Q
的东西,那将是非常感谢的。我正在自学Python和椭圆曲线,所以我相信任何入门级程序员都可以在2分钟内弄清楚我已经学了一段时间了。
2条答案
按热度按时间x0fgdtte1#
我不知道你在计算什么,但下面是可以绘制图形的代码:
图表:
djmepvbi2#
根据HYRY的回答,我只是更新了一些细节,使其更好:
我们得到了这个多边形的根:[3. 2.4444444 0. ] [5. 3.7037037 2. ]