分布函數(shù)分位點及pvalue計算及程序?qū)崿F(xiàn)

1、 分布函數(shù)分位點及p-value計算及程序?qū)崿F(xiàn)各分布函數(shù)分位點及p-value計算及程序?qū)崿F(xiàn)一、各分布函數(shù)計算1.標(biāo)準(zhǔn)正態(tài)分布函數(shù)的計算方法（1）：利用泰勒展示求解 由于分布函數(shù)f(t)=t-x22dx不能被求解成為簡單的初等函數(shù)，所以不能直接用定義初等復(fù)合函數(shù)的方法計算其值。通過把正態(tài)密度函數(shù)f(x)=展開為泰勒級數(shù),對泰勒級數(shù)的各項求定積分，然后相加各項數(shù)值即得分布函數(shù)f（t）的值。f(x)=t-x22-x22x2x4x6x2nn （1） 1-+-+.+(-1)+.）123n1!22!23!2n!2-x22由（1），（2）得出t3t5t7t2n+1n f(t)=0.5t-+-+.+(-1)

2、+.) （3）123n1!232!253!27n!2(2n+1)因為（3）式是一個無窮級數(shù)，計算時只能對前面有限項相加。選擇越多，結(jié)果越準(zhǔn)確。但數(shù)值精度并不要求無限高，可選取50項就足夠，因而這里n=0,1,250。流程圖如下： 0f(t)=dx=0.5+et-x22dx （2）代開c+軟件，輸入以下程序代碼：#include<stdio.h>#include<math.h>#define pi 3.141592653main( )float z,t,y=0,jsx,f,b=1 ; /jsx為每個級數(shù)項的值，f為所要求分布函數(shù)的值int n=1,a=1,m=1 ;pri

3、ntf("please enter t:");scanf("%f",&t); /輸入分位點twhile(n<101) /n<101,得出的級數(shù)項值共有50個z=pow(t,n);jsx=a*z/(b*n);y+=jsx;a=-1*a ; /正負符號標(biāo)志b=b*2*m ;m+;n=n+2;f=0.5+y/sqrt(2*pi);printf("the result you want is:n");printf("%f",f);輸入數(shù)據(jù)：比如t=0.02,x=-0.01,x=0輸出: 方法（2）：積分

4、的近似算法： 正態(tài)分布函數(shù)f(t)=t-x22dx的計算歸為積分計算。b-a,k=0,1,.，n在每個子區(qū)間xk,xk+1（k=0,1,.,n-1）上采用梯形公式bb-af(x)dxf(a)+f(b)則可得復(fù)合梯形公式a2n-1hn-1htn=f(xk)+f(xk+1)=f(a)+2f(xk)+f(b)。如果將求積區(qū)間再二分2k=02k=1一次，則可提高求積精度。此時分點增至2n+1個，則由原來每個子區(qū)間xk,xk+11經(jīng)過二分只增加了一個分點x1=(xk+xk+1)，再用復(fù)合梯形公式求得該子區(qū)k+22設(shè)將區(qū)間a,b分為n等份，共有n+1個分點。分點xk=a+kh,h=b-ah間上的積分值為f

5、(xk)+2f(x1)+f(xk+1)（其中h=代表二分前的步k+n42hn-1hn-1長）。將每個子區(qū)間上的積分值相加得：t2n=f(xk)+f(xk+1)+f(x1)k+4k=02k=02c+程序算法：#include<stdio.h>#include<math.h>double f(double x)return exp(-x*x/2);double f(double a,double b,double ep=1e-6)double h,s1=0,s2=(b-a)*(f(a)+f(b)/2;int n,k;for(int n=1;fabs(s1-s2)>ep

