形态滤波的matlab代码

clear all
x=0:1:1000;
Fs=1000; 
nfft=1024;
y=3*square（2*pi*0.01*x50）;       %矩形波
subplot（321）;plot（xy）;grid;axis（[01000-55]）;
title（‘模拟信号时域波形‘）xlabel（‘时间 t/s‘）ylabel（‘电压 V/v‘）;
yBW1=abs（fft（ynfft））;     %模拟信号的 频域波形
subplot（322）plot（（0:nfft/2-1）/nfft*FsyBW1（1:nfft/2））;
axis（[0250-50400]）;
title（‘模拟信号频谱‘）xlabel（‘频率 f/Hz‘）ylabel（‘振幅‘）;

yn=awgn（y15‘measured‘）;         %增加噪声信号
subplot（323）;plot（yn）;grid;axis（[01000-55]）;  
title（‘含噪信号时域波形‘）xlabel（‘时间 t/s‘）ylabel（‘电压 V/v‘）;
yBW2=abs（fft（ynnfft））;   %含噪信号的 频域波形   
subplot（324）plot（（0:nfft/2-1）/nfft*FsyBW2（1:nfft/2））; 
axis（[0250-50400]）;
title（‘含噪信号频谱‘）xlabel（‘频率 f/Hz‘）ylab