計算方法上機題doc

1、第八題:題目，基本方法從函數(shù)表X0.00.10.1950.30.4010.5f（x）：0.398940.39695 :0.391420.38138:0.368120.35206出發(fā)，用下列方法計算f(0.15),f(0.31),f(0.47)的近似值：(1)分段線性插值；(2)分段二次插值;(3)全區(qū)間上拉格朗日插值(要求一次性輸入整張函數(shù)表，并利用計算機選擇在插值計算中所需的結(jié)點)(1) 先選取兩個結(jié)點與X使xx，X然后在區(qū)間【Xi，Xi上作線性插值，即得 f(x) Li(x)= ypXl y?？赬iXiXiXd(2) 選取距點X最近的三個結(jié)點X4，Xi，X 1進行二次插值，即取i卅十 X

2、X.f(X) l(x)八H。k-L 丿 j j戎 Xk_Xj（3）n=5,nnLn（X）yklk（X）=7k =0（X _X0）（X _X2）（X _Xk i）（X _Xn）（Xk X0）（Xk Xk丿（Xk Xk i）（Xk Xn）三，計算結(jié)果（1）X0.150.310.47f（x）0.3940390.3800670.356927X0.150.310.47f（x）0.3944820.3802250.361399x0.150.310.47f(x)0.3944730.3802190.357222四，結(jié)果分析編程基本按照課本套公式，方法從(1)至(3),計算結(jié)果越來越精確 五，源程序(1)#i n

3、cludeint main(v oid)int i,j;double a6,b6,x,y;prin tf(E nter x:);for(i=0;i6;i+)scan f(%lf,&ai);prin tf(E nter f(x):);for(i=0;i6;i+)sca nf(%lf,&bi);for(j=0;jai;i+); y=bi-1*(x-ai)/(ai-1-ai)+bi*(x-ai-1)/(ai-ai-1); prin tf(x=%lf , ,x);prin tf(y=%lfn,y);return 0;(2)#i ncludeint main(v oid)int i,j;double a

4、6,b6,x,y;prin tf(E nter x:);for(i=0;i6;i+)scan f(%lf,&ai);prin tf(E nter f(x):);for(i=0;i6;i+)sea nf(%lf,&bi);for(j=0;jai;i+);if(ai+1-x)(x-ai-2) y=bi-1*(x-ai)/(ai-1-ai)*(x-ai+1)/(ai-1-ai+1); y=y+bi*(x-ai-1)/(ai-ai-1)*(x-ai+1)/(ai-ai+1); y=y+bi+1*(x-ai-1)/(ai+1-ai-1)*(x-ai)/(ai+1-ai);prin tf(x=%lf ,

5、,x); prin tf(y=%lfn,y);else y=bi-2*(x-ai-1)/(ai-2-ai-1)*(x-ai)/(ai-2-ai); y=y+bi-1*(x-ai-2)/(ai-1-ai-2)*(x-ai)/(ai-1-ai); y=y+bi*(x-ai-2)/(ai-ai-2)*(x-ai-1)/(ai-ai-1);prin tf(x=%lf , ,x); prin tf(y=%lfn,y);return 0;(3)#i ncludeint main(v oid)int i,j,k;double a6,b6,x,y,p;prin tf(E nter x:);for(i=0;i6

