數(shù)據(jù)預(yù)處理的MATLAB實(shí)現(xiàn)

1、% 調(diào)用 smooth 函數(shù)進(jìn)行加噪數(shù)據(jù)的平滑處理% %* 產(chǎn)生加噪正弦波信號(hào)，繪制加噪波形圖 * t = linspace(0,2*pi,500); % 產(chǎn)生一個(gè)從 0 到 2*pi 的向量，長(zhǎng)度為 500 y = 100*sin(t); % 產(chǎn)生正弦波信號(hào)% 產(chǎn)生 500 行 1 列的服從 N(0,152) 分布的隨機(jī)數(shù)，作為噪聲信號(hào)noise = normrnd(0,15,500,1);y = y + noise; % 將正弦波信號(hào)加入噪聲信號(hào) figure; % 新建一個(gè)圖形窗口 plot(t,y); % 繪制加噪波形圖 xlabel(t); % 為 X 軸加標(biāo)簽ylabel(y =

2、sin(t) + 噪聲);% 為 Y 軸加標(biāo)簽%* 利用移動(dòng)平均法對(duì)加噪信號(hào)進(jìn)行平滑處理， 繪制平滑波形圖 * yy1 = smooth(y,30); % 利用移動(dòng)平均法對(duì) y 進(jìn)行平滑處理figure; % 新建一個(gè)圖形窗口 plot(t,y,k:);% 繪制加噪波形圖hold on;plot(t,yy1,k,linewidth,3); % 繪制平滑后波形圖xlabel(t); % 為 X 軸加標(biāo)簽 ylabel(moving); % 為 Y 軸加標(biāo)簽 legend(加噪波形，平滑后波形);%* 利用 lowess 方法對(duì)加噪信號(hào)進(jìn)行平滑處理， 繪制平滑波形圖 * yy2 = smooth(

3、y,30,lowess); % 利用 lowess 方法對(duì) y 進(jìn)行平滑處理figure; % 新建一個(gè)圖形窗口plot(t,y,k:);% 繪制加噪波形圖hold on;plot(t,yy2,k,linewidth,3); % 繪制平滑后波形圖xlabel(t); % 為 X 軸加標(biāo)簽 ylabel(lowess); % 為 Y 軸加標(biāo)簽 legend(加噪波形，平滑后波形);%* 利用 rlowess 方法對(duì)加噪信號(hào)進(jìn)行平滑處理， 繪制平滑波形圖 *yy3 = smooth(y,30,rlowess); % 利用 rlowess 方法對(duì) y 進(jìn)行平滑處理figure; % 新建一個(gè)圖形窗口

