%--------------------------------------------------------------------------
% 调用 smooth 函数进行加噪数据的平滑处理
%--------------------------------------------------------------------------
%*****************产生加噪正弦波信号,绘制加噪波形图*************************
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 轴加标签
%*************利用移动平均法对加噪信号进行平滑处理,绘制平滑波形图************
*
yy1 = smooth(y,30); % 利用移动平均法对 y 进行平滑处理
figure; % 新建一个图形窗口
plot(t,y,'k:'); % 绘制加噪波形图
hold on;
plot(t,yy1,'k','linewidth',3); % 绘制平滑后波形图
xlabel('t'); % 为 X 轴加标签
ylabel('moving'); % 为 Y 轴加标签
legend('加噪波形','平滑后波形');
%*************利用 lowess 方法对加噪信号进行平滑处理,绘制平滑波形图************
*
yy2 = smooth(y,30,'lowess'); % 利用 lowess 方法对 y 进行平滑处理
figure; % 新建一个图形窗口
plot(t,y,'k:'); % 绘制加噪波形图
hold on;
plot(t,yy2,'k','linewidth',3); % 绘制平滑后波形图
xlabel('t'); % 为 X 轴加标签
ylabel('lowess'); % 为 Y 轴加标签
legend('加噪波形','平滑后波形');
%*************利用 rlowess 方法对加噪信号进行平滑处理,绘制平滑波形图***********
*
yy3 = smooth(y,30,'rlowess'); % 利用 rlowess 方法对 y 进行平滑处理
figure; % 新建一个图形窗口