6、;i+)scan f(%lf,&ai);prin tf(E nter f(x):);for(i=0;i6;i+)sea nf(%lf,&bi);for(j=0;j3;j+)y=0;prin tf(E nter x:);sea nf(%lf, &x);for(i=0;i6;i+)p=1;for(k=0;k3)a5r：BiaF(tliF,l.Hii)Sfnm )-； J I1-J|iJ*l*(M-3|1|-31-1); priintfrK-lIF “聲prfcntirry-tiFW-yb:rcurri a：第AJE (1) II rrrgr(9j, II turnlnglfLC口h Duild 1

7、 Dcbui： XFanl lii JxLfj XXFlu! jo Fale E二張aulq=二 KL ELvdiiv了(2)亡-口皿h典E血M問氓牛曲I EVd I曲HdpGTdbj|f pr-Lntf FnEirr Fx FM泊t出l.FwCJ-H;JlEl；i*li； FtJ(i+iJ-*K(a Jti-ijn曾仙(litar k 10.4?Rm any hny lu cphL Irwiubi i-ii *ca.-ari i/ 町 i-i -ji ii-c/c，卩t j-ji i+n)； y-u+bl 1 (!- i-l/iE4l -H 1-1 )*(- 11+I|/C*i -a( i+

8、mn J-y*b| i*1*|；3r-sl 1-1 )/( 1*1 |-a t-1J*si-a| i)/Ui);prints.mdln& QiS ! H iSrilriEa*eInt #linVQld)M-b( i-3 (i-a(l-1J/C4l-2 |-a| i-1 )*M-a( i3ZC4i-2 laO H J: y-U+b| i-13 JC*i-11-aG L- b*H-a| L j/(Mi-11-at iy-y+blilCM-jll-EJJ/t+ri -j| i-E-J|.x； prlntFCnFn*Entar k IV W.l B.l#& N.J 19.4113. ”l5iEntar

9、也肝RM B.M1595 .3a2B.M412 .沾昭EilnrUEW0 ., yB.WHflHH=f ,W(3)H /:弭hlgwli 別“ f+ = 丿1J I， 血 d U0田 idH 苗em Jnavrl PjroircJE H.uiH I-odaHnpM Ci iLnc：lutii!- int iMlnCwldli-r| |ANmtmbHf -1| tMto*xlfli 8.9If.3 H.6LEhtr mH.Htn6.J9142 H.3413A A.3&B12 .J52M.DrLJ- xEBi 15LseemE*t m=9.31*Mn-HLIWWH rutr xr0i?Pnsau

10、fUi hvyV4.：INUJTy超L沾？ .v coni iriurprltftnter k r); fDr(l*lai nJ - f1, Bd3 IJ4 / Ftbd ik Filsi 1 1 FJhd bn Fl*j 2 Ri：iil.? SL Jrtujsl i第九題:，題目已知直升機旋轉(zhuǎn)機翼外形曲線輪廓線上的某些型值點及端點處的一階導(dǎo)數(shù)值y(xo)=1.86548,y(x)=-O.O46115,試計算該曲線上橫坐標(biāo)為2,4,6,12,16,30,60,110,180,280,400,5處點的縱坐標(biāo)（要求該 曲線具有二階光滑度）。二，基本方法根據(jù)課本圖4-9流程圖及追趕法解方程。三，計

11、算結(jié)果x246121630y7.83433110.47688012.35848516.57965319.09317325.386597x60110180280400515y32.80436936.64788635.91714129.36842716.7991410.524459四，結(jié)果分析利用追趕法求 Mi，比高斯消去法步驟更簡單， 雖然結(jié)果有 微小偏差，尤其在s515處與課本結(jié)果偏差較大。 但在精度 要求不高的情況下，追趕法可以減少工作量， 提高工作效率。 五，源程序#i nclude int main(v oid)int i,j;doublea19,b19,h19,s18,t18,g19,

12、x,y,p,q,f19,m20, n20;printf(輸入左端點一階導(dǎo):);scan f(%lf,&p);printf(輸入右端點一階導(dǎo):);sca nf(%lf,&q);prin tf(E nter x:);for(i=0;i19;i+)scan f(%lf,&ai);prin tf(E nter f(x):);for(i=0;i19;i+)sca nf(%lf,&bi);for(i=0;i19;i+) hi=ai-ai-1;for(i=1;i18;i+)si=hi/(hi+hi+1);for(i=1;i18;i+)ti=1-si;for(i=1;i18;i+)gi=6/(hi+hi+1)

13、*(bi+1-bi)/hi+1-(bi-bi-1)/hi); g0=6/h1*(b1-b0)/h1-p);g18=6/h18*(q-(b18-b17)/h 18);f1=0.5;for(i=2;i19;i+)fi=ti-1/(2-si-1*fi-1);n1=g0/2; for(i=2;i0;i-)mi=ni-fi*mi+1;printf(輸出 M:n);for(i=1;i20;i+)prin tf(%lf ,mi);for(j=0;jai;i+);y=mi*(ai-x)*(ai-x)*(ai-x)/hi/6; y=y+mi+1*(x-ai-1)*(x-ai-1)*(x-ai-1)/hi/6;

14、y=y+(bi-1-mi*hi*hi/6)*(ai-x)/hi; y=y+(bi-mi+1*hi*hi/6)*(x-ai-1)/hi;prin tf(f(x):%lf,y);return 0;六，附圖1-91scanFCtlfpKp尸iriE.*l.HiE*r Ktai.； far(L-fl!：l!lf：i+3!printFfTHter 何切；卄TrzwnliiriK” Kh:l Jrtur41B.：EC1 ;!+hl-aiij-a|l-l4f |l;K(b(1i-* D2f ria1*CR-b| 1|/h| H);f|1-lu53fBF(LE311；i+JFi-tL-1|/(“iF| i I

15、js；farU-lSKMSl*)n n-M ri-1|-= 1-1 *n( l-1)3/CS-s(l-1J*Fl-igi|i1 屯鳳 EZ 如Er蟲Hd斷!kw 出旨| |GlDb<|Tl 回1亡冃創(chuàng)|PM”申 CJ * I n *iLnclud- Iwt第十二題一，題目某類疾病發(fā)生率y %o和年齡段x（每五年為一段，例如0-5歲為第一段， 6-10歲為第二段）之間有形如y中咕乂的關(guān)系。試根據(jù)觀測得到的如下 數(shù)據(jù)表，用最小二乘法確定式中的參數(shù) a和b,并計算相應(yīng)的均方誤 差與最大偏差。x123456789y0.8982.383.071.842.021.942.222.774.02x1011

16、1213141516171819y4.765.466.5310.916.522.535.750.661.681.8二，基本方法（1） 將非線性方程化為線性方程，對y二aebx等號兩邊取對數(shù)得到In y = ln a bx。（2）按照最小二乘法的原始定義，不對方程進行變換，而均方誤差的方程式，通過求解非線性方程組來解決問題，此時要用到牛頓迭代法。三，計算結(jié)果（1）a=0.681361,b=0.230620,e=77.307455,m=18.328522（2）第一次：a=0.721560,b=0.249856,e=17.925350,m=8.1115274 第二次：a=0.239368,b=0.3

17、10747,e=5.430872,m=3.467987.四，結(jié)果分析由以上結(jié)果對比，第二種方法結(jié)果更加精確，但是采用迭代法計算時, 又必須多次選擇初值才有可能取得較好的結(jié)果。五，源程序（1）#i nclude#in cludeint main（v oid）int i;double a,b,e,m,a0,a1,b0,b1,y1,y2,x19,y19;for(i=0;i19;i+)xi=i+1;prin tf(x:);for(i=0;i19;i+)prin tf(%lf ,xi);prin tf(n y:);for(i=0;i19;i+)sea nf(%lf, &yi);for(i=0;i19;i

18、+)yi=log(yi);a0=a仁b0=b仁y1=y2=0;for(i=0;i19;i+)a仁 a1+xi;y仁 y1+yi;b0=b0+xi;b仁 b1+pow(xi,2);y2=y2+xi*yi;a0=19;b=(y2-y1*(b0/a0)/(b1-a1*(b0/a0);a=(y1-b*a1)/a0;a=exp(a);for(i=0;i19;i+)yi=exp(yi);e=m=0;for(i=0;i19;i+) e=e+pow(a*exp(i+1)*b)-yi),2);for(i=0;im) m=pow(a*exp(i+1)*b)-yi),2);else;prin tf(a=%lf b=

19、%lf e=%lf m=%lfn,a,b,e/19,sqrt(m); return 0;(2)#i nclude#in clude int main(v oid)int i,j;double a,b,e,m,f,f1,f2,g,g1,g2,c,d,x19,y19;for(i=0;i19;i+)xi=i+1;prin tf(x:);for(i=0;i19;i+)prin tf(%lf ,xi);prin tf(n y:);for(i=0;i19;i+)sca nf(%lf, &yi);prin tf(E nter a0:);sca nf(%lf,&a);prin tf(E nter b0:);s

20、ca nf(%lf,&b);f=f仁 f2=g=g 仁 g2=0;for(j=0;j20;j+)for(i=0;i19;i+)f=f+2*(a*exp(i+1)*b)-yi)*exp(i+1)*b);f1=f1+2*exp(2*(i+1)*b);f2=f2+4*(i+1)*a*exp(2*(i+1)*b)-2*(i+1)*yi*exp(i+1)*b);g=g+2*(a*exp(i+1)*b)-yi)*(i+1)*a*exp(i+1)*b);g仁 g1+4*(i+1)*a*exp(2*(i+1)*b)-2*(i+1)*yi*exp(i+1)*b);g2=g2+4*pow(i+1)*a,2)*ex

