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

hpucq

clear
ncinfo('rhum.mon.mean.nc');
: }; J7 `1 k) }y_len=60;
t_len=854;
" l) b, J8 R6 P7 g: [' yz_index=7
x_index=120
/ h3 S! e% Q+ p, y0 _x = ncread('rhum.mon.mean.nc','lon');%x
  U, U6 }. ~6 F3 cz = ncread('rhum.mon.mean.nc', 'level');%z
% V+ D' V: b$ [- F/ t0 p! Uy = ncread('rhum.mon.mean.nc', 'lat');%y
t = ncread('rhum.mon.mean.nc', 'time');%t
rhum = ncread('rhum.mon.mean.nc', 'rhum');
figure; hold on;
! m5 Q6 n* L5 Y# X) j. Mfor i = 30:y_len
  for j=10:t_len
[ii,jj]=meshgrid(j,i)
array(i,j)=rhum(x_index,i,z_index,j);
1 S8 t, c2 c5 N5 E' q) _0 g4 I# ncontour(i,j,z,[rhum(x_index,i,z_index,j)])
  end
/ n4 _& _' d# \2 ]0 V/ M6 a& lend

. w5 F  x1 M" f  j3 s* V" e9 o0 o+ L错误使用 contour (line 84)
Z 的大小必须为 2x2 或更大。


这个contour函数里 i,j很好定义，主要这是这个rhum，也就是contour里的矩阵z，rhum现在是固定两维，两个自变量。麻烦大神指点一下如何定义这个rhum， 成为一个合适的矩阵z？非常感谢！

Colbie

for i = 30:y_len
2 P( b- a1 S  J6 V% n2 b0 K- T* ~  for j=10:t_len
! o5 P* D7 E7 _: B, `[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)])
  end
! _1 b8 |" |/ }0 b2 T" t; d0 hend

, F8 ]" j! o8 V2 T* {* P0 L这里不对，不应在循环里面嵌套contour，因为对循环每一步而言，i和j都只有一个数值，meshgrid不能把它们展开成矩阵，而且contour(i,j,z)里面的这个z也没赋过值。
4 v6 H0 a6 s6 s- |# j; y举个正确用法的例子，循环每一步计算Z的一个元素，全部计算完毕后，用meshgrid把横纵坐标向量展开成矩阵，最后在循环外面用contour绘图
draw_x = 1:1/4:7;
draw_y = 2:1/3:6;
  u5 b. B* a/ x5 ^f = @(x,y) x.^(exp( -y.^2/50 ));
for ii = 1:1:numel(draw_x)
    for jj = 1:1:numel(draw_y)
        Z(jj,ii) = f(draw_x(ii),draw_y(jj));
6 v) ^! ~5 O+ y7 r  w    end
6 N  h1 F: m9 x" k/ J# q$ Cend
[draw_x,draw_y] = meshgrid(draw_x,draw_y);
* a( ]! b$ j3 Y, `& \contour (draw_x,draw_y,Z);