干货!从示波器到MATLAB:如何实现电子信号的数学分析

Matlab是一款应用较广的数学软件,常用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境。

Matlab语言十分简洁,这使得它在数值计算、矩阵计算等方面相对于C、Fortran等传统语言有着无以伦比的优越性。而其集成了易于使用的视窗环境,可以非常方便的绘制函数和数据图层。Matlab常常引用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测等领域。

Matlab是对已有数据的分析,而数据的采集则交由外部的采集仪来完成,比如麦克风采集音频信号、摄像机采集图像光学信号、示波器采集电路信号等。示波器作为最常用的测量仪器可以直接或间接测量出大部分领域的可感知信号,包括音频、光强、电磁波等(当然前提得转化为电信号才行)。

下面我们以mobileScope示波器(www.fescale.com)为例讲解配合matlab对信号的分析。

首先采集一组三角波和正弦波信号。


采集一组三角波和正弦波信号

点击导出,可将图层数据导出来。


“导出”功能按钮


数据文件所在目录

打开软件所在的文件夹,可以看到里面有个analyse.m 的脚本文件和matlab的文件夹,导出的数据就放在里面命名方式为xxx.dat,将这两个文件拷贝到电脑。

Analyse.m文件中集成了几个比较经典的demo,我们来看下几个主要demo的代码示例。

我们先打开analyse.m文件,将里面需要载入波形数据文件修改为实际采集的。

[if !supportLists]§ [endif]

chDatas = load('wave17_01_03_17_42_34_824.dat');

[if !supportLists]§ [endif]

matlab载入波形数据还是相对简单的,一行语句就搞定,而且还自动生成相应的矩阵。当然matlab目前也支持使用C/C++,但显然那会是一个相当烦躁的编码过程。

运行:


matlab描绘的两组波形

可以看出这个和示波器显示的是一致的,比例有所压缩,这个我们可以随意用鼠标调整。

基本矩阵运算

Matlab支持丰富的矩阵运算,加减乘除、微分、倒数、绝对值等等,都是几行语句就搞定的事。

[if !supportLists]§ [endif]

addition = ch1 + ch2; % 加法

[if !supportLists]§ [endif]

[if !supportLists]§ [endif]

subtraction = ch1 - ch2; % 减法

[if !supportLists]§ [endif]

[if !supportLists]§ [endif]

abs2 = abs(ch2); % 绝对值

[if !supportLists]§ [endif]


波形加减后的图形

上图所示是波形加减后的图形

李莎茹曲线

李莎茹图形,又称XY模式,可以用来测量两路信号相位差,评定频率比参数等强大的功能。这个在matlab上只需要一句语句就能画出来的。

[if !supportLists]§ [endif]

plot(ch1, ch2, '-g');

[if !supportLists]§ [endif]


三角波和正弦波合成的李莎茹图形

上图所示既是三角波和正弦波合成的李莎茹图形

时频域运算

对于信号处理来说,时域和频域的转换几乎是用的最多的。离散傅里叶变换FFT可以方便的把一组时域信号转换频率的权重集合,直观的反应出一个信号的频率成分。

[if !supportLists]§ [endif]

ch1FFT = abs(fft(ch1, FFT_N));

[if !supportLists]§ [endif]


频率成份

上图显示出两组信号的频率成份

坐标系转换

直角坐标系转换为极坐标显示。

[if !supportLists]§ [endif]

polarTRefer = -pi:2*pi/length(ch1):pi-2*pi/length(ch1);

[if !supportLists]§ [endif]

[if !supportLists]§ [endif]

polar(polarTRefer, ch1, 'g');

[if !supportLists]§ [endif]


三角波信号的在极坐标系下的显示


正选波信号的在极坐标系下的显示

上面图显示两个通道信号的在极坐标系下的显示

滤波处理

滤波是数字信号处理的一个比较复杂的问题,一个好的滤波器可以真实还原源数据信号有效屏蔽噪声。Matlab集成低通滤波器。

[if !supportLists]§ [endif]

filter_out = filter(filter_bz, filter_az, filter_in);

[if !supportLists]§ [endif]


绿色是原始带噪声的信号,红色为经过滤波后的信号波形

上图绿色是原始带噪声的信号,红色为经过滤波后的信号波形

曲线拟合

示波器采集的是离散点信号,有时候想要得到其拟合的曲线,以用于建立模型或者预测数据。典型的例子为,电池在没电情况下的放电曲线,信号传输过程中的衰减特性曲线或者热传导曲线等。示波器采集的是离散带噪声的样本,拟合的意义就是剔除干扰找到最佳的曲线。


曲线拟合的效果

©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 221,695评论 6 515
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 94,569评论 3 399
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 168,130评论 0 360
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 59,648评论 1 297
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 68,655评论 6 397
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 52,268评论 1 309
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 40,835评论 3 421
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 39,740评论 0 276
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 46,286评论 1 318
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 38,375评论 3 340
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 40,505评论 1 352
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 36,185评论 5 350
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 41,873评论 3 333
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 32,357评论 0 24
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 33,466评论 1 272
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 48,921评论 3 376
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 45,515评论 2 359

推荐阅读更多精彩内容