matlab中smooth函数平滑处理数据实例.docxVIP

t=linspace(0,2*pi,500);  %产生一个从  0到2*pi  的向量，长度  500 y=100*sin(t);  %产生正弦波信号 %产生500行1列的听从N（0,15*15）散布的随机数，作为噪声信号 noise=normrnd(0,15,500,1); y=y+noise; figure; plot(t,y); xlabel(t); ylabel(y=sin(t)+噪声); yy1=smooth(y,30); figure; plot(t,y,k:); holdon; plot(t,yy1,k,linewidth,3);  %将噪声信号加入正弦波信号 %新建一个图形窗口 %绘制加噪波形图 %为X轴加标签 %为Y轴加标签 %利用挪动均匀法对y做光滑办理 %新建一个图形窗口 %绘制加噪波形图 %绘制光滑后波形图 xlable(t); xlabel(t); ylabel(moving); legend(加噪波形,光滑后波形); yy2=smooth(y,30,lowess);%利用lowess方法y做光滑办理 figure; plot(t,y,k); plot(t,y,k:); holdon; plot(t,yy2,k,linewidth,3); xlabel(t); ylabel(lowess); legend(加噪波形,光滑后波形); yy3=smooth(y,30,rlowess);%利用rlowess方法y做光滑办理 figure; plot(t,y,k:); holdon; plot(t,yy3,k,linewidth,3); xlabel(t); ylabel(rlowess); legend(加噪波形,光滑后波形); yy4=smooth(y,30,loess); figure; plot(t,y,k:); holdon; plot(t,yy4,k,linewidth,3);  %利用loess方法y做光滑办理 xlabel(t); ylabel(loess); legend(加噪波形,光滑后波形); yy5=smooth(y,30,sgolay,3);%利用sgolay方法y做光滑办理 figure; plot(t,y,k:); holdon; plot(t,yy5,k,linewidth,3); xlabel(t); ylabel(sgolay); legend(加噪波形,光滑后波形); 以下为成效图： (加噪波形图) （加噪波形图2） （光滑后） （加标签） lowess法办理）