Lisp-Stat 数据读取与绘图

2019独角兽企业重金招聘Python工程师标准>>>

    应研发要求，需要研究工频对产品的影响，理论不需赘述，直接上数据：

数据预处理

原始数据如下：

一般用宏处理成csv文件，分隔符为逗号和空格居多，Matlab和Lisp-Stat处理就和便捷了。

        

数据读取

Matlab的话，直接使用csvread函数即可处理。

Lisp-Stat的话，可以使用read-data-columns函数达到相同目的（以空格分隔符为例）。

> (def data (read-data-columns "e:/data.txt" 2))
((25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 13 12 12 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 19 19 19 19 19 19 19 19 19 20 20 20 20 20 19 19 19 19 19 19 19 19 19 19 18 18 17 17 17 16 16 16 15 15 15 14 14 13 13 13 12 12 11 11 11 11 11 11 11 11 11 11 11 12 12 13 13 14 14 15 15 15 16 17 17 18 18 19 19 19 20 21 21 22 22 23 23 23 24 24 25 25 26 26 27 27 27 28 28 29 29 29 30 30 30 31 31 31 31 31 31 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 31 31 31 31 31 31 31 31 31 30 30 29 29 29 28 28 27 27 27 26 25 25 24 24 23 22 21 20 19 19 18 16 15 14 13 11 11 11 12 14 15 17 18 19 20 21 22 23 24 24 25 25 26 27 27 28 28 29 29 30 30 31 31 31 31 31 31 32 32 32 32 32 32 33 32 33 32 32 32 32 32 32 32 32 32 31 32 31 32 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 32 31 31 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 31 31 31 31 31 31 31 30 30 30 30 30 29 29 29 29 28 28 28 27 27 26 26 25 25 25 24 24 23 22 22 21 20 20 19 18 17 16 15 14 13 12 11 11 11 12 13 15 16 17 18 19 20 21 22 23 23 24 25 25 26 27 27 28 28 29 29 30 30 31 31 31 31 32 32 32 32 32 32 32 32 32 32 32 32 31 31 31 31 30 30 29 29 29 28 27 27 26 25 24 23 23 21 21 19 18 17 15 14 12 11 11 12 13 15 16 18 19 20 21 22 23 24 25 25 27 27 27 28 29 29 29 30 30 31 31 31 31 31 32 32 32 32 32 32 32 32 32 32 32 31 31 31 31 30) (15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 9 9 9 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 12 12 12 12 12 11 11 11 11 11 11 10 10 9 9 9 9 9 8 8 8 8 8 8 8 8 8 8 8 9 9 9 9 10 10 11 11 11 11 12 12 12 12 13 13 13 13 13 13 14 14 14 14 14 15 15 15 15 16 16 16 16 16 16 16 17 17 17 17 17 17 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 17 17 17 17 17 16 16 16 16 16 16 15 15 15 15 14 14 13 13 13 13 12 11 11 10 9 8 8 8 9 10 11 12 12 13 13 13 14 14 15 15 15 15 16 16 16 16 16 17 17 17 17 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 17 17 17 17 17 17 17 17 17 16 16 16 16 16 16 16 15 15 15 15 15 14 14 14 13 13 13 13 12 12 11 11 10 9 9 8 8 8 9 9 11 11 12 12 13 13 13 14 14 14 15 15 15 16 16 16 16 16 17 17 17 17 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 17 17 17 17 17 16 16 16 16 15 15 14 14 13 13 13 12 12 11 10 9 8 8 9 9 11 11 12 13 13 13 14 14 15 15 15 16 16 16 16 17 17 17 17 17 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 17))

> (def ad (car data))
(25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 13 12 12 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 19 19 19 19 19 19 19 19 19 20 20 20 20 20 19 19 19 19 19 19 19 19 19 19 18 18 17 17 17 16 16 16 15 15 15 14 14 13 13 13 12 12 11 11 11 11 11 11 11 11 11 11 11 12 12 13 13 14 14 15 15 15 16 17 17 18 18 19 19 19 20 21 21 22 22 23 23 23 24 24 25 25 26 26 27 27 27 28 28 29 29 29 30 30 30 31 31 31 31 31 31 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 31 31 31 31 31 31 31 31 31 30 30 29 29 29 28 28 27 27 27 26 25 25 24 24 23 22 21 20 19 19 18 16 15 14 13 11 11 11 12 14 15 17 18 19 20 21 22 23 24 24 25 25 26 27 27 28 28 29 29 30 30 31 31 31 31 31 31 32 32 32 32 32 32 33 32 33 32 32 32 32 32 32 32 32 32 31 32 31 32 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 32 31 31 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 31 31 31 31 31 31 31 30 30 30 30 30 29 29 29 29 28 28 28 27 27 26 26 25 25 25 24 24 23 22 22 21 20 20 19 18 17 16 15 14 13 12 11 11 11 12 13 15 16 17 18 19 20 21 22 23 23 24 25 25 26 27 27 28 28 29 29 30 30 31 31 31 31 32 32 32 32 32 32 32 32 32 32 32 32 31 31 31 31 30 30 29 29 29 28 27 27 26 25 24 23 23 21 21 19 18 17 15 14 12 11 11 12 13 15 16 18 19 20 21 22 23 24 25 25 27 27 27 28 29 29 29 30 30 31 31 31 31 31 32 32 32 32 32 32 32 32 32 32 32 31 31 31 31 30)

> (def curr (car (cdr data)))
(15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 9 9 9 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 12 12 12 12 12 11 11 11 11 11 11 10 10 9 9 9 9 9 8 8 8 8 8 8 8 8 8 8 8 9 9 9 9 10 10 11 11 11 11 12 12 12 12 13 13 13 13 13 13 14 14 14 14 14 15 15 15 15 16 16 16 16 16 16 16 17 17 17 17 17 17 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 17 17 17 17 17 16 16 16 16 16 16 15 15 15 15 14 14 13 13 13 13 12 11 11 10 9 8 8 8 9 10 11 12 12 13 13 13 14 14 15 15 15 15 16 16 16 16 16 17 17 17 17 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 17 17 17 17 17 17 17 17 17 16 16 16 16 16 16 16 15 15 15 15 15 14 14 14 13 13 13 13 12 12 11 11 10 9 9 8 8 8 9 9 11 11 12 12 13 13 13 14 14 14 15 15 15 16 16 16 16 16 17 17 17 17 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 17 17 17 17 17 16 16 16 16 15 15 14 14 13 13 13 12 12 11 10 9 8 8 9 9 11 11 12 13 13 13 14 14 15 15 15 16 16 16 16 17 17 17 17 17 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 17)

绘图

> (plot-lines (iseq 1 (length ad)) ad)
#<Object: 14c63c4, prototype = SCATTERPLOT-PROTO>

从图形可以推断出：工频对采样有周期性影响，该影响又与磁路结构有关。分析可由此展开。

转载于:https://my.oschina.net/u/1011760/blog/475581