第一章案例計(jì)算示意Y=a+bx

假設(shè)已有5組實(shí)驗(yàn)數(shù)據(jù)（x,y)如上表，

若要計(jì)算a和b,只要打開(kāi)光盤(pán)，點(diǎn)擊“cip”;找到“各章程序及操作文件”，再打開(kāi)“第一章及一次擬合文件夾”，點(diǎn)擊“DEM一次擬合.vbp”，打開(kāi)計(jì)算程序如下：x12345y7891012Y=a+bx運(yùn)行該程序，出現(xiàn)界面如下：將實(shí)驗(yàn)次數(shù)改為“5”，點(diǎn)擊“開(kāi)始計(jì)算”出現(xiàn)以下界面：Y=a+bx根據(jù)提示，分別輸入5組實(shí)驗(yàn)數(shù)據(jù)，可得結(jié)果，如下圖：Y=a+bx根據(jù)計(jì)算結(jié)果可知，a=5.6,b=1.2,若x=1代入，則得y=6.8,有0.2的偏差。整體的平均絕對(duì)偏差為0.24。注意可以通過(guò)改變程序中以下語(yǔ)句（共有4句），可改變結(jié)果的小數(shù)點(diǎn)位數(shù)“A00=Int(a0*100000+0.5)/100000”原設(shè)置為“5”位小數(shù)點(diǎn)，若改為“6”，只需將該語(yǔ)句改為：“A00=Int(a0*1000000+0.5)/1000000”即可，其他3句也仿照。Y=a+bx若要計(jì)算P4中的例1-1，由于數(shù)據(jù)太多，如果仍像剛才一樣輸入，就可能輸錯(cuò)數(shù)據(jù)，也無(wú)法判斷數(shù)據(jù)輸入是否正確，可通過(guò)數(shù)據(jù)文件的形式輸入。具體方法如下：先將原程序中:Fori=1Tonx(i)=InputBox("x("&i&")=")y(i)=InputBox("y("&i&")=")Nexti改為：Open“shujv1.dat"ForInputAs#1Fori=1TonInput#1,ax,ay//需在前面增加：Dimax,ayx(i)=axy(i)=ayNextiClose#1Y=a+bx數(shù)據(jù)文件“shujv1.dat”可用“記事本”建立。Y=a+bx文件名為“shujv1.dat”，保存類(lèi)型為“所有文件”；可直接放在“vb”程序的目錄下，這樣可直接引用。Y=a+bx計(jì)算結(jié)果文件存放形式Y(jié)=a+bx^n

只要將原來(lái)一次擬合的程序中的一個(gè)語(yǔ)句加以修改即可，將原程序中：“x(i)=InputBox(”x(“&i&”)=“)””的后面增加一句：x(i)=(x(i))^n即可，(x(i)=ax^3);例1-5，P19，具體數(shù)據(jù)如下表：x-3-2-124y14.38.34.78.322.7Y=aebx

將擬合公式兩邊取對(duì)數(shù)，可得線性形式：lnY=lna+bx，程序修改如下：原y(i)=InputBox(“y(”&i&“)=”)后面，增加一句：

y(i)=log(y(i))原a0=(m*d-c*p)/(n*d-c^2),增加一句：

a0=exp(a0),即可。例1-6，P20Y=aXb將擬合公式兩邊取對(duì)數(shù)，可得線性形式：lnY=lna+blnX，程序修改如下：原X(i)=InputBox(“X(”&i&“)=”)后面，增加一句：

X(i)=log(X(i))

y(i)=InputBox(“y(”&i&“)=”)后面，增加一句：

y(i)=log(y(i))原a0=(m*d-c*p)/(n*d-c^2),增加一句：

a0=exp(a0),即可。P28,第8題。Y=aXbFori=1Ton

x(i)=Exp(x(i))y(i)=Exp(y(i))eer=eer+(A00*x(i)^A11-y(i))^2sd=sd+Abs(A00*x(i)^A11-y(i))Nexti二次擬合程序的應(yīng)用y=a0+a1x+a2x2

已知下表實(shí)驗(yàn)數(shù)據(jù)，用二次擬合。只要打開(kāi)光盤(pán)，點(diǎn)擊“cip”;找到“各章程序及操作文件”，再打開(kāi)“第一章及二次擬合文件夾”，點(diǎn)擊“DEM二次擬合.vbp”，打開(kāi)計(jì)算程序如下。

