计算方法上机题doc

1、第八题:题目，基本方法从函数表X0.00.10.1950.30.4010.5f（x）：0.398940.39695 :0.391420.38138:0.368120.35206出发，用下列方法计算f(0.15),f(0.31),f(0.47)的近似值：(1)分段线性插值；(2)分段二次插值;(3)全区间上拉格朗日插值(要求一次性输入整张函数表，并利用计算机选择在插值计算中所需的结点)(1) 先选取两个结点与X使xx，X然后在区间【Xi，Xi上作线性插值，即得 f(x) Li(x)= ypXl y。口XiXiXiXd(2) 选取距点X最近的三个结点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）三，计算结果（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四，结果分析编程基本按照课本套公式，方法从(1)至(3),计算结果越来越精确 五，源程序(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第九题:，题目已知直升机旋转机翼外形曲线轮廓线上的某些型值点及端点处的一阶导数值y(xo)=1.86548,y(x)=-O.O46115,试计算该曲线上横坐标为2,4,6,12,16,30,60,110,180,280,400,5处点的纵坐标（要求该 曲线具有二阶光滑度）。二，基本方法根据课本图4-9流程图及追赶法解方程。三，计

11、算结果x246121630y7.83433110.47688012.35848516.57965319.09317325.386597x60110180280400515y32.80436936.64788635.91714129.36842716.7991410.524459四，结果分析利用追赶法求 Mi，比高斯消去法步骤更简单， 虽然结果有 微小偏差，尤其在s515处与课本结果偏差较大。 但在精度 要求不高的情况下，追赶法可以减少工作量， 提高工作效率。 五，源程序#i nclude int main(v oid)int i,j;doublea19,b19,h19,s18,t18,g19,

12、x,y,p,q,f19,m20, n20;printf(输入左端点一阶导:);scan f(%lf,&p);printf(输入右端点一阶导:);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亡冃创|PM”申 CJ * I n *iLnclud- Iwt第十二题一，题目某类疾病发生率y %o和年龄段x（每五年为一段，例如0-5岁为第一段， 6-10岁为第二段）之间有形如y中咕乂的关系。试根据观测得到的如下 数据表，用最小二乘法确定式中的参数 a和b,并计算相应的均方误 差与最大偏差。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等号两边取对数得到In y = ln a bx。（2）按照最小二乘法的原始定义，不对方程进行变换，而均方误差的方程式，通过求解非线性方程组来解决问题，此时要用到牛顿迭代法。三，计算结果（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.四，结果分析由以上结果对比，第二种方法结果更加精确，但是采用迭代法计算时, 又必须多次选择初值才有可能取得较好的结果。五，源程序（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办公室卫生管理制度一、主要内容与适用范围1 本制度规定了办公室卫生管理的工作内容和要求及检查与考核。2 此管理制度适用于本公司所有办公室卫生的管理二、定义 1 公共区域：包括办公室走道、会议室、卫生间，每天由行政文员进行清扫；2个人区域：包括个人办公桌及办公区域由各部门工作人员每天自行清扫。1. 公共区域环境卫生应做到以下几点：1) 保持公共区域及个人区域地面干净清洁、无污物、污水、浮土，无死角。2) 保持门窗干净、