MATLAB解微分方程ODE初值问题的数值解PDE问题的数值解[精品].ppt

* * PDETOOL 有限元方法 * F=inline([y(2);-y(1)],t,y) ode23(F,[0 2*pi],[1;0]) opts=odeset(‘reltol’,1.e-4,’abstol’,1.e-6,’outputfcn’) Examples * ode23(@twobody,[0 2*pi],[1;0;0;1]); 0 1 2 3 4 5 6 7 -1.5 -1 -0.5 0 0.5 1 1.5 Examples * y0=[1;0;0;3]; ode23(@twobody,[0 2*pi],y0); 0 1 2 3 4 5 6 7 -2 0 2 4 6 8 10 12 14 16 18 * y0=[1;0;0;3]; [t,y]=ode23(@twobody,[0 2*pi],y0); plot(y(:,1),y(:,2)); axis equal * y0=[1;0;0;3]; [t,y]=ode23(@twobody,[0 2*pi],y0); plot(y(:,1)) plot(y(:,2)) * A problem is stiff if the solution being sought is varying slowly, but there are nearby solutions that very rapidly, so the numerical method must take small steps to obtain satisfactory results. A model of flame propagation（火焰燃烧）: y是球的半径，y^2和y^3与球的表面积和体积有关 想一下，点燃一根火柴，光球迅速增长，到达一个关键的大小，然后维持它的大小（由于进入球内氧气和消耗的氧气平衡） Stiff Problem(刚性问题) * 0 10 20 30 40 50 60 70 80 90 100 0 0.2 0.4 0.6 0.8 1 1.2 1.4 eta=0.02;sym y; F=inline(‘y^2-y^3’,’t’,’y’);ode23(F,[0 2/eta],eta); * eta=0.00002;sym y; F=inline(‘y^2-y^3’,’t’,’y’);ode23(F,[0 2/eta],eta); * 0 1 2 3 4 5 6 7 8 9 10 x 10 4 0 0.2 0.4 0.6 0.8 1 1.2 1.4 ode23s eta=0.00002; ode23s(inline(y^2-y^3,t,y),[0 2/eta],eta); * 例 蛇形曲线的常微分方程初值问题 MATLAB数值求解命令 F=inline(1./(1+x.^2)-2*y.^2)； ode23(F,[0,6],0) 输出结果为图形 * [T,y]=ode23(f,[0,6],0)将得到自变量和函数的离散数据 T =0 0.0001 0.0005 0.0025 0.0125 0.0496 0.1085 0.1863 0.2837 0.4091 0.5991 0.8513 1.0567 1.2680 1.5110 1.8050 2.1788 2.6842 3.2842 3.8842 4.4842 5.0842 5.6842 6.0000 y =0 0.0001 0.0005 0.0025 0.0125 0.0495 0.1073 0.1800 0.2626 0.3505 0.4411 0.4944 0.5000 0.4868 0.4607 0.4242 0.3793 0.3270 0.2783 0.2411 0.2122 0.1891 0.1705 0.1620 * 例 洛伦兹模型由如下常微分方程组描述 取 ?=8/3，?=10，?=28。 初值：x(0)=0，y(0)=0，z(0)=0.01。 利用MATLAB求解常微分方程数值解命令计算出t∈[0，80]内，三个未知函数的数据值，并绘出相空间在Y-X平面的投影曲线 * 气象学家Lorenz提出一篇论文，名叫「一只蝴蝶拍一下翅膀会不会在Taxas州引起龙卷风？」论述某系统如果初期条件差一点点，结果会很不