EDA365电子论坛网
标题:
contour函数报错怎么办？
[打印本页]
作者:
panphp
时间:
2020-6-2 16:07
标题:
contour函数报错怎么办？
clear
. c2 M( [5 ^! J
ncinfo('rhum.mon.mean.nc');
" ]. |# y; r; J+ Z! q2 z
y_len=60;
- U% g/ b& K/ y5 a! d6 O
t_len=854;
4 e3 o4 C# w- C' {
z_index=7
+ Y j# b, A+ F+ r$ |: q
x_index=120
' j0 b7 E4 a; {1 r
x = ncread('rhum.mon.mean.nc','lon');%x
$ m. [% |4 J& c
z = ncread('rhum.mon.mean.nc', 'level');%z
9 R I. S1 v, L9 T( @' [8 }
y = ncread('rhum.mon.mean.nc', 'lat');%y
: J( m3 P" z* N
t = ncread('rhum.mon.mean.nc', 'time');%t
* A2 U' Y* g, ^
rhum = ncread('rhum.mon.mean.nc', 'rhum');
$ h/ [1 Z) e8 q& m
figure; hold on;
- _" k" O+ L1 G5 |2 J) W9 Y- ]5 C- m) t
draw_y = 30:1:60;
7 ^ p3 a1 r2 Y1 r, S; y& t. Z$ W+ `8 w
draw_t = 1:1:854;
: X4 M& [( I9 ^; W7 e% Q5 P
f = rhum(x_index,draw_y,z_index,draw_t)
" T: P/ o- l) \( P8 E
for ii = 1:1:numel(draw_y)
& H5 M) A7 f: T% H$ m
for jj = 1:1:numel(draw_t)
9 _6 { Z% v# y( d
Z(jj,z_index,ii,x_index) = f(x_index,draw_y(ii),z_index,draw_t(jj));
6 N, a. R7 l4 ?
end
. W& {$ ]. N" E% W) k, s
end
3 U1 z& M& Q4 d1 ]! _
[draw_y,draw_t] = meshgrid(draw_y,draw_t);
* ?/ u, K `0 ]
contour (draw_y,draw_t,Z);
1 ?/ g% z X8 |# a& d" w7 z! H
0 A* W. Y- ^ D- S) U
错误显示
试图访问 f(120,30,7,1)；由于 size(f)=[1,31,1,854]，索引超出范围。
# u" F# G& _, R/ f( G
改了循环里面之后会显示 试图访问 f(1,32,1,1)；由于 size(f)=[1,31,1,854]，索引超出范围。
6 T- _1 M7 Y% ~; c: p+ A
请教大神，非常感谢！
: I5 \. U9 M6 o2 r! ^9 x3 w
作者:
Sky_wm
时间:
2020-6-2 16:16
可能参数超出了范围
欢迎光临 EDA365电子论坛网 (https://bbs.eda365.com/)
Powered by Discuz! X3.2