这是一个计算中断概率的程序,程序运行到q(i)=1-exp(-(Ncsma*Prs(1-PrR)))就报错Indexmustbeapositiveintegerorlogical.,q=zeros(size(t));C=zeros(size(t));symsvgm1[yyy1]=F();eq2=10^(-6)+yyy1;df2=solve(eq2);rv=R*(rs*df2)^0.25;

Question

这是一个计算中断概率的程序,程序运行到q(i)=1-exp(-(Ncsma*Prs(1-PrR)))就报错Index must be a positive integer or logical.,
q=zeros(size(t));
C=zeros(size(t));
syms v g m1
[yyy1] = F();
eq2=10^(-6)+yyy1;
df2=solve(eq2);
rv=R*(rs*df2)^0.25; %观察区域的计算 Fb是随机变量g(xj,0)/g(xi,0)的ccdf gb是随机变量g(xj,0)/g(xi,0)的pdf 参考文献29断概率的求解
Fg=((Bs*g)^2/n)/(rv^2);%随机变量z的cdf
fg=g^(a1-1)*gamma((a1+1)/b1)/((a1+1)^2*a1*gamma(a1)*gamma(1/b1));%修改之后的
Hg=Fg*fg;
rg=rv^n/Bs*R^n;
Prs=int(Hg,g,0,rg)+int(fg,g,rg,20);%在观察范围内一个发射机发射的信号就可以对接收信号带来破坏性的概率
yy2=g^(m1-1)*gamma(2*m1)/(g+1)^2*m1;
V=subs(fm,'m',1);%和g0的值是一样的
V1=subs(fm,'m','v');
V2=subs(yy2,'g','V1');
rv1=rv-R;
rv2=rv+R;
%if(0

Answer 1

q(i)=1-exp(-(Ncsma*Prs(1-PrR)))中的1-PrR估计有问题吧.PrR是积分得到的,积分的值如果不是整数这个式子就执行不了了.