正弦函數(shù)擬合計(jì)算_第1頁
正弦函數(shù)擬合計(jì)算_第2頁
正弦函數(shù)擬合計(jì)算_第3頁
正弦函數(shù)擬合計(jì)算_第4頁
正弦函數(shù)擬合計(jì)算_第5頁
已閱讀5頁,還剩1頁未讀 繼續(xù)免費(fèi)閱讀

下載本文檔

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

文檔簡(jiǎn)介

1、正弦函數(shù)擬合計(jì)算、正弦函數(shù)的一般表達(dá)式的建立正弦函數(shù)的一般表達(dá)式為:y = x0si n(%t x2) + X3(1)對(duì)于一系列的n個(gè)點(diǎn)(n _ 3):(ti,yj,i =0,1廠,n -1( 2)要用點(diǎn)(ti,yj,i =0,1, ,n -1擬合計(jì)算上述方程,則使:n 二S - ' XoSin(X1ti X2) X3 - yi 2i =0最小。PS要使得S最小,應(yīng)滿足:S =0, k =0,1,2,3%即::x0:S:X1%cS芒X3=x2X0Sin(X1tj X2) X3 - yi 】sin(X1tj x?)2ksin(X1tiX2) X3 - yi ktj cosgti x?)2

2、x°s in( X1ti X2) X3 yX0 cos(X1tiX2)2X0Sin(X1tjX2) X3 - yi 丨X0 =0二X°sin%ti+X2)+ X3 yibi nx(ti+X2) = 0瓦X°sinW+X2)+ X3 yi.ti cos(ti +x2)= 0瓦X0sinXI+x2)+ x3 yito s<ti+ x2) = 0ZX0sinX11+ x2)+ x3 -yi】=0解上述4元非線性方程組,即可得到正弦函數(shù)的一般表達(dá)式的系數(shù):X0,X1,X2,X3。、多元非線性方程組解法對(duì)于n元非線性方程組,記:F X - lf°(X), f

3、i(X), , fz(X)lT , X - lx°,Xi, ,Xn以及雅克比矩陣F'(X) =:%(X)CXo吋1(X)cXo吋2(X) I cXo 心)IL :'Xo拼0(X)計(jì)o(X)?。╔)図cX2陜吋1(X)吋1(X)拼1(X)cX1cX2CX3拼2(X)吋2(X)拼2(X)cXiCX2CX3吋3(X)吋3(X)質(zhì)3(X)cX1cX2CX3像單個(gè)方程的Newton迭代法-櫛,采用逐次線性化的方法構(gòu)造方程(1)的Newton迭代法*在某個(gè)近似解X.處,將向盤函數(shù)F(X) Taylor JR開,則仃:從而得到(刃的近似方程F(兀)+玖兀)(心廠血)=0令朋上二區(qū)田

4、-X- 將Newton法的迭代公式改寫為:rxk+l=xt+AxtpttXAX1£) = -FQCl)毎-步迭代均盅解Newdon方程俎;Fr(XtXAXt) = -F(XJ這是個(gè)線性方程組,可用高斯消去袪求解.即:fo(X)巒0(X)&(X)蠲(X)故0做區(qū)(X)並(X)血0血1龜(X)f3(X)-血0£f°(X)£fo(X)1CX2誦(X)誦(X)CX2Ef2(X)£f2(X)CX2西(X)西(X)cx2啟一3J fo(X)Az、-f1(X)Ax2-f2(X)0X3 )<-f3(X)(5)三、正弦函數(shù)的一般表達(dá)式系數(shù)求解要擬合

5、正弦函數(shù)的一般表達(dá)式(1)的系數(shù),線性方程組(5)中的表達(dá)式為:,o(X)=送 X°sin(Xiti +X2)+ X3 yi kin(x1ti +x?) f(X)=送 RoSinddi + x2) + x3 yi E cos(x1t x2) f2(X)八Xosin(x1tix2px3-yi lcos(x1tix2)f3(X)二 x°sin(xitiX2)X3-yi 1Wo(X)寸.2/*、=S sin (xiti +X2)&o旨o(jì)(X) =£ Ixoti sin 2(xiti +X2)】中瓦(冷 一yjti cos(xiti + X2)】"o(X

6、)=送 Ixosin2(xdi+X2)】+送 収3 W)cos(xdi+x?)】&2書o(X) 丁:=送 sin(Xiti +X2)&3.:fi(X) i-ti sin 2(xiti x2):Xo2fi(X)八 Xocos2(xiti X2) L ' (x3-yjt'si n(xiti X2):xifi(X)八 lxoticos2(xiti x2 (xyi)tisin(捲1 x2) 1:x2fi(X) = ti cos(xiti x2):X3二丄、sin2(Xiti X2)lx。2_ ;(X)= Roti cos2(xiti +X2)】一瓦 kyJtiSin(xiti+x2)】:xi一f2(X)' Xo cos2(Xiti X2)-' (X3-yjs in (x-I x?)】:x2-f2(X)= y cos(x-ti X2)f3(X)CXg f3(X) *做 f3(X)cX2f3(X).CX3X3-、sin(Mx2)-'J Xoti cos(x1ti x2)】八 Xo cos(x1tiX2)=n根據(jù)前面所述的Newton迭代法,先給出xo,x1,x2,x3的初值Xo,代入公式(5)求得:=Xk - L% LX1 LX2:LX3 fk = 0,1,2為

溫馨提示

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

最新文檔

評(píng)論

0/150

提交評(píng)論