for m=[301:330]
Sum=0;
for k=[0,m-301]
A=random('poisson',5,1,m-300);
Pk=A(k+1);
sum=Pk*(m-300-k);
Sum=sum;
end
B=[B Sum];
end
x=[301:330]
y=19x/3600-10/9-B/100;
plot(x,y,'r')
hold on
BB=[]
for m=[301:330]
Sum=0;
for k=[0,m-301]
A=random('poisson',10,1,m-300);
Pk=A(k+1);
sum=Pk(m-300-k);
Sum=sum;
end
BB=[BB Sum];
end
y2=x/200-13/16-BB/100;
plot(x,y2)
![img](
答思路和尝试过的方法
难道pk不是用泊松分布吗