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

uiabluqp

d=0:0.1:20;
% N9 w: l3 r# s: r2 B7 Zl1=d.^3/41.67;
2 F  Y- f0 N: Il2=(d.^4/1.62).^(1/3);
D=2.5^(1/3);
4 u1 \+ L4 N! Z4 eL=10;
l=0:0,2:20;
plot(d,l1,d,l2,D,l,d,L);
xlabel('d(cm)');
- r. x2 k* T: ~7 c8 q) f& U0 kylabel('l(cm)');
0 D* [8 V3 k& ^6 o. t* y! z" ?8 m# r[a,b]=meshgrid(d,l);
/ q! s( }$ t$ y  Ez=0.00613*(d.^2);
Q=z.*l;contour(a,b,Q,15);


/ L0 E  R- z- h3 [+ w, u- L
; Q7 `% b8 b$ k0 a( P求助，运行不起来，划红线的是l，怕跟数字1混淆



% p2 J* ~( Q9 h
5 N1 t6 I( C7 f, F6 y
, R- M" A9 U5 v* _

xiaogegepcb

帮顶

haidaowang

eshgrid函数已经将原本是一维向量的d和l生成对应的二维的矩阵a和b了，那么用于绘制等高线图的Q就应当采用这些矩阵进行运算，即将原本使用d的替换为a，原本使用l的替换成b，注意乘除幂都应当使用点运算。
写成
3 A, [$ @5 v& @: r4 m

• 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 );
  / N1 y8 z3 p+ ^+ i3 M  z

  U$ b: h* b2 {; T, U; T
, F) R4 K! K4 r4 I8 k
# J* w! G0 Y* _% C6 H

uiabluqp

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

$ }9 f' x& F, Y  q* W) q太感谢了
7 G8 }  `4 Z6 V2 Q8 v. N: p