matlab - Matlab：绘制平滑阴影的计数图

Question

在这个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))

而不是imagesc你也可以使用surf：

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