matlab數據平滑程序_第1頁
matlab數據平滑程序_第2頁
matlab數據平滑程序_第3頁
matlab數據平滑程序_第4頁
matlab數據平滑程序_第5頁
已閱讀5頁,還剩4頁未讀 繼續免費閱讀

下載本文檔

版權說明:本文檔由用戶提供并上傳,收益歸屬內容提供方,若內容存在侵權,請進行舉報或認領

文檔簡介

1、matlab 數據平滑程序一、matlab自帶平滑函數clear,clctic,t=linspace(0,2*pi,500)'%產生一個從0到2*pi的向量,長度500y=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=smoo

2、th(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');pl

3、ot(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&#

4、39;,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('加噪波形','平

5、滑后波形'); 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,

6、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

7、(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=y1;endplot(t,y1,

溫馨提示

  • 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯系上傳者。文件的所有權益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙,網頁內容里面會有圖紙預覽,若沒有圖紙預覽就沒有圖紙。
  • 4. 未經權益所有人同意不得將文件中的內容挪作商業或盈利用途。
  • 5. 人人文庫網僅提供信息存儲空間,僅對用戶上傳內容的表現方式做保護處理,對用戶上傳分享的文檔內容本身不做任何修改或編輯,并不能對任何下載內容負責。
  • 6. 下載文件中如有侵權或不適當內容,請與我們聯系,我們立即糾正。
  • 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

評論

0/150

提交評論