matlab数据平滑程序文件_第1页
matlab数据平滑程序文件_第2页
matlab数据平滑程序文件_第3页
matlab数据平滑程序文件_第4页
matlab数据平滑程序文件_第5页
已阅读5页,还剩6页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

1、. 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 轴加标签

2、 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); xlab

3、el(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:); .专业专注. .

4、 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 他们依次平滑后生成的曲线 .专业专注. . .专业专注. . .专业专注. . 二、七点平

5、滑法的 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)

6、/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)=(39.0*y(500)+8.0*y(499)-4.0*y(498)- 4.0*y(497)+1.0*y(496)+4*y(495)-2*y(494)/42.0; y=

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论