Lisp-Stat 数据读取与绘图

原创
2015/07/07 12:19
阅读数 202

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

数据预处理

原始数据如下:

一般用宏处理成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>

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

展开阅读全文
加载中

作者的其它热门文章

打赏
1
0 收藏
分享
打赏
2 评论
0 收藏
1
分享
返回顶部
顶部