数学建模扩展：2 数据预处理调用smooth多种平滑函数对加噪数据进行平滑处理.zip

%% 学习目标： 调用smooth多种平滑函数对加噪数据进行平滑处理 %% QQ：1960009019 %% 更多matlab精彩专题课程和案例，可以搜索微信公众号：大仙一品堂 %% 产生加噪正弦波信号，绘制加噪波形图*** 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; % 新建一个图形窗口 plot(t,y,'k:'); % 绘制加噪波形图 hold on; plot(t,yy3,'k','linewidth',3); % 绘制平滑后波形图 xlabel('t'); % 为X轴加标签 ylabel('rlowess'); % 为Y轴加标签 legend('加噪波形','平滑后波形'); %% 利用loess方法对加噪信号进行平滑处理，绘制平滑波形图 yy4 = smooth(y,30,'loess'); % 利用loess方法对y进行平滑处理 figure; % 新建一个图形窗口 plot(t,y,'k:'); % 绘制加噪波形图 hold on; plot(t,yy4,'k','linewidth',3); % 绘制平滑后波形图 xlabel('t'); % 为X轴加标签 ylabel('loess'); % 为Y轴加标签 legend('加噪波形','平滑后波形'); %% 利用sgolay方法对加噪信号进行平滑处理，绘制平滑波形图 yy5 = smooth(y,30,'sgolay',3); % 利用sgolay方法对y进行平滑处理 figure; % 新建一个图形窗口 plot(t,y,'k:'); % 绘制加噪波形图 hold on; plot(t,yy5,'k','linewidth',3); % 绘制平滑后波形图 xlabel('t'); % 为X轴加标签 ylabel('sgolay'); % 为Y轴加标签 legend('加噪波形','平滑后波形');