各位好
一點問題請教
我目前有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
不知道這樣指令是否正確?????
感謝