由於我在用SPH方法仿真時用的是FORTRAN語言，並且沒有找到直接輸出vtk文件的代碼，因此偷懶通過MATLAB轉換一下數據。
用到的Matlab子程序可通過一下鏈接找到。
Matlab處理數據導出Paraview可讀的vtk文件（一）
假設我每個粒子輸出的格式按照下面輸出的
x(i), y(i), z(i), vx(i), vy(i), vz(i), p(i), rho(i), phase(i)
含義為每個粒子的坐標，速度，壓強，密度，相。
vtkwrite(‘output.vtk’,'unstructured_grid',x(:),y(:),z(:),'vectors','velocity',vx(:),vy(:),vz(:),'scalars', 'pressure',p(:), 'scalars', 'rho',rho(:),'scalars', 'phase',phase(:));
解釋：
1，輸出文件名
2，定義非結構網格，還有其他選項
3，后面跟每個點的坐標
4，‘vectors’定義向量，后面‘velocity’是自定義的名字，比如速度，然后后面跟速度的三個分量
5，‘scalars’定義數值，后面引號中自定義名字，比如壓強和密度等
6，向量和數值最后在paraview中都有對應的顯示方式。
將生成的vtk文件拖進paraview就可以在Representation中選擇Point Gaussian。
在Coloring中選擇要顯示的數據
二、顯示向量
1、選擇菜單中filter->Alphabetical->Glyph
或者點下面
/2、設置你要顯示的向量
