想请问一下这个怎么改 ？

panphp

- S' Z' q4 a! w8 {  ?6 R8 t
clear;
x_len=144;
4 g6 j: Q4 v! Y! F1 @, y% r- n3 v  vy_len=60;
z_len=17;
t_len=854;
z_index=7;
x_index=120;
( e/ ~  N+ n0 n1 t: K" j" Z0 l
z = ncread('rhum.mon.mean.nc', 'level');%z
! V$ N5 F! }* P# R7 _* ^, Ey = ncread('rhum.mon.mean.nc', 'lat');%y
; r& I" U9 a( Cx = ncread('rhum.mon.mean.nc', 'lon');%x
- q! r3 @- Y1 F  W5 N( b% L; [) Yt = ncread('rhum.mon.mean.nc', 'time');%t
rhum = ncread('rhum.mon.mean.nc', 'rhum');
. T/ b0 v% ^% e
array = zeros(x_len,y_len);
for i = 30:y_len
    for j = 1:t_len
. X; ^/ R& v4 h" i0 W- p0 w' h; [        array(i,j)=rhum(x_index,i,z_index,j);
& ]8 u) n0 C/ ~$ u0 O1 |    end
- w, r3 l* H6 h- Zend
. d' B' Z6 o6 P5 d) {0 kmesh(t,y,array)
错误使用 mesh (line 79)
数据维度必须一致。
) X' H! G8 w( D. ~想请问一下这个怎么改 ？

House

需要先用meshgrid根据横、纵坐标重新生成对应的矩阵
另外注意循环嵌套的顺序，是Z(jj,ii)而不是Z(ii,jj)
/ L( ^2 `! S/ p/ m6 X& d' a/ Idraw_x = 1:1/4:7;
draw_y = 2:1/3:6;
9 r% x1 H5 m' df = @(x,y) x.^(exp( -y.^2/50 ));
2 l! Z4 o6 }/ a! h. _9 q2 Xfor ii = 1:1:numel(draw_x)
* _& X) b+ C" `2 x, g" v9 J    for jj = 1:1:numel(draw_y)
        Z(jj,ii) = f(draw_x(ii),draw_y(jj));
    end
  w+ s1 c, i2 s' jend
[draw_x,draw_y] = meshgrid(draw_x,draw_y);
+ d8 i- D2 [3 I( S7 O" E4 Y: }mesh(draw_x,draw_y,Z);