北航數(shù)值分析大作業(yè)(三)

1、北航數(shù)值分析大作業(yè)（三）（看到那個收10個財富值的就不舒服，共享的精髓就在于免費）一、題目關(guān)于x，y，t，u，v，w的下列方程組0.5cost+u+v+w-x=2.67t+0.5sinu+v+w-y=1.070.5t+u+cosv+w-x=3.74t+0.5u+v+sinw-y=0.79以及關(guān)于z，t，u的下列二維數(shù)表z ut00.40.81.21.620-0.5-0.340.140.942.063.50.2-0.42-0.5-0.260.31.182.380.4-0.18-0.5-0.5-0.180.461.420.60.22-0.34-0.58-0.5-0.10.620.80.78-0.0

2、2-0.5-0.66-0.5-0.021.01.50.46-0.26-0.66-0.74-0.5確定了一個二元函數(shù)z=f(x，y。1試用數(shù)值方法求出f(x，y在區(qū)域D=(x,y0x0.8,0.5y1.5上的一個近似表達式要求p(x,y一最小的k值達到以下的精度其中xi=0.08i，yj=0.5+0.05j。2計算f(xi*,yj*,p(xi*,yj*(i=1,2,8;j=1,2,5的值，以觀察p(x,y逼近f(x,y的效果，其中xi*=0.1i,yj*=0.5+0.2j。二、算法設(shè)計方案1.將和代入非線性方程組中，用牛頓法解出和；2.以采取分片二次插值，選擇（m，n）滿足如果或，則m=1或4；

3、如果或，則n=1或n=4。選擇為插值節(jié)點，相應(yīng)的Lagrange形式的插值多項式為其中 (k=m-1, m, m+1 (r=n-1, n, n+1并將和代入，便得到了數(shù)表。3.進行曲面擬和系數(shù)矩陣， 其中，k從0逐漸增大，直到，便得到了要求精度的系數(shù)。三、全部源程序/在visual Studio2010的c語言環(huán)境下編譯通過#include #include #include #define n 4static double xn,x_n,U1121,t1121,u1121,tt85,uu85,UU85,c99;int k;double MAX(double an/求數(shù)組中的最大值int i;

4、double max;max=fabs(a0;for(i=0;i if(fabs(ai>maxmax=fabs(ai; return(max;void DooLittle(double ann,double bn/選主元的DooLittle 分解法求線性方程組int i,j,k,t,i_k,mn;double unn,lnn,sn,yn;double u_x,l_u,max,temp;for(k=1;k<=n;k+ for(i=k;i<=n;i+ l_u=0;for(t=1;t<=k-1;t+ l_u=l_u+li-1t-1*ut-1k-1;si-1=ai-1k-1-l

5、_u;max=fabs(sk-1;i_k=k;mk-1=k;for(i=k;i<=n;i+if(fabs(si-1>maxmax=fabs(si-1;i_k=i;mk-1=i_k;if(i_k!=kfor(t=1;t<=k-1;t+ temp=lk-1t-1;lk-1t-1=li_k-1t-1;li_k-1t-1=temp;for(t=k;t<=n;t+ temp=ak-1t-1;ak-1t-1=ai_k-1t-1;ai_k-1t-1=temp;temp=sk-1;sk-1=si_k-1;si_k-1=temp;uk-1k-1=sk-1;for(j=k+1;j<=

6、n;j+ l_u=0;for(t=1;t<=k-1;t+ l_u=l_u+lk-1t-1*ut-1j-1;uk-1j-1=ak-1j-1-l_u;for(i=k+1;i<=n;i+ li-1k-1=si-1/uk-1k-1; for(k=1;k<=n-1;k+ temp=bk-1;bk-1=bmk-1-1;bmk-1-1=temp;y0=b0;for(i=2;i<=n;i+l_u=0;for(t=1;t<=i-1;t+ l_u=l_u+li-1t-1*yt-1; yi-1=bi-1-l_u;x_n-1=yn-1/un-1n-1;for(i=n-1;i>=1;

7、i-u_x=0;for(t=i+1;t<=n;t+u_x=u_x+ui-1t-1*x_t-1;x_i-1=(yi-1-u_x/ui-1i-1;void NEWTON(double x_x,double y_y/牛頓法求解非線性方程組子程序 double Fn,F_1nn;int i,j;for(i=0;i x_i=xi=1;for(i=0;i<2000;i+F0=(2.67-0.5*cos(x0-x1-x2-x3+x_x;F1=(1.07-x0-0.5*sin(x1-x2-x3+y_y;F2=(3.74-0.5*x0-x1-cos(x2-x3+x_x;F3=(0.79-x0-0.5

8、*x1-x2-sin(x3+y_y;F_100=-0.5*sin(x0;F_101=1; F_102= 1;F_103= 1;F_110=1; F_111=0.5*cos(x1;F_112= 1; F_113= 1;F_120=0.5;F_121=1;F_122=-sin(x2; F_123= 1;F_130=1;F_131=0.5;F_132= 1;F_133=cos(x3;DooLittle(F_1,F;if(MAX(x_/MAX(x<=1.0e-12return; for(j=0;j xj+=x_j; cout<<"迭代2000次沒有帶到精度要求"&

9、lt; return;void Solve_tu(/求解對應(yīng)的t、u子程序int i,j;double x_x11,y_y21; for(i=0;i<11;i+x_xi=0.08*i; for(j=0;j<21;j+y_yj=0.5+0.05*j;for(i=0;i<11;i+for(j=0;j<21;j+ NEWTON(x_xi,y_yj;tij=x0;uij=x1; for(i=0;i<8;i+x_xi=0.1*(i+1; for(j=0;j<5;j+y_yj=0.5+0.2*(j+1;for(i=0;i<8;i+for(j=0;j<5;j+

10、 NEWTON(x_xi,y_yj;ttij=x0;uuij=x1;void Lagrange(double *t,double *u,double *U,int la,int lb,int flag/分片二次代數(shù)插值int i,j,k,m,p,c,d,q;double a1121,b1121,Z66,temp1,temp2,L13,L23;Z00=-0.5 ; Z01=-0.34; Z02= 0.14; Z03= 0.94; Z04= 2.06; Z05= 3.5;Z10=-0.42; Z11=-0.5 ; Z12=-0.26; Z13= 0.3 ; Z14= 1.18; Z15= 2.38

11、;Z20=-0.18; Z21=-0.5 ; Z22=-0.5 ; Z23=-0.18; Z24= 0.46; Z25= 1.42;Z30= 0.22; Z31=-0.34; Z32=-0.58; Z33=-0.5 ; Z34=-0.1 ; Z35= 0.62;Z40= 0.78; Z41=-0.02; Z42=-0.5 ; Z43=-0.66; Z44=-0.5 ; Z45=-0.02;Z50= 1.5 ; Z51= 0.46; Z52=-0.26; Z53=-0.66; Z54=-0.74; Z55=-0.5; for(i=0;i for(j=0;j c=int(*(u+i*lb+j/0.

12、4;aij=(c+d*0.4; c=int(*(t+i*lb+j/0.2;bij=(c+d*0.2;if(aij<0.4 aij=0.4;if(aij>1.6 aij=1.6;if(bij<0.2 bij=0.2;if(bij>0.8 bij=0.8;for(k=0;k<3;k+ temp1=1;temp2=1;for(m=0;m<3;m+if(m!=ktemp1*=(*(t+i*lb+j-(bij+(m-1*0.2/(bij+(k-1*0.2-(bij+(m-1*0.2;temp2*=(*(u+i*lb+j-(aij+(m-1*0.4/(aij+(k-1*

13、0.4-(aij+(m-1*0.4;L1k=temp1;L2k=temp2;temp1=0;for(k=0;k<3;k+for(m=0;m<3;m+p=int(bij/0.2-1+k;q=int(aij/0.4-1+m; temp1+=L1k*L2m*Zpq;*(U+i*lb+j=temp1;if(flagprintf("%s%2d%s%.2f%s%2d%s%.2f%s%19.11e%s","x(",i,"=",0.08*i,", y(",j,"=",0.5+0.05*j,"

14、;, ",*(U+i*lb+j," "if(!fmod(j+1,2cout< if(j=lb-1cout< void Multi(double *a, double *b, double *c, int la, int lb, int lc, int r, int s, int t/求r*s階矩陣A與s*t階矩陣B的乘積矩陣Cint i, j, k;for (i=0; i for (j=0; j *(c+i*lc+j=0;for (k=0; k double Inverse(double *a, double *b, int la, int lb, i

15、nt N /求n階方陣A的逆矩陣Bint i, j, k;double temp;for(i=0; i for(j=0; j if (i=j*(b+i*lb+j=1;else*(b+i*lb+j=0;for (k=0; k j=k;for (i=k+1; i if (fabs(*(a+i*la+k>fabs(*(a+j*la+k j=i;if (j!=kfor (i=0; i temp=*(a+j*la+i;*(a+j*la+i=*(a+k*la+i;*(a+k*la+i=temp;temp=*(b+j*lb+i;*(b+j*lb+i=*(b+k*lb+i;*(b+k*lb+i=temp

16、;if (*(a+k*la+k=0return 0;if (temp=*(a+k*la+k!=1for (i=0; i *(a+k*la+i/=temp;*(b+k*lb+i/=temp;for (i=0; i if (*(a+i*la+k!=0 && (i!=ktemp=*(a+i*la+k;for (j=0; j *(a+i*la+j-=temp*(*(a+k*la+j;*(b+i*lb+j-=temp*(*(b+k*lb+j;return 0;void solve_C(int i,j,r,s;double t12121, t22121, t32121,d99,e99;do

17、uble B119, B_T911, G219, G_T921,P1121;double temp, FangCha;for(i=0;i<9;i+for(j=0;j<11;j+Bji=pow(0.08*j,i;B_Tij=pow(0.08*j,i;for(j=0;j<21;j+Gji=pow(0.5+0.05*j,i;G_Tij=pow(0.5+0.05*j,i; for (k=0; k<9; k+ FangCha=0;Multi(B_T0, B0, t10, 11, 9, 21, k+1, 11, k+1; Inverse(t10, c0, 21, 9, k+1;Mu

18、lti(e0, c0, d0, 9, 9, 9, k+1, k+1, k+1;Multi(c0, B_T0, t10, 9, 11, 21, k+1, k+1, 11;Multi(t10, U0, t20, 21, 21, 21, k+1, 11, 21;Multi(G_T0, G0, t10, 21, 9, 21, k+1, 21, k+1;Inverse(t10, c0, 21, 9, k+1;Multi(G0, c0, t30, 9, 9, 21, 21, k+1, k+1;Multi(t20, t30, c0, 21, 21, 9, k+1, 21, k+1;for(i=0;i<

19、11;i+for(j=0;j<21;j+temp=0;for(r=0;r for(s=0;s temp+=crs*Bir*Gjs;Pij=temp; FangCha+=(Uij-temp*(Uij-temp;printf("%s%d%s%.11e%s","k=",k," Sigma=",FangCa,"n"if(FangCha<=1.0e-7printf("%s%d%s%.11e%s","達到精度要求時: k=",k," sigma=",Fan

20、gCha,"n系數(shù)c(r,s如下：n"for(i=0;i for(j=0;j printf("%s%d%s%d%s%19.11e%s","c(",i,",",j,"=",cij," "if(j=(k-1/2cout< cout< cout< return;cout<<"經(jīng)過8次擬合沒有達到所需精度；"< return;void Check(/驗證子程序 int i,j,r,s;double B89,G59,temp;La

21、grange( tt0,uu0,UU0,8,5,0; for(i=0;i for(j=0;j<8;j+Bji=pow(0.1*(j+1,i;for(j=0;j<5;j+Gji=pow(0.5+0.2*(j+1,i; for(i=0;i<8;i+for(j=0;j<5;j+temp=0;for(r=0;r for(s=0;s temp+=crs*Bir*Gjs;printf("%s%d%s%.2f%s%d%s%.2f","x*(" , i+1, "=" , 0.1*(i+1 , ", y*("

22、,j+1,"=",0.5+0.2*(j+1;printf("%s%19.11e%s%19.11e%s" , ", f(x*,y*=",UUij,", p(x*,y*=",temp,"n"void main( void NEWTON(double,double;void DooLittle(double *,double *;void Solve_tu(;void solve_C(;void Check(;void Lagrange(double *t,double *u,double *U,in

23、t la,int lb,int ;Solve_tu(;Lagrange(t0,u0,U0,11,21,1;solve_C(;Check(;getchar(;四、程序運行結(jié)果：1. 數(shù)表：xi, yj, f(xi, yjx( 0=0.00, y( 0=0.50, 4.46504018481e-001; x( 0=0.00, y( 1=0.55, 3.24683262928e-001;x( 0=0.00, y( 4=0.70, 3.40189556268e-003; x( 0=0.00, y( 5=0.75, -8.87358136380e-002;x( 0=0.00, y( 6=0.80, -1

24、.73371632750e-001; x( 0=0.00, y( 7=0.85, -2.50534611467e-001;x( 0=0.00, y( 8=0.90, -3.20276506388e-001; x( 0=0.00, y( 9=0.95, -3.82668066110e-001;x( 0=0.00, y(10=1.00, -4.37795766738e-001; x( 0=0.00, y(11=1.05, -4.85758941444e-001;x( 0=0.00, y(12=1.10, -5.26667254884e-001; x( 0=0.00, y(13=1.15, -5.6

25、0638479797e-001;x( 0=0.00, y(14=1.20, -5.87796538768e-001; x( 0=0.00, y(15=1.25, -6.08269779090e-001;x( 0=0.00, y(16=1.30, -6.22189452876e-001; x( 0=0.00, y(17=1.35, -6.29688378186e-001;x( 0=0.00, y(18=1.40, -6.30899760003e-001; x( 0=0.00, y(19=1.45, -6.25956152545e-001;x( 1=0.08, y( 0=0.50, 6.38015

26、226511e-001; x( 1=0.08, y( 1=0.55, 5.06611755147e-001;x( 1=0.08, y( 2=0.60, 3.82176369277e-001; x( 1=0.08, y( 3=0.65, 2.64863491154e-001;x( 1=0.08, y( 6=0.80, -4.34680402049e-002; x( 1=0.08, y( 7=0.85, -1.31601056789e-001;x( 1=0.08, y(12=1.10, -4.63904911519e-001; x( 1=0.08, y(13=1.15, -5.0923672470

27、1e-001;x( 1=0.08, y(14=1.20, -5.47761117962e-001; x( 1=0.08, y(15=1.25, -5.79594388339e-001;x( 1=0.08, y(16=1.30, -6.04857258890e-001; x( 1=0.08, y(17=1.35, -6.23673421332e-001;x( 1=0.08, y(18=1.40, -6.36168248413e-001; x( 1=0.08, y(19=1.45, -6.42467656690e-001;x( 1=0.08, y(20=1.50, -6.42697102700e-

28、001;x( 2=0.16, y( 0=0.50, 8.40081395767e-001; x( 2=0.16, y( 1=0.55, 6.99764165673e-001;x( 2=0.16, y( 2=0.60, 5.66061442352e-001; x( 2=0.16, y( 3=0.65, 4.39171608118e-001;x( 2=0.16, y( 6=0.80, 1.00638523891e-001; x( 2=0.16, y( 7=0.85, 2.06074006784e-003;x( 2=0.16, y( 8=0.90, -8.93540247670e-002; x( 2

29、=0.16, y( 9=0.95, -1.73626968865e-001;x( 2=0.16, y(10=1.00, -2.50799956160e-001; x( 2=0.16, y(11=1.05, -3.20932269445e-001;x( 2=0.16, y(12=1.10, -3.84097735005e-001; x( 2=0.16, y(13=1.15, -4.40382175418e-001;x( 2=0.16, y(14=1.20, -4.89881152313e-001; x( 2=0.16, y(15=1.25, -5.32697965534e-001;x( 2=0.

30、16, y(16=1.30, -5.68941879292e-001; x( 2=0.16, y(17=1.35, -5.98726549515e-001;x( 2=0.16, y(18=1.40, -6.22168629750e-001; x( 2=0.16, y(19=1.45, -6.39386535697e-001;x( 2=0.16, y(20=1.50, -6.50499350788e-001;x( 3=0.24, y( 0=0.50, 1.05151509180e+000; x( 3=0.24, y( 1=0.55, 9.02927430831e-001;x( 3=0.24, y

31、( 2=0.60, 7.60580266860e-001; x( 3=0.24, y( 3=0.65, 6.24715198146e-001;x( 3=0.24, y( 4=0.70, 4.95519756001e-001; x( 3=0.24, y( 5=0.75, 3.73134042775e-001;x( 3=0.24, y( 6=0.80, 2.57656748872e-001; x( 3=0.24, y( 7=0.85, 1.49150559410e-001;x( 3=0.24, y( 8=0.90, 4.76469867734e-002; x( 3=0.24, y( 9=0.95,

32、 -4.68493232015e-002;x( 3=0.24, y(12=1.10, -2.88573700635e-001; x( 3=0.24, y(13=1.15, -3.55406364786e-001;x( 3=0.24, y(20=1.50, -6.39746839258e-001;x( 4=0.32, y( 2=0.60, 9.64607727216e-001; x( 4=0.32, y( 3=0.65, 8.20347369475e-001;x( 4=0.32, y( 4=0.70, 6.82447678179e-001; x( 4=0.32, y( 5=0.75, 5.510

33、85208598e-001;x( 4=0.32, y( 6=0.80, 4.26392385902e-001; x( 4=0.32, y( 7=0.85, 3.08462995633e-001;x( 4=0.32, y( 8=0.90, 1.97357129692e-001; x( 4=0.32, y( 9=0.95, 9.31056208594e-002;x( 4=0.32, y(10=1.00, -4.28599223403e-003; x( 4=0.32, y(11=1.05, -9.48339252969e-002;x( 4=0.32, y(12=1.10, -1.7857299036

34、4e-001; x( 4=0.32, y(13=1.15, -2.55553779055e-001;x( 4=0.32, y(14=1.20, -3.25840150158e-001; x( 4=0.32, y(15=1.25, -3.89506988363e-001;x( 4=0.32, y(16=1.30, -4.46638204599e-001; x( 4=0.32, y(17=1.35, -4.97324951768e-001;x( 4=0.32, y(18=1.40, -5.41664032699e-001; x( 4=0.32, y(19=1.45, -5.79756479795e

35、-001;x( 5=0.40, y( 0=0.50, 1.49832105248e+000; x( 5=0.40, y( 1=0.55, 1.33499863207e+000;x( 5=0.40, y( 4=0.70, 8.78960023174e-001; x( 5=0.40, y( 5=0.75, 7.39145108704e-001;x( 5=0.40, y( 6=0.80, 6.05744871487e-001; x( 5=0.40, y( 7=0.85, 4.78883861067e-001;x( 5=0.40, y( 8=0.90, 3.58650625882e-001; x( 5

36、=0.40, y( 9=0.95, 2.45102236196e-001;x( 5=0.40, y(10=1.00, 1.38268350928e-001; x( 5=0.40, y(11=1.05, 3.8154865070e-002;x( 5=0.40, y(12=1.10, -5.52528211681e-002; x( 5=0.40, y(13=1.15, -1.41986880814e-001;x( 5=0.40, y(14=1.20, -2.22094439096e-001; x( 5=0.40, y(15=1.25, -2.95635232460e-001;x( 5=0.40,

37、y(16=1.30, -3.62679511503e-001; x( 5=0.40, y(17=1.35, -4.23306164224e-001;x( 5=0.40, y(18=1.40, -4.77601036132e-001; x( 5=0.40, y(19=1.45, -5.25655426667e-001;x( 5=0.40, y(20=1.50, -5.67564743655e-001;x( 6=0.48, y( 0=0.50, 1.73189274038e+000; x( 6=0.48, y( 1=0.55, 1.56203457721e+000;x( 6=0.48, y( 2=

38、0.60, 1.39721691821e+000; x( 6=0.48, y( 3=0.65, 1.23780100674e+000;x( 6=0.48, y( 4=0.70, 1.08408753268e+000; x( 6=0.48, y( 5=0.75, 9.36322772315e-001;x( 6=0.48, y( 6=0.80, 7.94704449054e-001; x( 6=0.48, y( 7=0.85, 6.59387198028e-001;x( 6=0.48, y( 8=0.90, 5.30487586840e-001; x( 6=0.48, y( 9=0.95, 4.0

39、8088685454e-001;x( 6=0.48, y(10=1.00, 2.92244201230e-001; x( 6=0.48, y(11=1.05, 1.82982206854e-001;x( 6=0.48, y(12=1.10, 8.03084940354e-002; x( 6=0.48, y(13=1.15, -1.57904130516e-002;x( 6=0.48, y(14=1.20, -1.05344551621e-001; x( 6=0.48, y(15=1.25, -1.88398090610e-001;x( 6=0.48, y(16=1.30, -2.6500714

40、9319e-001; x( 6=0.48, y(17=1.35, -3.35237838904e-001;x( 6=0.48, y(18=1.40, -3.99164503887e-001; x( 6=0.48, y(19=1.45, -4.56868143302e-001;x( 6=0.48, y(20=1.50, -5.08434993278e-001;x( 7=0.56, y( 0=0.50, 1.97122178640e+000; x( 7=0.56, y( 1=0.55, 1.79532959950e+000;x( 7=0.56, y( 2=0.60, 1.62406711323e+

41、000; x( 7=0.56, y( 3=0.65, 1.45783058271e+000;x( 7=0.56, y( 6=0.80, 9.92349533324e-001; x( 7=0.56, y( 7=0.85, 8.49032663329e-001;x( 7=0.56, y(10=1.00, 4.56658513233e-001; x( 7=0.56, y(11=1.05, 3.3865449139e-001;x( 7=0.56, y(12=1.10, 2.27108255770e-001; x( 7=0.56, y(13=1.15, 1.22025089193e-001;x( 7=0

42、.56, y(14=1.20, 2.33922196376e-002; x( 7=0.56, y(15=1.25, -6.88187019710e-002;x( 7=0.56, y(16=1.30, -1.54649344213e-001; x( 7=0.56, y(17=1.35, -2.34152666459e-001;x( 7=0.56, y(18=1.40, -3.07391091913e-001; x( 7=0.56, y(19=1.45, -3.74434862348e-001;x( 7=0.56, y(20=1.50, -4.35360556536e-001;x( 8=0.64,

43、 y( 0=0.50, 2.21566786369e+000; x( 8=0.64, y( 1=0.55, 2.03420113361e+000;x( 8=0.64, y( 2=0.60, 1.85695514362e+000; x( 8=0.64, y( 3=0.65, 1.68435816416e+000;x( 8=0.64, y( 4=0.70, 1.51677635240e+000; x( 8=0.64, y( 5=0.75, 1.35451904115e+000;x( 8=0.64, y( 8=0.90, 9.02046083802e-001; x( 8=0.64, y( 9=0.9

44、5, 7.63226477663e-001;x( 8=0.64, y(10=1.00, 6.30604821954e-001; x( 8=0.64, y(11=1.05, 5.04252814597e-001;x( 8=0.64, y(12=1.10, 3.84216715546e-001; x( 8=0.64, y(13=1.15, 2.70520476641e-001;x( 8=0.64, y(14=1.20, 1.63168572400e-001; x( 8=0.64, y(15=1.25, 6.21485581168e-002;x( 8=0.64, y(16=1.30, -3.2566

45、6193968e-002; x( 8=0.64, y(17=1.35, -1.21016534844e-001;x( 8=0.64, y(18=1.40, -2.03251399623e-001; x( 8=0.64, y(19=1.45, -2.79330359558e-001;x( 8=0.64, y(20=1.50, -3.49319957540e-001;x( 9=0.72, y( 0=0.50, 2.46468422266e+000; x( 9=0.72, y( 1=0.55, 2.27805897940e+000;x( 9=0.72, y( 2=0.60, 2.0952512508

46、4e+000; x( 9=0.72, y( 3=0.65, 1.91671812800e+000;x( 9=0.72, y( 4=0.70, 1.74285462878e+000; x( 9=0.72, y( 5=0.75, 1.57399842733e+000;x( 9=0.72, y( 6=0.80, 1.41043483523e+000; x( 9=0.72, y( 7=0.85, 1.25240175061e+000;x( 9=0.72, y(12=1.10, 5.50774848504e-001; x( 9=0.72, y(13=1.15, 4.28825676973e-001;x(

47、 9=0.72, y(16=1.30, 1.00345478241e-001; x( 9=0.72, y(17=1.35, 3.26856518657e-003;x( 9=0.72, y(18=1.40, -8.76530659133e-002; x( 9=0.72, y(19=1.45, -1.72467247819e-001;x( 9=0.72, y(20=1.50, -2.51230220752e-001;x(10=0.80, y( 0=0.50, 2.71781110947e+000; x(10=0.80, y( 1=0.55, 2.52639950126e+000;x(10=0.80

48、, y( 4=0.70, 1.97457455665e+000; x(10=0.80, y( 5=0.75, 1.79951057910e+000;x(10=0.80, y( 6=0.80, 1.62944822055e+000; x(10=0.80, y( 7=0.85, 1.46465004375e+000;x(10=0.80, y(10=1.00, 1.00383741991e+000; x(10=0.80, y(11=1.05, 8.6191233228e-001;x(10=0.80, y(12=1.10, 7.25992371111e-001; x(10=0.80, y(13=1.1

49、5, 5.96137711520e-001;x(10=0.80, y(14=1.20, 4.72386627914e-001; x(10=0.80, y(15=1.25, 3.54758095898e-001;x(10=0.80, y(16=1.30, 2.43254184181e-001; x(10=0.80, y(17=1.35, 1.37862222525e-001;x(10=0.80, y(18=1.40, 3.85567703264e-002; x(10=0.80, y(19=1.45, -5.46985959345e-002;x(10=0.80, y(20=1.50, -1.419

50、49659709e-001;2. 選擇過程的k,sigma的值：k=0; Sigma=1.44288077184e+002;k=1; Sigma=3.22090897364e+000;k=3; Sigma=1.72117537914e-004;KHD:工藝基本知識,th=578用一只標準壓力表檢定甲乙里兩只壓力k,sigma 以及 p(x,y 中的系數(shù) c_rs(r = 0, 1,.,k; s = 0, 1,.,k k = 5, sigma = 2.541378464155e-008c(0,0= 2.02123051865e+000; c(0,1=-3.66842679633e+000; c(

51、0,2= 7.09248619136e-001;c(1,3= 1.63118340106e+000; c(1,4=-4.84719972982e-001; c(1,5= 6.06142812008e-002;c(2,0= 2.56889811789e-001; c(2,1= 1.57991867296e+000; c(2,2=-4.63408147906e-001;c(3,0=-2.69260334060e-001; c(3,1=-7.30247662993e-001; c(3,2= 1.07614508118e+000;c(3,3=-8.07012812750e-001; c(3,4= 3.

52、02872801922e-001; c(3,5=-4.59726349324e-002;c(4,3= 2.43330467049e-001; c(4,4=-1.41334736110e-001; c(4,5= 2.65102405350e-002;c(5,0=-5.59032635880e-002; c(5,1= 1.43199240531e-001; c(5,2=-1.36270369212e-001;c(5,3= 4.07196214869e-002; c(5,4= 3.77503223942e-003; c(5,5=-2.66770114760e-003;x*(1=0.10, y*(1=

53、0.70, f(x*,y*= 1.94720407918e-001, p(x*,y*= 1.94730358440e-001;表時，讀得標準表的指示值為，甲乙兩表的讀數(shù)各為.x*(1=0.10, y*(4=1.30, f(x*,y*=-5.97566707641e-001, p(x*,y*=-5.97558849575e-001;x*(1=0.10, y*(5=1.50, f(x*,y*=-6.46459593901e-001, p(x*,y*=-6.46446099003e-001;x*(2=0.20, y*(1=0.70, f(x*,y*= 4.05979189288e-001, p(x*

54、,y*= 4.05989540586e-001;x*(2=0.20, y*(2=0.90, f(x*,y*=-2.25159583746e-002, p(x*,y*=-2.25211142281e-002;x*(2=0.20, y*(3=1.10, f(x*,y*=-3.38220816040e-001, p(x*,y*=-3.38224018827e-001;x*(2=0.20, y*(4=1.30, f(x*,y*=-5.44437831522e-001, p(x*,y*=-5.44430444410e-001;x*(2=0.20, y*(5=1.50, f(x*,y*=-6.47361338568e-001, p(x*,y*=-6.47348000423e-001;x*(3=0.30, y*(1=0.70, f(x*,y*= 6.34777195151e-001, p(x*,y*= 6.34787453517e-001;x*(3=0.30, y*(2=0.90, f(x*,y*= 1.58801168839e-001, p(x