|
|
EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
我通过for循环求解一个带有变系数微分方程后，plot出各个方程的曲线，想将各个曲线对应legend打印出来，但是始终无法成功，不知道哪里出了问题，代码如下：
1 O, d8 c# }. F6 ]% z6 o, V [c = (4-19.2/2)*pi/180;8 ?2 K; X/ x8 ]- J4 Y5 m6 _- z( J
R1 = 38;1 w6 w$ Y, F0 p7 e4 I
t1 = 5*pi/180;/ |3 n& t; Z0 x M' ^( Y$ F7 M
t2 = pi/2;- h9 _9 |+ Y: ?. ?4 c& h
A=R1.^2 * sin(t1) * tan(t1) * tan(t2/2);- a" C6 U' Z$ r; u5 Z8 P
, V) _) g) h8 Q) ~
S={'aaa'};% 定义一个元胞数组,并随意存储一个字符串
& n8 _8 }; c# I+ @- W0 r; afor KB = [0, 0.2, 0.4, 0.6, 0.8, 1]
|7 t% ` g0 z( S' C a1 = (2*1500*9*0.8*10^6/20)/(1500*9*13888.2*(2*1.313+KB+1)/60);/ O) O- s |9 `- B) s" D0 Z) ?
b1 = 1500/(2*1.313+KB+1);$ G' v- k1 V# L, }8 k& t
y0 = 0.1;9 E; u: r5 b8 S3 P% l3 A9 o$ i- |
yp0 = b1*sin(c);
: _ ^5 T# T5 d1 } [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);' W2 L3 Q: u. H, d
h = plot(x,y);
- K1 N' I6 y8 W' L3 N kb = sprintf('KB=%.1f',KB);%将KB的转化为字符串存在一个变量中- x& u" `* Z+ k3 \; e
S(i) = cellstr(kb);%将变量kb存储在S元胞数组中3 E7 O* ]4 f/ P8 n5 H/ s$ o
hold on% ^* ]$ {# y7 i6 R/ ]9 H
end
; J( p4 [& ?1 S4 h3 d; ?4 s' Tlegend(S);
# N. Q/ T+ _7 jgrid1 X! p" |/ R, Y8 W9 v! J
x6 ` [( b9 V6 Z6 ^此外，能否在for循环中使每个plot出来的曲线的线型不一样，比如说KB=0时是虚线，KB=1时是实线等等# W/ L) [8 [/ q, q
请教大神，非常感谢！! ]6 u3 D! l$ E0 A: k0 Z9 r
出现的问题如下图7 _6 t# x4 u) @9 C$ a
& d2 }/ q0 k+ E) A- P; R& E
9 b! i: X3 c5 _ l4 J
3 o& X6 }# S, O( Q- ?& [* K( U3 J |
|
/1 
发表于 2020-5-27 16:24
收藏
淘帖
支持！
反对！