小弟想要把avg值給固定住
以下是我的程式碼
for j=1:1000
for i=1:1000
data=normrnd(mu,sigma,m,n);%隨機產生常態機率的值(-3~3間)
data_total=[data_total; data];
if (data>3) || (data<-3)%若在-3與3之外,就停止迴圈
break
end
end
d=[d;i];
end
avg=mean(d)%ARL值
以上avg輸出結果都會落在350~380左右,現在我想要把avg值給固定住
就是我這次avg值是375,下次在run一次也是375,run個n次還是375,
於是在我的程式碼中加入了rng,如下
for j=1:1000
for i=1:1000
rng('default');
rng(i)
data=normrnd(mu,sigma,m,n);%隨機產生常態機率的值(-3~3間)
data_total=[data_total; data];
if (data>3) || (data<-3)%若在-3與3之外,就停止迴圈
break
end
end
d=[d;i];
end
avg=mean(d)%ARL值
以上的avg輸出結果都一值落在72,但並不是我要的結果,
因為我點開d矩陣一看,1~1000個都是72,
又點了data矩陣一看,發現1~72筆數據、73~144、145~216以此類堆...
也就是每72筆亂數都一樣,所以才會造成我的avg都一值落在72
主要是因為我的rng(i)碰上了if與break,所以第72筆break後就從1開始
而小弟要的結果是d矩陣不能都相同、rng(i)不受if迴圈控制、avg值也跟原始一樣落在
350~380左右,但下次再run時要跟上次一樣
謝謝!!麻煩大大了!!