尝试在 Python 中运行使用 Sage Math 的代码时出现 NameError 和 AttributeError

Question

我有以下 SageMath 代码，在 CoCalc 中完美运行：

def mean_x(factor, values):
    return sum([cos(2*pi*v/factor) for v in values])/len(values)

def mean_y(factor, values):
    return sum([sin(2*pi*v/factor) for v in values])/len(values)

def calculatePeriodAppeal(factor, values):
    mx = mean_x(factor, values)
    my = mean_y(factor, values)
    appeal = sqrt(mx^2+my^2)
    return appeal

def calculateBestLinear(factor, values):
    mx = mean_x(factor, values).n()
    my = mean_y(factor, values).n()
    y0 = factor*atan2(my,mx)/(2*pi).n()
    err = 1-sqrt(mx^2+my^2).n()
    return [factor*x + y0, err]

def calculateGCDAppeal(factor, values):
    mx = mean_x(factor, values)
    my = mean_y(factor, values)
    appeal = 1 - sqrt((mx-1)^2+my^2)/2
    return appeal

testSeq = [0,125,211,287,408,520,650,735,816,942,1060] 
gcd = calculateGCDAppeal(x, testSeq)
agcd = find_local_maximum(gcd,2,100)
print(agcd)
plot(gcd,(x, 2, 100))

输出是来自testSeq 的数字的最佳近似最大公约数，以及一个绘图。

如何在 Python 中使用此代码？

这是当前的 Python 版本，它还不能工作：

import numpy as np
import sage as sm

def mean_x(factor, values):
    return sum([np.cos(2*np.pi*v/factor) for v in values])/len(values)

def mean_y(factor, values):
    return sum([np.sin(2*np.pi*v/factor) for v in values])/len(values)

def calculatePeriodAppeal(factor, values):
    mx = mean_x(factor, values)
    my = mean_y(factor, values)
    appeal = np.sqrt(mx**2+my**2)
    return appeal

def calculateBestLinear(factor, values):
    mx = mean_x(factor, values).n()
    my = mean_y(factor, values).n()
    y0 = factor*np.atan2(my,mx)/(2*np.pi).n()
    err = 1-np.sqrt(mx**2+my**2).n()
    return [factor*x + y0, err]

def calculateGCDAppeal(factor, values):
    mx = mean_x(factor, values)
    my = mean_y(factor, values)
    appeal = 1 - np.sqrt((mx-1)**2+my**2)/2
    return appeal

testSeq = [0,125,211,287,408,520,650,735,816,942,1060] 
gcd = calculateGCDAppeal(x, testSeq)
agcd = sm.find_local_maximum(gcd,2,100)
print(agcd)

我得到的错误是：

Traceback (most recent call last):

  File "<ipython-input-805-d2a8b405fd43>", line 30, in <module>
    gcd = calculateGCDAppeal(x, testSeq)

NameError: name 'x' is not defined

我不明白，因为这段代码在 CoCalc 中有效。

Traceback (most recent call last):
    
      File "<ipython-input-803-80cdeb0485a5>", line 33, in <module>
        agcd = sm.find_local_maximum(gcd,2,100)
    
    AttributeError: module 'sage' has no attribute 'find_local_maximum'

...但我知道 SageMath 有函数find_local_maximum。

如果我改用numerical.optimize.find_local_maximum，我会得到：

Traceback (most recent call last):

  File "<ipython-input-842-949de8b03df5>", line 34, in <module>
    agcd = sm.numerical.optimize.find_local_maximum(gcd,2,100)

AttributeError: module 'sage' has no attribute 'numerical'

我不知道如何添加“数字”属性。

Answer 1

首先，这段代码中的任何地方都没有 x 的定义。另外find_local_maximum的语法是sage.numerical.optimize.find_local_maximum不是sage.find_local_maximum（来自the SageMath documentation）

Answer 2

当 Sage 启动时，它会导入很多东西。

使用 Python 时，需要自己导入。

特别是，Sage 在启动时定义了x。

如果使用 Python 则导入它：

from sage.calculus.predefined import x

使用import_statements 函数找出答案 在 Sage 会话中。

快速解决方法是万能的

from sage.all import *

但最好只导入需要的内容。

查看此相关问题的答案

• How can I translate this SageMath code to run in Python?