数学建模与数学实验重点分析.ppt

西安建筑科技大学 数学建模与数学实验 1. 求微分方程组 解：先写出微分方程的函数文件： function dx=example1(t,x) dx=zeros(2,1); dx(1)=x(2)+cos(t); dx(2)=sin(2*t); 课后作业： 西安建筑科技大学 数学建模与数学实验 再编写求解该微分方程的MATLAB程序： [t X]=ode45(@example1,[0,50],[0.5,-0.5]); figure(1) subplot(211),plot(t,X(:,1),b-),title(x_1),grid on subplot(212),plot(t,X(:,2),r-),title(x_2),grid on figure(2) plot(X(:,1),X(:,2),g-),title(轨迹),grid on xlabel(x_1),ylabel(x_2) 西安建筑科技大学 数学建模与数学实验 西安建筑科技大学 数学建模与数学实验 西安建筑科技大学 数学建模与数学实验 湖水污染物质量的变化量=流入湖水的污染物质量-流出湖水的污染物质量。其数学表达式为： 或 西安建筑科技大学 数学建模与数学实验 西安建筑科技大学 数学建模与数学实验 r=dsolve(Dw=r*(w1w)/V,w(0)=0.1,t) 求解得到： r =w1 - exp(-(r*t)/V)*(w1 - 1/10) 即 将题目数值带入求解： r =0.05 - exp(-((4.4534e+10)*100)/(5.941e+12))*(0.05 - 1/10) 得到：r =0.0736  西安建筑科技大学 数学建模与数学实验 * * * * * * * * * * * * * * * * * * * 西安建筑科技大学 数学建模与数学实验 微分方程的MATLAB求解 微分方程的解析解 在MATLAB中，求解微分方程(组)的函数是desolve，调用格式如下： r = dsolve(eq1,eq2,..., cond1,cond2,... , v) 或 r = dsolve(eq1,eq2,...,cond1,cond2,... , v) 微分方程模型 西安建筑科技大学 数学建模与数学实验 解：命令为 r=dsolve(Dy=a*x,x) 结果为r=1/2*a*x^2+C1。 解：命令为 r=dsolve(Dx=x+sin(t),x(0)=1,t) 结果为r =-1/2*cos(t)-1/2*sin(t)+3/2*exp(t)。 其中eq1,eq2,…用来表示常微分方程(组)，cond1, cond2,…表示初始或边界条件，v表示自变量，缺省时默认为t。 西安建筑科技大学 数学建模与数学实验 解：命令为 r=dsolve(D2y-Dy^2/y=0,y(0)=1,Dy(0)=2, x) 结果为r =exp(2*x)。 解：命令为[x,y]=dsolve(Dx+5*x+y=exp(t),Dy-x-3*y=exp(2*t),t) 西安建筑科技大学 数学建模与数学实验 2.2 微分方程的数值解 解含初始条件的常微分方程的MATLAB格式 [t,Y] = solver(‘odefun’,tspan,y0) 其中，odefun表示以函数文件存储的微分方程(不含初始条件)，tspan为二维向量，用来表示区间的起点和终点，y0表示自变量取tspan区间起点时的初始条件。solver指ode45、ode23、 ode113、ode15s、ode23s、ode23t和ode23tb七种求解方法中的一种。 西安建筑科技大学 数学建模与数学实验 西安建筑科技大学 数学建模与数学实验 例8. 求下列微分方程 解：先建立函数文件来表示微分方程： function dy = rigid(t,y) dy = zeros(3,1); % dy应为列向量 dy(1) = y(2) * y(3); dy(2) = -y(1) * y(3); dy(3) = -0.51 * y(1) * y(2); 西安建筑科技大学 数学建模与数学实验 再在脚本函数或命令窗口中输入： [t,Y] = ode45(@rigid,[0 12],[0 1 1]); plot(t,Y(:,1),ro-,t,Y(:,2),g*:,t,Y(:,3),b