本帖最后由 Freeman_2016 于 2019-9-18 11:22 编辑 * Z% N5 G" F5 h1 n P0 M# U
5 k$ C( W+ n" o- rMatlab是一款应用较广的数学软件,常用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境。 Matlab语言十分简洁,这使得它在数值计算、矩阵计算等方面相对于C、Fortran等传统语言有着无以伦比的优越性。而其集成了易于使用的视窗环境,可以非常方便的绘制函数和数据图层。Matlab常常引用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测等领域。 Matlab是对已有数据的分析,而数据的采集则交由外部的采集仪来完成,比如麦克风采集音频信号、摄像机采集图像光学信号、示波器采集电路信号等。示波器作为最常用的测量仪器可以直接或间接测量出大部分领域的可感知信号,包括音频、光强、电磁波等(当然前提得转化为电信号才行)。 下面我们以mobileScope示波器(www.fescale.com)为例讲解配合matlab对信号的分析。 首先采集一组三角波和正弦波信号。
/ c% X+ n! Q1 S9 c3 X1 I
' J; a' n, ~1 ?! c/ q3 s3 D
采集一组三角波和正弦波信号 * P' p7 N* a [4 m( \0 {
7 l3 ~/ Y( K; q" M2 J
点击导出,可将图层数据导出来。 “导出”功能按钮 / m! S+ i& Y4 G* d
! ?3 |9 n) z, R) @$ F
数据文件所在目录 , V0 K5 W6 c) W% @8 l
& z' b1 W' J) T( j
打开软件所在的文件夹,可以看到里面有个 analyse.m 的脚本文件和matlab的文件夹,导出的数据就放在里面命名方式为xxx.dat,将这两个文件拷贝到电脑。 Analyse.m文件中集成了几个比较经典的demo,我们来看下几个主要demo的代码示例。 我们先打开analyse.m文件,将里面需要载入波形数据文件修改为实际采集的。 § chDatas = load('wave17_01_03_17_42_34_824.dat'); § matlab载入波形数据还是相对简单的,一行语句就搞定,而且还自动生成相应的矩阵。当然matlab目前也支持使用C/C++,但显然那会是一个相当烦躁的编码过程。 运行: ' P% T. q$ k6 A; G3 H0 {+ m
$ S6 w( U6 N7 w
matlab描绘的两组波形
3 K& ^2 W+ q3 a v" N4 y$ e, y; y
# e& _# T$ ?& N) ^- P可以看出这个和示波器显示的是一致的,比例有所压缩,这个我们可以随意用鼠标调整。 基本矩阵运算Matlab支持丰富的矩阵运算,加减乘除、微分、倒数、绝对值等等,都是几行语句就搞定的事。 § addition = ch1 + ch2; % 加法 § § subtraction = ch1 - ch2; % 减法 § § abs2 = abs(ch2); % 绝对值 § " i& `$ L. M3 _! `5 C" M
( {! p2 w' c$ Q: g/ y6 N
波形加减后的图形
/ j; R. }' ~3 ^ j5 [8 P8 i
' ~/ n8 G2 B- {; i% O9 H
上图所示是波形加减后的图形 李莎茹曲线李莎茹图形,又称XY模式,可以用来测量两路信号相位差,评定频率比参数等强大的功能。这个在matlab上只需要一句语句就能画出来的。 § plot(ch1, ch2, '-g'); §
3 b; Q1 T8 [( `, ~) p+ a5 U
& p2 I2 ~8 \& h# i* Z t; U# K3 o三角波和正弦波合成的李莎茹图形
( u @5 |0 @0 a" ]$ A+ M: d
) }" G1 S7 \9 m3 h7 c! ?
上图所示既是三角波和正弦波合成的李莎茹图形 时频域运算对于信号处理来说,时域和频域的转换几乎是用的最多的。离散傅里叶变换FFT可以方便的把一组时域信号转换频率的权重集合,直观的反应出一个信号的频率成分。 § ch1FFT = abs(fft(ch1, FFT_N)); §
& A$ y2 m1 G2 y @/ s) v; F7 h! D0 o. c6 C9 {! J6 a5 J! r
频率成份 8 |( q# j+ k& }* j4 P
. w9 \" J) z- B/ r! N) b8 J( u' u
上图显示出两组信号的频率成份 坐标系转换直角坐标系转换为极坐标显示。 § polarTRefer = -pi:2*pi/length(ch1):pi-2*pi/length(ch1); § § polar(polarTRefer, ch1, 'g'); §
) L b; B+ u8 a/ i1 c, i& _# U3 u" ~9 A* \
三角波信号的在极坐标系下的显示 * M; a8 h1 b/ n( W4 s# r7 k
3 n! D' m1 l' k- H$ o7 ]
正选波信号的在极坐标系下的显示 & _, N. O7 @; d3 j. D3 b7 }
+ R- K8 z$ n9 w; r+ C5 X& o z
上面图显示两个通道信号的在极坐标系下的显示 滤波处理滤波是数字信号处理的一个比较复杂的问题,一个好的滤波器可以真实还原源数据信号有效屏蔽噪声。Matlab集成低通滤波器。 § filter_out = filter(filter_bz, filter_az, filter_in); § ' j: m% T* w( g u7 j
& t( s7 j( b% e
绿色是原始带噪声的信号,红色为经过滤波后的信号波形 8 K& b7 G$ Y/ Y+ I$ K
8 J: Y: `- T( U/ C" _6 p! {, N9 c
; W! \* L+ ]: g/ p4 Z) @' g" `* ~
上图绿色是原始带噪声的信号,红色为经过滤波后的信号波形 曲线拟合示波器采集的是离散点信号,有时候想要得到其拟合的曲线,以用于建立模型或者预测数据。典型的例子为,电池在没电情况下的放电曲线,信号传输过程中的衰减特性曲线或者热传导曲线等。示波器采集的是离散带噪声的样本,拟合的意义就是剔除干扰找到最佳的曲线。 , l* V2 O7 y2 B
/ U0 V. D! l" ?8 A$ `* u- r
, s8 ^2 S% E4 r( c; I7 F9 O
3 Y7 y1 s5 G+ D# S2 S/ d
# P% `' p1 B$ s6 {
/ J1 W$ X6 U* M) _+ k2 s4 h3 `# Q+ T# O& v$ X" j, k
% m, e" O4 G5 g4 R* [1 V2 _& t
/ k8 G9 `/ T4 a& W0 Y! d7 F |