matlab数学实验实验报告欧拉公式ROSSLER微分方程.docVIP

数学实验—实验报告 一、实验项目： 二、实验目的和要求 1、本章将对人口变化、动物种群变迁、网络系统的可靠性分析，介绍微分方程（组）的模型建立、数值解和图形解等方法，并用MATLAB几何直观地展示各种求解方法的求解结果。 利用欧拉公式求解方程 三、实验题目 问题一：求微分方程的解析解，并画出它们的图形， y’=y+2x, y(0)=1, 0x1; 问题二：用向前欧拉公式和改进的欧拉公式求方程 y’=y-2x/y, y(0)=1的数值解（0≤x ≤1 , h=0.1） 要求编写程序。 问题三：Rossler微分方程组 当固定参数b=2, c=4时，试讨论随参数a由小到大变化（如a∈(0,0.65]）而方程解的变化情况，并且画出空间曲线图形，观察空间曲线是否形成混沌状？ 问题四：水的流出时间 一横截面积为常数A，高为H的水池内盛满水，由池底一横截面积为B的小孔放水。设水从小孔流出的速度为v=(2gh)0.5，求在任意时刻的水面高度和将水放空所需的时间。 时间t→高度h。 问题五：考虑相互竞争模型 两种相似的群体之间为了争夺有限的同一种事物来源和生存空间而进行生死存亡竞争时，往往是竞争力较弱的种群灭亡，而竞争力较强的种群达到环境容许的最大数量 假设有甲、乙两个生物种群，当它们各自生存于一个自然环境中，均服从Logistics规律。 三、实验过程 问题一： 用matlab编写代码： x=[0,1] y=dsolve(Dy=y+2*x) y=dsolve(Dy=y+2*x, y(0)=1, x) ezplot(x,y) 输出：y =-2*x+exp(t)*C1 （通解）y =-2*x-2+3*exp(x) 画图：x=0:0.01:1;y =-2*x-2+3*exp(x);plot(x,y) 问题二： 1、分析： 解：（1）解析解法得到其精确解： （2） 向前欧拉法： 迭代公式为?，其中 （3）改进欧拉法： 迭代公式为?，其中 2、Matlab编码 x1(1)=0;y1(1)=1;y2(1)=1;h=0.1; for k=1:10 x1(k+1)=x1(k)+h; y1(k+1)=(1-h)*y1(k)+2*h*x1(k)/y1(k); y2(k+1)=(1+h)*y2(k)+(h*h)/2-h*x1(k)/y2(k)-h*(x1(k)+h)/(y2(k)+h); end x=0:0.1:1; y=(2*x+1).^(1/2); x1=x1(1:11),y=y(1:11),y1=y1(1:11),y2=y2(1:11), plot(x,y,x1,y1,k:,x1,y2,r--) 显示图像及结果： x1 = 0 0.1000 0.2000 0.3000 0.4000 0.5000 0.6000 0.7000 0.8000 0.9000 1.0000 y = 1.0000 1.0954 1.1832 1.2649 1.3416 1.4142 1.4832 1.5492 1.6125 1.6733 1.7321 y1 = 1.0000 0.9000 0.8322 0.7971 0.7926 0.8143 0.8557 0.9103 0.9731 1.0402 1.1092 y2 = 1.0000 1.0959 1.1847 1.2679 1.3468 1.4222 1.4948 1.5653 1.6340 1.7016 1.7683 图中，蓝色曲线是精确解，红色曲线是向前欧拉法曲线，黑色曲线是改进后欧拉法曲线 问题三 1、matlab编程 function r=rossler(t,x) global a; global b; global c; r=[-x(2)-x(3);x(1)+a*x(2);b+x(3)*(x(1)-c)]; global a; global b; global c; b=2; c=4; t0=[0,200]; for a=0:0.02:0.65 [t,x]=ode45(rossler,t0,[0,0,0]); subplot(1,2,1); plot(t,x(:,1),r,t,x(:,2),g,t,x(:,3),b); title(x(oìé?),y(?ìé?),z(àoé?)??t±??ˉ?é??);xlabel(t); subplot(1,2,2