x123456y616325482116二次擬合程序的應(yīng)用將Open“dem.dat”ForInputAs#1改為Open“erci1.dat”ForInputAs#1,并將表上的數(shù)據(jù)利用上面介紹的記事本保存為erci1.dat，運(yùn)算可得結(jié)果。二次擬合程序的應(yīng)用點(diǎn)擊”開(kāi)始計(jì)算”，輸入”6”，回車(chē)，得下圖結(jié)果，有關(guān)小數(shù)點(diǎn)的問(wèn)題同一次擬合。二次擬合程序的應(yīng)用

利用二次擬合程序，通過(guò)適當(dāng)?shù)淖儞Q，可以擬合許多類(lèi)型的實(shí)驗(yàn)數(shù)據(jù)，在程序中的修改工作量也非常小，可以非常方便使用，當(dāng)然目前也有許多軟件可以勝任該工作，但作者覺(jué)得，與其化很大的精力去找到擬合函數(shù)，還不如自己構(gòu)建一個(gè)函數(shù)來(lái)得快，下面我們通過(guò)幾個(gè)案例來(lái)說(shuō)明。y=a0+a1xn+a2xm

上面這個(gè)擬合公式，只要實(shí)驗(yàn)數(shù)據(jù)在不同的次方下有意義，就可以非常方便地利用光盤(pán)上的程序。如擬合成：(erci2.dat)

y=a0+a1x0.5+a2x1.8數(shù)據(jù)為：只要“x1(i)=xxx2(i)=xx^2x123456y613.825415980

x1(i)=xx^0.5x2(i)=xx^1.8y=a0+a1x0.5+a2x1.8Nu=c1Rec2Prc3

只要將上式兩邊進(jìn)行對(duì)數(shù)運(yùn)算，就可以得到線性表達(dá)式：lnNu=lnc1+c2lnRe+c3lnPr以P11例1-3的數(shù)據(jù)為例，建立erci4.dat數(shù)據(jù)文件，程序的修改工作如下：（dimxx1)Input#1,xx,YYx1(i)=xxx2(i)=xx^2y(i)=YY

Input#1,xx,xx1,YYx1(i)=Log(xx)x2(i)=Log(xx1)y(i)=Log(YY)Nu=c1Rec2Prc3C1=exp(-3.77255)=0.023;C2=0.80003C3=0.30003注意必須對(duì)常數(shù)項(xiàng)數(shù)據(jù)進(jìn)行變換處理。y=a0+a1x1

+a2x2

已知如下表實(shí)驗(yàn)數(shù)據(jù)：Open"erci5.dat"ForInputAs#1Fori=1TomInput#1,xx1,xx2,YYx1(i)=xx1x2(i)=xx2y(i)=YYNexti建立erci5.dat數(shù)據(jù)文件，修改部分見(jiàn)右邊，計(jì)算結(jié)果如下：x1123456x24861221y122118311211單變量任意次擬合擬合公式：計(jì)算公式單變量任意次擬合

已知下列數(shù)據(jù)，請(qǐng)分別用2次，3次，4次，5次擬合。(假設(shè)已建立renci1.dat）調(diào)用光盤(pán)第一章的”任意次擬合中DEMN.frm”,將”m=7”,改為”m=6”,將Open“g:dem.dat”ForInputAs#1改為Open“renci1.dat"ForInputAs#1x4861221y57.87221.1124.8526.117.096.01

數(shù)據(jù)按下式構(gòu)建：1+2*x+3*x^2+0.01x^3+0.0001*x^4+0.0002*x^5單變量任意次擬合點(diǎn)擊運(yùn)行，分別輸入2，輸入3，輸入4，輸入5，可以得到左圖的計(jì)算結(jié)果，當(dāng)輸入5時(shí)，得到的結(jié)果和構(gòu)建的公式相同。1+2*x+3*x^2+0.01x^3+0.0001*x^4+0.0002*x^5單變量任意次擬合

該程序沒(méi)有太多的變型，只能將x變成sinx或cosx或logx等簡(jiǎn)單的變化，在化工應(yīng)用中，某些特殊的場(chǎng)合，強(qiáng)制要求A