不知道哪里出了问题

shelby

我通过for循环求解一个带有变系数微分方程后，plot出各个方程的曲线，想将各个曲线对应legend打印出来，但是始终无法成功，不知道哪里出了问题，代码如下：
c = (4-19.2/2)*pi/180;
2 @6 e( v4 l6 ~& ~( `R1 = 38;
t1 = 5*pi/180;
t2 = pi/2;
; O; Z0 Y1 f+ T2 F5 v, |A=R1.^2 * sin(t1) * tan(t1) * tan(t2/2);

S={'aaa'};% 定义一个元胞数组,并随意存储一个字符串
for KB = [0, 0.2, 0.4, 0.6, 0.8, 1]
    a1 = (2*1500*9*0.8*10^6/20)/(1500*9*13888.2*(2*1.313+KB+1)/60);
    b1 = 1500/(2*1.313+KB+1);
* _/ p' u# P; u- j    y0 = 0.1;
    yp0 = b1*sin(c);
& H  a  |2 M1 t2 e    [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);
# M( X0 d2 s; T/ d% X    h = plot(x,y);
. w. S) b+ m1 ]7 t    kb = sprintf('KB=%.1f',KB);%将KB的转化为字符串存在一个变量中
: @- t! }7 G: @5 J$ j0 ^    S(i) = cellstr(kb);%将变量kb存储在S元胞数组中
5 h1 U( C. T% A; f    hold on
5 A4 n/ M; p  z! C2 c) I. D& iend
$ N; V! a: F6 F9 y3 R' W* E: T) `legend(S);
! l/ G1 }. d* H) \1 o- j' Vgrid
# }1 D3 z! T% l! t. H
此外，能否在for循环中使每个plot出来的曲线的线型不一样，比如说KB=0时是虚线，KB=1时是实线等等
( u0 L  d0 x0 _8 w9 V. m; C, l请教大神，非常感谢！
出现的问题如下图

, i, P: p- F; j7 K, q
/ {* p& H  c! {2 F4 L) j
4 v( B0 u( i" r+ D4 T9 v/ n

Taio

$ J5 V) L+ y7 p6 k3 t8 B3 ~; zS(i) = cellstr(kb);这句话有问题吧
) n' u+ \/ N  e: ^不应该是
$ g" {: ^, b. W% }S{i} = cellstr(kb);吗？我没调试过，只是第一眼看见这个代码的感觉是这样的。