t(1)=0;
vs(1)=0;
vx(1)=0;
xr(1)=0;
nr(1)=0;
tp(1)=0;
mr(1)=0;
sr(1)=0;
st(1)=0;
xt(1)=0;
tc(1)=0;
tp(1)=0;
tp1(1)=0;
sr1(1)=0;
for n=2:1:26
t(n)=unidrnd(-3,3)
vs(n)=unidrnd(1/60,7/60)
vx(n)=unidrnd(1/15,1/5)
xr(n)=unidrnd(3,7)
nr(n)=unidrnd(20,50)
mr(n)=50-nr(n)+xr(n)
sr(n)=mr(n)
st(n)=sr(n).*vs(n)
xt(n)=xr(n).*vx(n)
tc(n)=st(n)+xt(n)
tp(n)=30-t(n-1)+t(n)-tc(n-1)+t(n)
tp1(n)=tp(n).*0.4
dr=poissrnd(tp1(n))
if dr<=mr
sr1(n)=0
else sr1(n)=dr(n)-mr(n)
end
end
想输出DR,SR,SR1这三个，估计各种错误，大神来帮助啊