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

uiabluqp

d=0:0.1:20;
% ?2 L% V6 o! V! `. pl1=d.^3/41.67;
l2=(d.^4/1.62).^(1/3);
: o% }5 Q+ \/ n/ D" @" B$ Z% dD=2.5^(1/3);
" r9 Y2 K7 e# n0 m" FL=10;
l=0:0,2:20;
plot(d,l1,d,l2,D,l,d,L);
xlabel('d(cm)');
ylabel('l(cm)');
2 A' A# ~0 M6 t[a,b]=meshgrid(d,l);
$ C, z3 y& D* b" Az=0.00613*(d.^2);
Q=z.*l;contour(a,b,Q,15);

& _: m: r' |7 m4 n
# g; ]# }4 Y: ^$ M
求助，运行不起来，划红线的是l，怕跟数字1混淆
: I3 i* L# Q% a& S+ R" \


2 E& G5 f8 V* D! X, }- s
  J6 G9 Z8 D# m2 i3 ?! b6 q

xiaogegepcb

帮顶

haidaowang

eshgrid函数已经将原本是一维向量的d和l生成对应的二维的矩阵a和b了，那么用于绘制等高线图的Q就应当采用这些矩阵进行运算，即将原本使用d的替换为a，原本使用l的替换成b，注意乘除幂都应当使用点运算。
- [$ i$ Z8 m4 n4 H6 w写成

• 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 );
  / A5 k0 W) v+ k( n( t

  {0 d# x# T. m# B. g
0 C* h$ w# h) Y1 _& K

uiabluqp

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

+ x9 |. X$ b- c0 x( I9 X6 W0 q! T/ m太感谢了