我和楼主写法与遇到的问题一致，说一下我的解决方法，希望有所帮助～
传入的n,f=0.1,50先不用动，用它们计算t,b,r,l时保证t>0,b=-t,r>0,l=-r就好。之后再将n,f各自取反，即n,f=-0.1,-50，这样就能得到正确的P矩阵。