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

uiabluqp

d=0:0.1:20;
/ e( v/ m: B2 @" E% ^. Ul1=d.^3/41.67;
/ N, k; p: w3 e2 b/ {: S# K; Il2=(d.^4/1.62).^(1/3);
7 A! s" }* k1 t6 ~D=2.5^(1/3);
3 Q7 ?4 T9 ]* \) h% ?- G- ]- g- p8 kL=10;
l=0:0,2:20;
plot(d,l1,d,l2,D,l,d,L);
xlabel('d(cm)');
ylabel('l(cm)');
% p, r  L3 o* Y0 w1 j* g[a,b]=meshgrid(d,l);
z=0.00613*(d.^2);
Q=z.*l;contour(a,b,Q,15);
4 M/ g( b% R9 y! `' c( d- O; V" W
4 X- I( Z. i; z9 b3 ?

$ o, b# ~5 B# d: X0 {/ Q求助，运行不起来，划红线的是l，怕跟数字1混淆


6 x- Y. T6 P1 \- v

& X# t, f& a& D- Y7 J5 j

' p+ g  |- A( `, G- H; Y$ ?6 h$ N# i2 P
( f3 [2 b) q) W. Z/ L7 {1 ]) V$ b

uiabluqp

haidaowang 发表于 2020-7-9 15:54
. S. A- S% W/ y; Q( L* ]: s7 Y0 Ueshgrid函数已经将原本是一维向量的d和l生成对应的二维的矩阵a和b了，那么用于绘制等高线图的Q就应当采用这 ...

' F9 m; Y" t. z4 k太感谢了

haidaowang

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

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

- l5 T* O5 z# l+ c" y5 X! q
, ]: N) P" N6 t
+ ?) B7 e4 q9 m

xiaogegepcb

帮顶