数据预处理的MATLAB实现

1、%-% 调用smooth函数进行加噪数据的平滑处理%-%*产生加噪正弦波信号,绘制加噪波形图* t = linspace(0,2*pi,500' % 产生一个从0到2*pi的向量,长度为500y = 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 +

2、 噪声' % 为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('加噪波形','平滑后波形'%*利用lowe

3、ss方法对加噪信号进行平滑处理,绘制平滑波形图* 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

4、方法对加噪信号进行平滑处理,绘制平滑波形图* 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方

5、法对加噪信号进行平滑处理,绘制平滑波形图* 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方法对加噪信号

6、进行平滑处理,绘制平滑波形图* 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('加噪波形','平滑后波形'%-% 调用smoothts函数对

7、金融时间序列数据进行平滑处理%-%*读取金融时间序列数据,绘制波形图* x = xlsread('examp03_02.xls' % 从文件examp03_02.xls中读取数据price = x(:,4' % 提取矩阵x的第4列数据,即收盘价数据figure; % 新建一个图形窗口plot(price,'k','LineWidth',2; % 绘制日收盘价曲线图,黑色实线,线宽为2xlabel('观测序号' ylabel('上海股市日收盘价' % 为X轴和Y轴加标签%*利用盒子法对数据进行平滑处理,绘制平

8、滑波形图* output1 = smoothts(price,'b',30; % 用盒子法平滑数据,窗宽为30output2 = smoothts(price,'b',100; % 用盒子法平滑数据,窗宽为100figure; % 新建一个图形窗口plot(price,'.' % 绘制日收盘价散点图hold onplot(output1,'k','LineWidth',2; % 绘制平滑后曲线图,黑色实线,线宽为2plot(output2,'k-.','LineWidth',2; %

9、 绘制平滑后曲线图,黑色点划线,线宽为2xlabel('观测序号' ylabel('Box method' % 为X轴和Y轴加标签% 为图形加标注框legend('原始散点','平滑曲线(窗宽30','平滑曲线(窗宽100','location','northwest'%*利用高斯窗方法对数据进行平滑处理,绘制平滑波形图* output3 = smoothts(price,'g',30; % 窗宽为30,标准差为默认值0.65output4 = smoothts(pr

10、ice,'g',100,100; % 窗宽为100,标准差为100figure; % 新建一个图形窗口plot(price,'.' % 绘制日收盘价散点图hold onplot(output3,'k','LineWidth',2; % 绘制平滑后曲线图,黑色实线,线宽为2plot(output4,'k-.','LineWidth',2; % 绘制平滑后曲线图,黑色点划线,线宽为2xlabel('观测序号' ylabel('Gaussian window method'

11、 % 为X轴和Y轴加标签legend('原始散点','平滑曲线(窗宽30,标准差0.65',.'平滑曲线(窗宽100,标准差100','location','northwest'%*利用指数法对数据进行平滑处理,绘制平滑波形图* output5 = smoothts(price,'e',30; % 用指数法平滑数据,窗宽为30output6 = smoothts(price,'e',100; % 用指数法平滑数据,窗宽为100figure; % 新建一个图形窗口plot(price,

12、'.' % 绘制日收盘价散点图hold onplot(output5,'k','LineWidth',2; % 绘制平滑后曲线图,黑色实线,线宽为2plot(output6,'k-.','LineWidth',2; % 绘制平滑后曲线图,黑色点划线,线宽为2xlabel('观测序号' ylabel('Exponential method' % 为X轴和Y轴加标签legend('原始散点','平滑曲线(窗宽30','平滑曲线(窗宽100'

13、,'location','northwest'%-% 调用medfilt1函数对加噪正弦波信号进行平滑处理%-%*产生加噪正弦波信号,绘制加噪波形图* t = linspace(0,2*pi,500' % 产生一个从0到2*pi的向量,长度为500y = 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进行中值滤波fig