4、plot(t,y,k:);% 繪制加噪波形圖hold on;plot(t,yy3,k,linewidth,3); % 繪制平滑后波形圖xlabel(t); % 為 X 軸加標(biāo)簽ylabel(rlowess); lege nd (加噪波形% 為 Y 軸加標(biāo)簽，平滑后波形 ）;%*利用 loess 方法對(duì)加噪信號(hào)進(jìn)行平滑處理，繪制平滑波形圖*yy4 = smooth(y,30,loess); % 利用 loess 方法對(duì) y 進(jìn)行平滑處理 figure; % 新建一個(gè)圖形窗口 plot(t,y,k:); % 繪制加噪波形圖 hold on;plot(t,yy4,k,linewidth,3); %

5、繪制平滑后波形圖 xlabel(t); % 為 X 軸加標(biāo)簽 ylabel(loess); % 為 Y 軸加標(biāo)簽lege nd （加噪波形，平滑后波形 ）;%*利用 sgolay 方法對(duì)加噪信號(hào)進(jìn)行平滑處理， 繪制平滑波形圖 *yy5 = smooth(y,30,sgolay,3); % 利用 sgolay 方法對(duì) y 進(jìn)行平滑處理 figure; % 新建一個(gè)圖形窗口 plot(t,y,k:); % 繪制加噪波形圖hold on; plot(t,yy5,k,linewidth,3); % 繪制平滑后波形圖 xlabel(t); % 為 X 軸加標(biāo)簽 ylabel(sgolay); % 為 Y

6、 軸加標(biāo)簽legend(加噪波形，平滑后波形);%調(diào)用 smoothts 函數(shù)對(duì)金融時(shí)間序列數(shù)據(jù)進(jìn)行平滑處理%*讀取金融時(shí)間序列數(shù)據(jù)，繪制波形圖 *x = xlsread（examp03_02.xls）; % 從文件 examp03_02.xls 中讀取數(shù)據(jù) price = x（:，4）;% 提取矩陣 x 的第 4 列數(shù)據(jù)，即收盤(pán)價(jià)數(shù)據(jù)figure; % 新建一個(gè)圖形窗口*plot（price，k，LineWidth，2）; % 繪制日收盤(pán)價(jià)曲線(xiàn)圖，黑色實(shí)線(xiàn)，線(xiàn)寬為 2 xlabel（觀測(cè)序號(hào)）;ylabel（上海股市日收盤(pán)價(jià)）;%為X軸和Y軸加標(biāo)簽%* 利用盒子法對(duì)數(shù)據(jù)進(jìn)行平滑處理，繪制平滑

7、波形圖output1 = smoothts（price，b，30）; % 用盒子法平滑數(shù)據(jù)，窗寬為 30output2 = smoothts（price，b，100）; % 用盒子法平滑數(shù)據(jù)，窗寬為 100figure; % 新建一個(gè)圖形窗口 plot(price,.); % 繪制日收盤(pán)價(jià)散點(diǎn)圖 hold on plot(output1,k,LineWidth,2); % 繪制平滑后曲線(xiàn)圖，黑色實(shí)線(xiàn)，線(xiàn)寬為 2plot(output2,k-.,LineWidth,2); % 繪制平滑后曲線(xiàn)圖，黑色點(diǎn)劃線(xiàn)，線(xiàn)寬為 2xlabel(觀測(cè)序號(hào));ylabel(Box method); % 為 X 軸

8、和 Y 軸加標(biāo)簽% 為圖形加標(biāo)注框legend(原始散點(diǎn)，平滑曲線(xiàn)(窗寬30),平滑曲線(xiàn)(窗寬100),location,northwest);%* 利用高斯窗方法對(duì)數(shù)據(jù)進(jìn)行平滑處理， 繪制平滑波形圖 *output3 = smoothts(price,g,30); % 窗寬為 30，標(biāo)準(zhǔn)差為默認(rèn)值 0.65output4 = smoothts(price,g,100,100); % 窗寬為 100，標(biāo)準(zhǔn)差為 100figure; % 新建一個(gè)圖形窗口plot(price,.); % 繪制日收盤(pán)價(jià)散點(diǎn)圖hold onplot(output3,k,LineWidth,2); % 繪制平滑后曲線(xiàn)圖

9、，黑色實(shí)線(xiàn)，線(xiàn)寬為 2 plot(output4,k-.,LineWidth,2); % 繪制平滑后曲線(xiàn)圖，黑色點(diǎn)劃線(xiàn)，線(xiàn)寬為 2 xlabel( 觀測(cè)序號(hào) ); ylabel(Gaussian window method); % 為 X 軸和 Y 軸加標(biāo)簽 lege nd(原始散點(diǎn)，平滑曲線(xiàn)(窗寬30,標(biāo)準(zhǔn)差0.65),*平滑曲線(xiàn) (窗寬 100，標(biāo)準(zhǔn)差 100),location,northwest);%* 利用指數(shù)法對(duì)數(shù)據(jù)進(jìn)行平滑處理，繪制平滑波形圖output5 = smoothts(price,e,30); % 用指數(shù)法平滑數(shù)據(jù)，窗寬為 30 output6 = smoothts(p

10、rice,e,100); % 用指數(shù)法平滑數(shù)據(jù)，窗寬為 100 figure; % 新建一個(gè)圖形窗口plot(price,.); % 繪制日收盤(pán)價(jià)散點(diǎn)圖hold on plot(output5,k,LineWidth,2); % 繪制平滑后曲線(xiàn)圖，黑色實(shí)線(xiàn)，線(xiàn)寬為 2 plot(output6,k-.,LineWidth,2); % 繪制平滑后曲線(xiàn)圖，黑色點(diǎn)劃線(xiàn)，線(xiàn)寬為 2 xlabel（ 觀測(cè)序號(hào) ）; ylabel（Exponential method）; % 為 X 軸和 Y 軸加標(biāo)簽 legend（原始散點(diǎn)，平滑曲線(xiàn)（窗寬30）,平滑曲線(xiàn)（窗寬100）,location,northwe

11、st）;%調(diào)用 medfilt1 函數(shù)對(duì)加噪正弦波信號(hào)進(jìn)行平滑處理% %* 產(chǎn)生加噪正弦波信號(hào)，繪制加噪波形圖 *t = linspace(0,2*pi,500);% 產(chǎn)生一個(gè)從 0 到 2*pi 的向量，長(zhǎng)度為 500y = 100*sin(t);% 產(chǎn)生正弦波信號(hào)%產(chǎn)生500行1列的服從N(0,152)分布的隨機(jī)數(shù)，作為噪聲信號(hào)noise = normrnd(0,15,500,1);y = y + noise; % 將正弦波信號(hào)加入噪聲信號(hào) figure; % 新建一個(gè)圖形窗口 plot(t,y); % 繪制加噪波形圖 xlabel(t); % 為 X 軸加標(biāo)簽ylabel(y = sin(t) + 噪聲);% 為 Y 軸加標(biāo)簽%* 調(diào)用 medfilt1 對(duì)加噪正弦波信號(hào) y 進(jìn)行中值濾波， 并繪制波形圖 * yy = medfilt1(y,30)