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

问题描述

我有以下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).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)
    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,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)
    appeal = np.sqrt(mx**2+my**2)
    return appeal

def calculateBestLinear(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,values)
    appeal = 1 - np.sqrt((mx-1)**2+my**2)/2
    return appeal

testSeq = [0,testSeq)
agcd = sm.find_local_maximum(gcd,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,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,100)

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

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

解决方法

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

当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?

attributes python sage