clear
clc
fie=fopen（‘D:\25hz?1.dat‘‘r‘）;????%讀入脈搏數(shù)據(jù)
[a1count]=fread（fieinf‘float‘‘b‘）;
col=fclose（fie）
n=length（a1）??%轉(zhuǎn)換頻率
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的轉(zhuǎn)置
?subplot（611）;plot（a）;grid;title（‘原信號(hào)‘）;
?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（‘去噪后信號(hào)‘）;
?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（‘算子‘）;

?屬性????????????大小?????日期????時(shí)間???名稱(chēng)
-----------?---------??----------?-----??----

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

-----------?---------??----------?-----??----

??????????????????804????????????????????1