微分方程省公开课一等奖全国示范课微课金奖课件.pptx

1第9章

MATLAB微分方程求解第1页

29.1常微分方程初值问题常微分方程初值问题是科学计算中常见问题，MATLAB提供了求解常微分方程（ODE）初值问题一系列ode函数。ode函数主要采取Runge-Kutta法求解常微分方程。惯用ode函数函数名含义函数名含义ode23普通2-3阶法解ODEode23s变阶法解刚性odeode45普通4-5阶法解ODEode23t解适度刚性odeode113普通变阶法解ODEode23tb低阶法解刚性ode第2页

39.1.1ODE机理将高阶常微分方程初值问题，转变为1阶常微分方程组初值问题。9.1.2普通2-3阶法解ODE例9-1编制常微分方程函数文件dfun.mfunctiondy=dfun(x,y)dy=y-2*x/y编制以下计算程序clear;clc;X12=linspace(0,5,50);[X,Y]=ode23(dfun,X12,2);plot(X,Y)运行结果如图9-1所表示图9-1第3页

49.1.3普通4-5阶法解ODE例9-2编制微分方程组函数文件dfun.mfunctiondy=dfun(x,y)dy(1,1)=-2*y(1)+y(2)+2*sin(x);dy(2,1)=10*y(1)-9*y(2)+9*(cos(x)-sin(x));编制以下程序clear;clc;X12=linspace(0,9,30);ode45(dfun,X12,[2,3]);运行结果如图9-2所表示图9-2第4页

59.2常微分方程边值问题9.2.1边值问题介绍普通来说，微分方程边值问题可能有解、也可能无解，可能有唯一解、也可能有没有数解。在假定有唯一解前提下，边值问题有三种基本解法：迭加法：假如微分方程和边界条件是线性，边值问题能够转化为初值问

题。试射法：将问题转化为对漏缺初值搜索，一旦漏缺初值确定，就能够按

初值问题求解。松弛法：先猜测满足边界条件区间网点上解值，然后利用微分方程进

行迭代改进。第5页

69.2.2bvp4c求解思绪为求解常微分方程边值问题，MATLAB提供了bvp4c函数文件。该文件是依据有限元法中配置法编写。其解题基本思绪以下：把待解问题转化为标准边值问题为期望解指定初始猜测利用原微分方程，结构残差函数用原微分方程和边界条件，借助迭代使残差不停减小为了bvp4c使用方便，MATLAB提供了能产生很好猜测解bvpinit指令，减轻了用户进行初始猜测负担。另外解算指令bvp4c也能处理包含未知参数微分方程边值问题。第6页

79.2.3bvp4c配套指令例9-3编制微分方程组函数文件dfun.mfunctiondy=dfun(x,y,c)dy(1,1)=y(2);dy(2,1)=-c*abs(y(1));编制残差方程组函数文件bfun.mfunctionres=bfun(ya,yb,c)res(1,1)=ya(1);res(2,1)=yb(1)+2编制以下计算程序sol=bvp4c(dfun,bfun,sinit,[],c);plot(sol.x,sol.y(1,:),r--)holdonc=1.5;sol=bvp4c(dfun,bfun,sinit,[],c);plot(sol.x,sol.y(1,:),b-)legend(c=1,c=1.5)运行结果如图9-3所表示clear;clc;x=linspace(0,4,5);v=[1;0];sinit=bvpinit(x,v);c=1;图9-3第7页

89.3常微分方程解析解9.3.1微分方程通解例9-4编制以下程序clear;clc;fun=Dx=-a*tx=dsolve(fun,t)运行结果为fun=Dx=-a*tx=-1/2*a*t^2+C1例9-5编制以下程序clear;clc;fun1=Df=f+g;fun2=Dg=g-f;[f,g]=dsolve(fun1,fun2,t)运行结果为f=-exp(t)*(C1*cos(t)-C2*sin(t))g=exp(t)*(C1*sin(t)+C2*cos(t))第8页

99.3.2微分方程特解例9-6编制以下程序clear;clc;fun=D2y=4x=dsolve(fun,y(0)=1,y(1)=3,x)运行结果为fun=D2y=4x=2*x^2+1例9-7编制以下程序clear;clc;fun1=Df=f+g;fun2=Dg=g-f;v1=f(0)=2;v2=g(0)=2;[f