>>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);
评论2