随机误差的统计规律程序代码.docVIP

x=[8.56 8.53 8.50 8.56 8.50 8.53 8.59 8.50 8.56 8.60 8.53 8.53 8.62 8.53 8.59 8.63 8.63 8.59 8.50 8.44 8.56 8.59 8.53 8.50 8.68 8.38 8.53 8.59 8.54 8.50 8.62 8.59 8.63 8.61 8.47 8.62 8.65 8.56 8.62 8.66 8.59 8.55 8.63 8.65 8.75 8.53 8.44 8.50 8.50 8.62 8.47 8.53 8.56 8.53 8.53 8.47 8.62 8.66 8.56 8.59 8.56 8.54 8.56 8.62 8.56 8.66 8.50 8.68 8.50 8.47 8.66 8.59 8.50 8.65 8.56 8.53 8.59 8.53 8.63 8.53 8.53 8.66 8.50 8.66 8.66 8.69 8.62 8.50 8.53 8.75 8.63 8.53 8.44 8.56 8.50 8.69 8.69 8.57 8.53 8.47 8.53 8.56 8.66 8.53 8.53]; X=x/5; P=sum(X)/105; s=sqrt((sum(X.^2)-P*sum(X))/104); A=min(X); B=max(X); a=P-3.62*s; b=P+3.62*s; X1=[X(1:10)]; X2=[X(1:20)]; X3=[X(1:30)]; X4=[X(1:40)]; X5=[X(1:50)]; X6=[X(1:60)]; X7=[X(1:70)]; X8=[X(1:80)]; X9=[X(1:90)]; X10=[X(1:100)]; X11=[X(1:105)]; P1=sum(X1)/10; P2=sum(X2)/20; P3=sum(X3)/30; P4=sum(X4)/40; P5=sum(X5)/50; P6=sum(X6)/60; P7=sum(X7)/70; P8=sum(X8)/80; P9=sum(X9)/90; P10=sum(X10)/100; P11=sum(X11)/105; s1=sqrt((sum(X1.^2)-P1*sum(X1))/9); s2=sqrt((sum(X2.^2)-P2*sum(X2))/19); s3=sqrt((sum(X3.^2)-P3*sum(X3))/29); s4=sqrt((sum(X4.^2)-P4*sum(X4))/39); s5=sqrt((sum(X5.^2)-P5*sum(X5))/49); s6=sqrt((sum(X6.^2)-P6*sum(X6))/59); s7=sqrt((sum(X7.^2)-P7*sum(X7))/69); s8=sqrt((sum(X8.^2)-P8*sum(X8))/79); s9=sqrt((sum(X9.^2)-P9*sum(X9))/89); s10=sqrt((sum(X10.^2)-P10*sum(X10))/99); s11=sqrt((sum(X11.^2)-P11*sum(X11))/104); P0=[P1 P2 P3 P4 P5 P6 P7 P8 P9 P10 P11]; s0=[s1 s2 s3 s4 s5 s6 s7 s8 s9 s10 s11]; N=[10 20 30 40 50 60 70 80 90 100 105]; plotyy(N,P0,N,s0); xlabel(n); text(5,P+0.001,P/s); text(110,P+0.001,s(x)/s); title(统计随测量次数的增加P及s（x）的变化图); E=(B-A)/7; L0=A; L1=A+E; L2=A+2*E; L3=A+3*E; L4=A+4*E; L5=A+5*E; L6=A+6*E; L7=A+7*E; n1=0; for i=1:105 n=((X(i)=L0)(X(i)L1)) n1=n1+n end n2=0; for i=1:105 n=((X(i)=L1)(X(i)L2)) n2=n2+n end n3=0; for i=1:105 n=((X(i)=L2)(X(i)L3)) n3=n3+n end n4=0; for i=1:105 n=((X(i)=L3)(X(i)L4)) n4=n4+n end n5=0; for i=1:105 n=((X(i)=L4)(X(i)L5)) n5=n5+