一维加噪信号的平滑处理（2）

%*****************利用盒子法对数据进行平滑处理，绘制平滑波形图****************
output1 = smoothts(price,'b',30);  % 用盒子法平滑数据，窗宽为30
output2 = smoothts(price,'b',100);  % 用盒子法平滑数据，窗宽为100
figure;  % 新建一个图形窗口
plot(price,'.');  % 绘制散点图
hold on
plot(output1,'k','LineWidth',2);  % 绘制平滑后曲线图，黑色实线，线宽为2
plot(output2,'k-.','LineWidth',2);  % 绘制平滑后曲线图，黑色点划线，线宽为2
xlabel('观测序号'); ylabel('Box method'); % 为X轴和Y轴加标签
% 为图形加标注框
legend('原始散点','平滑曲线(窗宽30)','平滑曲线(窗宽100)','location','northwest');

%****************利用高斯窗方法对数据进行平滑处理，绘制平滑波形图**************
output3 = smoothts(price,'g',30);  % 窗宽为30，标准差为默认值0.65
output4 = smoothts(price,'g',100,100);  % 窗宽为100，标准差为100
figure;  % 新建一个图形窗口
plot(price,'.');  % 绘制散点图
hold on
plot(output3,'k','LineWidth',2);  % 绘制平滑后曲线图，黑色实线，线宽为2
plot(output4,'k-.','LineWidth',2);  % 绘制平滑后曲线图，黑色点划线，线宽为2
xlabel('观测序号'); ylabel('Gaussian window method'); % 为X轴和Y轴加标签
legend('原始散点','平滑曲线(窗宽30，标准差0.65)',...
          '平滑曲线(窗宽100，标准差100)','location','northwest');

%*****************利用指数法对数据进行平滑处理，绘制平滑波形图****************
output5 = smoothts(price,'e',30);  % 用指数法平滑数据，窗宽为30
output6 = smoothts(price,'e',100);  % 用指数法平滑数据，窗宽为100
figure;  % 新建一个图形窗口
plot(price,'.');  % 绘制散点图
hold on
plot(output5,'k','LineWidth',2);  % 绘制平滑后曲线图，黑色实线，线宽为2
plot(output6,'k-.','LineWidth',2);  % 绘制平滑后曲线图，黑色点划线，线宽为2
xlabel('观测序号'); ylabel('Exponential method'); % 为X轴和Y轴加标签
legend('原始散点','平滑曲线(窗宽30)','平滑曲线(窗宽100)','location','northwest');