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

hpucq

0 V+ `3 r6 g5 s( v. p; c
2 O; X% q) {  d2 C# y) p: xclear
9 j3 ?5 Y* r" K2 Tncinfo('rhum.mon.mean.nc');
- K( ~) L. d/ z7 \6 Uy_len=60;
/ u) p2 Q: Q5 W4 O% x+ o3 x+ gt_len=854;
z_index=7
x_index=120
x = ncread('rhum.mon.mean.nc','lon');%x
z = ncread('rhum.mon.mean.nc', 'level');%z
# j/ P, @* q6 T) sy = ncread('rhum.mon.mean.nc', 'lat');%y
t = ncread('rhum.mon.mean.nc', 'time');%t
, @2 N1 [" \! p& l( Xrhum = ncread('rhum.mon.mean.nc', 'rhum');
. w" ~$ J7 j0 s  _- T0 e. Bfigure; hold on;
7 z, }6 Z( N- J* s2 t/ w) U- ifor i = 30:y_len
  for j=10:t_len
" C2 a( q4 w! y[ii,jj]=meshgrid(j,i)
array(i,j)=rhum(x_index,i,z_index,j);
+ G: o. S1 F0 t8 f) ?' Jcontour(i,j,z,[rhum(x_index,i,z_index,j)])
  end
, g$ A. u! Z+ m$ bend

错误使用 contour (line 84)
Z 的大小必须为 2x2 或更大。
7 d& c) Z# P: u

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

Colbie

for i = 30:y_len
  for j=10:t_len
[ii,jj]=meshgrid(j,i)
# y' X; Y( B% {' S1 |array(i,j)=rhum(x_index,i,z_index,j);
: B  p3 _2 n8 a2 s3 W9 kcontour(i,j,z,[rhum(x_index,i,z_index,j)])
- b9 z: \, A0 i( s  end
end

这里不对，不应在循环里面嵌套contour，因为对循环每一步而言，i和j都只有一个数值，meshgrid不能把它们展开成矩阵，而且contour(i,j,z)里面的这个z也没赋过值。
6 ^) `7 L" N" X举个正确用法的例子，循环每一步计算Z的一个元素，全部计算完毕后，用meshgrid把横纵坐标向量展开成矩阵，最后在循环外面用contour绘图
draw_x = 1:1/4:7;
  ?9 q- y2 F& D9 L) M* H/ Rdraw_y = 2:1/3:6;
  Y- {6 L# a6 Z9 ef = @(x,y) x.^(exp( -y.^2/50 ));
for ii = 1:1:numel(draw_x)
    for jj = 1:1:numel(draw_y)
- ]1 u( I& u0 m        Z(jj,ii) = f(draw_x(ii),draw_y(jj));
5 m/ X  ^: I% W    end
end
[draw_x,draw_y] = meshgrid(draw_x,draw_y);
3 f* v" M- D4 }4 ?7 y" d; ~contour (draw_x,draw_y,Z);