版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

>虫友互识 (5764)
>论坛更新日志 (3287)
>留学DIY (3103)
>留学生活 (2461)
>海外博后 (1486)
>育儿交流 (1200)
>海归之家 (1000)
>医药类考试 (1000)
>材料综合 (799)
>人文社科 (798)
>信息科学 (760)
>复合材料 (658)
>地学 (640)
>文学芳草园 (601)
>转基因 (600)
>我与小木虫的故事 (600)

返回列表

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

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

回复此楼

» 猜你喜欢

哪位老哥知道今年的国自然具体哪一天放榜？已经有6人回复
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急已经有6人回复
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急已经有7人回复
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急已经有5人回复
欢迎发来filecode的Mz6后的代码验证其规律已经有50人回复
filecode 已经有11人回复
售SCI文章，我:8O.5.5.1O.54,科目全,可十急已经有5人回复
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急已经有5人回复
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急已经有4人回复
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急已经有4人回复

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

安装配置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
专业: 原子和分子物理

★ ★ ★ ★ ★
jjdg: 金币+5, 感谢分享经验 2017-04-12 00:20:02

楼上怎么出乱码了！？！？

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()

赞一下(1人)

回复此楼

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

10楼2017-04-11 21:27:23

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

查看全部 10 个回答

ConAntares

捐助贵宾 (初入文坛)

混沌的学渣

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

链接处的线该如何消除？

我做了差不多两百个这样的数据，如果再做一次那工程量我实在是承受不起。

附件里面有示范数据，求好心人帮帮忙，老司机带带我。

TvT

figure_2.png

赞一下

回复此楼

» 本帖附件资源列表

欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
附件 1 : CurrentLTVC.dat

2017-04-07 04:19:15, 5.83 M

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

2楼2017-04-07 04:19:54

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

ConAntares

捐助贵宾 (初入文坛)

混沌的学渣

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

再补充一下，我的Python版本是3.5，Matplotlib是最新的版本，希望有大佬带我飞一波。
我的论文其他配图都用Matplotlib做好了，就差几个三维的数据可视化图形了。

赞一下

回复此楼

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

3楼2017-04-07 05:01:59

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

约束陶瓷装甲

金虫 (正式写手)

应助: 2 (幼儿园)
金币: 928.6
红花: 3
帖子: 305
在线: 197.4小时
虫号: 4027532
注册: 2015-08-18
专业: 医学图像数据处理与分析

★
jjdg: 金币+1, 感谢参与 2017-04-12 00:18:41

x=x1,x2,x3,...,xn要用循环分成n次来plot，而不是一次plot，一次plot的结果就是连在一起。

[ 发自手机版 http://muchong.com/3g ]

赞一下(2人)

回复此楼

4楼2017-04-07 05:38:50

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

查看全部 10 个回答

最具人气热帖推荐 [查看全部]		作者	回/看	最后发表

[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜？ +6	Ldrop2023 2026-08-13	6/300	2026-08-15 10:09 by tuanggou
[公派出国] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3	k0dTPqJtl0jt 2026-08-14	5/250	2026-08-15 07:09 by 4wMiSEwB6436
[公派出国] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3	7lpolszZVXgi 2026-08-14	5/250	2026-08-15 04:41 by 4wMiSEwB6436
[考研] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3	HFw0lei2R37i 2026-08-14	7/350	2026-08-15 00:40 by 4wMiSEwB6436
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +4	HFw0lei2R37i 2026-08-14	8/400	2026-08-15 00:16 by 4wMiSEwB6436
[论文投稿] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3	HFw0lei2R37i 2026-08-14	5/250	2026-08-15 00:16 by 4wMiSEwB6436
[基金申请] 奇怪，两个人的filecode固定段从头到尾一模一样 +8	布布和一二 2026-08-10	11/550	2026-08-14 14:58 by Equinoxhua
[基金申请] 有时候，自然基金真的不能太认真（我的申报经验） +6	majunge000 2026-08-11	8/400	2026-08-14 14:44 by 医学老男孩
[基金申请] FileCode能看出啥？ +10	要乐观耀哥 2026-08-10	32/1600	2026-08-14 09:37 by 要乐观耀哥
[基金申请] 我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人 +10	家与远方 2026-08-10	15/750	2026-08-14 02:08 by 绵羊哥哥
[基金申请] 关于Filecode分析方法 +9	majunge000 2026-08-10	12/600	2026-08-13 23:20 by iwuli
[文学芳草园] 阿姨 +4	汪汪锅 2026-08-09	4/200	2026-08-13 19:43 by arzu_hma
[基金申请] 静等基金结果 +8	gjjjzhong 2026-08-10	21/1050	2026-08-13 17:56 by 且听虎啸
[基金申请] 2019年青年基金涵评意见，大家看看几个A，几个B? +11	Tide man 2026-08-11	11/550	2026-08-13 07:35 by 撸猫猫
[基金申请] 为什么网上很多人说本周 12号出结果 +6	瞬息宇宙 2026-08-10	7/350	2026-08-11 19:25 by Tide man
[基金申请] 什么时候出结果，有咨询渠道？？？ +3	Tide man 2026-08-11	3/150	2026-08-11 17:54 by kudofaye
[基金申请] 确定了，国自然21号放榜 +6	布布和一二 2026-08-10	7/350	2026-08-10 19:15 by 2000zf36392
[基金申请] 国自然结果 +4	Vierhys 2026-08-10	8/400	2026-08-10 15:06 by Vierhys
[基金申请] 据悉今年马上要出结果了 +7	瞬息宇宙 2026-08-10	8/400	2026-08-10 12:42 by Vivilian
[基金申请] 面上项目filecode邪修 +5	西山十月 2026-08-09	7/350	2026-08-10 07:32 by 仁砚薪传

24小时热门版块排行榜

ConAntares

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

» 猜你喜欢

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

ConAntares

ConAntares

» 本帖附件资源列表

ConAntares

约束陶瓷装甲

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