24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

ConAntares

捐助贵宾 (初入文坛)

混沌的学渣

应助: 0 (幼儿园)
金币: 1.1
散金: 10
帖子: 24
在线: 7.5小时
虫号: 3894049
注册: 2015-05-24
专业: 原子和分子物理

[求助] Matplotlib 三维绘图问题已有1人参与

首先，我用Fortran程序计算出我想要的数据，为三列，保存在一个 .dat的文件里。
格式如下：（x, y, z代表各自列上的数据，z和y之间存在关系 z=f(y)）
x1 y1 z1=f(y1)
x1 y2 z2=f(y2)
x1 y3 z3=f(y3)
……
x2 y1 z'1=f'(y1)
x2 y2 z'2=f'(y2)
x2 y3 z'3=f'(y3)
……
……
x2 y1 z''1=f'(y1)
x2 y2 z''2=f'(y2)
x2 y3 z''3=f'(y3)
……
大致意思是x每做一个循环，都算一遍 z=f(y) 的函数图像，其中x和f是有关系的。

而且y轴很长，我想只画出前一的部分。

想画出如 mat.jpg 所示的图

后来我的Matplotlib语言如下：

CODE:

import matplotlib

import mpl_toolkits.mplot3d

import numpy as np

import pylab

from matplotlib import cm

from matplotlib import rc

from mpl_toolkits.mplot3d import axes3d

fig = plt.figure()

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

#X, Y, Z = axes3d.get_test_data(0.05)

with open('CurrentLTVC.dat', 'r') as fig:

    x = []

    y = []

    z = []

    i = 0

    for line in fig:

        data = line.split()

        x.append((float(data[0])))

        y.append((float(data[1])))

        z.append((float(data[2])))

ax.plot_wireframe(x, y, z, rstride=100, cstride=100, color='red')

plt.ylim(0.1, 20)

plt.show()

但是绘出图像: figure_1.png 的效果却不尽人意，不仅没法控制绘出图像的范围，而且还有连接处的斜线。
求教如何修改成我想要的结果。

mat.jpg

figure_1.png

回复此楼

» 猜你喜欢

有没有人能给点建议已经有5人回复
假如你的研究生提出不合理要求已经有12人回复
实验室接单子已经有7人回复
全日制（定向）博士已经有5人回复
萌生出自己或许不适合搞科研的想法，现在跑or等等看？已经有4人回复
Materials Today Chemistry审稿周期已经有4人回复
参与限项已经有3人回复
对氯苯硼酸纯化已经有3人回复
所感已经有4人回复
要不要辞职读博？已经有7人回复

» 本主题相关价值贴推荐，对您同样有帮助:

安装配置Python+Numpy+Scipy+Matplotlib环境（win7）已经有20人回复
python/matplotlib 图例如何放在外面？已经有8人回复
python/matplotlib 如何在默认目录下打开上一级目录的文件？已经有5人回复
【新书】Python Data Analysis（作者 Ivan Idris）ISBN: 9781783553358 已经有65人回复
mayavi以及matplotlib 已经有1人回复
推荐几个vasp处理计算结果的python脚本已经有66人回复
急急急，在ReadHat上安装phononpy-1.5.1,numpy不成功，感激！已经有11人回复
2012版机器学习的好书 Machine Learning in Action【转载】已经有63人回复
有没有人用过phonoy，老师把任务交给我，各位大神救命啊! 已经有4人回复
python之蟒开启理论计算化学的新时代已经有22人回复
用phonopy计算声子谱已得到FORCE_SETS，之后我怎么得到声子谱？已经有7人回复
QECalc 已经有1人回复
【原创】python 写的计算 PCA 的小脚本已经有4人回复
【转帖】强烈推荐使用phonopy计算声子谱已经有44人回复
【原创】科学计算与绘图--python解决方案已经有27人回复

我很好玩哦，不信你玩玩看~

1楼 2017-04-07 04:13:51

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

ConAntares

捐助贵宾 (初入文坛)

混沌的学渣

应助: 0 (幼儿园)
金币: 1.1
散金: 10
帖子: 24
在线: 7.5小时
虫号: 3894049
注册: 2015-05-24
专业: 原子和分子物理

为了帮助有同样问题的人，我把绘图程序代码贴上上来了。

处理数据文件的 bat 执行程序是我上网找的，其实我也不会编，好在资源很丰富。

