当我运行失败的热方程代码时，我在FEniCS中遇到编码错误，error=np.abs(u_e.vector().array()-u.vector().array())AttributeError:'dolfin.cpp.la.PETScVector'对象没有属性'array'我相信这个想法是计算u和每个时间级别的精确解之间的差异，并将它们放入数组中，然后我们将取最大值。(如果不是，请纠正我!)计算误差的原始代码是u_e=interpolate(u_D,V)error=np.abs(u_e.vector().array()-u.vector().array()).max()print('
我有一个字符串，它是一个数学方程式，但带有一些自定义函数。我需要找到所有这些函数并用一些代码替换它们。例如我有一个字符串:a+b+f1(f2(x,y),x)我想要将f2(x,y)替换为x+y^2和f1(x,y)的代码code>与sin(x+y)。如果像示例中一样支持嵌套函数，那就太理想了。但是，如果不支持嵌套，它仍然有用。据我从类似主题中了解到，这可以使用编译器模块(如compiler.parse(eq))来完成。我如何使用compiler.parse(eq)创建的AST对象重建我的字符串，替换所有找到的函数？我只需要执行替换，然后字符串就会在其他程序中使用。不需要评估。
基于我发现的一个以前的例子here，我试图找出如何生成与我需要坚持的特定形式相对应的Fortran代码。所需的FORTRAN代码如下所示(itisbasedontheFitzHugh–Nagumomodel):SUBROUTINEFF(NE,U,PAR,F)!------------!DefinethenonlineartermIMPLICITNONEINTEGER,INTENT(IN)::NEDOUBLEPRECISION,INTENT(IN)::U(NE),PAR(*)DOUBLEPRECISION,INTENT(OUT)::F(NE)DOUBLEPRECISIONu,v,e,a1
有没有我可以在网上找到专门针对python的二分法？例如，给定这些方程，我如何使用二分法求解它们？x^3=93*x^3+x^2=x+5cos^2x+6=x 最佳答案 使用scipy.optimize.bisect:importscipy.optimizeasoptimizeimportnumpyasnpdeffunc(x):returnnp.cos(x)**2+6-x#0optimize.bisect调用_zeros._bisect，它是用C实现的。 关于python-使用二分法求解方程
以下公式用于对二维空间中的点进行分类:f(x1,x2)=np.sign(x1^2+x2^2-.6)所有点都在空间X=[-1,1]x[-1,1]中，每个点都有均匀的概率被选中。现在我想可视化等于的圆:0=x1^2+x2^2-.6x1的值应该在x轴上，x2的值应该在y轴上。这一定是可能的，但我很难将方程式转换为绘图。 最佳答案 您可以使用等高线图，如下所示(基于http://matplotlib.org/examples/pylab_examples/contour_demo.html中的示例):importnumpyasnpimpor
2.2.2一维热传导方程热传导方程是描述热量在介质中传导的数学模型。在许多实际应用中，我们需要预测物体的温度随时间和空间的演化情况，这就需要用到热传导方程。热传导方程的背景可以追溯到18世纪，当时科学家们对热的本质和热量如何传递产生了浓厚的兴趣。傅里叶在他的《热理论》一书中，提出了一种新的方法，可以用一组三角函数来表示任何连续函数。这种方法后来被称为傅里叶级数。傅里叶发现，可以用傅里叶级数展开热传导问题的解，从而得到了热传导方程。热传导在三维的等方向均匀介质里的传播可用以下方程表达:∂u∂t=div(Uu)=k
我想开发一个GUI应用程序来显示给定的数学方程式。当您单击方程式中的特定变量以表示它是未知变量(即要计算的变量)时，方程式会自行转换以评估所需的未知变量。例如:a=(b+c*d)/e让我们假设我点击“d”表示它是未知变量。那么等式应该重组为:d=(a*e-b)/c到目前为止，我只想知道如何根据用户输入重新排列给定的方程式。我从我兄弟那里得到的一个建议是在后端使用前缀/后缀符号表示来评估它。这是唯一的方法还是有更简单的建议？此外，我不仅会使用基本的数学函数，还会使用三角函数和微积分(我认为是基本的。没有偏微积分等等)。我认为前/后修复符号评估可能对评估更高的数学函数没有帮助。这只是我个人
我有一个类Point，由一个具有x和y坐标的点组成，我必须编写一个方法来计算并返回连接Point对象和作为参数传递的另一个Point对象的直线方程(my_point.get_straight_line(my_point2)。我知道如何在纸上用y-y1=m(x-x1)计算，我已经有一个方法my_point.slope(my_point2)来计算m，但我不能我真的全神贯注于如何将等式转换为Python。这是整个类(class):classPoint:def__init__(self,initx,inity):self.x=initxself.y=initydefgetx(self):ret
我有一个类Point，由一个具有x和y坐标的点组成，我必须编写一个方法来计算并返回连接Point对象和作为参数传递的另一个Point对象的直线方程(my_point.get_straight_line(my_point2)。我知道如何在纸上用y-y1=m(x-x1)计算，我已经有一个方法my_point.slope(my_point2)来计算m，但我不能我真的全神贯注于如何将等式转换为Python。这是整个类(class):classPoint:def__init__(self,initx,inity):self.x=initxself.y=initydefgetx(self):ret
是否可以从特定网站导出文本、图像和LaTeX方程式，以便您可以直接自定义您自己的PDF而不会模糊对象？只有图像具有固定分辨率。我知道有几种间接生成PDF的方法。试图在RiemannZetaFunction上呈现来自WolframMathWorld的PDF，例如，可以通过Chrome将其打印并另存为PDF，但是当您放大得更近时，LaTeX方程式和文本自然会变得模糊。我尝试下载“Wolfram的CDF播放器”，但它仅包含Mathematica库的语法-而不是WolframMathWorld提供的有用解释。我需要什么才能提取PDF文件中的文本、图像和LaTeX方程而不使它们变得模糊？
