数据预处理的MATLAB实现.doc

%-% 调用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) + 噪声); % 为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(加噪波形,平滑后波形);%-% 调用smoothts函数对金融时间序列数据进行平滑处理%-%*读取金融时间序列数据，绘制波形图*x = xlsread(examp03_02.xls); % 从文件examp03_02.xls中读取数据price = x(:,4); % 提取矩阵x的第4列数据，即收盘价数据figure; % 新建一个图形窗口plot(price,k,LineWidth,2); % 绘制日收盘价曲线图，黑色实线，线宽为2xlabel(观测序号); ylabel(上海股市日收盘价); % 为X轴和Y轴加标签%*利用盒子法对数据进行平滑处理，绘制平滑波形图*output1 = smoothts(price,b,30); % 用盒子法平滑数据，窗宽为30output2 = smoothts(price,b,100); % 用盒子法平滑数据，窗宽为100figure; % 新建一个图形窗口plot(price,.); % 绘制日收盘价散点图hold onplot(output1,k,LineWidth,2); % 绘制平滑后曲线图，黑色实线，线宽为2plot(output2,k-.,LineWidth,2); % 绘制平滑后曲线图，黑色点划线，线宽为2xlabel(观测序号); ylabel(Box method); % 为X轴和Y轴加标签% 为图形加标注框legend(原始散点,平滑曲线(窗宽30),平滑曲线(窗宽100),location,northwest);%*利用高斯窗方法对数据进行平滑处理，绘制平滑波形图*output3 = smoothts(price,g,30); % 窗宽为30，标准差为默认值0.65output4 = smoothts(price,g,100,100); % 窗宽为100，标准差为100figure; % 新建一个图形窗口plot(price,.); % 绘制日收盘价散点图hold onplot(output3,k,LineWidth,2); % 绘制平滑后曲线图，黑色实线，线宽为2plot(output4,k-.,LineWidth,2); % 绘制平滑后曲线图，黑色点划线，线宽为2xlabel(观测序号); ylabel(Gaussian window method); % 为X轴和Y轴加标签legend(原始散点,平滑曲线(窗宽30，标准差0.65),. 平滑曲线(窗宽100，标准差100),location,northwest);%*利用指数法对数据进行平滑处理，绘制平滑波形图*output5 = smoothts(price,e,30); % 用指数法平滑数据，窗宽为30output6 = smoothts(price,e,100); % 用指数法平滑数据，窗宽为100figure; % 新建一个图形窗口plot(price,.); % 绘制日收盘价散点图hold onplot(output5,k,LineWidth,2); % 绘制平滑后曲线图，黑色实线，线宽为2plot(output6,k-.,LineWidth,2); % 绘制平滑后曲线图，黑色点划线，线宽为2xlabel(观测序号); ylabel(Exponential method); % 为X轴和Y轴加标签legend(原始散点,平滑曲线(窗宽30),平滑曲线(窗宽100),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); % 指