qq_41660047
2019-10-03 08:52
采纳率: 34.2%
浏览 651

请教各位大神如何用Python画出如下函数图形

(x^2+9/4*y^2+z^2-1)^3-x^2*z^3-9/80*y^2*z^3==0

  • 写回答
  • 好问题 提建议
  • 关注问题
  • 收藏
  • 邀请回答

2条回答 默认 最新

  • threenewbee 2019-10-03 09:16

    你有三个变量,可以在一定区域内,以某个间隔精度,遍历x y,利用根据方程算出z,这个方程式你可以自己转换下。
    然后调用Matplotlib绘三维图
    https://blog.csdn.net/jasonzhoujx/article/details/81780774


    补充一下,数值算法理论上可行,代码如下

    import numpy as np
    import matplotlib.pyplot as plt
    import math
    import mpl_toolkits.mplot3d
    import sympy as sympy
    
    
    def solvez(x, y):
        z = sympy.symbols('z')
        result = 0.0    
        try:
            result = float(sympy.solve((x**2+9/4*y**2+z**2-1)**3-x**2*z**3-9/80*y**2*z**3, z)[0].evalf())
        except TypeError as e:
            pass
        return result
    
    fig = plt.figure()  #定义新的三维坐标轴
    ax3 = plt.axes(projection='3d')
    
    #定义三维数据
    xx = np.arange(-2,2,0.1)
    yy = np.arange(-2,2,0.1)
    X, Y = np.meshgrid(xx, yy)
    
    Z = np.zeros((len(X), len(Y)), dtype=float)
    
    
    for i in range(0,len(X)):
        for j in range(0,len(Y)):
            Z[i][j] = solvez(X[i][j],Y[i][j])
    
    #作图
    ax3.plot_surface(X,Y,Z,cmap='rainbow')
    #ax3.contour(X,Y,Z, zdim='z',offset=-2,cmap='rainbow)   #等高线图,要设置offset,为Z的最小值
    plt.show()
    

    实际上比较困难,在于没有办法知道解的可行范围,算出来区域之外无解都是复数解

    评论
    解决 无用
    打赏 举报
  • dabocaiqq 2019-10-03 09:32
    import numpy as np
    import matplotlib.pyplot as plt
    import math
    import mpl_toolkits.mplot3d
    
    x, y = np.mgrid[-2:2:0.01, -2:2:0.01]
    
    z=根据xy算出z的表达式,根据你的自己可以改写的
    ax = plt.subplot(111, projection='3d') 
    ax.plot_surface(x, y, z, rstride=1, cstride=1, cmap='rainbow', alpha=0.9)#绘面
    ax.set_xlabel('x')
    ax.set_ylabel('y')
    ax.set_zlabel('z')
    plt.show()
    
    评论
    解决 无用
    打赏 举报