首页 > 其他分享 >VMD计算原子到面的最短距离,TCL脚本

VMD计算原子到面的最短距离,TCL脚本

时间:2024-11-07 11:16:26浏览次数:5  
标签:atomselect set dist index VMD 原子 到面 TCL top

#calculate the distance between atom and planes in the trajectory# 

 set result [open dist-h-mental.txt w]
 set nframe [molinfo top get numframes]
 #目标原子
 set p [atomselect top "index 48"]
 #组成平面的三个原子
 set a [atomselect top "index 44"]
 set b [atomselect top "index 45"]
 set c [atomselect top "index 46"]
 set sum 0
 #遍历每一帧
 for {set i 0} {$i<$nframe} {incr i 1} {
 $p frame $i
 $a frame $i
 $b frame $i
 $c frame $i
 $p update
 $a update
 $b update
 $c update
 #计算向量 ba 和 bc
 set vecba [vecsub [measure center $b] [measure center $a]]
 set vecbc [vecsub [measure center $b] [measure center $c]]
 #计算法向量及其长度
 set normal [veccross $vecba $vecbc ]
 set normal_length [veclength $normal]
 #计算向量pa
 set target_to_plane [vecsub [measure center $p] [measure center $a]]
 #计算距离
 set distance [ expr abs([vecdot $target_to_plane $normal]) / $normal_length]
 puts $result "$i $distance"
 set sum [expr $sum+$distance]
 }
 close $result
 puts "the average diatance is : [expr double($sum)/$nframe]"

标签:atomselect,set,dist,index,VMD,原子,到面,TCL,top
From: https://blog.csdn.net/weixin_44620047/article/details/143589422

相关文章

  • Tcl脚本介绍(一)
      芯冰乐知识星球入口:芯冰乐TCL就是ToolCommandLanguage的简称,广泛应用在各大eda工具中。EDA工具能够高度自动化的运行得益于TCL。用户在编写完相关的tcl脚本后,让eda工具自动运行,便能够很大程度地解放双手了。其实,网络上介绍TCL脚本的课程层出不穷了。应很多粉丝的要求......
  • 20240712 TCL 笔试
    岗位:华星-电子设计类(提前批)笔试题型:20道不定项选择题1、不定项选择1.1【单选】下列关于理想二极体的说法,错误的是?(D)顺向时视为短路,逆向时视为开路顺向电阻等于零,逆向电阻无限大无顺向电压降,无逆向电流顺向电压等于零,逆向电流无限大解答:理想二极管:顺向导通时视为......
  • 冲一区!双重分解+粒子群优化+深度学习多元时序预测!CEEMDAN-Kmeans-VMD-PSO-Transformer
    目录效果一览基本介绍程序设计参考资料效果一览基本介绍1.Matlab实现CEEMDAN-Kmeans-VMD-PSO-Transformer融合K均值聚类的数据双重分解+粒子群优化+Transformer多元时间序列预测(完整源码和数据)运行环境Matlab2023b及以上。2.CEEMDAN分解,计算样本......
  • UEC++中的GetClass和StaticClass函数
    GetClass()用途:GetClass() 是 UObject 类的一个实例方法,用于获取调用它的对象的类信息。返回类型:返回 UClass*,即指向调用对象的类的 UClass 对象的指针。使用场景:当你有一个 UObject 或其子类的实例,并且想要获取这个实例所属类的信息时,你会使用 GetClass()。例......
  • KVM虚拟机下利用宝塔面板中的 Docker 配置 Nextcloud + onlyoffice
    第一部分:安装配置nextcloud准备(1)启动一个KVM虚拟机,见下图。该虚拟机为anlisos8…0.2虚拟机的ssh、hostname、IP地址都已配置好。(2)宝塔面板也已安装好docker一、环境docker版本:29.0.4DockerCompose版本:2.27.0mariadb-docker版本:10.6.17......
  • TCL中环开工率下滑,员工集体要求解约赔偿
    “尽管中环的市占率有所提高,但是高开工率也带来了巨量硅片库存,严重拖累了公司业绩。”转载:@科技新知 原创作者丨依蔓编辑丨蕨影因大幅下调开工率,光伏硅片龙头TCL中环疑似遭遇员工“离职潮”? 近日,“财联社”注意到,TCL中环子公司天津环智和天津环欧出现员工集体......
  • 轴承故障诊断 (12)基于交叉注意力特征融合的VMD+CNN-BiLSTM-CrossAttention故障识别模
    往期精彩内容:Python-凯斯西储大学(CWRU)轴承数据解读与分类处理Pytorch-LSTM轴承故障一维信号分类(一)-CSDN博客Pytorch-CNN轴承故障一维信号分类(二)-CSDN博客Pytorch-Transformer轴承故障一维信号分类(三)-CSDN博客三十多个开源数据集|故障诊断再也不用担心数据集了!P......
  • scran doubletCluster
      Identifypotentialclustersofdoubletcellsbasedonwhethertheyhaveintermediateexpressionprofiles,i.e.,theirprofilesliebetweentwoother“source”clusters.Thisfunctionisnowdeprecated,usefindDoubletClustersfromscDblFinderinstead.......
  • VMD-DBO-CNN-BiLSTM四模型多变量时间序列光伏功率预测一键对比 Matlab代码
    基于VMD-DBO-CNN-BiLSTM、VMD-CNN-BiLSTM、VMD-BiLSTM、BiLSTM四模型多变量时间序列光伏功率预测一键对比(仅运行一个main即可)[原创未发表]Matlab代码每个模型的预测结果和组合对比结果都有!运行步骤:1.先运行main1进行VMD分解2.在运行main2进行四模型一键对比代码......
  • 【风电功率预测】【多变量输入单步预测】基于VMD-CNN-LSTM的风电功率预测研究(Matlab代
                      ......