想请问一下这个怎么改 ？

panphp

/ S1 E% @; @0 {( @  @/ u# @0 l2 B, wclear;
$ P0 @. V# @% }$ p, Px_len=144;
y_len=60;
z_len=17;
  i+ I9 e  O7 I& S2 et_len=854;
z_index=7;
x_index=120;

z = ncread('rhum.mon.mean.nc', 'level');%z
y = ncread('rhum.mon.mean.nc', 'lat');%y
9 r8 L: q, M9 }6 v' R! ax = ncread('rhum.mon.mean.nc', 'lon');%x
  x5 U8 h' g) X! v$ ^9 t/ L, Qt = ncread('rhum.mon.mean.nc', 'time');%t
5 q' Z0 M3 H2 z6 Vrhum = ncread('rhum.mon.mean.nc', 'rhum');

7 o0 g: i4 g* w6 Larray = zeros(x_len,y_len);
: t& f: ?& h  kfor i = 30:y_len
: Y# W& _( r# z. e/ H    for j = 1:t_len
        array(i,j)=rhum(x_index,i,z_index,j);
- c  d" e* ?3 ~/ W    end
end
4 ^! R2 M, |0 Dmesh(t,y,array)
  `$ L8 g4 V" n% I错误使用 mesh (line 79)
. \- f9 J: E8 J5 B数据维度必须一致。
想请问一下这个怎么改 ？

House

需要先用meshgrid根据横、纵坐标重新生成对应的矩阵
另外注意循环嵌套的顺序，是Z(jj,ii)而不是Z(ii,jj)
draw_x = 1:1/4:7;
5 b5 X3 H: ~* \draw_y = 2:1/3:6;
f = @(x,y) x.^(exp( -y.^2/50 ));
" Q( ?2 z* l2 Gfor ii = 1:1:numel(draw_x)
% @5 N* T6 l8 E: E4 Y  S! v    for jj = 1:1:numel(draw_y)
        Z(jj,ii) = f(draw_x(ii),draw_y(jj));
9 f$ _; R7 ?2 \7 I% N" _/ Q    end
5 R4 k: [" A9 Z6 v* }end
[draw_x,draw_y] = meshgrid(draw_x,draw_y);
, c- `# i- S9 g& A' imesh(draw_x,draw_y,Z);