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

想请问一下这个怎么改 ?

[复制链接]

该用户从未签到

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

EDA365欢迎您登录!

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

x
: c" W% J; L/ J; C: |+ Z

" g5 s- O! J0 y0 j0 m& E" ^clear;
$ w& O1 n' g( p- [* I; rx_len=144;
) a+ Y  ~: ^0 j' p' d0 Dy_len=60;
' j( u0 @  `6 |! gz_len=17;
6 `% S" ^/ K" D" j4 ?0 k: Ot_len=854;% g3 x$ c" A7 X1 w
z_index=7;* b, u9 p5 j; K" a
x_index=120;7 d, a- n4 n  f$ c8 W  D: q2 T+ U# u

  U! U  B1 ~, w6 e* @9 h9 `% Y3 {z = ncread('rhum.mon.mean.nc', 'level');%z
; H- B( t9 p/ P0 }$ z3 E2 y4 x3 c- m6 By = ncread('rhum.mon.mean.nc', 'lat');%y4 I. [$ o3 a( v+ L( U' }
x = ncread('rhum.mon.mean.nc', 'lon');%x
! j, @( Q, q3 _! a) it = ncread('rhum.mon.mean.nc', 'time');%t
- G+ I" ]- r  h# ~$ vrhum = ncread('rhum.mon.mean.nc', 'rhum');& d7 z5 c- m2 B" N) S5 H8 n

) t6 W6 x6 m' _8 iarray = zeros(x_len,y_len);
" {5 J3 t1 Q. vfor i = 30:y_len( L) Y. a1 e- u  q9 H
    for j = 1:t_len( _8 ]( s1 q% n3 w6 q+ U  R
        array(i,j)=rhum(x_index,i,z_index,j);* q* K$ [& _* X9 Z
    end
7 ?+ ^) v+ P& z% l5 i8 ]# ^end
5 i4 \$ L: L& U7 D; L" q; |mesh(t,y,array)
* C: Q4 W6 L' H9 b  L错误使用 mesh (line 79): Z) T' `) o) \/ n* P0 ^
数据维度必须一致。
) b1 _9 y' K. k想请问一下这个怎么改 ?

该用户从未签到

2#
发表于 2020-6-5 13:28 | 只看该作者
需要先用meshgrid根据横、纵坐标重新生成对应的矩阵
% b9 {$ [% G5 s3 h: `/ J另外注意循环嵌套的顺序,是Z(jj,ii)而不是Z(ii,jj)$ u) _. p! l* X
draw_x = 1:1/4:7;
* d; `% p: T, i. E4 s! c% `4 jdraw_y = 2:1/3:6;4 O& q7 S5 I5 d+ Q& g! r4 f2 N
f = @(x,y) x.^(exp( -y.^2/50 ));
7 c8 K) T2 y" C, H; x. sfor ii = 1:1:numel(draw_x); g! V1 m  {; D6 v/ z/ o
    for jj = 1:1:numel(draw_y)
1 Q" c- v( z1 B/ T        Z(jj,ii) = f(draw_x(ii),draw_y(jj));& `0 G" p' R) {  q2 B1 T
    end
! O6 P# N, ^) f7 ~0 s! ~end
9 [* Z, |  {% a& ]2 p[draw_x,draw_y] = meshgrid(draw_x,draw_y);
6 R  [5 J* A1 R; d# o% r) }. ~! Vmesh(draw_x,draw_y,Z);
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

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

EDA365公众号

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

GMT+8, 2025-6-24 19:36 , Processed in 0.078125 second(s), 23 queries , Gzip On.

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

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

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