21、p(2*(i+1)*b)-2*pow(i+1),2)*a*yi*exp(i+1)*b);c=f1*a+f2*b-f;d=g1*a+g2*b-g;b=(d-c*(g1/f1)/(g2-f2*(g1/f1);a=(c-f2*b)/f1;e=0;m=0;for(i=0;i19;i+)e=e+pow(a*exp(i+1)*b)-yi),2);for(i=0;im)m=pow(a*exp(i+1)*b)-yi),2);prin tf(a=%lfb=%lf 均方 誤差 e=%lf 最大偏 差m=%lfn,a,b,e/19,sqrt(m); return 0;六，附圖六,附圖茴 H+zK - Mkraufr

22、 ViuuJ C+ - lii+zlxppl茴 H+zK - Mkraufr ViuuJ C+ - lii+zlxpplb-fl.2S2337b-fl月曲BE*hfl.317257hH_22JaT7b-e.2K3&4bih1H.UUU4V-6.4986-fl.156713-fl.aiasa?fl-142ni3I 方 i；亍iTl.BWME j律1E1*541794P2SJW ZFMH* F2U9-a3B.?34WF.A3.4|n-7H.%n-4L4bIE* dCihVJhivfra -RH71-14 活拒 eiiI-A115E71I-19.97E4L1i-39u73?ma 1嚴( | n | 6 3C斷B - JI1If?i-Sla/ITSBTtr6ia59&?3*SttnSttnLWt辦公室衛(wèi)生管理制度一、主要內(nèi)容與適用范圍1 本制度規(guī)定了辦公室衛(wèi)生管理的工作內(nèi)容和要求及檢查與考核。2 此管理制度適用于本公司所有辦公室衛(wèi)生的管理二、定義 1 公共區(qū)域：包括辦公室走道、會議室、衛(wèi)生間，每天由行政文員進行清掃；2個人區(qū)域：包括個人辦公桌及辦公區(qū)域由各部門工作人員每天自行清掃。1. 公共區(qū)域環(huán)境衛(wèi)生應(yīng)做到以下幾點：1) 保持公共區(qū)域及個人區(qū)域地面干凈清潔、無污物、污水、浮土，無死角。2) 保持門窗干凈、