数学形态滤波

包括打开，滤波，处理，去噪等一系列功能。代码实验无误

clear
clc
fie=fopen（‘D:\25hz 1.dat‘‘r‘）;    %读入脉搏数据
[a1count]=fread（fieinf‘float‘‘b‘）;
col=fclose（fie）
n=length（a1）  %转换频率
j=1
for i=1:n;
  if  rem（i50）==0;
       b（j:）=a1（i:）;
       j=j+1
  end
end

a=b./30;
a=a‘;%a的转置
 subplot（611）;plot（a）;grid;title（‘原信号‘）;
 se=strel（‘line‘70）;
 se1=strel（‘line‘300）;
 a1=imopen（imclose（ase）se）;
 a2=imclose（imopen（ase）se）;
 a3=（a1+a2）./2;
 subplot（612）;plot（a3）;grid;title（‘去噪后信号‘）;
 a4=a3-imopen（a3se1）;
 subplot（613）;plot（a4）;grid;title（‘白TOP-HAT‘）;
 a5=imclose（a3se1）-a3;
 subplot（614）;plot（a5）;grid;title（‘黑TOP-HAT‘）;
 a6=imdilate（a3se1）-imerode（a3se1）;
 subplot（615）;plot（a6）;grid;title（‘GRAD‘）;
 a7=2.*a3-imopen（a3se1）-imclose（a3se1）;
 subplot（616）;plot（a7）;grid;title（‘算子‘）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        804  2012-05-10 15:04  yuanxinhao.m

----------- ---------  ---------- -----  ----

                  804                    1