程序显然有问题,在循环里只对s(100,100)赋值,其他元素当然为0.
猜测程序的用意,可能是要把下面这句
s(100,100)=sum(find(z<=0.05));
改成
s(round(x*100+1),round(y*100+1))=sum(find(z<=0.05));
但由于没有研究原题,这样改是否符合题意我就说不准了。
按上面的程序,一开始s没有初始化,第一次执行
s(100,100)=sum(find(z<=0.05));
这句话时会生成一个100*100的矩阵s,其他元素都是0,(100,100)处是上式算出来的值;然后后面的循环过程中还会一遍遍的覆盖这个值。。这是要干什么额