原網(wǎng)址如下:
https://ww2.mathworks.cn/help/matlab/creating_plots/highlight-specific-contour-levels-1.html
代碼如下。在MATLAB可以運(yùn)行,在天元不能運(yùn)行。
Z = peaks(100);
zmin = floor(min(Z(:)));
zmax = ceil(max(Z(:)));
zinc = (zmax - zmin) / 40;
zlevs = zmin:zinc:zmax;
figure
contour(Z,zlevs)
zindex = zmin:2:zmax;
hold on
contour(Z,zindex,'LineWidth',2)
hold off
根據(jù)北太天元中contour函數(shù)的功能情況做一些調(diào)整即可正確繪制等高線了:第一處zlevs包含超出Z實(shí)際數(shù)據(jù)范圍的層級(因floor(min(Z))可能小于實(shí)際最小值),需使用min(Z(:))和max(Z(:))直接獲取數(shù)據(jù)范圍;第二處zinc計(jì)算可能使zlevs末位略微超過zmax(如zmax=8.0但實(shí)際數(shù)據(jù)最大為7.999),需要修改求zlevs的方式使得確保不越界,比如