خب با اين حساب من كل برنامتون رو اينطوري خلاصه كردم
clear
clc
kd1=eye(1,100);
k=[-1,0,1];
kd=k(randi([1,3],1,100 ))
Prate=0.85+(0.9-0.85)*rand(1,100);
ereq1=kd.*Prate;
Ereq.d=sum(ereq1);
Id=zeros(1,100);
Id(kd==0)=1
هرجا kd صفر بود Id ميشه يك و در غير اين صورت همون صفر باقي ميمونه