等高线图未显示结果

Question

我只是想知道我有以下代码：

[r,R] = meshgrid(2:0.01:10,-4:0.01:4);
Z=(((R).^2)-1 );
   A=2+ (R.*(r+r^-1));
   B=sqrt(Z).*(r-r^-1);
   C=((r+r^-1)./(sqrt(Z))).*acoth(A./B);
    figure
contourf(r,R,C);
    h=colorbar;
    set(get(h,'ylabel'),'string','\gamma P_P L')
    xlabel('$r$','Interpreter','latex','FontSize',18)
    ylabel('$D \over 2\sqrt{M}$','Interpreter','latex','FontSize',18)

我只是得到A和B的所有答案都是无限的，但是我不确定为什么吗？

Answer 1

你忘了点。 r^-1与inv(r)相同。 并且由于r的所有行都与第一个相同，因此它的行减少的梯形形式几乎具有全零，而逆数将具有无限值。 我相信你的意思

A=2+ (R.*(r+r.^-1));

或者，更明确的是，

A=2+ (R.*(r+ 1./r));

对于其余的r^-1 ，依此类推。