>>t=linspace(0,2*pi,500)'; %产生一个从0到2*pi的向量,长度500 >> y=100*sin(t); %产生正弦波信号
%产生500行1列的服从N(0,15*15)分布的随机数,作为噪声信号
>> 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); %绘制平滑后波形图
>> xlable('t');
>> xlabel('t');
>> ylabel('moving');
>> legend('加噪波形','平滑后波形');
>>
>> yy2=smooth(y,30,'lowess'); %利用lowess方法y做平滑处理
>> figure;
>> plot(t,y,'k');
>> plot(t,y,'k:');
>> hold on;
>> plot(t,yy2,'k','linewidth',3);
>> xlabel('t');
>> ylabel('lowess');
>> legend('加噪波形','平滑后波形');
>>
>> yy3=smooth(y,30,'rlowess'); %利用rlowess方法y做平滑处理
>> figure;
>> plot(t,y,'k:');
>> hold on;
>> plot(t,yy3,'k','linewidth',3);
>> xlabel('t');
>> ylabel('rlowess');
>> legend('加噪波形','平滑后波形');
>>
>> yy4=smooth(y,30,'loess'); %利用loess方法y做平滑处理
>>figure;
>>plot(t,y,'k:');
>>hold on;
>>plot(t,yy4,'k','linewidth',3);