Lisp-Stat 數據讀取與繪圖

    應研發要求，須要研究工頻對產品的影響，理論不需贅述，直接上數據：

數據預處理

原始數據以下：

通常用宏處理成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>

從圖形能夠推斷出：工頻對採樣有周期性影響，該影響又與磁路結構有關。分析可由此展開。