各位好
一点问题请教
我目前有20个分子座标的XYZ档案,资料储存在一个档案1.xyz之下
1_X1 1_Y1 1_Z1
1_X2 1_Y2 1_Z2
1_X3 1_Y3 1_Z3
2_X1 2_Y1 2_Z1
2_X2 2_Y2 2_Z2
2_X3 2_Y3 2_Z3
2_X4 2_Y4 2_Z4
2_X5 2_Y5 2_Z5
这20个分子的大小皆不同 但有个最大的size, ex N=20
我想做一件事情那就是 把这20个分子label成 1~20
然后
我希望储存成 这样的阵列
假设分子一 他的座标就是X(i,1) Y(i,1),Z(i,1), i=1~N(1)
假设分子二 他的座标就是X(i,2) Y(i,2),Z(i,2) , i=1~N(2)
index 1 ,2,...20 标示第几个分子
请问我要怎么写我的指令???
我想过这样
for k=1:20
x[:,k],y[:,k],z[:,k]=importdata('1.xyz')
end
不知道这样指令是否正确?????
感谢