基于matlab的中值滤波算法浅析

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

•