多元函数连续

最新推荐文章于 2024-09-17 23:15:58 发布

jennydale

最新推荐文章于 2024-09-17 23:15:58 发布

阅读量46

点赞数

分类专栏：多元函数立体图像文章标签： python numpy 开发语言

本文链接：https://blog.csdn.net/jennydale/article/details/129768627

版权

多元函数立体图像专栏收录该内容

3 篇文章 0 订阅

订阅专栏

这段代码展示了如何利用Python的numpy和matplotlib库创建3D图形。用户可以输入x轴和y轴的长度以及x和y的幂指数，程序会生成一个基于给定数学模型（z=sqrt(x^n1+y^n2)）的3D图形，并允许自定义坐标轴标签和标题。

摘要由CSDN通过智能技术生成

在这里插入图片描述
下面展示一些 内联代码片。

import numpy as np

import matplotlib.pyplot as plt


fig=plt.figure()

ax = fig.add_subplot(111, projection='3d')

print("Welcometo use 3D graphic generator, the formula model is z=sqrt(x^n1+y^n2).\n")

a1=input("pleaseenter x axis length value(x),actual display range is from - x to x:\n")

a1=float(a1)

print("Thex axis actual display range is from - "+str(a1)+" to"+str(a1)+":\n")

a2=input("pleaseenter y axis length value(y),actual display range is from - y to y:\n")

a2=float(a2)

print("They axis actual display range is from - "+str(a2)+" to"+str(a2)+":\n")

a3=input("pleaseenter n power number of x value:(please enter an integer)\n")

a3=int(a3)

print("The power number of x value is "+str(a3)+".\n")

a4=input("pleaseenter n power number of y value:(please enter an integer)\n")

a4=int(a4)

print("Thepower number of y value is "+str(a4)+".\n")

X=np.arange(-a1,a1,0.5)

Y=np.arange(-a2,a2,0.5)

X,Y=np.meshgrid(X,Y)

epsilon = 1e-15  # 设置一个极小值
Z = (X**2*Y)/(X**2+Y**2+epsilon)
ax.plot_surface(X,Y,Z,rstride=2,cstride=2,alpha=0.5,cmap=plt.cm.rainbow)

x1=input("pleaseenter x label name:\n")

x2=input("pleaseenter y label name:\n")

x3=input("pleaseenter z label name:\n")

x4=input("pleaseenter title label name:\n")

ax.set_xlabel(x1,fontsize=20)

ax.set_ylabel(x2,fontsize=20)

ax.set_zlabel(x3,fontsize=20)

ax.set_title(x4,fontsize=25)

ax.view_init(30,35)

b1=a1+4

b2=a2+4

cset=ax.contour(X,Y, Z, zdir = 'x', offset = -b1, cmap = plt.cm.hot)

cset=ax.contour(X,Y, Z, zdir = 'y', offset = -b2, cmap = plt.cm.hot)

cset=ax.contour(X,Y, Z, zdir = 'z', offset = -7, cmap = plt.cm.hot)

ax.set_xlim(-b1,b1)

ax.set_ylim(-b2,b2)

ax.set_zlim(-8,8)

plt.show()