小编Wil*_*l C的帖子

在C中跳过读取文本文件的标题

我正在从文件中读取数据orderedfile.txt。有时该文件具有以下形式的标题:

BEGIN header

       Real Lattice(A)               Lattice parameters(A)    Cell Angles
   2.4675850   0.0000000   0.0000000     a =    2.467585  alpha =   90.000000
   0.0000000  30.0000000   0.0000000     b =   30.000000  beta  =   90.000000
   0.0000000   0.0000000  30.0000000     c =   30.000000  gamma =   90.000000

 1                            ! nspins
25   300   300                ! fine FFT grid along <a,b,c>
END header: data is "<a b c> pot" in units of Hartrees

 1     1     1            0.042580
 1     1     2            0.049331
 1     1     3            0.038605
 1     1     4            0.049181
Run Code Online (Sandbox Code Playgroud)

有时没有标题,数据从第一行开始。我的读取数据的代码如下所示。当数据从第一行开始但不存在标题时它会起作用。有办法解决这个问题吗?

int …
Run Code Online (Sandbox Code Playgroud)

c file-io

5
推荐指数
1
解决办法
2841
查看次数

从(x,y)坐标中找到行进的距离

我目前有一个python脚本,该脚本读取一个3列的文本文件,其中包含助行器的x和y坐标以及行人走路的时间。

我已读取此数据并将其分配给numpy数组,如以下代码所示:

import numpy as np
import matplotlib.pyplot as plt

data = np.loadtxt("info.txt", delimiter = ',')

x = data[:,0]
y = data[:,1]
t = data[:,2]
Run Code Online (Sandbox Code Playgroud)

文件采用以下格式(x,y,t):

5907364.2371    -447070.881709  2193094
5907338.306978  -447058.019176  2193116
5907317.260891  -447042.192668  2193130
Run Code Online (Sandbox Code Playgroud)

现在,我想找到步行者随时间变化的距离。我可以想到的一种方法是在循环中将x坐标的差异和y坐标的所有差异相加。但是,这似乎是一个长篇大论的方法,我认为可以用一种数值积分方法来解决。有人对我能做什么有任何想法吗?

python scipy numerical-methods numerical-integration

1
推荐指数
1
解决办法
1259
查看次数