麻烦大神指点一下如何定义这个rhum， 成为一个合适的矩阵z？非常感谢！

hpucq

) R$ `3 K; z7 P1 s) `9 V4 A
; d+ p  ]; F4 b/ T9 ^, oclear
  u. f5 f' S' I/ w( G' I, wncinfo('rhum.mon.mean.nc');
y_len=60;
. S3 D# P( V6 Z% h( t7 ~t_len=854;
6 c# b6 G" Z& |7 [  l/ Yz_index=7
x_index=120
$ s# f, y7 c9 _0 ^* [  t' Hx = ncread('rhum.mon.mean.nc','lon');%x
z = ncread('rhum.mon.mean.nc', 'level');%z
y = ncread('rhum.mon.mean.nc', 'lat');%y
t = ncread('rhum.mon.mean.nc', 'time');%t
$ j# I* A% B0 O: t' @: Grhum = ncread('rhum.mon.mean.nc', 'rhum');
figure; hold on;
for i = 30:y_len
  for j=10:t_len
[ii,jj]=meshgrid(j,i)
array(i,j)=rhum(x_index,i,z_index,j);
contour(i,j,z,[rhum(x_index,i,z_index,j)])
5 j+ D( n- w. \( v% v# j* O  end
end
$ f6 g  K/ U0 w
错误使用 contour (line 84)
Z 的大小必须为 2x2 或更大。
5 e* R9 N) m) G. I. C
0 e1 P" }/ n( e8 i/ u& e5 h& r
这个contour函数里 i,j很好定义，主要这是这个rhum，也就是contour里的矩阵z，rhum现在是固定两维，两个自变量。麻烦大神指点一下如何定义这个rhum， 成为一个合适的矩阵z？非常感谢！

Colbie

for i = 30:y_len
) I7 ?8 q. ?6 x& ~3 X6 K* E) K  for j=10:t_len
[ii,jj]=meshgrid(j,i)
6 c6 Y$ U) }5 x' u+ l! H$ g# Yarray(i,j)=rhum(x_index,i,z_index,j);
contour(i,j,z,[rhum(x_index,i,z_index,j)])
  end
8 ?' ~7 g  R6 Z9 L- M' E$ Iend

这里不对，不应在循环里面嵌套contour，因为对循环每一步而言，i和j都只有一个数值，meshgrid不能把它们展开成矩阵，而且contour(i,j,z)里面的这个z也没赋过值。
* u7 F2 T9 m8 }4 v& A0 X5 n) ?& i举个正确用法的例子，循环每一步计算Z的一个元素，全部计算完毕后，用meshgrid把横纵坐标向量展开成矩阵，最后在循环外面用contour绘图
% X9 B  {; ^- `" Gdraw_x = 1:1/4:7;
8 }  w; i9 k6 `: x3 F/ ?draw_y = 2:1/3:6;
f = @(x,y) x.^(exp( -y.^2/50 ));
# B. N( z+ K$ a, i/ {4 cfor ii = 1:1:numel(draw_x)
    for jj = 1:1:numel(draw_y)
3 [6 z- \  q+ B        Z(jj,ii) = f(draw_x(ii),draw_y(jj));
. d5 Y/ P7 t: F% b    end
end
[draw_x,draw_y] = meshgrid(draw_x,draw_y);
contour (draw_x,draw_y,Z);