vpasolve 总是返回一个空的解决方案结构

如何解决vpasolve 总是返回一个空的解决方案结构

我需要解决这个可行性问题。

代码如下：

pc = zeros([1,n]);
ac = zeros([1,n]);
al = zeros([m,n]);
pc(:)= 10 + rand([1,n])*10; 
pl = zeros([m,n]);
pl = 1 + rand([m,n])* 10; 
ac(1,:) = 0.9;
al(1,:) = 0.5;
al(2,:) = 0.6;
al(3,:) = 0.7; 
A = [al ; ac];
P = [pl ; pc];
T = 10000
syms x [k,n]
assume(x >= 0)
eqn1 = sum(sum(P(1:m,:).*x(1:m,:))) <= T;
eqn2 = sum(P(k,:).*x(k,:))<= T;

for i = 1:n
     val(i) = sum(x(:,i));
     eqn(i) = val(i)==1;
end 
eqns = [eqn1,eqn2,eqn]
sol = vpasolve(eqns,x )

由于T足够大，我认为这个问题可能有一个可行的解决方案，但显然似乎没有，因为sol是由struct [0x1 sym]组成的

有人可以帮助我了解问题出在哪里吗？

谢谢