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

uiabluqp

d=0:0.1:20;
% H5 a: L* R: `( W  [1 Ql1=d.^3/41.67;
5 j* Q& |% l0 e8 _2 S6 E$ a% |  J. El2=(d.^4/1.62).^(1/3);
3 E# x. z  N& K6 lD=2.5^(1/3);
: G. T1 |" M4 s  x& Q5 {L=10;
l=0:0,2:20;
plot(d,l1,d,l2,D,l,d,L);
# o7 q6 n" x' v4 d# q& qxlabel('d(cm)');
3 O9 `0 b7 t2 d: u% q4 e& K2 K+ ?ylabel('l(cm)');
[a,b]=meshgrid(d,l);
1 H1 V1 t( U3 C$ k: }z=0.00613*(d.^2);
, E( T. U' N7 g3 H1 iQ=z.*l;contour(a,b,Q,15);
1 P6 Y4 V. e# h( d: N


$ ?5 y- y$ z/ U3 }2 x- K求助，运行不起来，划红线的是l，怕跟数字1混淆
$ k8 q3 X3 a7 p( [
8 v5 w4 S6 I* W" b8 K/ R2 l6 V
8 t  M' K5 v( ~% d
3 f5 m8 y- P( C: e


6 o) L4 W; K% }; n, O9 v

xiaogegepcb

帮顶

haidaowang

eshgrid函数已经将原本是一维向量的d和l生成对应的二维的矩阵a和b了，那么用于绘制等高线图的Q就应当采用这些矩阵进行运算，即将原本使用d的替换为a，原本使用l的替换成b，注意乘除幂都应当使用点运算。
1 ?1 D3 h  ?* D) S. B+ 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 );
    S* `4 c# e, R6 e5 D, [0 h; y$ {" E0 @

9 @' h5 A, V: i# R5 H
% F# x$ J5 h* T$ A  u, U6 {$ R

uiabluqp

haidaowang 发表于 2020-7-9 15:54
7 D$ D' T! B0 Q/ d4 t: {6 i# T, Yeshgrid函数已经将原本是一维向量的d和l生成对应的二维的矩阵a和b了，那么用于绘制等高线图的Q就应当采用这 ...

0 z- E5 A9 z7 ^0 ~) m4 d太感谢了
2 ?  H7 `6 f5 u; A. x% @3 a