参考链接：http://bbs.pinggu.org/thread-3777396-1-1.html

smoothts函数

调用格式：

output = smoothts(input)

output = smoothts(input, ‘b’, wsize) % 盒子法

output = smoothts(input, ‘g’, wsize, stdev) % 高斯窗方法

output = smoothts(input, ‘e’, n) % 指数法

【例7.1-2】现有上海股市日开盘价、最高价、最低价、收盘价、收益率等数据，时间跨度为2005年1月4日至2007年4月3日，共510组数据。完整数据保存在文件examp7_1_2.xls中，其中部分数据如下图所示。试调用smoothts函数对日收盘价数据进行平滑处理

绘制日收盘价曲线图：

% 从文件examp7_1_2.xls中读取数据

>> x = xlsread('examp7_1_2.xls');

>> price = x(:,4)'; % 提取矩阵x的第4列数据，即收盘价数据

>> figure; % 新建一个图形窗口% 绘制日收盘价曲线图，黑色实线，线宽为2

>> plot(price,'k','LineWidth',2); % 为X轴和Y轴加标签

>> xlabel('观测序号');

>> ylabel('上海股市日收盘价');

盒子法：

% 用盒子法平滑数据，窗宽为30

>> output1 = smoothts(price,'b',30);

% 用盒子法平滑数据，窗宽为100

>> output2 = smoothts(price,'b',100);

>> figure; % 新建一个图形窗口

>> plot(price,'.'); % 绘制日收盘价散点图

>> hold on

% 绘制平滑后曲线图，黑色实线，线宽为2

>> plot(output1,'k','LineWidth',2);

% 绘制平滑后曲线图，黑色点划线，线宽为2

>> plot(output2,'k-.','LineWidth',2);

% 为X轴和Y轴加标签

>> xlabel('观测序号'); ylabel('Box method');

% 为图形加标注框

>> 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 % 绘制平滑后曲线图，黑色实线，线宽为2

>> plot(output3,'k','LineWidth',2); % 绘制平滑后曲线图，黑色点划线，线宽为2

>> plot(output4,'k-.','LineWidth',2); % 为X轴和Y轴加标签>> xlabel('观测序号'); ylabel('Gaussian window method');

>> 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 % 绘制平滑后曲线图，黑色实线，线宽为2

>> plot(output5,'k','LineWidth',2); % 绘制平滑后曲线图，黑色点划线，线宽为2

>> plot(output6,'k-.','LineWidth',2); % 为X轴和Y轴加标签>> xlabel('观测序号'); ylabel('Exponential method');

>> legend('原始散点','平滑曲线(窗宽30)','平滑曲线(窗宽100)','location','northwest');

三、medfilt1函数(一维中值滤波)

调用格式：

y = medfilt1(x,n)

y = medfilt1(x,n,blksz)

y = medfilt1(x,n,blksz,dim)

【例7.1-3】产生一列正弦波信号，加入噪声信号，然后调用medfilt1函数对加入噪声的正弦波进行滤波(平滑处理)

% 产生一个从0到2*pi的向量，长度为500

>> t = linspace(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轴加标签

中值滤波：

% 调用medfilt1对加噪正弦波信号y进行中值滤波，并绘制波形图

>> yy = medfilt1(y,30); % 指定窗宽为30，对y进行中值滤波

>> figure; % 新建一个图形窗口>> plot(t,y,'k:'); % 绘制加噪波形图

>> hold on % 绘制平滑后曲线图，黑色实线，线宽为3

>> plot(t,yy,'k','LineWidth',3);

>> xlabel('t'); % 为X轴加标签

>> ylabel('中值滤波'); % 为Y轴加标签

>> legend('加噪波形','平滑后波形');