概述
一个小朋友问到这个一维信号的平滑处理方法的时候,我就用了中值平滑做的效果不是很好,师兄推荐的模板方法也不错,后来自己学习的过程中,有一段代码很好,就加以贴出来,跟大家共享一下。
%*****************产生加噪正弦波信号,绘制加噪波形图*************************
t = linspace(0,2*pi,500)'; % 产生一个从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轴加标签

%*************利用移动平均法对加噪信号进行平滑处理,绘制平滑波形图*************
yy1 = smooth(y,30); % 利用移动平均法对y进行平滑处理
figure; % 新建一个图形窗口
plot(t,y,'k:'); % 绘制加噪波形图
hold on;
plot(t,yy1,'k','linewidth',3); % 绘制平滑后波形图
xlabel('t'); % 为X轴加标签
ylabel('moving'); % 为Y轴加标签
legend('加噪波形','平滑后波形');

%*************利用lowess方法对加噪信号进行平滑处理,绘制平滑波形图*************
yy2 = smooth(y,30,'lowess'); % 利用lowess方法对y进行平滑处理
figure; % 新建一个图形窗口
plot(t,y,'k:'); % 绘制加噪波形图
hold on;
plot(t,yy2,'k','linewidth',3); % 绘制平滑后波形图
xlabel('t'); % 为X轴加标签
ylabel('lowess'); % 为Y轴加标签
legend('加噪波形','平滑后波形');

%*************利用rlowess方法对加噪信号进行平滑处理,绘制平滑波形图************
yy3 = smooth(y,30,'rlowess'); % 利用rlowess方法对y进行平滑处理
figure; % 新建一个图形窗口
plot(t,y,'k:'); % 绘制加噪波形图
hold on;
plot(t,yy3,'k','linewidth',3); % 绘制平滑后波形图
xlabel('t'); % 为X轴加标签
ylabel('rlowess'); % 为Y轴加标签
legend('加噪波形','平滑后波形');

%*************利用loess方法对加噪信号进行平滑处理,绘制平滑波形图*************
yy4 = smooth(y,30,'loess'); % 利用loess方法对y进行平滑处理
figure; % 新建一个图形窗口
plot(t,y,'k:'); % 绘制加噪波形图
hold on;
plot(t,yy4,'k','linewidth',3); % 绘制平滑后波形图
xlabel('t'); % 为X轴加标签
ylabel('loess'); % 为Y轴加标签
legend('加噪波形','平滑后波形');

%*************利用sgolay方法对加噪信号进行平滑处理,绘制平滑波形图*************
yy5 = smooth(y,30,'sgolay',3); % 利用sgolay方法对y进行平滑处理
figure; % 新建一个图形窗口
plot(t,y,'k:'); % 绘制加噪波形图
hold on;
plot(t,yy5,'k','linewidth',3); % 绘制平滑后波形图
xlabel('t'); % 为X轴加标签
ylabel('sgolay'); % 为Y轴加标签
legend('加噪波形','平滑后波形');

最后
以上就是傲娇萝莉为你收集整理的一维加噪信号的平滑处理(1)的全部内容,希望文章能够帮你解决一维加噪信号的平滑处理(1)所遇到的程序开发问题。
如果觉得靠谱客网站的内容还不错,欢迎将靠谱客网站推荐给程序员好友。
本图文内容来源于网友提供,作为学习参考使用,或来自网络收集整理,版权属于原作者所有。
发表评论 取消回复