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

想请问一下这个怎么改 ?

[复制链接]

该用户从未签到

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

EDA365欢迎您登录!

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

x

9 T6 |+ u  ^3 F# [) ~
- g( ?; ^0 ^: d3 Jclear;
2 B+ \; W+ ]$ v' u4 vx_len=144;1 c5 o% p. p# e3 x# J3 [
y_len=60;8 F3 R+ h) L' e( t! U/ R3 h2 k
z_len=17;) p& x: j2 Z5 G! ~
t_len=854;% I9 Y7 b- K* K" L
z_index=7;& C9 ^* ~7 c3 c  G, s# u
x_index=120;
5 O9 _; o# z9 _& e( i! f/ z; Q/ q4 b8 E1 I
z = ncread('rhum.mon.mean.nc', 'level');%z) k1 e6 ?' l$ a6 X6 s# J8 {
y = ncread('rhum.mon.mean.nc', 'lat');%y1 y. |! ?( ~$ E5 B5 {6 N- h
x = ncread('rhum.mon.mean.nc', 'lon');%x# Q6 X3 [6 x8 A% D& p
t = ncread('rhum.mon.mean.nc', 'time');%t3 D$ m4 p3 a0 X( J9 K( K
rhum = ncread('rhum.mon.mean.nc', 'rhum');3 [- \  n% H5 v  D5 C2 C4 t. i+ j

, _, ~, z: S2 g0 F$ i9 A$ C, ]array = zeros(x_len,y_len);
4 O, A4 N" T/ o: @6 E8 X: k+ Jfor i = 30:y_len
1 F# C0 [4 d  t/ C+ u    for j = 1:t_len
. e; o4 n# \: ]+ r5 }9 m5 o0 J        array(i,j)=rhum(x_index,i,z_index,j);
  ~. S% y, z, X5 h& f% o    end. _3 J8 W0 w) d8 V- g
end
. c+ \# W2 W  R/ H0 A+ U. z! ]2 ~4 S5 Xmesh(t,y,array)
" \) l+ M5 j& Z2 D5 d错误使用 mesh (line 79)
1 a: Z* Y2 i2 a% Q7 t( z数据维度必须一致。; G( _) L! u5 N% \( F, |, N
想请问一下这个怎么改 ?

该用户从未签到

2#
发表于 2020-6-5 13:28 | 只看该作者
需要先用meshgrid根据横、纵坐标重新生成对应的矩阵) d8 y" v# O( I  L- C3 }0 k* ]
另外注意循环嵌套的顺序,是Z(jj,ii)而不是Z(ii,jj)) U$ ^. l! q1 j6 @" K% a
draw_x = 1:1/4:7;
, B: |& [* B6 D; }' m4 I. Odraw_y = 2:1/3:6;
) B8 {1 h, n6 m* p) |  `f = @(x,y) x.^(exp( -y.^2/50 ));
! q8 e( c6 |, [) r) `! U) d3 Y: zfor ii = 1:1:numel(draw_x)
$ p& l1 }/ h' ?2 Y    for jj = 1:1:numel(draw_y)+ w$ m+ U4 \0 J: V+ v7 f# ]; Z' j
        Z(jj,ii) = f(draw_x(ii),draw_y(jj));) T' T+ m* F% z6 M& Q8 k
    end
0 E' n" j& r% n0 N( k6 y9 pend
7 v; f$ J2 L/ u; e4 ][draw_x,draw_y] = meshgrid(draw_x,draw_y);
- p  d+ x& p$ E8 h+ imesh(draw_x,draw_y,Z);
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

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

EDA365公众号

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

GMT+8, 2025-11-6 00:30 , Processed in 0.140625 second(s), 24 queries , Gzip On.

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

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

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