推荐于2016-02-12 · 知道合伙人互联网行家
程序
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clf;t=6*pi*(0:100)/100;
y=1-exp(-0.3*t).*cos(0.7*t);
plot(t,y,'r-','LineWidth',3)%<3>
hold on
tt=t(find(abs(y-1)>0.05));ts=max(tt);%<5>
plot(ts,0.95,'bo','MarkerSize',10)%<6>
hold off
axis([-inf,6*pi,0.6,inf])
set(gca,'Xtick',[2*pi,4*pi,6*pi],'Ytick',[0.95,1,1.05,max(y)])%<9>
set(gca,'XtickLabel',{'2*pi';'4*pi';'6*pi'})%<10>
set(gca,'YtickLabel',{'0.95';'1';'1.05';'max(y)'})%<11>
grid on
text(13.5,1.2,'\fontsize{12}{\alpha}=0.3')%<13>
text(13.5,1.1,'\fontsize{12}{\omega}=0.7')%<14>
cell_string{1}='\fontsize{12}\uparrow';%<15>
cell_string{2}='\fontsize{16} \fontname{隶书}镇定时间';
cell_string{3}='\fontsize{6}';
cell_string{4}=['\fontsize{14}\rmt_{s} = 'num2str(ts)];%<18>
text(ts,0.85,cell_string,'Color','b','HorizontalAlignment','Center')%<19>
title('\fontsize{14}\it y = 1 - e^{ -\alpha t}cos{\omegat}')%<20>
xlabel('\fontsize{14} \bft \rightarrow')
ylabel('\fontsize{14} \bfy\rightarrow')%<22>
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%