在Octave中使用textscan（）...如何正确格式化？

Question

我在磁盘上有一个数据集，正在逐行读取，我想将数据列之一转换为浮点向量（范围为0-23.99999）（当天）。

数据如下所示：

2010/01/01,00:00:00.979131, 27.4485,  51.9362, 14.8,  6
2010/01/01,00:00:01.021977, 27.5149,  51.9375, 16.0,  6
2010/01/01,00:00:01.074032, 27.4797,  51.9446, 14.5, 10
2010/01/01,00:00:01.663689, 25.8441,-152.8141, 14.6,  6
2010/01/01,00:00:01.639541, 25.8744,-152.6122,  1.5,  5
2010/01/01,00:00:02.232099, -2.2447,  11.5023, 18.8,  6
2010/01/01,00:00:02.256351, -0.8135,  27.3139, 17.7,  5
2010/01/01,00:00:02.306734, -2.7797,  28.5109, 26.0,  5
2010/01/01,00:00:02.620765, 25.6656,-154.2029, 26.2,  9
2010/01/01,00:00:02.658495, 25.6698,-154.2157, 23.0,  6
2010/01/01,00:00:02.731266, -5.7106, 126.4517,  3.6,  5
2010/01/01,00:00:02.787495, -5.7138, 126.5210, 24.4,  8
2010/01/01,00:00:02.811636, -3.2453, 124.6919, 21.1,  8

第2栏（例如00：00：00.979131）很有趣，我想做类似的事情

setenv GNUTERM 'x11';
fid = fopen('myfile.txt', 'r');
m = textscan(fid, '%d%d%d%d/%d%d/%d%d, %d%d:%d%d:%d%f, %f, %f, %f, %d');
mx = m(:, 5); %here, I would expect to grab 14.8, 16.0, etc
my = m(:, 2) / 24.0; %here, all data from timestamp column (00:00:00.979131, for ex)
plot(mx, my);

问题是我传递给的字符串的textscan数据格式错误。该数字的格式为“ hrs：minutes：seconds”（军事时间）。如何访问/转换var mx或的值my？

谢谢，jml

Answer 1

的输出textscan是一个单元数组。如果您在答案中使用命令：

m = textscan(fid, '%d/%d/%d %d:%d:%f %f %f %f %d', 'delimiter', ',');

然后获得14.9、16.0、14.5的垂直向量：

MyNinthField = m{9};

MyNinthField =
14.8000
16.0000
14.5000
14.6000
 1.5000
18.8000
17.7000
26.0000
26.2000
23.0000
 3.6000
24.4000
21.1000

然后，获取时间戳（从一天开始算起的秒数）：

Hours = double(m{4});
Minutes = double(m{5});
Seconds = m{6};

对于Seconds double不需要，因为m{6}已经是double。然而，m{4}和m{5}都是int32。

要获得以秒为单位的一天中的时间，您需要做的是：

TimeOfDayInSeconds = 3600*Hours+60*Minutes+Seconds;

TimeOfDayInSeconds =
0.97913
1.02198
1.07403
1.66369
1.63954
2.23210
2.25635
2.30673
2.62077
2.65849
2.73127
2.78749
2.81164

如果您不进行从int32到类型的转换double，则Octave会将值截断为整数。但是，MATLAB甚至不允许整数和双精度数组之和。

希望这可以帮助。