|
|
EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
我通过for循环求解一个带有变系数微分方程后，plot出各个方程的曲线，想将各个曲线对应legend打印出来，但是始终无法成功，不知道哪里出了问题，代码如下：( }; ~, e$ n: B. J
c = (4-19.2/2)*pi/180;
0 \ |- N4 h! G& G% vR1 = 38;+ v5 q# [! o1 A2 ?5 p
t1 = 5*pi/180;2 o( u3 C3 M2 f1 J2 N- u
t2 = pi/2;
/ L4 U) g" o2 V4 X y1 tA=R1.^2 * sin(t1) * tan(t1) * tan(t2/2);" U$ P$ p( I7 Y `9 M7 O
. @3 T1 }( \! c6 b J E2 TS={'aaa'};% 定义一个元胞数组,并随意存储一个字符串8 K7 q" S8 P4 f3 k$ g+ q
for KB = [0, 0.2, 0.4, 0.6, 0.8, 1]2 T. e0 d% \5 X0 l7 z4 A$ Y
a1 = (2*1500*9*0.8*10^6/20)/(1500*9*13888.2*(2*1.313+KB+1)/60);4 B7 Z$ a- G- c6 U6 k9 U
b1 = 1500/(2*1.313+KB+1);/ H- l6 O& W, D) E2 I
y0 = 0.1;
& X& ~8 R1 @( ]+ r% j2 S8 C, B yp0 = b1*sin(c);) g/ G% M3 r+ S& F8 p; V' l
[x,y] = ode15i(@(x,y,yp)yp - a1*(A*x.^2)*sqrt(abs(35-y))*sign(35-y) - b1*sin(c+x) , [0, 19.2*pi/180], y0, yp0);
# e+ P' A. d/ d& i. \4 E h = plot(x,y);' b1 I2 ^2 Y7 v2 J- `" U
kb = sprintf('KB=%.1f',KB);%将KB的转化为字符串存在一个变量中
% V f# V& @' d* e2 l S(i) = cellstr(kb);%将变量kb存储在S元胞数组中$ h9 p0 i4 i& j) |5 z
hold on% I {. p% o1 p" N! C+ I, c% d# U
end
: H* S' h1 c( |legend(S);
5 s' Q" h9 g4 f, C; c5 @7 Lgrid
1 U' _2 Z2 M/ i1 ~# l
' ~% h! e- c; P此外，能否在for循环中使每个plot出来的曲线的线型不一样，比如说KB=0时是虚线，KB=1时是实线等等
4 N& u7 a& b/ u! D1 y请教大神，非常感谢！; Y% Z. P$ N+ W; p; C5 q, E
出现的问题如下图6 T' u6 L& `0 i) m( o2 u
9 U; W4 h( U# F8 T# F- ?7 S- y/ x" V7 [# `$ L5 s
& q5 E+ K7 ` g0 W/ } |
|
/1 
发表于 2020-5-27 16:24
收藏
淘帖
支持！
反对！