>>t=linspace(0,2*pi,500)';
%产生 500 行 1 列的服从 N(0,15*15)分布的随机数,作为噪声信号
>> noise=normrnd(0,15,500,1);
%新建一个图形窗口
%绘制加噪波形图
%为 X 轴加标签
>> ylabel('y=sin(t)+噪声');
>>
>> yy1=smooth(y,30);
>> figure;
%利用移动平均法对 y 做平滑处理
%新建一个图形窗口
>> plot(t,y,'k:');
>> hold on;
>> plot(t,yy1,'k','linewidth',3);
>> xlable('t');
>> ylabel('moving');
>> legend('加噪波形','平滑后波形');
>>
>> yy2=smooth(y,30,'lowess');
>> 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');
>> figure;
>> plot(t,y,'k:');
>> hold on;
>> plot(t,yy3,'k','linewidth',3);
>> xlabel('t');
>> ylabel('rlowess');
>> legend('加噪波形','平滑后波形');
>>
>> yy4=smooth(y,30,'loess');
>>figure;
>>plot(t,y,'k:');
>>hold on;
>>plot(t,yy4,'k','linewidth',3);
评论0
最新资源