I1=qfuncinv(BER0)*(2./(sqrt(2.^(1-1).*SNR.*1)));
I2=qfuncinv(BER0)*(2./(sqrt(2.^(2-1).*SNR.*2)));
syms x;
K=besselk(a-b,2*sqrt(a*b*sqrt(x./SNR)));
f=((a.*b)^((a+b)./2))./(gamma(a).*gamma(b).*SNR.^((a+b)./4))*x.^((a+b)./4-1).*K;
F1=int(f,x,0,I1);
F2=int(f,x,0,I2);
F=[F1;F2];
会出现这样的情况是为什么?求指教