做动力学模拟分析时,xvg文件画图常用origin,当分析的文件多时,每个xvg文件都需要画图,或者不同模拟需要比较时也通过图比较,再手动画图让人比较焦躁,所以就在网上查询利用python画图。

其中xvg文件经过 grep 命令处理,处理后的disatance/54HN-363010.xvg可以直接被python读取

grep -v '@' 54HN-363010.xvg | grep -v '#' > disatance/54HN-363010.xvg
import matplotlib.pyplot as pltX, Y = [], []for line in open('54HN-363O10.xvg','r'):values = [float(s) for s in line.split()]X.append(values[0])Y.append(values[1])
plt.plot(X, Y, color='k', linestyle='dashed',marker='o')
plt.title("54VAL-UNK distance", fontsize=14)
plt.xlabel("Time (ps)")
plt.ylabel("54VAL-HN_UNK-O10 (nm)")
plt.show()

得到下面这样的图:

import matplotlib.pyplot as pltfilename = '54HN-363O10.xvg'
filename1 = '223HH-363O8.xvg'
X, Y = [], []
X1, Y1 = [], []with open(filename, 'r') as f:lines = f.readlines()for line in lines:value = [float(s) for s in line.split()]X.append(value[0])Y.append(value[1])with open(filename1, 'r') as f1:lines = f1.readlines()for line in lines:value = [float(s) for s in line.split()]X1.append(value[0])Y1.append(value[1])plt.plot(X, Y, color='black', linestyle='dashed', marker='o', label='54VAL-UNK')
plt.plot(X, Y1, color='red', linestyle='dashed', marker='o', label='223HH-UNK')
plt.legend(loc= 'upper right')
plt.xlabel("Time (ps)")
plt.ylabel("Distance (nm)")plt.show()

得到的图是:

将两个xvg进行比较的python脚本比较笨重,还没想出怎么可以读取两个文件,所以就直接粘贴复制前面的命令,执行两次。后面有优化后的命令在补充到里面。

[参考博文]
(https://blog.csdn.net/Mr_Cat123/article/details/78828634)(https://blog.csdn.net/qq_35077107/article/details/97619685)

python_matplotlib画图相关推荐

  1. KNN之约会网站测试

    感谢Jack-Cui大佬的知识分享 机器学习专栏点击这里 目录 感谢Jack-Cui大佬的知识分享 项目背景介绍 案例步骤 1.准备数据:解析数据 代码块 运行结果 相关方法介绍 with 语句 re ...

  2. python把标签图例放图外的下面 自动放在外侧

    参数loc: plt.legend(bbox_to_anchor=(0.5, -0.2),loc=8,ncol=10) # , borderaxespad=0 bbox_to_anchor=(0.5, ...

  3. 3 用python进行OpenCV实战之画图(直线,矩形,圆形)

    前言 在上一节我们通过使用NumPy的数组分割成功的在我们的图像上画了一个绿色的方块,但是如果我们想画一个单一的线条或者圆圈该怎么办呢?NumPy没有提供相关的功能,但是OpenCV提供了相关的函数, ...

  4. java画笔覆盖在界面_Java实现画图程序和重绘

    上次聊了一下事件监听机制,今天就来聊一下怎么实现一个画图程序并且实现重绘. 一.实现画图程序 1.实现一个画图程序所需的API类? JFrame窗体容器组件类 JPanel 面板元素组件类 JButt ...

  5. matlab在曲线给命名,matlab 利用xlsread画图,怎么将一组excel数据导入,通过matlab作图...

    Matlab 循环 for 语句 xlsread EXCEL表格数据导入 画图 Matlab的 xlsread() 函数可以将Excel数据到matlab工作空间,然后就可以根据读入据作图.下面给出操 ...

  6. 计算机画图工具怎么缩小图片,Win10如何放大或缩小图片?利用win10画图工具放大、缩小图片教程...

    在日常使用电脑过程中,我们经常会碰到需要放大或缩小照片(图片)的情况.那么,win10系统下该如何扩大或者缩小照片(图片)呢?其实,我们可以通过使用win10系统自带的画图工具来实现.下面小编就向大家 ...

  7. java画出斜椭圆_【转】画图java源代码,只画直线,矩形,椭圆

    /* *只画直线,矩形,椭圆,只能向右下角画 * *PainterPanel extends JPanel implements MouseListener *addMouseListener(thi ...

  8. r语言 断轴 画图_R语言基础画图/绘图/作图

    R语言基础画图 R语言免费且开源,其强大和自由的画图功能,深受广大学生和可视化工作人员喜爱,这篇文章对如何使用R语言作基本的图形,如直方图,点图,饼状图以及箱线图进行简单介绍. 0 结构 每种图形构成 ...

  9. CCF系列之画图(201409-2)

    试题编号: 201409-2 试题名称: 画图 时间限制: 1.0s 内存限制: 256.0MB 问题描述: 问题描述 在一个定义了直角坐标系的纸上,画一个(x1,y1)到(x2,y2)的矩形指将横坐 ...

最新文章

  1. 校准曲线(calibration curve)是什么?如何绘制校准曲线(calibration curve)?如何通过过校准曲线进行分析?什么是高估?什么是低估?
  2. 天地图专题五:在天地图上绘制电子区域并保存数据
  3. 13篇京东CVPR 2019论文!你值得一读~
  4. SQL优化的一些总结
  5. 博士后斯坦福大学计算机学院,美国斯坦福大学博士后职位
  6. QQ自动播放动态视频json卡片怎么换地址
  7. mysql deadlock found_MySQL遇到Deadlock found when trying to get lock,解决方案
  8. 网络服务器主机和普通家用电脑主机的区别在哪里?
  9. 没有 本地计算机策略组,本地组策略编辑器没有mmc
  10. coron 小娜_Win10小娜近乎报废!教你卸载毫无用处的Cortana
  11. 庄子·内篇·逍遥游第一(郭象注)
  12. linux 中du和df的使用
  13. 安装Win10 Ubuntu20.04双系统
  14. HTML调用Java函数或语句,在动态THML语句中调用JS函数传递带空格参数的问题
  15. 自动化测试练习项目环境搭建
  16. 求质数(Prime Number 素数)的方法——厄拉多塞筛法
  17. 手机LBS位置应用服务盘点
  18. 建模大佬们都不外传的3D max游戏建模神操作技巧
  19. 共有24款C/C++ 地理信息系统GIS开源软件
  20. 工业智能网关BL110应用之37:实现台达 PLC Delta_DVP 接入华为云平台

热门文章

  1. 【语音识别】基于BP神经网络实现语音特征信号分类附matlab代码
  2. runat=server的根本意义是什么?
  3. dell4528s linux_戴尔台式机和笔记本电脑上的 Linux
  4. 苹果公司的专利战其实是营销战略?
  5. PHP语言打一生肖,花言巧语打一生肖 是什么动物
  6. 微信瑞文智力测试1分_瑞文智力测试
  7. 42表盘直径是从哪测量_表盘直径多大才合适?如何测量手表表盘的直径?
  8. 拼多多电商玩家快速采集平台数据
  9. YY一下上网本的新变体:个人数码伴侣
  10. 人卫第七版-流行病学期末复习重点(超全)