请问Z怎样设置才是矩阵呢？

uiabluqp

d=0:0.1:20;
7 i& H7 _: O0 n3 z/ K  I' _l1=d.^3/41.67;
4 L- x7 \( b" ?% ?( ul2=(d.^4/1.62).^(1/3);
# g9 v/ ?+ Q$ RD=2.5^(1/3);
; \, ], k* |$ F. D+ Q' iL=10;
$ n$ d& _% V% ]7 |) cl=0:0,2:20;
plot(d,l1,d,l2,D,l,d,L);
$ E' _/ O: c3 c2 C" e9 jxlabel('d(cm)');
ylabel('l(cm)');
& o% A) H" j! \* p$ s# K* m- a[a,b]=meshgrid(d,l);
1 ]6 d/ {( v0 ?. n: ^' qz=0.00613*(d.^2);
' t* ?% k' I: e0 OQ=z.*l;contour(a,b,Q,15);
5 @& G3 \) z: j8 u3 K- r
- b0 Z& Y* U' c* J+ ~

0 i0 ]2 ^/ b$ G) V; B& P求助，运行不起来，划红线的是l，怕跟数字1混淆
' c& [# i- H* C  l% t
( Z$ W, f4 b7 I# M. ^1 R


' R' L- v1 G; j9 S6 V


8 w2 ]( z. s, E: h9 ?7 K

xiaogegepcb

帮顶

haidaowang

eshgrid函数已经将原本是一维向量的d和l生成对应的二维的矩阵a和b了，那么用于绘制等高线图的Q就应当采用这些矩阵进行运算，即将原本使用d的替换为a，原本使用l的替换成b，注意乘除幂都应当使用点运算。
8 P7 J' }7 f' v& J7 O写成

• d = 0 : 0.1 : 20;
• l1 = d.^3/41.67;
• l2 = (d.^4/1.62).^(1/3);
• D = 2.5^(1/3);
• L = 10;
• l = 0 : 0.2 : 20;
• [ a, b ] = meshgrid( d, l );
• z = 0.00613*( a.^2 );
• Q = z .* b;
• subplot( 121 )
• plot( d, l1, 'ko-', d, l2, 'bs-.', D, l, 'rp--', d, L, 'ch-' );
• xlabel('d(cm)');
• ylabel('l(cm)');
• subplot( 122 )
• contour( a, b, Q, 15 );
  2 E5 X9 S& }( Z

0 b& `* j8 X7 V) a- _/ N( R) i* l

uiabluqp

haidaowang 发表于 2020-7-9 15:54
eshgrid函数已经将原本是一维向量的d和l生成对应的二维的矩阵a和b了，那么用于绘制等高线图的Q就应当采用这 ...

' ]& t0 s  D; r& ?( Q太感谢了
9 V; m" u% N# ^" Y- U