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

hpucq

6 [9 O: A4 ~- c0 @5 m
clear
ncinfo('rhum.mon.mean.nc');
9 ~1 `% r% W8 c! z' Wy_len=60;
t_len=854;
z_index=7
6 g; t7 v  F( p; B8 A5 d8 t0 w1 Ux_index=120
x = ncread('rhum.mon.mean.nc','lon');%x
# d; y7 H3 l% {. E& d* L: K3 Sz = ncread('rhum.mon.mean.nc', 'level');%z
y = ncread('rhum.mon.mean.nc', 'lat');%y
t = ncread('rhum.mon.mean.nc', 'time');%t
rhum = ncread('rhum.mon.mean.nc', 'rhum');
: \% j4 ^' j) y0 s4 l3 ~figure; hold on;
for i = 30:y_len
2 [$ D0 L9 Y! J* B2 x  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)])
6 m9 K4 J8 v% F6 O6 n5 `  end
0 ]. m* l: |1 ~0 N) k  D/ v2 m8 Gend

错误使用 contour (line 84)
Z 的大小必须为 2x2 或更大。


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

Colbie

for i = 30:y_len
* X+ E& K) R- p! n5 w  for j=10:t_len
; v9 ?* F2 u; P6 a/ r+ w[ii,jj]=meshgrid(j,i)
  z: l( L* |- p  e! Z& rarray(i,j)=rhum(x_index,i,z_index,j);
( Y' Y! W- Z6 Hcontour(i,j,z,[rhum(x_index,i,z_index,j)])
  end
end

这里不对，不应在循环里面嵌套contour，因为对循环每一步而言，i和j都只有一个数值，meshgrid不能把它们展开成矩阵，而且contour(i,j,z)里面的这个z也没赋过值。
5 P+ l4 w9 U- b0 E5 H7 [' r举个正确用法的例子，循环每一步计算Z的一个元素，全部计算完毕后，用meshgrid把横纵坐标向量展开成矩阵，最后在循环外面用contour绘图
1 |0 R/ f* D4 [: \+ |draw_x = 1:1/4:7;
( ?0 [2 [; L+ T6 ydraw_y = 2:1/3:6;
f = @(x,y) x.^(exp( -y.^2/50 ));
for ii = 1:1:numel(draw_x)
0 [. f! p! I9 J    for jj = 1:1:numel(draw_y)
5 l; P" k# x3 Q; u2 e        Z(jj,ii) = f(draw_x(ii),draw_y(jj));
    end
end
( P0 S7 ?' l0 M4 a[draw_x,draw_y] = meshgrid(draw_x,draw_y);
contour (draw_x,draw_y,Z);