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

想请问一下这个怎么改 ?

[复制链接]

该用户从未签到

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

EDA365欢迎您登录!

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

x
; ^: Y. Z! j8 A& Y0 v: @2 {. V1 j5 z6 l
* \3 Q' c+ E! K, s: s+ c% ?
clear;; w6 X  I% o( H8 S9 ]
x_len=144;; ?1 y* q- m" z; m
y_len=60;% y! r$ g) K2 \+ k- b  O3 P
z_len=17;0 A( ^0 j( n2 P. O" ]& Q
t_len=854;; _: r# U" |" m* o6 C/ k
z_index=7;
+ _, W! L; B; t( q+ ?$ {* Q! jx_index=120;; g1 ~# l9 @8 ^2 i6 L

3 K. ]5 D$ X  V' h' `, m! mz = ncread('rhum.mon.mean.nc', 'level');%z
0 q% M: y* h. A$ H) v! Qy = ncread('rhum.mon.mean.nc', 'lat');%y0 _) z7 k! l0 {* M: M6 V
x = ncread('rhum.mon.mean.nc', 'lon');%x. h" H& c% O) Y9 O0 Y
t = ncread('rhum.mon.mean.nc', 'time');%t
4 P0 @- s5 ?4 a4 F1 hrhum = ncread('rhum.mon.mean.nc', 'rhum');8 Q0 f7 p% H5 M0 y3 b( T) c
$ r. N& k6 d8 f: J! r8 T
array = zeros(x_len,y_len);9 l) N3 D# Y: j5 K
for i = 30:y_len& e" G% R: n. l8 I9 d# c0 [' w
    for j = 1:t_len
5 x4 `  Y5 b9 e9 G0 [        array(i,j)=rhum(x_index,i,z_index,j);
2 `. Z/ K+ M/ B! U8 }    end
7 d8 v% |' c0 {' o% i! Q: Gend
! \1 j: w: D: @% F/ N4 q9 Tmesh(t,y,array)
4 F! \$ W1 d, W0 [5 Q5 H$ o" `错误使用 mesh (line 79)5 X% U2 c2 l- Q. U
数据维度必须一致。/ A+ h7 ?! k. A0 O
想请问一下这个怎么改 ?

该用户从未签到

2#
发表于 2020-6-5 13:28 | 只看该作者
需要先用meshgrid根据横、纵坐标重新生成对应的矩阵7 I" N2 |5 k, q6 P) j6 `$ [
另外注意循环嵌套的顺序,是Z(jj,ii)而不是Z(ii,jj)- o* X( K  ~/ U+ u' O
draw_x = 1:1/4:7;0 |* f: A( l* c. I! |4 g! f1 j
draw_y = 2:1/3:6;) M, j0 y& _- f4 \' I/ L
f = @(x,y) x.^(exp( -y.^2/50 ));
. M' F$ I- o  z; N+ \for ii = 1:1:numel(draw_x)% C3 \8 z3 E: ?
    for jj = 1:1:numel(draw_y)
5 h  j0 y+ W1 }. S9 H3 w( s        Z(jj,ii) = f(draw_x(ii),draw_y(jj));
$ E/ X8 `5 c8 P: G2 s4 e1 w& v    end
$ }1 W- l% T$ e( b' O. i+ s- \& Uend) c9 n: E6 [3 u( Z
[draw_x,draw_y] = meshgrid(draw_x,draw_y);
( v7 N' O. J) l: v: O2 x; a- nmesh(draw_x,draw_y,Z);
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

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

EDA365公众号

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

GMT+8, 2025-11-24 14:36 , Processed in 0.156250 second(s), 23 queries , Gzip On.

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

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

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