weifen.m

计算微分方程程序

clf;%清屏
t=6*pi*(0:100)/100;%矢量化赋值
y=1-exp(-0.3*t).*cos(0.7*t);%函数表示
tt=t(find(abs(y-1)>0.05));%
ts=max(tt);		%求最大值					
plot(t,y,'r-','LineWidth',3)%绘制图形
axis([-inf,6*pi,0.6,inf])%inf 表示无穷大
set(gca,'Xtick',[2*pi,4*pi,6*pi],'Ytick',[0.95,1,1.05,max(y)])	%<5>
grid on																	%<6>
title('\it y = 1 - e^{ -\alphat}cos{\omegat}')	%注意输入方法
text(13.5,1.2,'\fontsize{12}{\alpha}=0.3')
text(13.5,1.1,'\fontsize{12}{\omega}=0.7')							%<9>
hold on;
plot(ts,0.95,'bo','MarkerSize',10);
hold off				%<10>
cell_string{1}='\fontsize{12}\uparrow';								%<11>
cell_string{2}='\fontsize{16} \fontname{隶书}镇定时间';			%<12>
cell_string{3}='\fontsize{6}  ';										%<13>
cell_string{4}=['\fontsize{14}\rmt_{s} = ' num2str(ts)];			%<14>
text(ts,0.85,cell_string)												%<15>
xlabel('\fontsize{14} \bft \rightarrow')							%<16>
ylabel('\fontsize{14} \bfy \rightarrow')				 			%<17>