数据预处理的MATLAB实现.docxVIP

  • 24
  • 0
  • 约4.26千字
  • 约 4页
  • 2021-01-31 发布于天津
  • 举报
% % 调用 smooth 函数进行加噪数据的平滑处理 % %***************** 产生加噪正弦波信号,绘制加噪波形图 ************************* 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); lege nd (加噪波形 % 为 Y 轴加标签 ,平滑后波形 ); %************* 利用 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 轴加标签 lege nd (加噪波形 ,平滑后波形 ); %************* 利用 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

文档评论(0)

1亿VIP精品文档

相关文档