最小二乘法數(shù)值計(jì)算實(shí)驗(yàn)報(bào)告.doc

數(shù)學(xué)與計(jì)算科學(xué)學(xué)院實(shí) 驗(yàn) 報(bào) 告實(shí)驗(yàn)項(xiàng)目名稱(chēng) 最小二乘多項(xiàng)式擬合 所屬課程名稱(chēng) 數(shù)值計(jì)算 實(shí) 驗(yàn) 類(lèi) 型 驗(yàn)證型 實(shí) 驗(yàn) 日 期 5.8.2012 班 級(jí) 隧道1002班 學(xué) 號(hào) 201008020233 姓 名 李彬彬 成 績(jī) 一、實(shí)驗(yàn)概述：【實(shí)驗(yàn)?zāi)康摹?通過(guò)上機(jī)計(jì)算，對(duì)曲線的最小二乘法的擬合有進(jìn)一步的掌握，并且能夠熟練的運(yùn)用這種方法。【實(shí)驗(yàn)原理】在科學(xué)實(shí)驗(yàn)數(shù)據(jù)處理中，往往要根據(jù)一組給定的實(shí)驗(yàn)數(shù)據(jù)，求出自變量x與因變量y的函數(shù)關(guān)系，這是為待定參數(shù)，由于觀測(cè)數(shù)據(jù)總有誤差，且待定參數(shù)ai的數(shù)量比給定數(shù)據(jù)點(diǎn)的數(shù)量少(即nm)，因此它不同于插值問(wèn)題.這類(lèi)問(wèn)題不要求通過(guò)點(diǎn)，而只要求在給定點(diǎn)上的誤差的平方和最小.當(dāng)時(shí)，即(4.4.1)這里是線性無(wú)關(guān)的函數(shù)族，假定在上給出一組數(shù)據(jù)，以及對(duì)應(yīng)的一組權(quán)，這里為權(quán)系數(shù)，要求使最小，其中這就是最小二乘逼近，得到的擬合曲線為y=s(x)，這種方法稱(chēng)為曲線擬合的最小二乘法.(4.4.2)中實(shí)際上是關(guān)于的多元函數(shù)，求I的最小值就是求多元函數(shù)I的極值，由極值必要條件，可得(4.4.3)根據(jù)內(nèi)積定義(見(jiàn)第三章)引入相應(yīng)帶權(quán)內(nèi)積記號(hào)(4.4.4)則(4.4.3)可改寫(xiě)為這是關(guān)于參數(shù)的線性方程組，用矩陣表示為(4.4.5)(4.4.5)稱(chēng)為法方程.當(dāng)線性無(wú)關(guān)，且在點(diǎn)集上至多只有n個(gè)不同零點(diǎn)，則稱(chēng)在X上滿(mǎn)足Haar條件，此時(shí)(4.4.5)的解存在唯一(證明見(jiàn)3).記(4.4.5)的解為 從而得到最小二乘擬合曲線(4.4.6)可以證明對(duì)，有故(4.4.6)得到的即為所求的最小二乘解.它的平方誤差為(4.4.7)均方誤差為在最小二乘逼近中，若取，則，表示為(4.4.8) 【實(shí)驗(yàn)環(huán)境】Microsoft visual c+ 二、實(shí)驗(yàn)內(nèi)容：【實(shí)驗(yàn)方案】測(cè)得銅導(dǎo)線在溫度Ti（）時(shí)的電阻Ri（）如表，求電阻R與溫度T的近似函數(shù)關(guān)系。i0123456Ti（）19.125.030.136.040.045.150.0Ri（）76.3077.8079.2580.8082.3583.9085.10用計(jì)算機(jī)程序進(jìn)行多項(xiàng)式擬合上述T與R的近似函數(shù)關(guān)系?！緦?shí)驗(yàn)過(guò)程】（實(shí)驗(yàn)步驟、記錄、數(shù)據(jù)、分析）1， 確定數(shù)據(jù)，易知測(cè)得的數(shù)據(jù)接近一條直線，故擬合函數(shù)為R=a0+a1T；2， 編寫(xiě)計(jì)算機(jī)程序；3， 運(yùn)行程序，得出結(jié)果；4， 計(jì)算出的擬合多項(xiàng)式函數(shù)計(jì)算出的數(shù)據(jù)與原數(shù)據(jù)進(jìn)行比較【實(shí)驗(yàn)結(jié)論】（結(jié)果）【實(shí)驗(yàn)小結(jié)】（收獲體會(huì)）通過(guò)本實(shí)驗(yàn)使用多項(xiàng)式對(duì)數(shù)據(jù)進(jìn)行擬合，我掌握數(shù)據(jù)擬合的基本原理，并且掌握最小二乘法的計(jì)算方法，同時(shí)學(xué)會(huì)使用數(shù)學(xué)的方法對(duì)數(shù)據(jù)擬合的情況進(jìn)行判斷。對(duì)我們以后對(duì)數(shù)據(jù)進(jìn)行分析很有幫助三、指導(dǎo)教師評(píng)語(yǔ)及成績(jī)：評(píng) 語(yǔ)評(píng)語(yǔ)等級(jí)優(yōu)良中及格不及格1.實(shí)驗(yàn)報(bào)告按時(shí)完成,字跡清楚,文字?jǐn)⑹隽鲿?邏輯性強(qiáng)2.實(shí)驗(yàn)方案設(shè)計(jì)合理3.實(shí)驗(yàn)過(guò)程（實(shí)驗(yàn)步驟詳細(xì),記錄完整,數(shù)據(jù)合理,分析透徹）4實(shí)驗(yàn)結(jié)論正確. 成 績(jī)： 指導(dǎo)教師簽名： 批閱日期：附錄1：源 程 序#includeint main() double a0,a1,sumTi=0,sumTi2=0,sumRi=0,TiRi=0; double T100,R100; int n,m,i,j; printf(輸入數(shù)據(jù)的個(gè)數(shù)n:); scanf(%d,&n); printf(輸入各組數(shù)據(jù)Ti與Ri:n); m=n; i=0;j=0; while(n-) scanf(%lf%lf,&Ti,&Ri+); printf(%dn,i); for(j=0;ji;j+) sumTi+=Tj; sumRi+=Rj; sumTi2+=Tj*Tj; TiRi+=Tj*Rj; a0=-(sumTi*TiRi-sumTi2*sumRi)/(m*sumTi2-sumTi*sumTi); a1=(m*TiRi-sumTi*sumRi)/(m*sumTi2-sumTi*sumTi); p