| 1 |
cnh |
1.1 |
function [c,h]=mycoutourf(x,y,a,ci,varargin) |
| 2 |
|
|
% mycontourf(x,y,a,ci) |
| 3 |
|
|
|
| 4 |
|
|
|
| 5 |
|
|
% Choose contour levels |
| 6 |
|
|
if prod(size(ci))==1 |
| 7 |
|
|
amin=min(min(a)); |
| 8 |
|
|
amax=max(max(a)); |
| 9 |
|
|
cmin=floor(amin/ci)*ci; |
| 10 |
|
|
cmax=ceil(amax/ci)*ci; |
| 11 |
|
|
if cmin==cmax |
| 12 |
|
|
error(sprintf('cmin=cmax=%g',cmin)); |
| 13 |
|
|
end |
| 14 |
|
|
if cmin<0 & cmax>0 |
| 15 |
|
|
cmax=max(abs([cmin cmax])); |
| 16 |
|
|
cmin=-cmax; |
| 17 |
|
|
end |
| 18 |
|
|
CI=cmin:ci:cmax; |
| 19 |
|
|
else |
| 20 |
|
|
CI=ci; |
| 21 |
|
|
end |
| 22 |
|
|
cmin=min(CI); |
| 23 |
|
|
cmax=max(CI); |
| 24 |
|
|
|
| 25 |
|
|
[c,h]=contourf(x,y,a,CI,varargin{:}); |
| 26 |
|
|
caxis([cmin cmax]); |