主要内容

重新取样和过滤非均匀采样信号

打开生活的脚本

一个人记录自己的体重在2012磅在闰年。这个人没有每天记录自己的体重,所以数据不均匀。使用信号分析仪应用进行预处理,研究记录体重。应用程序使您可以填写缺失的数据点通过插值信号均匀网格。(这个过程给最好的结果如果信号只有很小的差距。)

加载数据,并将测量值转换为公斤。数据文件已丢失的数据集南。有27个数据点失踪,其中大部分是在8月份为期两周的延伸。

wt = datetime (2012、1、1:366)”;负载weight2012.dat重量= weight2012 (:, 2) / 2.20462;validpoints = ~ isnan(重量);失踪= wt (~ validpoints);失踪(15:26)
ans =12 x1 datetime09 - 8月- 2012年10 - 8月- 2012年11 - 8月- 2012年12 - 2012年8月- 2012年8月15 - 16日- 2012年8月17日2012年8月- 2012年8月18 - 19日- 8月- 2012年8月20 - 22 - 2012 - 2012年8月23日- 8月- 2012年

将数据存储在一个MATLAB®的时间表。删除丢失的点。删除直流值集中在波动。转换的时间信息持续时间数组被减去第一个时间点。更多细节,请参阅数据类型支持的信号分析金宝app仪。

重量=重量(validpoints);重量=重量-意思是(重量);wt = wt (validpoints);wt = wt - wt (1);wg =时间表(wt、重量);

开放信号分析仪并拖动的时间表显示。在显示选项卡上,单击光谱打开一个光谱的观点。在时间选项卡上,选择显示标记。放大失踪的延伸通过设置时间限制来200年和250年天。

右键单击并选择信号的信号表重复的。重命名复制的预处理通过右击并选择信号的信号表重命名。选择预处理信号的信号表,在分析仪选项卡上,单击进行预处理。在功能画廊,选择重新取样。在函数参数面板出现,指定这些参数:

  • 重采样方法- - - - - -采样率

  • 频率的单位- - - - - -周期/天

  • 采样率- - - - - -1

  • 插值法- - - - - -形状保持立方方法

点击应用然后点击接受所有保存并退出预处理结果模式。重新取样信号叠加显示的选择其名称旁边的复选框。

缩小显示全年的数据。在光谱选项卡,设置泄漏到最大值。原始的光谱和重新取样信号对大多数频率吻合较好。频谱显示了两个明显的峰值,大约在0.14周期/一天,另一个在非常低的频率。更好地定位山峰,显示选项卡上,单击数据游标并选择两个。山峰上的游标。每个光标悬停在频率域更精确值的位置。

  • 中频峰在0.143 = 1/7周期/天,这对应于一个一周的周期。

  • 低频峰在0.005周期/天,这对应于一个210天的周期。

通过点击删除游标数据游标。删除原始信号的显示。过滤预处理信号去除的影响周期。预处理信号在信号中选择表,在分析仪选项卡上,单击进行预处理。在预处理模式:

  1. 消除低频周期,高通滤波器的信号。选择高通滤波从功能画廊。在函数参数出现面板中,输入一个通频带的频率0.05周期/天。使用其他参数的默认值。点击应用。

  2. 把每周循环,bandstop-filter信号。选择Bandstop从功能画廊。在函数参数面板中,输入一个低通频带的频率0.135周期/天和通带频率更高0.15周期/天。使用其他参数的默认值。点击应用。

检查结果和点击接受所有。情节都在显示信号。

预处理信号显示波动比原来的少。信号的形状表明,人的体重变化在夏季比冬季,但这可能是一个工件的重采样。单击图标信息列的信号表条目预处理执行信号的预处理步骤。

signalAnalyzer_weights5_22a.png

看到一个完整的总结的预处理步骤,包括所有你选择的设置,点击生成函数在分析仪选项卡。生成函数出现在MATLAB®编辑器。

函数(y, ty) =预处理(x, tx)%预处理输入x%这个函数需要一个输入向量x和一个向量的时间值% tx. tx数值向量以秒为单位。%按照时间表文档(类型“医生的时间表”%命令行)学习如何索引表变量及其时间%值,这样您就可以将它们传递到这个函数。%由MATLAB 9.13 (R)和信号处理工具箱9.0。%上生成:03 - 1月- 2022 15:47:27targetSampleRate = 1.1574074074074073 e-05;(y, ty) =重新取样(x, tx targetSampleRate,“pchip”);Fs = 1 /意味着(diff (ty));%的平均采样率y =高通滤波(y, 5.787 e-07, Fs,“Steepnes”,0.85,“StopbandAttenuation”、60);y = bandstop (y, [1.5625 e-06 1.73611111111111 e-06], Fs,陡峭的,0.85,“StopbandAttenuation”、60);结束

另请参阅

应用程序

功能

相关的例子

更多关于