exm080714m_1.m

matlab mfiles实例分析

%exm080714m_1.m
clc
disp('按任意键，开始运行。')
disp('(1) 用符号计算求此方程解析解')
pause
xsym=dsolve('D2x+100*Dx+0.9999*x=0','x(0)=1,Dx(0)=0','t')
dxsym=diff(xsym,'t')

disp(' ')
disp('按任意键，继续运行。')
disp('(2) 三种计算结果的图形比较')
pause
t=(0:5000)/10;x=eval(char(dxsym));	%产生绘解析解图形的数据
Tspan=500;
opts=simset('Solver','ode45');	
[tt1,xx1,s]=sim('exm080714_1',Tspan,opts);%ODE45解Exm080714_1模型
opts=simset('Solver','ode15s');
[tt2,xx2,s]=sim('exm080714_1',Tspan,opts);%ODE15S解Exm080714_1模型
plot(t,x,'k',tt1,xx1(:,2),'b:',tt2,xx2(:,2),'r-.')	%绘比较图
axis([246 247 -8.55e-4  -8.35e-4])						%局部放大
legend('Symbolic','ODE45','ODE15S',0)
ns1=length(xx1)							%ODE45解点数
ns2=length(xx2)							%ODE15S解点数