想请问一下这个怎么改 ？

panphp

0 b6 ]# P* E; e3 bclear;
: J. ?8 P$ B1 D8 kx_len=144;
y_len=60;
1 P7 Q6 v9 k8 x. v, Y1 lz_len=17;
t_len=854;
: i) k( z3 g' m2 j8 j; Ez_index=7;
x_index=120;
& ]& e, F2 o3 x* i0 I
z = ncread('rhum.mon.mean.nc', 'level');%z
* K0 W/ Y# {+ R  qy = ncread('rhum.mon.mean.nc', 'lat');%y
! [+ Y1 w+ f8 o$ l1 |. Ix = ncread('rhum.mon.mean.nc', 'lon');%x
t = ncread('rhum.mon.mean.nc', 'time');%t
rhum = ncread('rhum.mon.mean.nc', 'rhum');
4 G) J) \- U2 w" u# R
array = zeros(x_len,y_len);
% P' _8 ^2 |# e% ]' C' Qfor i = 30:y_len
    for j = 1:t_len
        array(i,j)=rhum(x_index,i,z_index,j);
( F# h9 Z5 s1 I3 H" N3 c    end
- Q0 M, T3 g9 W8 K8 ~. a. Jend
mesh(t,y,array)
% W+ f: y- Y) g8 W+ G+ ]7 e错误使用 mesh (line 79)
8 |. z. F+ G3 P2 T$ D/ G数据维度必须一致。
) w: j4 Z0 a- j) C; Q4 V) T5 M想请问一下这个怎么改 ？

House

需要先用meshgrid根据横、纵坐标重新生成对应的矩阵
另外注意循环嵌套的顺序，是Z(jj,ii)而不是Z(ii,jj)
7 c: @( K: j; Y- \- c; g$ \: P5 O# pdraw_x = 1:1/4:7;
3 Q6 y2 @0 B1 Z$ _1 Idraw_y = 2:1/3:6;
f = @(x,y) x.^(exp( -y.^2/50 ));
for ii = 1:1:numel(draw_x)
' C6 n, z  N) ^% Y    for jj = 1:1:numel(draw_y)
" {2 ^' l8 D& C! }" L2 C# J        Z(jj,ii) = f(draw_x(ii),draw_y(jj));
    end
; R" N' {/ J" ~end
[draw_x,draw_y] = meshgrid(draw_x,draw_y);
1 Z6 u6 f7 G& u4 A" x9 Lmesh(draw_x,draw_y,Z);