找回密码
 注册
关于网站域名变更的通知
查看: 438|回复: 1
打印 上一主题 下一主题

contour函数报错怎么办?

[复制链接]

该用户从未签到

跳转到指定楼层
1#
发表于 2020-6-2 16:07 | 只看该作者 回帖奖励 |倒序浏览 |阅读模式

EDA365欢迎您登录!

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

x
clear$ N) j' I$ r7 l' @; t4 U
ncinfo('rhum.mon.mean.nc');! e1 Z8 U  o9 ]
y_len=60;% ]( ]0 Q  c0 h6 D/ V1 u7 r
t_len=854;( x% `! ^- x- [: ]& d$ }7 a
z_index=74 {- s3 q! o7 C3 E6 {' s! e0 q" B5 o: K
x_index=120
, S% F. b, i( p+ Vx = ncread('rhum.mon.mean.nc','lon');%x
; N# g* z! C$ iz = ncread('rhum.mon.mean.nc', 'level');%z) X( q6 o+ ^/ E. A) |1 @' w
y = ncread('rhum.mon.mean.nc', 'lat');%y: c+ n) a( a2 L9 Z$ s& K
t = ncread('rhum.mon.mean.nc', 'time');%t
7 G& K! `8 }3 w- Q  d) Qrhum = ncread('rhum.mon.mean.nc', 'rhum');
6 _% C* o$ d" e% Q% @4 P0 y" @figure; hold on;$ i, R9 L( }; a9 k: ], T
draw_y = 30:1:60;( F3 `# i( Z$ ]: J, d
draw_t = 1:1:854;8 N/ l& I- t4 ?; M. E$ E
f = rhum(x_index,draw_y,z_index,draw_t): M$ t9 \+ m- M  r" `6 H  l
for ii = 1:1:numel(draw_y)( b; `% p# z+ N
    for jj = 1:1:numel(draw_t)
  X* p2 v! |: Z/ L        Z(jj,z_index,ii,x_index) = f(x_index,draw_y(ii),z_index,draw_t(jj));0 ?# F) g% a* J1 v$ Z3 A
    end1 v, h  b# {* i
end+ v# W2 b) z* e
[draw_y,draw_t] = meshgrid(draw_y,draw_t);6 U4 K7 h! {, E/ E- \
contour (draw_y,draw_t,Z);
3 f1 u, d' w* Z5 |. @. u* J) c) U3 ]4 ~* O, ?
错误显示 试图访问 f(120,30,7,1);由于 size(f)=[1,31,1,854],索引超出范围。
6 f- b+ R4 O) g4 j改了循环里面之后会显示 试图访问 f(1,32,1,1);由于 size(f)=[1,31,1,854],索引超出范围。
* I/ ^+ H1 P/ y* G# \请教大神,非常感谢!
/ W; L, O1 s7 u; w

该用户从未签到

2#
发表于 2020-6-2 16:16 | 只看该作者
可能参数超出了范围
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

推荐内容上一条 /1 下一条

EDA365公众号

关于我们|手机版|EDA365电子论坛网 ( 粤ICP备18020198号-1 )

GMT+8, 2025-11-5 15:25 , Processed in 0.156250 second(s), 24 queries , Gzip On.

深圳市墨知创新科技有限公司

地址:深圳市南山区科技生态园2栋A座805 电话:19926409050

快速回复 返回顶部 返回列表