智能仪器作业教程.docx

智能仪器作业 姓 名：贺申羽 学 号 班 级：测控三班 指导老师：李艳萍 智能仪器作业 作业一 请完成决策滤波器的程序设计，并分别对温度、加速度数据进行滤波处理，比较不同滤波参数对数据滤波效果的影响。 要求： 设计程序，完成报告，报告内容包括以下部分： 决策滤波器算法说明； 算法实现流程图； 数据原始波形图； 采用不同滤波器参数（至少三种）的滤波后数据波形图； 滤波效果的对比分析说明； 推荐的滤波器参数设置。 报告格式：word 通过校网提交报告电子版即可 基于中值数绝对偏差的决策滤波器算法说明； 中值绝对偏差估计的决策滤波器能够判别出奇异数据，并以有效性的数值来取代。采用一个移动窗口，， 利用m个数据来确定的有效性。如果滤波器判定该数据有效，则输出，否则，如果判定该数据为奇异数据，用中值来取代。 一个序列的中值对奇异数据的灵敏度远无小于序列的平均值，用中值构造一个尺度序列，设{xi(k) }中值为Z，则 给出了每个数据点偏离参照值的尺度 令{d(k)}的中值为d，著名的统计学家FR.Hampel提出并证明了中值数绝对偏差MAD＝1.4826*d，MAD可以代替标准偏差σ。对3σ法则的这一修正有时称为“Hampel标识符”。 读取数据 计算数据大小； 确定m和L的值； 初始化滤波后数据前几项 建立移动窗口； 取窗口数据中值z； 计算尺度序列中值d; 滤波后数据等于原数据 qQ*L? 计算Q、q N 滤波后数据等于窗口中值 Y 窗口移至最后？ 窗口继续向后移动 结束 N Y 算法所用程序，以温度为例 load wendu.txt X=wendu; size=length(X); m=8;L=5; for j=1:m Y(j)=X(j); end for i=1:size-m b=X(i):X(i+m); z=median(b); d=median(b-z); Q=1.4826*d; q=abs(X(i+m)-z); if q(L*Q) Y(i+m)=X(i+m); else Y(i+m)=z; end end Y=Y; figure(1);plot(X); title(温度数据原始波形图); xlabel(个数);ylabel(幅值); figure(2);plot(Y); title(温度数据滤波后波形图); legend(m=8,L=5) xlabel(个数);ylabel(幅值); 数据原始波形图； 采用不同滤波器参数（至少三种）的滤波后数据波形图； 温度第一组 M=8 L=5 温度第二组 M=8 L=10 温度第三组 M=15 L=10 加速度数据 第一组 M=8 L=10 第二组数据 M=10 L=10 第三组 M=15 L=10 滤波效果的对比分析说明: 从上面六组数据中，我们可以分析得知 在温度的三组数据中，在未经滤波器处理之前，图中有较多明显毛刺，与温度的实际变化情况不符。在经过滤波器处理之后，毛刺明显减少，并且从图中可知M,L对温度的幅值变化没有明显影响。 在加速度的三组数据中，在经滤波器处理之后，幅值范围明显变小。而且我们可以得知，M越大，它的幅值范围越大，L的值对幅值范围没有明显影响。