matlab - 等高线图显示没有结果

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。