首先简述一下中值滤波算法:以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,其返回值就是一个中位数。