數(shù)據(jù)擬合實驗舉例程序_第1頁
數(shù)據(jù)擬合實驗舉例程序_第2頁
數(shù)據(jù)擬合實驗舉例程序_第3頁
數(shù)據(jù)擬合實驗舉例程序_第4頁
全文預(yù)覽已結(jié)束

下載本文檔

版權(quán)說明：本文檔由用戶提供并上傳，收益歸屬內(nèi)容提供方，若內(nèi)容存在侵權(quán)，請進(jìn)行舉報或認(rèn)領(lǐng)

文檔簡介

1、實驗二：曲線擬合目的與要求：認(rèn)識最小二乘法的基來源理，用最小二乘法求擬合數(shù)據(jù)的多項式，做出失散函數(shù)和擬合函數(shù)的圖形，掌握利用最小二乘法進(jìn)行數(shù)據(jù)擬合的基本思想，熟習(xí)找尋最正確方法擬合曲線的方法，經(jīng)過計算機(jī)解決實驗問題例題1已知觀察數(shù)據(jù)x01y對實驗一中的例題2進(jìn)行曲線擬合x=0:1;y=;plot(x,y,k.,markersize,25)axis(0-216)p=polyfit(x,y,3);p1=polyfit(x,y,6);t=0:;s=polyval(p,t);s1=polyval(p1,t);holdonplot(t,s,k-,linewidth,2)plot(t,s,k-,linew

2、idth,2)grid;例題2已知觀察數(shù)據(jù)x1345678910y1054211234求一代數(shù)多項式曲線，使其最好地擬合這組給定數(shù)據(jù)。（1）畫出數(shù)據(jù)散布趨向圖xi=1345678910;yi=1054211234;plot(xi,yi,o)（2）成立數(shù)學(xué)模型y=a2x2+a1x+a0成立超定方程組系數(shù)矩陣A=ones(size(xi);xi;xi.2A=111139141615251636174918641981101003）求超定方程組的最小二乘解a=Ayia=（4）求擬合曲線方程b=b=y=poly2str(b,x)y=x2-x+f2=polyval(flipud(a),xi);plot(

3、xi,yi,bo,xi,f2,r-)5）用方程y=axb擬合x=ones(size(xi);log(xi);aa=xlog(yi)aa=yy=exp*xi.;yy=exp*xi.;plot(xi,yi,bo,xi,yy,r-,xi,f2,b-)例題3已知觀察數(shù)據(jù)xy用polyfit命令作擬合直線并計算誤差。fprintf(ix(i)y(i)y=c_1x+c_2誤差n);fori=1:6x=;y=;c=polyfit(x,y,1);error=y(i)-(c(1)*x(i)+c(2);fprintf(%2d%2.1f%5.2f%10.4f%10.5fn,i,x(i)end,y(i),c(1)*x(i)+c(2),error);ix(i)y(i)y=c_1x+c_2誤差123456任務(wù)：由化學(xué)實驗獲得某物質(zhì)濃度與時間的關(guān)系以下：時間t12345678濃度y時間t910111213141516濃度y求濃度

溫馨提示

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

最新文檔

評論

0/150

提交評論