想请问一下这个怎么改 ？

panphp

/ g! g& s- O7 e2 u
; f. U5 X: X" T' dclear;
x_len=144;
$ \7 v2 Y% p' Q0 D5 \% |( ?6 [( F$ Ry_len=60;
3 r- V; d0 G* F) `2 J2 Uz_len=17;
2 p/ E# c& P* J  _  T8 H5 }$ ht_len=854;
z_index=7;
x_index=120;
# T% a) w0 h, }" }5 {
z = ncread('rhum.mon.mean.nc', 'level');%z
4 W- s. M8 m/ _1 B% T/ f1 ey = ncread('rhum.mon.mean.nc', 'lat');%y
x = ncread('rhum.mon.mean.nc', 'lon');%x
4 Y# T  |* a' H: m, `# A/ Pt = ncread('rhum.mon.mean.nc', 'time');%t
rhum = ncread('rhum.mon.mean.nc', 'rhum');

array = zeros(x_len,y_len);
$ u$ g! g5 `& \; ~- n8 hfor i = 30:y_len
( K6 \! X+ u& V3 l) ~+ [    for j = 1:t_len
        array(i,j)=rhum(x_index,i,z_index,j);
4 B; E$ {$ I! @$ x; L    end
end
mesh(t,y,array)
5 C8 }. y9 B2 ~* n- E错误使用 mesh (line 79)
数据维度必须一致。
想请问一下这个怎么改 ？

House

需要先用meshgrid根据横、纵坐标重新生成对应的矩阵
另外注意循环嵌套的顺序，是Z(jj,ii)而不是Z(ii,jj)
draw_x = 1:1/4:7;
$ m" r1 y5 Y/ @draw_y = 2:1/3:6;
f = @(x,y) x.^(exp( -y.^2/50 ));
- ~9 y; j- a8 b+ Afor ii = 1:1:numel(draw_x)
' G. d: G- a1 S! r    for jj = 1:1:numel(draw_y)
' v3 V8 j  n9 t1 c: k. i        Z(jj,ii) = f(draw_x(ii),draw_y(jj));
- |. C( c# {+ S0 M    end
- V. C* u9 H* _end
) M3 G# W. b4 U: M3 N[draw_x,draw_y] = meshgrid(draw_x,draw_y);
mesh(draw_x,draw_y,Z);