你可以用matlab求解，代码如下：
syms ts tb ns1 nb1 c fs fb pb1 pb2 ps1 ps2
[pb1,pb2,ps1,ps2] = solve(2*ts*ns1-pb1*nb1-ps1*nb1+2*c*nb1+fs,2*tb*nb1-pb1*ns1-ps1*ns1+2*c*ns1+fb,2*ts*ns1-pb2*nb1-ps2*nb1+2*c*nb1+pb2+ps2-2*ts-2*c-fs,2*tb*nb1-pb2*ns1-ps2*ns1+2*c*ns1+pb2+ps2-2*tb-2*c-fb)
结果如下：
pb1 =

-(fb*ns1 - fs*nb1 - 2*c*nb1^2 + 2*c*ns1^2 + nb1^2*ps1 - ns1^2*ps1 + 2*fb*nb1^2*ns1 - 2*fs*nb1*ns1^2 + 2*nb1*ns1^2*ps1 - 2*nb1^2*ns1*ps1 - 3*fb*nb1*ns1 + 3*fs*nb1*ns1 - 4*c*nb1*ns1^2 + 4*c*nb1^2*ns1)/(- 2*nb1^2*ns1 + nb1^2 + 2*nb1*ns1^2 - ns1^2)


pb2 =

-(fb*nb1 - fs*ns1 - 2*c*nb1^2 + 2*c*ns1^2 - 2*fb*nb1^2 + 2*fs*ns1^2 + nb1^2*ps2 - ns1^2*ps2 + 2*fb*nb1^2*ns1 - 2*fs*nb1*ns1^2 + 2*nb1*ns1^2*ps2 - 2*nb1^2*ns1*ps2 - fb*nb1*ns1 + fs*nb1*ns1 - 4*c*nb1*ns1^2 + 4*c*nb1^2*ns1)/(- 2*nb1^2*ns1 + nb1^2 + 2*nb1*ns1^2 - ns1^2)


ps1 =

(fs*ns1 - fb*nb1 + fb*ns1^2 - 3*fs*ns1^2 + 2*fs*ns1^3 - 2*fb*nb1*ns1^2 + 2*fb*nb1*ns1)/(2*(- 2*nb1^2*ns1 + nb1^2 + 2*nb1*ns1^2 - ns1^2))


ps2 =

-(fb*nb1 - fs*ns1 - 3*fb*nb1^2 + 2*fb*nb1^3 + fs*nb1^2 - 2*fs*nb1^2*ns1 + 2*fs*nb1*ns1)/(2*(- 2*nb1^2*ns1 + nb1^2 + 2*nb1*ns1^2 - ns1^2))

show me the equations pls

2*ts*ns1-pb1*nb1-ps1*nb1+2*c*nb1+fs=0
2*tb*nb1-pb1*ns1-ps1*ns1+2*c*ns1+fb=0
2*ts*ns1-pb2*nb1-ps2*nb1+2*c*nb1+pb2+ps2-2*ts-2*c-fs=0
2*tb*nb1-pb2*ns1-ps2*ns1+2*c*ns1+pb2+ps2-2*tb-2*c-fb=0
求解pb1,pb2,ps1,ps2

given these 4 equations, Explicit solution could not be found.

这里有一个问题：解里面并不全是参数，而是也出现了变量。比如ps1的解出现了ps2，而ps2也是需要求的，这个有什么好办法么