想请问一下这个怎么改？

panphp · 发表于 2020-6-5 11:37

EDA365欢迎您登录！

您需要登录才可以下载或查看，没有帐号？注册

x

clear;
x_len=144;
y_len=60;
z_len=17;
t_len=854;
z_index=7;
x_index=120;

z = ncread('rhum.mon.mean.nc', 'level');%z
y = ncread('rhum.mon.mean.nc', 'lat');%y
x = ncread('rhum.mon.mean.nc', 'lon');%x
t = ncread('rhum.mon.mean.nc', 'time');%t
rhum = ncread('rhum.mon.mean.nc', 'rhum');

array = zeros(x_len,y_len);
for i = 30:y_len
for j = 1:t_len
array(i,j)=rhum(x_index,i,z_index,j);
end
end
mesh(t,y,array)
错误使用 mesh (line 79)
数据维度必须一致。
想请问一下这个怎么改？

House · 发表于 2020-6-5 13:28

需要先用meshgrid根据横、纵坐标重新生成对应的矩阵
另外注意循环嵌套的顺序，是Z(jj,ii)而不是Z(ii,jj)
draw_x = 1:1/4:7;
draw_y = 2:1/3:6;
f = @(x,y) x.^(exp( -y.^2/50 ));
for ii = 1:1:numel(draw_x)
for jj = 1:1:numel(draw_y)
Z(jj,ii) = f(draw_x(ii),draw_y(jj));
end
end
[draw_x,draw_y] = meshgrid(draw_x,draw_y);
mesh(draw_x,draw_y,Z);

帐号		自动登录	找回密码
密码			注册

想请问一下这个怎么改？

EDA365欢迎您登录！

浏览过的版块

推荐内容 /1

想请问一下这个怎么改 ？

EDA365欢迎您登录！

浏览过的版块

推荐内容 /1

想请问一下这个怎么改？