0

我画了 4D 情节。

我的情节有一些部分（表面）。

我在不同的数组中生成每个部分（每个表面）。

图的最后一部分（最后一个表面）有两个部分。

1个在上面，第二个在后面。

但在图中pylab这是不正确的。

最后一个表面始终位于图表的顶部。（如果你旋转mat了，那么表面会显示在顶部）

混沌图 在此处输入图像描述

EDIT： 编码：

#!/usr/bin/python3
from scipy.integrate import odeint
from numpy import var
import numpy as np
from mpl_toolkits.mplot3d import Axes3D
from matplotlib.ticker import LinearLocator, FormatStrFormatter
import matplotlib.pyplot as plt
import matplotlib.cm as cm
import time as tim
import matplotlib as mpl
epsel=np.array(1.*10**(-5),float)
ru=np.array(.264,float) #free parameter
ddd=100
ddd2=(ddd+1)**2
bi =  np.zeros((4,ddd2,4),float)
bbb=1.
def boundary1(ru,ddd=100):#x=y
    namnam=0
    for zz in np.linspace(-bbb,bbb+1.,ddd):
        for tt in np.linspace(-bbb,bbb,ddd):
            bi[0,namnam,:]=[tt+epsel,tt+epsel,zz+epsel,0]
            namnam=namnam+1
    return(namnam)
def boundary2(ru,ddd=100):#x=0############################################
    namnam=0
    for zz in np.linspace(-bbb,bbb+1,ddd):
        for yy in np.linspace(-bbb,bbb,ddd):
            bi[1,namnam,:]=[0.+epsel,yy+epsel,zz+epsel,1]
            namnam=namnam+1
    return(namnam)
def boundary3(ru,ddd=100):#z=0
    namnam=0
    for xx in np.linspace(-bbb,bbb,ddd):
        for yy in np.linspace(-bbb,bbb,ddd):
            bi[2,namnam,:]=[xx+epsel,yy+epsel,0.+epsel,2]
            namnam=namnam+1
    return(namnam)
def boundary4(ru,ddd=100):
    namnam=0
    for zz in np.linspace(-bbb,bbb+1.,ddd):
        for tt in np.linspace(ru**(.5),bbb,ddd/4):
            bi[3,namnam,:]=[tt+epsel,1./tt+epsel,zz+epsel,3]
            namnam=namnam+1
        for tt in np.linspace(ru**(.5),bbb,ddd/4):
            bi[3,namnam,:]=[1./tt+epsel,tt+epsel,zz+epsel,3]
            namnam=namnam+1
        for tt in np.linspace(ru**(.5),bbb,ddd/4):
            bi[3,namnam,:]=[-tt+epsel,-1./tt+epsel,zz+epsel,3]
            namnam=namnam+1
        for tt in np.linspace(ru**(.5),bbb,ddd/4):
            bi[3,namnam,:]=[-1./tt+epsel,-tt+epsel,zz+epsel,3]
            namnam=namnam+1
    return(0)
def bound(ru,ddd=100):
    bbb=2.*ru
    boundary4(ru,ddd=100)
    boundary3(ru,ddd=100)
    boundary2(ru,ddd=100)
    boundary1(ru,ddd=100)
    return(0)

bound(ru,ddd=100)
fig = plt.figure()
ax = fig.add_subplot(111, projection='3d')
for i in range(4):
    ax.scatter(bi[i,:,0],bi[i,:,1],bi[i,:,2],c=bi[i,:,3], cmap=plt.hot())
plt.show()
4

0 回答 0