首先简述一下中值滤波算法:以1000个数据为例,取奇数个数为一组(如5个数为一组),则共有995组数据。
为什么是995组?假如我们取第1,2,3,4,5个数为第一组,那么第二组是第2,3,4,5,6这5位数。若第2组是6,7,8,9,10,那么总的数据会少了将近800组,严重损坏数据。
简单简述下代码
首先是数据提取。为了简单起见,我从一个二维txt文件中提取数据,并取其中的一列(代码显示是取第7列)
data_sum_A=importdata('E:\MATLAB代码\toulan.txt');
data_B=data_sum_A(:,7);
importdata函数是matlab中非常好用的从文件中提取数据的算法。具体这里不多说了,可以在命令窗口中用 help importdata 查一下
接下来先把原图画出来
subplot(2,1,1);
address_B = plot(data_B);
title('原始数据');
get(address_B);
set(address_B,'LineStyle','-.');
用get以及set函数设定图像的样式。
用上面说的算法处理数据:
i=1;
N=4;
L=length(data_B);
for i=1:L-N
data_C(i)=median(data_B(i:i+N));
end
subplot(2,1,2);
plot(data_C);
title('中值滤波后');
这里用到了取中位数的函数median,其返回值就是一个中位数。