我的程序是 Python3.5 中的 Matplotlib2.0 数据可视化程序包

CODE:

################################################################ ####    A 3D Drawing Program ####    Version of April 11, 2017 ####    Lu Niu ################################################################  import matplotlib.pyplot as plt import numpy as np  from matplotlib import pyplot as plt from matplotlib import rc from mpl_toolkits.mplot3d import Axes3D  fig = plt.figure(figsize=(8, 6)) ax = Axes3D(fig)  plt.rc('font', family='serif') ax.text2D(0.20, 0.72, 'Left Current', fontsize = 16, transform=ax.transAxes) ax.text2D(0.22, 0.68, '$\\tau_s=0.1ps$', fontsize = 16, transform=ax.transAxes) #### J10 with open('G-E01-Ts001-J10-e01-P10-TVC/CurrentLTVC00.dat', 'r') as fig:     x = 0.0     y = []     z = []     for line in fig:         data = line.split()         y.append((float(data[1])))         z.append((float(data[2]))) ax.plot_wireframe(x, y, z, rstride=100, cstride=100, color='#D62728', lw='1.6', antialiased=True) with open('G-E01-Ts001-J10-e01-P10-TVC/CurrentLTVC02.dat', 'r') as fig:     x = 0.2     y = []     z = []     for line in fig:         data = line.split()         y.append((float(data[1])))         z.append((float(data[2]))) ax.plot_wireframe(x, y, z, rstride=100, cstride=100, color='#D62728', lw='1.6', antialiased=True) with open('G-E01-Ts001-J10-e01-P10-TVC/CurrentLTVC04.dat', 'r') as fig:     x = 0.4     y = []     z = []     for line in fig:         data = line.split()         y.append((float(data[1])))         z.append((float(data[2]))) ax.plot_wireframe(x, y, z, rstride=100, cstride=100, color='#D62728', lw='1.6', antialiased=True) with open('G-E01-Ts001-J10-e01-P10-TVC/CurrentLTVC06.dat', 'r') as fig:     x = 0.6     y = []     z = []     for line in fig:         data = line.split()         y.append((float(data[1])))         z.append((float(data[2]))) ax.plot_wireframe(x, y, z, rstride=100, cstride=100, color='#D62728', lw='1.6', antialiased=True) with open('G-E01-Ts001-J10-e01-P10-TVC/CurrentLTVC08.dat', 'r') as fig:     x = 0.8     y = []     z = []     for line in fig:         data = line.split()         y.append((float(data[1])))         z.append((float(data[2]))) ax.plot_wireframe(x, y, z, rstride=100, cstride=100, color='#D62728', lw='1.6', antialiased=True) with open('G-E01-Ts001-J10-e01-P10-TVC/CurrentLTVC10.dat', 'r') as fig:     x = 1.0     y = []     z = []     for line in fig:         data = line.split()         y.append((float(data[1])))         z.append((float(data[2]))) ax.plot_wireframe(x, y, z, rstride=100, cstride=100, color='#D62728', lw='1.6', antialiased=True)  #### J01 with open('G-E01-Ts001-J01-e01-P10-TVC/CurrentLTVC00.dat', 'r') as fig:     x = 0.0     y = []     z = []     for line in fig:         data = line.split()         y.append((float(data[1])))         z.append((float(data[2]))) ax.plot_wireframe(x, y, z, rstride=100, cstride=100, color='#FF7F0E', lw='1.6', antialiased=True) with open('G-E01-Ts001-J01-e01-P10-TVC/CurrentLTVC02.dat', 'r') as fig:     x = 0.2     y = []     z = []     for line in fig:         data = line.split()         y.append((float(data[1])))         z.append((float(data[2]))) ax.plot_wireframe(x, y, z, rstride=100, cstride=100, color='#FF7F0E', lw='1.6', antialiased=True) with open('G-E01-Ts001-J01-e01-P10-TVC/CurrentLTVC04.dat', 'r') as fig:     x = 0.4     y = []     z = []     for line in fig:         data = line.split()         y.append((float(data[1])))         z.append((float(data[2]))) ax.plot_wireframe(x, y, z, rstride=100, cstride=100, color='#FF7F0E', lw='1.6', antialiased=True) with open('G-E01-Ts001-J01-e01-P10-TVC/CurrentLTVC06.dat', 'r') as fig:     x = 0.6     y = []     z = []     for line in fig:         data = line.split()         y.append((float(data[1])))         z.append((float(data[2]))) ax.plot_wireframe(x, y, z, rstride=100, cstride=100, color='#FF7F0E', lw='1.6', antialiased=True) with open('G-E01-Ts001-J01-e01-P10-TVC/CurrentLTVC08.dat', 'r') as fig:     x = 0.8     y = []     z = []     for line in fig:         data = line.split()         y.append((float(data[1])))         z.append((float(data[2]))) ax.plot_wireframe(x, y, z, rstride=100, cstride=100, color='#FF7F0E', lw='1.6', antialiased=True) with open('G-E01-Ts001-J01-e01-P10-TVC/CurrentLTVC10.dat', 'r') as fig:     x = 1.0     y = []     z = []     for line in fig:         data = line.split()         y.append((float(data[1])))         z.append((float(data[2]))) ax.plot_wireframe(x, y, z, rstride=100, cstride=100, color='#FF7F0E', lw='1.6', antialiased=True)  #### J00 with open('G-E01-Ts001-J00-e01-P10-TVC/CurrentLTVC00.dat', 'r') as fig:     x = 0.0     y = []     z = []     for line in fig:         data = line.split()         y.append((float(data[1])))         z.append((float(data[2]))) ax.plot_wireframe(x, y, z, rstride=100, cstride=100, color='#1F77B4', lw='1.6', antialiased=True) with open('G-E01-Ts001-J00-e01-P10-TVC/CurrentLTVC02.dat', 'r') as fig:     x = 0.2     y = []     z = []     for line in fig:         data = line.split()         y.append((float(data[1])))         z.append((float(data[2]))) ax.plot_wireframe(x, y, z, rstride=100, cstride=100, color='#1F77B4', lw='1.6', antialiased=True) with open('G-E01-Ts001-J00-e01-P10-TVC/CurrentLTVC04.dat', 'r') as fig:     x = 0.4     y = []     z = []     for line in fig:         data = line.split()         y.append((float(data[1])))         z.append((float(data[2]))) ax.plot_wireframe(x, y, z, rstride=100, cstride=100, color='#1F77B4', lw='1.6', antialiased=True) with open('G-E01-Ts001-J00-e01-P10-TVC/CurrentLTVC06.dat', 'r') as fig:     x = 0.6     y = []     z = []     for line in fig:         data = line.split()         y.append((float(data[1])))         z.append((float(data[2]))) ax.plot_wireframe(x, y, z, rstride=100, cstride=100, color='#1F77B4', lw='1.6', antialiased=True) with open('G-E01-Ts001-J00-e01-P10-TVC/CurrentLTVC08.dat', 'r') as fig:     x = 0.8     y = []     z = []     for line in fig:         data = line.split()         y.append((float(data[1])))         z.append((float(data[2]))) ax.plot_wireframe(x, y, z, rstride=100, cstride=100, color='#1F77B4', lw='1.6', antialiased=True) with open('G-E01-Ts001-J00-e01-P10-TVC/CurrentLTVC10.dat', 'r') as fig:     x = 1.0     y = []     z = []     for line in fig:         data = line.split()         y.append((float(data[1])))         z.append((float(data[2]))) ax.plot_wireframe(x, y, z, rstride=100, cstride=100, color='#1F77B4', lw='1.6', antialiased=True)  ax.view_init(30, 30) ax.set_xlim3d(0, 1) ax.set_ylim3d(0, 4) ax.set_zlim3d(-0.10, 0.25) ax.set_xlabel(r'Voltage $(V)$', family='serif', fontsize=14) ax.set_ylabel(r'Time $(ps)$', family='serif', fontsize=14) ax.set_zlabel(r'Current $(nA)$', family='serif', fontsize=14)  plt.xticks([0, 0.2, 0.4, 0.6, 0.8, 1.0], fontsize=12) plt.yticks([0.0, 0.5, 1.0, 1,5, 2.0, 2.5, 3.0, 3.5, 4.0], fontsize=12) ax.set_zticks([-0.10,-0.05,0,0.05,0.10,0.15,0.20,0.25]) #ax.set_zticks(np.arange(0,6,1.0)) plt.gca().invert_xaxis()  plt.savefig('GV-Ts001-L.png', dpi=512) plt.savefig('GV-Ts001-L.eps', dpi=512) plt.show()