%--------------------------------------------------------------------------
% 调用medfilt1函数对加噪正弦波信号进行平滑处理
%--------------------------------------------------------------------------
%*****************产生加噪正弦波信号,绘制加噪波形图*************************
t = linspace(0,2*pi,500)'; % 产生一个从0到2*pi的向量,长度为500
y = 100*sin(t); % 产生正弦波信号
% 产生500行1列的服从N(0,152)分布的随机数,作为噪声信号
noise = normrnd(0,15,500,1);
y = y + noise; % 将正弦波信号加入噪声信号
figure; % 新建一个图形窗口
plot(t,y); % 绘制加噪波形图
xlabel('t'); % 为X轴加标签
ylabel('y = sin(t) + 噪声'); % 为Y轴加标签
%***********调用medfilt1对加噪正弦波信号y进行中值滤波,并绘制波形图***********
yy = medfilt1(y,30); % 指定窗宽为30,对y进行中值滤波
figure; % 新建一个图形窗口
plot(t,y,'k:'); % 绘制加噪波形图
hold on
plot(t,yy,'k','LineWidth',3); % 绘制平滑后曲线图,黑色实线,线宽为3
xlabel('t'); % 为X轴加标签
ylabel('中值滤波'); % 为Y轴加标签
legend('加噪波形','平滑后波形');