exm09714_0.m

来自「这是本人收集的一些程序源代码」· M 代码 · 共 22 行

%exm09714_0.m			例9.7.1.4_1的主运行程序
%						它调用exm9714_1.mdl
clc
xsym=dsolve('D2x+100*Dx+0.9999*x=0','x(0)=1,Dx(0)=0','t')
dxsym=diff(xsym,'t')  

t=(0:5000)/10;x=eval(char(dxsym));						%产生绘解析解图形的数据
Tspan=500;
opts=simset('Solver','ode45');	
[tt1,xx1,s]=sim('exm09714_1',Tspan,opts);				%ODE45解Exm09714_1模型
opts=simset('Solver','ode15s');
[tt2,xx2,s]=sim('exm09714_1',Tspan,opts);				%ODE15S解Exm09714_1模型
clf,shg
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)
title('对于Stiff方程，三种不同解法的比较。')
ns1=length(xx1);											%ODE45解点数
ns2=length(xx2);											%ODE15S解点数  
disp('ODE45解点数'),disp(ns1)
disp(' ')
disp('ODE15S解点数'),disp(ns2)