资源描述:
《数据预处理的MATLAB实现.doc》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库。
1、%--------------------------------------------------------------------------%调用smooth函数进行加噪数据的平滑处理%--------------------------------------------------------------------------%*****************产生加噪正弦波信号,绘制加噪波形图*************************t=linspace(0,2*pi,500)';%产生一个从0到2*pi的向量,长度为500y=100*sin(t);%产
2、生正弦波信号%产生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:');%绘制加噪波
3、形图holdon;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:');%绘制加噪波形图holdon;plot(t,yy2,'k','linewidth',3);%绘
4、制平滑后波形图xlabel('t');%为X轴加标签ylabel('lowess');%为Y轴加标签legend('加噪波形','平滑后波形');%*************利用rlowess方法对加噪信号进行平滑处理,绘制平滑波形图************yy3=smooth(y,30,'rlowess');%利用rlowess方法对y进行平滑处理figure;%新建一个图形窗口plot(t,y,'k:');%绘制加噪波形图holdon;plot(t,yy3,'k','linewidth',3);%绘制平滑后波形图xlabel('t');%为X轴加标签ylabel('rlowe
5、ss');%为Y轴加标签legend('加噪波形','平滑后波形');%*************利用loess方法对加噪信号进行平滑处理,绘制平滑波形图*************yy4=smooth(y,30,'loess');%利用loess方法对y进行平滑处理figure;%新建一个图形窗口plot(t,y,'k:');%绘制加噪波形图holdon;plot(t,yy4,'k','linewidth',3);%绘制平滑后波形图xlabel('t');%为X轴加标签ylabel('loess');%为Y轴加标签legend('加噪波形','平滑后波形');%*********
6、****利用sgolay方法对加噪信号进行平滑处理,绘制平滑波形图*************yy5=smooth(y,30,'sgolay',3);%利用sgolay方法对y进行平滑处理figure;%新建一个图形窗口plot(t,y,'k:');%绘制加噪波形图holdon;plot(t,yy5,'k','linewidth',3);%绘制平滑后波形图xlabel('t');%为X轴加标签ylabel('sgolay');%为Y轴加标签legend('加噪波形','平滑后波形');%------------------------------------------------
7、--------------------------%调用smoothts函数对金融时间序列数据进行平滑处理%--------------------------------------------------------------------------%*********************读取金融时间序列数据,绘制波形图************************x=xlsread('examp03_02.xls');%从文件examp03_02.xls中读取数据