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

uiabluqp

d=0:0.1:20;
6 ?) Y! G! a% pl1=d.^3/41.67;
l2=(d.^4/1.62).^(1/3);
D=2.5^(1/3);
L=10;
l=0:0,2:20;
6 s+ p) i: `# M" {, K$ I& Vplot(d,l1,d,l2,D,l,d,L);
) V' o! y: O7 r1 m: Q7 bxlabel('d(cm)');
- Y0 e5 n6 T1 _( _9 Yylabel('l(cm)');
[a,b]=meshgrid(d,l);
' G0 V/ v$ z/ }' gz=0.00613*(d.^2);
( B# P1 j' u5 ~, X$ |Q=z.*l;contour(a,b,Q,15);


/ x9 E1 r; T4 i! i3 P+ @2 r! R5 A
" a6 s2 d( Z% S% _5 M% Z求助，运行不起来，划红线的是l，怕跟数字1混淆
$ x* U6 k3 q, ?$ k7 E
* t5 J+ a( Z$ s, t5 X% h4 a
4 e6 ]9 v0 o$ v4 i. d


( Q+ {$ G7 ]( g1 }5 D( G9 t5 g/ h1 n

xiaogegepcb

帮顶

haidaowang

eshgrid函数已经将原本是一维向量的d和l生成对应的二维的矩阵a和b了，那么用于绘制等高线图的Q就应当采用这些矩阵进行运算，即将原本使用d的替换为a，原本使用l的替换成b，注意乘除幂都应当使用点运算。
% S) ?7 o" z! k写成
% M- g6 J2 @: n; o: K5 U& P2 s

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

8 e& o# a1 F* K

uiabluqp

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

太感谢了