偶然误差的定义
对误差的处理
由于评阅试卷的人为因素很大，因此我们着重分析偶然误差。假设批阅试卷的精准度受人的疲劳度情况，并且不同年龄段的人工作相同的时间疲劳状况不一样。又因为疲劳状况是一个模糊概念，因此我们引入模糊数学概念。因素集为：年龄段Ui={25—35,35—45，45—55，55—65}.判断集为：Si={不疲劳，不太疲劳，比较疲劳，疲劳}。查阅资料得这几个年段的人在工作的时候对应的不同疲劳度的概率如下：
我们将改组数据设置为权重，存放于一个4X4中。设阅卷组的年龄分布为A。对应不同的年龄组成可以得到不同的疲劳度的分布。最终通过不同的年龄组成确定每个阅卷老师的阅卷数。计算出批阅试卷失误的概率的期望值a作为本题的偶然误差标准
clear;clc; m=0;
for l=1:1000 t=0; while(t==0)
R = normrnd(70,10,1,100); round(R); for i=1:100
if R(i)>0 && R(i)
end end if(t~=100) t=0; else
[h,p]=sort(R,'descend'); end end
d = unifrnd(0,3,1,8); e = ceil(unifrnd(-2,1,1,8)); e=5*e;
n(1,:)=[12,13,12,13,12,13,12,13]; for i=1:8 for j=2:8
n(j,i)=n(j-1,i)-ceil(n(j-1,i)*0.3); if(n(j-1,i)==2) break; end end end m=0; for i=1:8 for j=1:8 m=m+n(i,j); end end m=max(m) t=0;
s=zeros(8,13,8); for k=1:8 for i=1:8
for j=1:n(k,i) if(k==1)
while(s(i,j,k)100) s(i,j,k)=normrnd(R(j+sum(n(1,1:k+i-2)))+e(i),d(i),1);
t=t+1; if(t>50)
s(i,j,k)=100; t=0; break; end end
elseif (mod(k+i-1,8)==0)
while(s(i,p1(i,j,k-1),k)100)
s(i,p1(i,j,k-1),k)=normrnd(R(p1(i,j,k-1)+sum(n(1,1:i-1)))+e(8),d(8),1); t=t+1; if(t>50)
s(i,j,k)=100; t=0; break; end end else
while(s(i,p1(i,j,k-1),k)100)
s(i,p1(i,j,k-1),k)=normrnd(R(p1(i,j,k-1)+sum(n(1,1:i-1)))+e(mod(k+i-1,8)),d(mod(k+i-1,8)),1); t=t+1; if(t>50)
s(i,j,k)=100; t=0; break; end end end end
[h1(i,:,k),p1(i,:,k)]=sort(s(i,:,k),'descend'); end end for i=1:8
w(i)=mean(s(i,p1(i,1,8),:)); end
[w,p2]=sort(w,'descend');
for i=1:3
p2(i)=p1(p2(i),1,8)+sum(n(1,1:p2(i)-1)); end t=0; for i=1:3 for j=1:6
if(p2(i)==p(j)) t=t+1; break; end end end if(t~=3) m=m+1; end end c=m/1000 q=1-c
clear;clc; m=0;
for l=1:1000 t=0; while(t==0)
R = normrnd(70,10,1,100); round(R); for i=1:100
if R(i)>0 && R(i)
[h,p]=sort(R,'descend');
end end
d = unifrnd(0,3,1,8); e = ceil(unifrnd(-2,1,1,8)); e=5*e;
n(1,:)=[12,13,12,13,12,13,12,13]; for i=1:8
n(2,i)=n(1,i)-ceil(n(1,i)*0.4); n(3,i)=n(2,i)-ceil(n(2,i)*0.3); for j=4:8
n(j,i)=n(j-1,i)-ceil(n(j-1,i)*0.2); if(n(j,i)==0) n(j,i)=1; end end end a=0; for i=1:8 for j=1:8 a=a+n(j,i); end end t=0; a=max(a);
s=zeros(8,13,8); for k=1:8 for i=1:8
for j=1:n(k,i) if(k==1)
while(s(i,j,k)100) s(i,j,k)=normrnd(R(j+sum(n(1,1:k+i-2)))+e(i),d(i),1);
t=t+1; if(t>50)
s(i,j,k)=100; t=0; break; end
end
elseif (mod(k+i-1,8)==0)
while(s(i,p1(i,j,k-1),k)100)
s(i,p1(i,j,k-1),k)=normrnd(R(p1(i,j,k-1)+sum(n(1,1:i-1)))+e(8),d(8),1); t=t+1; if(t>50)
s(i,j,k)=100; t=0; break; end end else
while(s(i,p1(i,j,k-1),k)100)
s(i,p1(i,j,k-1),k)=normrnd(R(p1(i,j,k-1)+sum(n(1,1:i-1)))+e(mod(k+i-1,8)),d(mod(k+i-1,8)),1); t=t+1; if(t>50)
s(i,j,k)=100; t=0; break; end end end end
[h1(i,:,k),p1(i,:,k)]=sort(s(i,:,k),'descend'); end end for i=1:8
w(i)=mean(s(i,p1(i,1,8),:)); end
[w,p2]=sort(w,'descend'); for i=1:3
p2(i)=p1(p2(i),1,8)+sum(n(1,1:p2(i)-1)); end t=0;
for i=1:3 for j=1:6
if(p2(i)==p(j)) t=t+1; break; end end end if(t~=3) m=m+1; end end c=m/1000 q=1-c
n(1,:)=[12,13,12,13,12,13,12,13]; for i=1:8
for j=2:8
n(j,i)=n(j-1,i)-ceil(n(j-1,i)*0.3); if(n(j-1,i)==2) break; end end end for i=1:8
for k=2:j
n(k,i)=ceil(n(k-1,i)*0.3); end end a=0; for i=1:8 for k=1:j-3 a=a+n(k,i); end
a=a+2*n(j-2,i); a=a+3*n(j-1,i);
a=a+4*n(j,i); end a=max(a)