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

uiabluqp

d=0:0.1:20;
, j  F! _" ~: V3 j5 x% J% M+ o( K& wl1=d.^3/41.67;
9 m) O) y" r, }7 k$ r# ul2=(d.^4/1.62).^(1/3);
: f) S) G0 t7 d: v' u  o4 h1 D6 \D=2.5^(1/3);
' G' T$ I) e5 k. G0 XL=10;
7 N: x- Z: X! z( Y% ^# v" {: gl=0:0,2:20;
plot(d,l1,d,l2,D,l,d,L);
xlabel('d(cm)');
+ r; T1 n0 m8 h" X8 V! A: |- dylabel('l(cm)');
& r% m/ I2 x- }  C7 q: L, ?/ J[a,b]=meshgrid(d,l);
z=0.00613*(d.^2);
Q=z.*l;contour(a,b,Q,15);
+ Z7 W, a% O) J  Y, H$ A
4 S4 E! u8 H# ?$ c5 E% C6 E' y
2 P1 i' A! w# y
求助，运行不起来，划红线的是l，怕跟数字1混淆
1 g! r1 ~4 |" V4 z( N5 C, t
1 S( C9 t. H. P7 s

7 A6 s. g. ~) i. t, n

. g! Z8 [8 z# H* u' b+ \2 {
1 b' p( q3 w1 O0 G

xiaogegepcb

帮顶

haidaowang

eshgrid函数已经将原本是一维向量的d和l生成对应的二维的矩阵a和b了，那么用于绘制等高线图的Q就应当采用这些矩阵进行运算，即将原本使用d的替换为a，原本使用l的替换成b，注意乘除幂都应当使用点运算。
写成
4 C$ R$ r# b* b% ]6 g: c" K

• 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 );
  

uiabluqp

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

太感谢了
, _9 [; Y6 M8 a- O