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

uiabluqp

d=0:0.1:20;
# p1 F, z9 c. `: z6 s0 k3 El1=d.^3/41.67;
l2=(d.^4/1.62).^(1/3);
* u' k, T; o1 eD=2.5^(1/3);
# D6 _+ Z$ [, I' P7 i! Q# FL=10;
l=0:0,2:20;
( c+ `' k8 t' j# m3 yplot(d,l1,d,l2,D,l,d,L);
8 D8 G# y% ^7 U5 R3 o2 kxlabel('d(cm)');
' v" ^' d/ F- _% t8 K; Z0 [/ aylabel('l(cm)');
/ d' ]) @" @& L! q6 `[a,b]=meshgrid(d,l);
5 Q1 A7 I" H4 U- |  A) Vz=0.00613*(d.^2);
Q=z.*l;contour(a,b,Q,15);
. J( i. V! S4 t" Q# I1 E& H4 p
. w- ?, a- E; V. p
5 y6 m. f7 A6 B  x  K4 Z+ K( a
) A3 |% |- I7 Y; D求助，运行不起来，划红线的是l，怕跟数字1混淆
) U, b  y; T( [6 e% W+ f
3 v1 w! x- @# h  v4 _) L+ v


  g  m& n6 w* A0 S

xiaogegepcb

帮顶

haidaowang

eshgrid函数已经将原本是一维向量的d和l生成对应的二维的矩阵a和b了，那么用于绘制等高线图的Q就应当采用这些矩阵进行运算，即将原本使用d的替换为a，原本使用l的替换成b，注意乘除幂都应当使用点运算。
+ l) t8 Y4 q. X( r, @" i) G) ]! a写成
! i$ t4 M# l: E

• 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 );
  9 z0 m( C7 j. Y

. l' r* P1 ^0 W! Y2 I

uiabluqp

haidaowang 发表于 2020-7-9 15:54
' `; l" j$ f+ W+ q5 Z( _, Neshgrid函数已经将原本是一维向量的d和l生成对应的二维的矩阵a和b了，那么用于绘制等高线图的Q就应当采用这 ...

/ x1 S1 J" H4 d+ v& o* ]太感谢了