lige826603209 发表于 2022-10-26 16:36
如果想统计的键数量比较少的话，可以手动把它设置进Set。

然后对数组@bondlength里记录的键长进行操作。

好的，感谢

如果想统计的键数量比较少的话，可以手动把它设置进Set。

1. my @bondlength;
   
2. my $trajectory = $Documents{".xtd"}->Trajectory;
   
3. for (my $i=1; $i<=$trajectory->NumFrames; ++$i) {
   
4.     my $frame = $trajectory->CurrentFrame = $i;
   
5.     @bondlength[$i-1] = $trajectory->Set("yourbond")->Item->Length;
   
6. }

复制代码

然后对数组@bondlength里记录的键长进行操作。

菜还不想动 发表于 2022-10-25 14:04
您好，看到文献中分子动力学部分中有 maximum bond length deviation（MAXD）的图，我看了之前 ...

不会写就学
简单的分析脚本不会写的话做MD寸步难行，不可能遇到需要专一性的分析的时候总是恰好就能找到个现成的脚本
VMD轨迹分析脚本的编写在北京科音分子动力学与GROMACS培训班（http://www.keinsci.com/workshop/KGMX_content.html）里有非常全面系统的讲解

sobereva 发表于 2022-10-25 01:00
提问时不要用不常见的缩写，MAXD别人都不知道你指什么

老师，是maximum bond length deviation (MAXD)，这个该怎么输出？

提问时不要用不常见的缩写，MAXD别人都不知道你指什么