|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
- Z$ w7 C9 O- K: ]
鉴于以后的目标主要是利用现有的Matlab数据(.mat或者.txt),主要考虑python导入Matlab数据的问题。以下代码可以解决python读取.mat文件的问题。主要使用sicpy.io即可。sicpy.io提供了两个函数loadmat和savemat,非常方便。 y/ q- e, u* Z) K8 b/ i; d6 O
# adapted from4 d( C: q' K& R5 r, b7 Y7 u2 W' n
import scipy.io as sio $ ~, _/ w O9 M4 `0 R
#import matplotlib.pyplot as plt
5 O- X3 H& d4 Ffrom pylab import *
0 _ o9 ~& a" P" ~import numpy as np
% B( k* S( r u ( P9 Z; A7 e$ B# L. u. I
matfn='E:\Pythonrun\myuse\matdata.mat' # the path of .mat data
& Q) U2 r4 R5 ^$ X" `data=sio.loadmat(matfn) 1 o$ N5 _" S7 X2 L, Q
7 Q9 A, c: N% v2 y
xx=data['matdata']* Q; b0 J) v7 \& d% K
figure(1)
* _4 k: e7 ^6 N# i+ J* tplot(xx)9 |( o6 y6 v; d5 e7 U/ N) G
show()% p/ g- f' t# k" y# y# ^
( s1 j3 i; p" i) S$ j$ s
以下代码是读入txt数据并转换成数组,方法比较笨,更有效的方法待研究。
; U, L+ V* ]$ x$ R, Ufrom numpy import *
+ I1 P, o) S- p) m& n' |
: u, x6 J& N' e3 w" _* \def file2list(filename): 1 D7 n- d! l/ R! u7 A
fr = open(filename) ) ~* f" e7 i: [' O2 _ U
array = fr.readlines() #以文件中的每行为一个元素,形成一个list列表
) C2 p, u2 B4 e# v* f num = len(array)
$ c2 l( F2 m# N7 _$ P) p) M returnMat = zeros((num,3))#初始化元素为0的,行号数个列表,其中每个元素仍是列表,元素数是3,在此表示矩阵 + E, ]: |0 }1 j3 E6 d8 |. r
index = 0
9 r: a& j& v6 \8 Y4 K' W# ~5 b for line in array: 5 J4 d3 R6 l( Y! T$ s4 Z3 k/ E+ r \# h
line = line.strip()#去掉一行后的回车符号
) Z8 e5 H3 i+ }0 _' K3 F. P. P3 k linelist = line.split(' ')#将一行根据分割符,划分成多个元素的列表
( J# C( d* O' \6 F& @: r returnMat[index,:] = linelist[0:3]#向矩阵赋值,注意这种赋值方式比较笨拙 1 P6 O1 }" G. @! d W: \
index +=1
2 W( i5 ?! |) E; \# x return returnMat
/ r7 n+ K/ v. i- o' i3 B9 x- v- E. H' M6 k' @& X# y
fname = 'E:\Pythonrun\myuse\num_data.txt'% K0 @+ y: P/ |: _
data= file2list(fname) |
|