Matlab：绘制平滑阴影的等高线图

Question

在this pdf的第61页，有图10.4。我正在尝试复制右下角的情节。但是，在他们的代码中，他们使用了 contourplot，这在 Matlab R2013a 中没有。

我尝试过使用contour、contourf、surf、surfc 等，但没有任何东西能复制它们的结果。

%This plots the stability region for RK-4

[X, Y] = meshgrid(-3:.1:1, -3:.1:3);
Mu = X + i*Y;
R = 1 + Mu + .5*Mu.^2 + (1/6)*Mu.^3 + (1/24)*Mu.^4;
Rhat = abs(R);
%contour(X, Y, Rhat)
%contourf(X, Y, Rhat)
%surfc(X, Y, Rhat)
surf(X, Y, Rhat)

在他们的代码中是否有模仿coutourplot 的绘图命令？我猜他们使用了不同版本的 Matlab，该命令可用。我在 R2013a 上。

这是他们情节中的图片：

Answer 1

从 PDF 来看，tt 似乎在值 1 处被截断。所以这看起来很相似：

[X, Y] = meshgrid(-3:.01:1, -3:.01:3);
Mu = X + i*Y;
R = 1 + Mu + .5*Mu.^2 + (1/6)*Mu.^3 + (1/24)*Mu.^4;
Rhat = abs(R);
Rhat = Rhat.*(Rhat<1);  %# here I truncate
imagesc([min(X(:)) max(X(:))],[min(Y(:)) max(Y(:))], Rhat)
colormap(flipud(gray))

您也可以使用surf 代替imagesc：

surf(X, Y, Rhat,'LineStyle','none')
view(0, 90);
colormap(flipud(gray))
xlim([-3 0.5]);
ylim([-3.2 3.2]);