6、;n*=2)h=(b-a)/n;s1 = s2;s2 = 0;for(int k=0;k<n;+k)s2 += h*f(a+(k+0.5)*h);s2 = (s1+s2)/2;return s2*sqrt(1/(8*atan(1.0);int main()double c=0,x;printf("please enter x:");scanf("%lf",&x);printf("%lf",0.5+f(c,x);return 0; 結(jié)果顯示：（輸入x的值為2）輸出： 方法（3）：連分式展開（算法程序見下面最后部分的總程序編

7、寫）（x）為標(biāo)準(zhǔn)正態(tài)分布函數(shù)，其是對稱的，只要求出x>0時（x）的值也就可以求得x<0時（x）=1-（-x）。由于（x）的近似式：1f(x)x-x22x2(-1)nnx22+1+3+5+.+(2n+1),當(dāng)0x3時（x）f(x)12n1-+.+,當(dāng)x3時xxxx其中f(x)=-x22采用遞推式表達，an+1=01f(x)xk2當(dāng)0x3時，令則（x）= +(-1)kxa=,k=n,n-1,.121+a1k(2k+1)+ak+1an+1=0f(x)當(dāng)x>3時，令則（x）= 1- ka=,k=n,n-1,.1x+a1kx+ak+1一般當(dāng)n>28時，精度可達10-122.bet

8、a分布函數(shù)的計算beta分布的分布函數(shù)ix(a,b) 遞推如下：1i(a+1,b)=i(a,b)-ux(a,b)xxai(a,b+1)=i(a,b)+1u(a,b)xxx1bu(a,b)=xa(1-x)b 其中xb(a,b)u(a+1,b)=a+bxu(a,b)xxaa+bux(a,b+1)=(1-x)ux(a,b)a由于利用beta分布的分布函數(shù)計算t分布，f分布，二項分布時，參數(shù)a,b的值要么是正整數(shù)，要么就是1/2的倍數(shù)。所以考慮參數(shù)a,b的值是正整數(shù)或者是1/2倍數(shù)情況下ix(a,b) 的計算。此時遞推公式初值選取有以下4種： 1111112當(dāng)a=,b=時，ux(,)=x(,)=1-

9、222222p111當(dāng)a=,b=1時，ux(,1)=-x),ix( 2221111當(dāng)a=1,b=時，ux(1,)=ix(1, 2222當(dāng)a=1,b=1時，ux(1,1)=x(1-x),ix(1,1)=x 3.t分布函數(shù)的計算ix(a,b) 與 t分布的關(guān)系n111-i(,)n2x22,t0t(t|n)=,其中x=。由于n為整數(shù)，所以n/2可以遞歸21n1,t0n+ti(,)x222到，即可以使用beta分布的遞推公式。 4.f分布函數(shù)計算ix(a,b) 與 f分布的關(guān)系mnmxf(x|m,n)=iy(,),其中y=。由于m,n均為整數(shù)，所以m/2,n/2均可以22n+mx遞歸到，即可以使用be

10、ta分布的遞推公式。 5二項分布函數(shù)的計算0,x0xkn-k由二項分布函數(shù)的定義：b（x|n,p)=(n,0xn k)p(1-p)k=01,xn再利用beta分布函數(shù)ix(a,b) 計算：即b（x|n,p)=i1-p(n-x,x+1),0xn6卡方分布函數(shù)的計算卡方分布函數(shù)h（x|n）的遞推計算公式如下：h(x|n)=h(x|n-2)-2f(x|n)nx1x2-1-2()e,(n=3,4,.)其中f(x|n)=xn2f(x|n)=f(x|n-2)2g()n-22xx-22h(x|2)=1-e,f(x|1)=則遞推初值為 x1-2h(x|1)=2f-1,f(x|2)=e27poisson分布函數(shù)

11、的計算,x0 p（x |l）=1-h（2l|2（x+1）二、各分布函數(shù)分位數(shù)的計算1.標(biāo)準(zhǔn)正態(tài)分布分位數(shù)的計算采用精度達到1.2x10-8的toda近似公式：ua(ybiy)0101i2,其中y=-ln4-(a1 )abi的各數(shù)值用數(shù)組b11表示如下：b11=(1.570796288,0.03706987906,-0.8364353589e-3,-0.2250947176e-3, 0.6841218299e-5,0.5824238515e-5,-0.1045274970e-5,0.8360937017e-7, -0.3231081277e-8,0.3657763036e-10,0.693623

12、3982e-12) 2. t分布函數(shù)的計算由t分布函數(shù)與ix(a,b) 的關(guān)系可得：11n1n由當(dāng)0a0,且ta滿足（tta| n）=1-ix(,)=1-a,其中x=2222n+ta2n1n1n,ta(n)=此得ix(,)，則b1-2a(,)=22222n+tp1當(dāng)a時,ta0，ta(n)=-t1-a(n) 23. f分布函數(shù)分位數(shù)的計算f分布的分位數(shù)fa(m,n)滿足f(fa|m,n)=1-a。由f分布函數(shù)與beta分布函數(shù)的關(guān)系得fa滿足f(fa|nm=,iymfamn(=-,a其中)1=22n+mfa，y則mfamn=ba(，) n+mfa22mnnba(,)所以fa(m,n)= mnm

13、(1-ba(,)22 4. 卡方分布函數(shù)分位數(shù)ca2(n)的計算當(dāng)n=1時，由h（x|1）=2-1=1-,得=1-/2，=ua，ca2(1)=(ua)222x-11當(dāng)n=2時，c(2)分布就是l=的指數(shù)分布e(),由h(x|2)=1-e2=1-a得22ca2(2)=-2lna 2當(dāng)n3時，ca2(n)n1-223其中ua是標(biāo)準(zhǔn)正態(tài)分布的分位數(shù)。 +ua，9n程序算法如下：#include<stdio.h>#include<stdlib.h>#include<math.h>#include<time.h>#define pi 3.141592653

14、#define maxtime 500 /定義最大迭代步數(shù)#define eps 1.0e-10 /定義精度double kafangfx(float x,int freedom);double poissonfx(float x,double p);double kafangfx1(float x,int freedom);double gaossfx1(float x);double kafangua0(float af,int freedom);double kafangua(float af,int freedom);double gama(int n);double kafangpx

15、(float x,int freedom);double gaossfx(float x);double gaossua(float af);double gaossua1(float af);double gaosspx(float x);double betafx(float x,double a,double b);double tdistfx(float x,int freedom);double fdistfx(float x,int freedom_m,int freedom_n); double betaua(float af,double a,double b);double

16、tdistua(float af,int freedom);double fdistua(float af,int freedom_m,int freedom_n); double binominalfx(float x,float p,int n); double betafx(float x,double a,double b)/貝塔分布函數(shù) int m,n;double i,u;double ta,tb;m=(int)(2*a);n=(int)(2*b);if(m%2=1&&n%2=1)ta=0.5;tb=0.5;u=sqrt(x*(1.0-x)/pi;i=1.0-2.0

17、/pi*atan(sqrt(1.0-x)/x);else if(m%2=1&&n%2=0)ta=0.5;tb=1;u=0.5*sqrt(x)*(1.0-x);i=sqrt(x);else if(m%2=0&&n%2=1)ta=1;tb=0.5;u=0.5*x*sqrt(1.0-x);i=1.0-sqrt(1.0-x);else if(m%2=0&&n%2=0)ta=1;tb=1;u=x*(1.0-x);i=x;while(ta<a)i=i-u/ta;u=(ta+tb)/ta*x*u;ta+=1;while(tb<b)i=i+u/tb;

18、u=(ta+tb)/tb*(1.0-x)*u;tb+=1;return i; double tdistfx(float x,int freedom_t) /t分布t(n)在x處的分布函數(shù)計算double t,prob;printf("輸入t分布的自由度freedom_t:n");scanf("%d",&freedom_t);printf("輸入t分布的x:n");scanf("%f",&x);t=freedom_t/(freedom_t+x*x);if(x>0)prob=1.0-0.5*bet

19、afx(t,freedom_t/2.0,0.5);elseprob=0.5*betafx(t,freedom_t/2.0,0.5);printf("t分布t（%f,%d）=%fn",x,freedom_t,prob); double fdistfx(float x,int freedom_m,int freedom_n) /f分布f（m,n）在x處分布函數(shù)計算double y,prob;printf("f分布的自由度freedom_m,freedom_n:n");scanf("%d%d",&freedom_m,&fre

20、edom_n);printf("輸入p(x<=x)分布值p中的x:n");scanf("%f",&x);if(x<=0) return 0.0;y=freedom_m*x/(freedom_n+freedom_m*x);prob=betafx(y,freedom_m/2.0,freedom_n/2.0);printf("f分布f（%f,%d,%d)=%fn",x,freedom_m,freedom_n,prob); double binominalfx(float x,float p,int n)/二項分布的分布函

21、數(shù)計算 double prob;printf("輸入二項分布參數(shù)p:n");scanf("%f",&p);printf("輸入二項分布x:n");scanf("%f",&x);printf("輸入二項分布n:n");scanf("%d",&n);if(x<0) prob=0.0;else if(x>=n) prob=1.0;else prob=betafx(1.0-p,n-int(x),int(x)+1);printf("二項分布

22、b(%f，%f，%d)=%fn",x,p,n,prob); double betaua(float af,double a,double b)/貝塔分布的上側(cè)分位數(shù)計算 int times=0;double x1,x2,xn;double f1,f2,fn,ua;x1=0.0;x2=1.0;f1=betafx(x1,a,b)-(1.0-af);f2=betafx(x2,a,b)-(1.0-af);while(fabs(x2-x1)/2.0)>eps)xn=(x1+x2)/2.0;fn=betafx(xn,a,b)-(1.0-af);if(f1*fn<0) x2=xn;el

23、se if(fn*f2<0) x1=xn;f1=betafx(x1,a,b)-(1.0-af);f2=betafx(x2,a,b)-(1.0-af);times+;if(times>maxtime) break;printf("times=%5dn",times);ua=xn;return ua; double tdistua(float af,int freedom)/t分布的上側(cè)分位數(shù)計算 double ua,tbp,bf;printf("t分布的上側(cè)分位概率af:n");scanf("%f",&af);pri

24、ntf("t分布的自由度freedomn");scanf("%d",&freedom);bf=1-af;printf("f分布的下側(cè)分位數(shù)bf=%fn",bf);if(af<=0.5)tbp=betaua(1-2*af,freedom/2.0,0.5);ua=sqrt(freedom/tbp-freedom);else if(af>0.5)tbp=betaua(1-2*(1.0-af),freedom/2.0,0.5);ua=-sqrt(freedom/tbp-freedom);printf("t分布下側(cè)

25、分位數(shù)t(%f,%d)=%fn",bf,freedom,ua); double fdistua(float af,int freedom_m,int freedom_n)/f分布的上側(cè)分位數(shù)計算float ua,tbp,bf;printf("f分布的上側(cè)概率分位數(shù)af:n");scanf("%f",&af);printf("f分布的自由度freedom_m,freedom_n:n");scanf("%d%d",&freedom_m,&freedom_n);bf=1-af;print

26、f("f分布的下側(cè)分位數(shù)bf=%fn",bf);tbp=betaua(af,freedom_m/2.0,freedom_n/2.0);ua=freedom_n*tbp/(freedom_m*(1.0-tbp);printf("f分布f(%f,%d,%d)=%fn",bf,freedom_m,freedom_n,ua); double kafangfx1(float x,int freedom)int k,n;double f,h,prob;k=freedom%2;if(k=1)f=exp(-x/2.0)/sqrt(2*pi*x);h=2.0*gaossf

27、x1(sqrt(x)-1.0;n=1;while(n<freedom)n=n+2;f=x/(n-2.0)*f;h=h-2.0*f;elsef=exp(-x/2.0)/2.0;h=1.0-exp(-x/2.0);n=2;while(n<freedom)n=n+2;f=x/(n-2.0)*f;h=h-2.0*f;prob=h;return prob; double kafangfx(float x,int freedom) /卡方分布函數(shù)計算 int k,n;double f,h,prob;printf("卡方分布的自由度freedom,n");scanf(&quo

28、t;%d",&freedom);printf("輸入p(x<=x)分布值p中的x:n");scanf("%f",&x);k=freedom%2;if(k=1)f=exp(-x/2.0)/sqrt(2*pi*x);h=2.0*gaossfx1(sqrt(x)-1.0;n=1;while(n<freedom)n=n+2;f=x/(n-2.0)*f;h=h-2.0*f;elsef=exp(-x/2.0)/2.0;h=1.0-exp(-x/2.0);n=2;while(n<freedom)n=n+2;f=x/(n-2.

29、0)*f;h=h-2.0*f;prob=h;printf("卡方分布k(%f,%d)=%fn",x,freedom,prob); double poissonfx(float x,float p) /poisson分布函數(shù)的計算 double prob;printf("輸入poisson分布參數(shù)p:n");scanf("%f",&p);printf("輸入poisson分布x:n");scanf("%f",&x);prob=1.0-kafangfx1(2*p,2*(int(x)+

30、1);printf("poisson分布p（%f,%f）=%fn",x,p,prob); double kafangua0(float af,int freedom)double ua,p,temp;if(freedom=1)p=1.0-(1.0-af+1.0)/2.0;temp=gaossua(p);ua=temp*temp;else if(freedom=2)ua=-2.0*log(af);elsetemp=1.0-2.0/(9.0*freedom)+sqrt(2.0/(9.0*freedom)*gaossua(af); ua=freedom*(temp*temp*te

31、mp);return ua; double kafangua(float af,int freedom)/卡方分布的上側(cè)分位數(shù)的計算 int times;double ua,x0,xn,bf;printf("輸入卡方分布的自由度freedom:n");scanf("%d",&freedom);printf("輸入卡方分布af:n");scanf("%f",&af);bf=1-af;printf("f分布的下側(cè)分位數(shù)bf=%fn",bf);x0=kafangua0(af,freed

32、om);printf("x0=%12.8fn",x0);if(freedom=1|freedom=2)ua=x0;else/采用牛頓迭代法times=1;xn=x0-(kafangfx1(x0,freedom)-1+af)/kafangpx(x0,freedom); while(fabs(xn-x0)>eps)x0=xn;xn=x0-(kafangfx1(x0,freedom)-1+af)/kafangpx(x0,freedom); times+;if(times>maxtime) break;printf("times=%5dn",time

33、s);ua=xn;printf("卡方分布下側(cè)分位數(shù)k(%f，%d)=%fn",bf,freedom,ua); double kafangpx(float x,int freedom)/卡方分布的密度函數(shù) double p,g;if(x<=0) return 0.0;g=gama(freedom);p=1.0/pow(2.0,freedom/2.0)/g*exp(-x/2.0)*pow(x,freedom/2.0-1.0); return p; double gama(int n)/伽馬分布函數(shù)gama(n/2)double g;int i,k;k=n/2;if(n%

34、2=1)g=sqrt(pi)*0.5;for(i=1;i<k;i+)g*=(i+0.5);elseg=1.0;for(i=1;i<k;i+)g*=i;return g; double gaossfx1(float x) double prob,t,temp;int i,n,symbol;temp=x;if(x<0)x=-x;n=28;if(x>=0&&x<=3.0)t=0.0;for(i=n;i>=1;i-)if(i%2=1) symbol=-1; else symbol=1;t=symbol*i*x*x/(2.0*i+1.0+t); pro

35、b=0.5+gaosspx(x)*x/(1.0+t); else if(x>3.0)t=0.0;for(i=n;i>=1;i-)t=1.0*i/(x+t);prob=1-gaosspx(x)/(x+t); x=temp;if(x<0)prob=1.0-prob;return prob; double gaossfx(float x)/正態(tài)分布函數(shù)的計算 double prob,t,temp;int i,n,symbol;printf("輸入p(x<=x)分布值p中的x:n"); scanf("%f",&x);temp=x;

36、if(x<0)x=-x;n=28;/連分式展開的階數(shù)if(x>=0&&x<=3.0)/連分式展開方法t=0.0;for(i=n;i>=1;i-)if(i%2=1) symbol=-1;else symbol=1;t=symbol*i*x*x/(2.0*i+1.0+t);prob=0.5+gaosspx(x)*x/(1.0+t);else if(x>3.0)t=0.0;for(i=n;i>=1;i-)t=1.0*i/(x+t);prob=1-gaosspx(x)/(x+t);x=temp;if(x<0)prob=1.0-prob;prin

37、tf("正態(tài)分布f(%f)=%fn",x,prob); double gaossua(float af)double b11=0.1570796288*10,0.3706987906*1.0e-1,-0.8364353589*1.0e-3,-0.2250947176*1.0e-3,0.6841218299*1.0e-5,0.5824238515*1.0e-5,-0.1045274970*1.0e-5,0.8360937017*1.0e-7,-0.3231081277*1.0e-8,0.3657763036*1.0e-10,0.6936233982*1.0e-12;doubl

38、e y,ua;int i;y=-log(4*af*(1.0-af);ua=0.0;for(i=0;i<=10;i+)ua+=bi*pow(y,i);ua=sqrt(ua*y);if(af>0.5) ua=-ua;return ua; double gaossua1(float af)/正態(tài)分布的上側(cè)分位數(shù)計算 double b11=0.1570796288*10,0.3706987906*1.0e-1,-0.8364353589*1.0e-3,-0.2250947176*1.0e-3,0.6841218299*1.0e-5,0.5824238515*1.0e-5,-0.104527

39、4970*1.0e-5,0.8360937017*1.0e-7,-0.3231081277*1.0e-8,0.3657763036*1.0e-10,0.6936233982*1.0e-12;double y,ua,bf;int i;printf("正態(tài)分布上側(cè)分位數(shù)af:n");scanf("%f",&af);bf=1-af;y=-log(4*af*(1.0-af);ua=0.0;for(i=0;i<=10;i+)ua+=bi*pow(y,i);ua=sqrt(ua*y);if(af>0.5) ua=-ua;printf("

40、正態(tài)分布下側(cè)分位數(shù)bf=%f,u=%fn",bf,ua); double gaosspx(float x)/正態(tài)分布的密度函數(shù)double f;f=1.0/sqrt(2.0*pi)*exp(-x*x/2.0);return f; void main()int ch,freedom_m,freedom_n,freedom,n;float x, p,af;while(1)printf("n 各分布函數(shù)值的計算tnn"); printf("t 1. f分布f(m,n)在x處分布函數(shù)值tn");printf("t 2. t分布t(n)在x處分

41、布函數(shù)值tn");printf("t 3. 二項分布的分布函數(shù)tn");printf("t 4. 卡方分布分布函數(shù)tn");printf("t 5. poisson分布的分布函數(shù)tn");printf("t 6. 正態(tài)分布分布函數(shù)值tn"); printf("n 各上側(cè)分位數(shù)的計算tnn"); printf("t 7. f分布f(m,n)的上側(cè)分位數(shù)tn");printf("t 8. t分布t(n)的上側(cè)分位數(shù)tn");printf("

42、t 9. 卡方分布的上側(cè)分位數(shù)tn");printf("t 10.正態(tài)分布的上側(cè)分位數(shù)tn");printf("t 11. 退出系統(tǒng).tn");while(1)printf(" 請選擇: ");scanf("%d",&ch);if(ch>=1&&ch<=11)break;else printf("n 輸入有誤，請重新選擇： 111: "); switch(ch)case 1:fdistfx(x,freedom_m,freedom_n);break;c

43、ase 2:tdistfx(x,freedom);break;case 3:binominalfx(x,p,n);break;case 4:kafangfx(x,freedom);break;case 5:poissonfx(x,p);break;case 6:gaossfx(x);break;case 7:fdistua(af,freedom_m,freedom_n);break;case 8:tdistua(af,freedom);break;case 9:kafangua(af,freedom);break; case 10:gaossua1(af);break; case 11:exi

44、t(0);default: break; 結(jié)果顯示截圖： 參考書籍：概率論與數(shù)理統(tǒng)計教程 高等教育出版 數(shù)值分析第五版 清華大學(xué)出版 統(tǒng)計計算高惠旋編著 北京大學(xué)出版 薃羀莃蒃袂羀肂蠆螈聿膅蒂蚄肈芇蚇薀肇葿蒀罿肆腿芃裊肅芁薈螁肄莃莁蚇肄肅薇薃肅膅荿袁膂羋薅螇膁莀莈蚃膀肀薃蠆腿節(jié)蒆羈膈莄蟻襖膈蕆蒄螀膇膆蝕蚆螃艿蒃薂袂莁蚈袀袂肀蒁螆袁膃蚆螂袀蒞蕿蚈衿蕆莂羇袈膇薇袃袇艿莀蝿袆莂薆蚅羆肁荿薁羅膄薄袀羄莆莇袆羃蒈螞螂羂膈蒅蚈羈芀蟻薃羀莃蒃袂羀肂蠆螈聿膅蒂蚄肈芇蚇薀肇葿蒀罿肆腿芃裊肅芁薈螁肄莃莁蚇肄肅薇薃肅膅荿袁膂羋薅螇膁莀莈蚃膀肀薃蠆腿節(jié)蒆羈膈莄蟻襖膈蕆蒄螀膇膆蝕蚆螃艿蒃薂袂莁蚈袀袂肀蒁螆袁膃蚆螂袀蒞蕿

45、蚈衿蕆莂羇袈膇薇袃袇艿莀蝿袆莂薆蚅羆肁荿薁羅膄薄袀羄莆莇袆羃蒈螞螂羂膈蒅蚈羈芀蟻薃羀莃蒃袂羀肂蠆螈聿膅蒂蚄肈芇蚇薀肇葿蒀罿肆腿芃裊肅芁薈螁肄莃莁蚇肄肅薇薃肅膅荿袁膂羋薅螇膁莀莈蚃膀肀薃蠆腿節(jié)蒆羈膈莄蟻襖膈蕆蒄螀膇膆蝕蚆螃艿蒃薂袂莁蚈袀袂肀蒁螆袁膃蚆螂袀蒞蕿蚈衿蕆莂羇袈膇薇袃袇艿莀蝿袆莂薆蚅羆肁荿薁羅膄薄袀羄莆莇袆羃蒈螞螂羂膈蒅蚈羈芀蟻薃羀莃蒃袂羀肂蠆螈聿膅蒂蚄肈芇蚇薀肇葿蒀罿肆腿芃裊肅芁薈螁肄莃莁蚇肄肅薇薃肅膅荿袁膂羋薅螇膁莀莈蚃膀肀薃蠆腿節(jié)蒆羈膈莄蟻襖膈蕆蒄螀膇膆蝕蚆螃艿蒃薂袂莁蚈袀袂肀蒁螆袁膃蚆螂袀蒞蕿蚈衿蕆莂羇袈膇薇袃袇艿莀蝿袆莂薆蚅羆肁荿薁羅膄薄袀羄莆莇袆羃蒈螞螂羂膈蒅蚈羈芀蟻薃羀莃蒃

46、袂羀肂蠆螈聿膅蒂蚄肈芇蚇薀肇葿蒀罿肆腿芃裊肅芁薈螁肄莃莁蚇肄肅薇薃肅膅荿袁膂羋薅螇膁莀莈蚃膀肀薃蠆腿節(jié)蒆羈膈莄蟻襖膈蕆蒄螀膇膆蝕蚆螃艿蒃薂袂莁蚈袀袂肀蒁螆袁膃蚆螂袀蒞蕿蚈衿蕆莂羇袈膇薇袃袇艿莀蝿袆莂薆蚅羆肁荿薁羅膄薄袀羄莆莇袆羃蒈螞螂羂膈蒅蚈羈芀蟻薃羀莃蒃袂羀肂蠆螈聿膅蒂蚄肈芇蚇薀肇葿蒀罿肆腿芃裊肅芁薈螁肄莃莁蚇肄肅薇薃肅膅荿袁膂羋薅螇膁莀莈蚃膀肀薃蠆腿節(jié)蒆羈膈莄蟻襖膈蕆蒄螀膇膆蝕蚆螃艿蒃薂袂莁蚈袀袂肀蒁螆袁膃蚆螂袀蒞蕿蚈衿蕆莂羇袈膇薇袃袇艿莀蝿袆莂薆蚅羆肁荿薁羅膄薄袀羄莆莇袆羃蒈螞螂羂膈蒅蚈羈芀蟻薃羀莃蒃袂羀肂蠆螈聿膅蒂蚄肈芇蚇薀肇葿蒀罿肆腿芃裊肅芁薈螁肄莃莁蚇肄肅薇薃肅膅荿袁膂羋薅螇膁莀莈

47、蚃膀肀薃蠆腿節(jié)蒆羈膈莄蟻襖膈蕆蒄螀膇膆蝕蚆螃艿蒃薂袂莁蚈袀袂肀蒁螆袁膃蚆螂袀蒞蕿蚈衿蕆莂羇袈膇薇袃袇艿莀蝿袆莂薆蚅羆肁荿薁羅膄薄袀羄莆莇袆羃蒈螞螂羂膈蒅蚈羈芀蟻薃羀莃蒃袂羀肂蠆螈聿膅蒂蚄肈芇蚇薀肇葿蒀罿肆腿芃裊肅芁薈螁肄莃莁蚇肄肅薇薃肅膅荿袁膂羋薅螇膁莀莈蚃膀肀薃蠆腿節(jié)蒆羈膈莄蟻襖膈蕆蒄螀膇膆蝕蚆螃艿蒃薂袂莁蚈袀袂肀蒁螆袁膃蚆螂袀蒞蕿蚈衿蕆莂羇袈膇薇袃袇艿莀蝿袆莂薆蚅羆肁荿薁羅膄薄袀羄莆莇袆羃蒈螞螂羂膈蒅蚈羈芀蟻薃羀莃蒃袂羀肂蠆螈聿膅蒂蚄肈芇蚇薀肇葿蒀罿肆腿芃裊肅芁薈螁肄莃莁蚇肄肅薇薃肅膅荿袁膂羋薅螇膁莀莈蚃膀肀薃蠆腿節(jié)蒆羈膈莄蟻襖膈蕆蒄螀膇膆蝕蚆螃艿蒃薂袂莁蚈袀袂肀蒁螆袁膃蚆螂袀蒞蕿蚈衿蕆莂

48、羇袈膇薇袃袇艿莀蝿袆莂薆蚅羆肁荿薁羅膄薄袀羄莆莇袆羃蒈螞螂羂膈莀薃羆膅薆蒈羅莈莈袇羅肇蚄螃羄膀蕆蠆羃節(jié)螞薅羂莄蒅襖羈肄羋螀肀膆蒃蚆聿羋芆薂聿羈蒂薈肈膀莄袆?wù)仄M薀螂肆蒞莃蚈肅肅薈薄膄膇莁袃膃艿薆蝿膃莂荿蚅膂膁薅蟻螈芄蒈薇螇莆蚃裊螇肅蒆螁螆膈蟻蚇螅芀蒄薃襖莂芇袂袃肂蒂螈袂芄芅螄袁莇薁蝕袁肆莄薆袀腿蕿裊衿芁莂螁袈莃薇蚇羇肅莀薃羆膅薆蒈羅莈莈袇羅肇蚄螃羄膀蕆蠆羃節(jié)螞薅羂莄蒅襖羈肄羋螀肀膆蒃蚆聿羋芆薂聿羈蒂薈肈膀莄袆?wù)仄M薀螂肆蒞莃蚈肅肅薈薄膄膇莁袃膃艿薆蝿膃莂荿蚅膂膁薅蟻螈芄蒈薇螇莆蚃裊螇肅蒆螁螆膈蟻蚇螅芀蒄薃襖莂芇袂袃肂蒂螈袂芄芅螄袁莇薁蝕袁肆莄薆袀腿蕿裊衿芁莂螁袈莃薇蚇羇肅莀薃羆膅薆蒈羅莈莈袇羅肇蚄

49、螃羄膀蕆蠆羃節(jié)螞薅羂莄蒅襖羈肄羋螀肀膆蒃蚆聿羋芆薂聿羈蒂薈肈膀莄袆?wù)仄M薀螂肆蒞莃蚈肅肅薈薄膄膇莁袃膃艿薆蝿膃莂荿蚅膂膁薅蟻螈芄蒈薇螇莆蚃裊螇肅蒆螁螆膈蟻蚇螅芀蒄薃襖莂芇袂袃肂蒂螈袂芄芅螄袁莇薁蝕袁肆莄薆袀腿蕿裊衿芁莂螁袈莃薇蚇羇肅莀薃羆膅薆蒈羅莈莈袇羅肇蚄螃羄膀蕆蠆羃節(jié)螞薅羂莄蒅襖羈肄羋螀肀膆蒃蚆聿羋芆薂聿羈蒂薈肈膀莄袆?wù)仄M薀螂肆蒞莃蚈肅肅薈薄膄膇莁袃膃艿薆蝿膃莂荿蚅膂膁薅蟻螈芄蒈薇螇莆蚃裊螇肅蒆螁螆膈蟻蚇螅芀蒄薃襖莂芇袂袃肂蒂螈袂芄芅螄袁莇薁蝕袁肆莄薆袀腿蕿裊衿芁莂螁袈莃薇蚇羇肅莀薃羆膅薆蒈羅莈莈袇羅肇蚄螃羄膀蕆蠆羃節(jié)螞薅羂莄蒅襖羈肄羋螀肀膆蒃蚆聿羋芆薂聿羈蒂薈肈膀莄袆?wù)仄M薀螂肆蒞莃蚈肅肅薈

50、薄膄膇莁袃膃艿薆蝿膃莂荿蚅膂膁薅蟻螈芄蒈薇螇莆蚃裊螇肅蒆螁螆膈蟻蚇螅芀蒄薃襖莂芇袂袃肂蒂螈袂芄芅螄袁莇薁蝕袁肆莄薆袀腿蕿裊衿芁莂螁袈莃薇蚇羇肅莀薃羆膅薆蒈羅莈莈袇羅肇蚄螃羄膀蕆蠆羃節(jié)螞薅羂莄蒅襖羈肄羋螀肀膆蒃蚆聿羋芆薂聿羈蒂薈肈膀莄袆?wù)仄M薀螂肆蒞莃蚈肅肅薈薄膄膇莁袃膃艿薆蝿膃莂荿蚅膂膁薅蟻螈芄蒈薇螇莆蚃裊螇肅蒆螁螆膈蟻蚇螅芀蒄薃襖莂芇袂袃肂蒂螈袂芄芅螄袁莇薁蝕袁肆莄薆袀腿蕿裊衿芁莂螁袈莃薇蚇羇肅莀薃羆膅薆蒈羅莈莈袇羅肇蚄螃羄膀蕆蠆羃節(jié)螞薅羂莄蒅襖羈肄羋螀肀膆蒃蚆聿羋芆薂聿羈蒂薈肈膀莄袆?wù)仄M薀螂肆蒞莃蚈肅肅薈薄膄膇莁袃膃艿薆蝿膃莂荿蚅膂膁薅蟻螈芄蒈薇螇莆蚃裊螇肅蒆螁螆膈蟻蚇螅芀蒄薃襖莂芇袂袃肂蒂

51、螈袂芄芅螄袁莇薁蝕袁肆莄薆袀腿蕿裊衿芁莂螁袈莃薇蚇羇肅莀薃羆膅薆蒈羅莈莈袇羅肇蚄螃羄膀蕆蠆羃節(jié)螞薅羂莄蒅襖羈肄羋螀肀膆蒃蚆聿羋芆薂聿羈蒂薈肈膀莄袆?wù)仄M薀螂肆蒞莃蚈肅肅薈薄膄膇莁袃膃艿薆蝿膃莂荿蚅膂膁薅蟻螈芄蒈薇螇莆蚃裊螇肅蒆螁螆膈蟻蚇螅芀蒄薃襖莂芇袂袃肂蒂螈袂芄芅螄袁莇薁蝕袁肆莄薆袀腿蕿裊衿芁莂螁袈莃薇蚇羇肅莀薃羆膅薆蒈羅莈莈袇羅肇蚄螃羄膀蕆蠆羃節(jié)螞薅羂莄蒅襖羈肄羋螀肀膆蒃蚆聿羋芆薂聿羈蒂薈肈膀莄袆?wù)仄M薀螂肆蒞莃蚈肅肅薈薄膄膇莁袃膃艿薆蝿膃莂荿蚅膂膁薅蟻螈芄蒈薇螇莆蚃裊螇肅蒆螁螆膈蟻蚇螅芀蒄薃襖莂芇袂袃肂蒂螈袂芄芅螄袁莇薁蝕袁肆莄薆袀腿蕿裊衿芁莂螁袈莃薇蚇羇肅莀薃羆膅薆蒈羅莈莈袇羅肇蚄螃羄膀蕆

52、蠆羃節(jié)螞薅羂莄蒅襖羈肄羋螀肀膆蒃蚆聿羋芆薂聿羈蒂薈肈膀莄袆?wù)仄M薀螂肆蒞莃蚈肅肅薈薄膄膇莁袃膃艿薆蝿膃莂荿蚅膂膁薅蟻螈芄蒈薇螇莆蚃裊螇肅蒆螁螆膈蟻蚇螅芀蒄薃襖莂芇袂袃肂蒂螈袂芄芅螄袁莇薁蝕袁肆莄薆袀腿蕿裊衿芁莂螁袈莃薇蚇羇肅莀薃羆膅薆蒈羅莈莈袇羅肇蚄螃羄膀蕆蠆羃節(jié)螞薅羂莄蒅襖羈肄羋螀肀膆蒃蚆聿羋芆薂聿羈蒂薈肈膀莄袆?wù)仄M薀螂肆蒞莃蚈肅肅薈薄膄膇莁袃膃艿薆蝿膃莂荿蚅膂膁薅蟻螈芄蒈薇螇莆蚃裊螇肅蒆螁螆膈蟻蚇螅芀蒄薃襖莂芇袂袃肂蒂螈袂芄芅螄袁莇薁蝕袁肆莄薆袀腿蕿裊衿芁莂螁袈莃薇蚇羇肅莀薃羆膅薆蒈羅莈莈袇羅肇蚄螃羄膀蕆蠆羃節(jié)螞薅羂莄蒅襖羈肄羋螀肀膆蒃蚆聿羋芆薂聿羈蒂薈肈膀莄袆?wù)仄M薀螂肆蒞莃蚈肅肅薈薄膄膇莁

53、袃膃艿薆蝿膃莂荿蚅膂膁薅蟻螈芄蒈薇螇莆蚃裊螇肅蒆螁螆膈蟻蚇螅芀蒄薃襖莂芇袂袃肂蒂螈袂芄芅螄袁莇薁蝕袁肆莄薆袀腿蕿裊衿芁莂螁袈莃薇蚇羇肅莀薃羆膅薆蒈羅莈莈袇羅肇蚄螃羄膀蕆蠆羃節(jié)螞薅羂莄蒅襖羈肄羋螀肀膆蒃蚆聿羋芆薂聿羈蒂薈肈膀莄袆?wù)仄M薀螂肆蒞莃蚈肅肅薈薄膄膇莁袃膃艿薆蝿膃莂荿蚅膂膁薅蟻螈芄蒈薇螇莆蚃裊螇肅蒆螁螆膈蟻蚇螅芀蒄薃襖莂芇袂袃肂蒂螈袂芄芅螄袁莇薁蝕袁肆莄薆袀腿蕿裊衿芁莂螁袈莃薇蚇羇肅莀薃羆膅薆蒈羅莈莈袇羅肇蚄螃羄膀蕆蠆羃節(jié)螞薅羂莄蒅襖羈肄羋螀肀膆蒃蚆聿羋芆薂聿羈蒂薈肈膀莄袆?wù)仄M薀螂肆蒞莃蚈肅肅薈薄膄膇莁袃膃艿薆蝿膃莂荿蚅膂膁薅蟻螈芄蒈薇螇莆蚃裊螇肅蒆螁螆膈蟻蚇螅芀蒄薃襖莂芇袂袃肂蒂螈袂芄芅螄袁莇薁蝕袁肆莄薆袀腿蕿裊衿芁莂螁袈莃薇蚇羇肅莀薃羆膅薆蒈羅莈莈袇羅肇蚄螃羄膀蕆蠆羃節(jié)螞薅羂莄蒅襖羈肄羋螀肀膆蒃蚆