matlab数据平滑程序文件.docxVIP

.. .. .. Matlab 数据平滑程序 一、matlab 自带平滑函数 clear,clc tic,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,50); %利用移动平均法对 y 做平 滑处理 . 专业专注 . .. .. .. figure; %新建一个图形窗口 plot(t,y,k:); %绘制加噪波形图 hold on; plot(t,yy1,k,linewidth,3); %绘制平滑后波形图 xlabel(t); ylabel(moving); legend( 加噪波形 ,平滑后波形 ); yy2=smooth(y,50,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,50,rlowess); %利用 rlowess 方法 y 做平 滑处理 figure; plot(t,y,k:); hold on; plot(t,yy3,k,linewidth,3); xlabel(t); ylabel(rlowess); legend( 加噪波形 , 平滑后波形 ); yy4=smooth(y,50,loess); %利用 loess 方法 y 做平滑 处理 figure; plot(t,y,k:); . 专业专注 . .. .. .. hold on; plot(t,yy4,k,linewidth,3); xlabel(t); ylabel(loess); legend( 加噪波形 ,平滑后波形 ); yy5=smooth(y,50,sgolay,3); %利用 sgolay 方法 y 做平滑处理 figure; plot(t,y,k:); hold on; plot(t,yy5,k,linewidth,3); xlabel(t); ylabel(sgolay); legend( 加噪波形 ,平滑后波形 ); . 专业专注 . .. .. .. toc 他们依次平滑后生成的曲线 . 专业专注 . .. .. .. . 专业专注 . .. .. .. . 专业专注 . .. .. .. 二、七点平滑法的 matlab 数据平滑 clear,clc tic,t=linspace(0,2*pi,500); y=100*sin(t); noise=normrnd(0,15,500,1); y=y+noise; . 专业专注 . .. .. .. figure; plot(t,y,k:); hold on; for j=0:10000 y1(1)=(39.0*y(1)+8.0*y(2)-4.0*y(3)-4.0*y(4)+1.0*y(5)+4*y(6)-2*y(7))/42.0; y1(2)=(8.0*y(2)+19.0*y(3)+16.0*y(4)+6.0*y(5)+4.0*y(6)-7*y(7)+4*y(8))/42.0; y1(3)=(-4.0*y(3)+16.0*y(4)+19.0*y(5)+12.0*y(6)+2.0*y(7)- 4*y(8)+1*y(9))/42.0; for i=4:497 y1(i)=(-2*(y(i-3)+y(i+3))+3.0*(y(i-2)+1*y(i+2))+6.0*(y(i- 1)+1*y(i+1))+7*y(i))/21.0; end y1(498)=(-4.0*y(498)+16.0*y(497)+19.0*y(496)+12.0*y(495)+2.0*y(494)- 4*y(493)+1*y(492))/42.0; y1(499)=(8.0*y(499)+19.0*y(498)+16.0*y(497)+6.0*y(496)+4.0*y(495)- 7*y(494)+4*y(493))/42.0; . 专业专注 . .. .. .. y1(500)