for i=1:length(c)
h_opt=odeset;y0=[x0;v0];
[t,x]=ode45(@dzqode,[0,tf],y0,h_opt,F0,w,c,k,m,i);
if c(i)/(2*m)==sqrt(k/m)
plot(t,x(:,1),'r');hold on;plot(t,x(:,2),'b');
elseif c(i)/(2*m)>sqrt(k/m)
plot(t,x(:,1),'g');hold on;
plot(t,x(:,2),'c');
else c(i)/(2*m)<sqrt(k/m)
plot(t,x(:,1),'k');hold on;plot(t,x(:,2),'m');
end;
end;
if length(c)>1
hold on;
end;
end;
xlabel('时间(s)');
ylabel('位移(m)和速度(m/s)');
title('位移和速度响应曲线');
hold off;
这个程序怎么用legend加图例啊?因为c是不定的,我搞了半天也搞不定,好象legend刷新什么的.