下载文档原格式
第五章 控制系统仿真§5.2 微分方程求解方法以一个自由振动系统实例为例进行讨论。
如下图1所示弹簧-阻尼系统，参数如下： M=5 kg, b=1 N.s/m, k=2 N/m, F=1NF图1 弹簧-阻尼系统假设初始条件为：00=t 时，将m 拉向右方，忽略小车的摩擦阻力，m x 0)0(= s m x /0)0(=•求系统的响应。
)用常微分方程的数值求解函数求解包括ode45、ode23、ode113、ode15s 、ode23s 等。
wffc1.m myfun1.m一、常微分方程的数值求解函数ode45求解 解：系统方程为 F kx x b x m =++•••这是一个单变量二阶常微分方程。
将上式写成一个一阶方程组的形式，这是函数ode45调用规定的格式。
令： x x =)1( （位移）)1()2(••==x x x （速度） 上式可表示成：⎥⎦⎤⎢⎣⎡--=⎥⎦⎤⎢⎣⎡=⎥⎥⎦⎤⎢⎢⎣⎡••)1(*20)2(*101)2()2()2()1(x x x x x x x 下面就可以进行程序的编制。
%写出函数文件myfun1.mfunction xdot=myfun1(t,x)xdot=[x(2);1-10*x(2)-20*x(1)];% 主程序wffc1.mt=[0 30];x0=[0;0];[tt,xx]=ode45(@myfun1,t,x0);plot(tt,yy(:,1),':b',tt,yy(:,2),'-r') legend('位移','速度')title('微分方程的解 x(t)')二、方法2：F kx x b x m =++•••251)()()(2++==s s s F s X s G%用传递函数编程求解ksys1.mnum=1;den=[5 1 2];%printsys(num,den)%t=0:0.1:10;sys=tf(num,den);figure(1)step(sys)figure(2)impulse(sys)figure(3)t=[0:0.1:10]';ramp=t;lsim(sys,ramp,t);figure(4)tt=size(t);noise=rand(tt,1);lsim(sys,noise,t)figure(5)yy=0.1*t.^2;lsim(num,den,yy,t)w=logspace(-1,1,100)';[m p]=bode(num,den,w);figure(6)subplot(211);semilogx(w,20*log10(m)); grid onsubplot(212);semilogx(w,p)grid on[gm,pm,wpc,wgc]=margin(sys)figure(7)margin(sys)figure(8)nyquist(sys)figure(9)nichols(sys)方法3：F kx x b x m =++•••125=++•••x x xx x x 4.02.02.0--=•••% 主程序wffc1.mt=[0 30];x0=[0;0];[tt,yy]=ode45(@myfun1,t,x0);figure(1)plot(tt,yy(:,1),':b',tt,yy(:,2),'-r') hold onplot(tt,0.2-0.2*yy(:,2)-0.4*yy(:,1),'-.k ')legend('位移','速度','加速度') title('微分方程的解')figure(2)plot(yy(:,1),yy(:,2))title('平面相轨迹')%写出函数文件myfun1.mfunction xdot=myfun1(t,x)xdot=[x(2);0.2-0.2*x(2)-0.4*x(1)];。
一、概述微分方程是自然科学和工程技术中常见的数学模型，它描述了连续系统的变化规律。
在实际应用中，求解微分方程是一项重要且复杂的工作。
而matlab是一种常用的科学计算软件，它提供了丰富的数学函数和工具，能够辅助工程师和科学家在求解微分方程方面取得良好的效果。
二、matlab差分法求解微分方程的基本原理差分法是一种常见的数值求解微分方程的方法。
它基于微分的定义，将微分方程中的微分运算用差分逼近来进行计算。
在matlab中，可以利用内置的数学函数和工具，通过差分法求解微分方程，得到数值解或者近似解。
三、matlab中使用差分法求解常微分方程的步骤1. 确定微分方程的类型和边界条件需要明确所要求解的微分方程是什么类型的，以及其所对应的边界条件是什么。
这对于后续的数值求解过程非常重要。
在matlab中，可以利用符号变量和函数来表示微分方程和边界条件。
2. 将微分方程离散化接下来，需要将微分方程进行离散化处理，将微分方程中的微分运算用差分逼近来进行计算。
这一步需要根据微分方程的具体形式和求解精度选择合适的差分方法，常见的有前向差分、后向差分和中心差分等方法。
3. 构建代数方程组将离散化后的微分方程转化为代数方程组。
这一步需要根据微分方程的离散化表达式和边界条件，利用matlab的矩阵和向量运算功能，构建代数方程组。
4. 求解代数方程组利用matlab的求解函数，求解构建得到的代数方程组，得到微分方程的数值解或者近似解。
在求解过程中，需要注意数值稳定性和收敛性，以及选择合适的数值积分方法和迭代算法。
四、实例：使用matlab差分法求解一阶常微分方程为了更好地理解matlab中使用差分法求解微分方程的过程，以下将通过一个具体的实例来演示。
假设要求解如下的一阶常微分方程：dy/dx = -2x + 1, y(0) = 11. 确定微分方程的类型和边界条件根据给定的方程，可以确定它是一阶常微分方程，且给定了初始条件y(0) = 1。
实验二： 微分方程模型Matlab 求解与分析一、实验目的[1] 掌握解析、数值解法，并学会用图形观察解的形态和进行解的定性分析； [2] 熟悉MATLAB 软件关于微分方程求解的各种命令；[3] 通过范例学习建立微分方程方面的数学模型以及求解全过程； [4] 熟悉离散 Logistic 模型的求解与混沌的产生过程。
二、实验原理1. 微分方程模型与MATLAB 求解解析解用MATLAB 命令dsolve(‘eqn1’,’eqn2’, ...) 求常微分方程（组）的解析解。
其中‘eqni'表示第i 个微分方程，Dny 表示y 的n 阶导数,默认的自变量为t 。
（1） 微分方程 例1 求解一阶微分方程 21y dxdy+= (1) 求通解 输入：dsolve('Dy=1+y^2')输出：ans =tan(t+C1)（2）求特解 输入：dsolve('Dy=1+y^2','y(0)=1','x')指定初值为1，自变量为x 输出：ans =tan(x+1/4*pi)例2 求解二阶微分方程 221()04(/2)2(/2)2/x y xy x y y y πππ'''++-=='=-原方程两边都除以2x ，得211(1)04y y y x x'''++-= 输入:dsolve('D2y+(1/x)*Dy+(1-1/4/x^2)*y=0','y(pi/2)=2,Dy(pi/2)=-2/pi','x')ans =- (exp(x*i)*(pi/2)^(1/2)*i)/x^(1/2) +(exp(x*i)*exp(-x*2*i)*(pi/2)^(3/2)*2*i)/(pi*x^(1/2))试试能不用用simplify 函数化简 输入: simplify(ans)ans =2^(1/2)*pi^(1/2)/x^(1/2)*sin(x) （2）微分方程组例3 求解 d f /d x =3f +4g ; d g /d x =-4f +3g 。
文章主题：探索数学求解软件Matlab在微分代数方程求解中的应用1. 引言微分代数方程（DAE）是描述物理系统中的相互依赖性和复杂性的数学模型。
解决这类方程对于现代科学和工程领域至关重要。
Matlab作为一种强大的数学计算软件，在微分代数方程求解中具有独特的优势。
本文将从简单到复杂的方式，探讨Matlab在DAE求解中的应用，并共享个人见解。
2. DAE的基本概念微分代数方程是描述包含未知函数及其导数或导数与未知函数的组合的方程。
通常的形式为F(x, x', t) = 0，其中x为未知函数，x'为其导数，t为自变量。
在实际应用中，这些方程往往伴随着初始条件和边界条件。
3. Matlab在解常微分方程（ODE）中的应用Matlab拥有丰富的ODE求解函数，如ode45、ode23等，可用于求解各种常微分方程。
这些函数可以自动选择适当的数值积分方法，并提供了方便的接口和参数设置，极大地简化了求解过程。
4. DAE求解方法的挑战与ODE相比，DAE的求解更具挑战性。
由于包含了代数变量和微分变量，常规的数值积分方法难以直接应用。
而且，方程的初始条件和边界条件也增加了求解的复杂性。
5. Matlab在DAE求解中的工具Matlab提供了一系列专门用于DAE求解的函数和工具包，如dare和ddesd等。
这些工具在模型建立、数值解法选择、收敛性分析等方面都具有独特的优势。
6. 案例分析：用Matlab求解电路模型的DAE以电路模型的DAE为例，通过Matlab可以快速建立系统的数学模型，并进行数值求解。
通过对参数的调节和模型的分析，可以更好地理解电路的动态特性，帮助优化设计和故障诊断。
7. 总结与展望通过本文的探讨，我们更深入地了解了Matlab在微分代数方程求解中的重要性和应用。
在未来，随着科学技术的发展，Matlab在此领域的功能和性能将得到进一步的提升，为工程科学领域提供更强大的支持。
个人观点：Matlab作为一种综合性的科学计算软件，对微分代数方程的求解起着至关重要的作用。
matlab微分方程模型Matlab微分方程模型是一种基于Matlab软件的数学建模方法，用于解决微分方程相关的问题。
微分方程是描述物理、工程和数学问题的重要工具，通过建立微分方程模型，可以对各种现象进行定量分析和预测。
在Matlab中，可以使用ode45函数求解常微分方程（ODE）或者ode15s函数求解刚性ODE。
这些函数可以通过数值方法近似求解微分方程的解析解，从而得到问题的数值解。
具体来说，可以通过在Matlab中定义微分方程的右侧函数，然后使用相应的ode函数进行求解。
例如，考虑一个简单的一阶线性微分方程模型：dy/dx = -ky，其中k为常数。
我们可以通过在Matlab中定义这个微分方程的右侧函数，并使用ode45函数求解。
具体步骤如下：1. 在Matlab中定义微分方程的右侧函数：function dydx = myODE(x,y)k = 0.1; % 设定常数k的值dydx = -k*y;end2. 使用ode45函数求解微分方程：xspan = [0 10]; % 设定求解区间y0 = 1; % 设定初始条件[x,y] = ode45(@myODE, xspan, y0);3. 绘制得到的数值解：plot(x,y);xlabel('x');ylabel('y');title('Solution of dy/dx = -ky');通过以上步骤，我们可以得到微分方程dy/dx = -ky的数值解，并绘制出解的图像。
这个简单的例子展示了如何使用Matlab微分方程模型求解微分方程。
除了一阶线性微分方程，Matlab微分方程模型还可以用于解决更复杂的微分方程问题，包括高阶线性微分方程、非线性微分方程、偏微分方程等。
通过定义相应的微分方程函数和合适的求解方法，可以在Matlab中进行数值求解。
此外，Matlab还提供了丰富的绘图和分析工具，可以对微分方程的解进行可视化和进一步分析。
解微分⽅程有两种解，⼀种是解析解，⼀种是数值解，这两种分别对应不同的解法利⽤dsolve 函数进⾏求解syms x;s = dsolve('eq1,eq2,...', ’cond1,cond2,...', 'v');%eq ：微分⽅程%cond ：条件%v ：独⽴变量%形如：⽅程：y'= f(t,y),初值：y(t0) = y0dsolve('Du = 1+ u^2','t')ans =tan(C2 + t)1i-1i求的解析解s = dsolve('D2y=3*y+2*x','x');% D2y ⽤以表⽰y 的⼆阶导数，默认是以t 为⾃变量的，所以最好指明⾃变量为x.syms y(x);s = dsolve([diff(y,x,2) == 3*y+2*x], [y(0) == 5])% diff 内依次是函数、⾃变量、微分阶数，⽅程⽤==表⽰相等⽽不是赋值求初值问题s = dsolve('Dy = y - 2*t / y','y(0) =1');求边界问题s = dsolve('x*D2y - 3*Dy =x^2','y(1)=0','y(5) = 0','x');求解⽅程s=dsolve('D2y =cos(2*x) - y','y(0) =1','Dy(0) = 0','x');simplify(s);(eqn,cond,‘IgnoreAnalyticConstraints’,false) %设置不化简结果求解⽅程组[f,g]= dsolve('Df = f + g','Dg = -f + g','f(0)=1','g(0) = 2','x');⽤Matlab 求解微分⽅程⽤Matlab 求解微分⽅程解析解1.求解析解2.初值问题3.边界问题4.⾼阶⽅程5.⽅程组问题⼀些例⼦dsolve('D2y+4*Dy+29*y = 0','y(0) = 0','Dy(0)= 15 ','x')ans =3*sin(5*x)*exp(-2*x)[x y z] = dsolve('Dx = 2*x-3*y+3*z','Dy = 4*x-5*y+3*z','Dz = 4*x-4*y+2*z')x =C7*exp(2*t) + C8*exp(-t)y =C7*exp(2*t) + C8*exp(-t) + C9*exp(-2*t)z =C7*exp(2*t) + C9*exp(-2*t)%可以对其进⾏简化操作x = simplify(x)x = C7*exp(2*t) + C8*exp(-t)y = simplify(y)y =exp(-2*t)*(C9 + C8*exp(t) + C7*exp(4*t))%龙格库塔法(Runge-Kutta 法)xfun=@(t,x)0.3.*x.*(1-x/8); %定义赋值函数r=0.3,k=8[tout,xout]=ode45(fun,[0,40],0.1) %⽅程数值解，四五阶RK 法[tout,xout]=ode23(xfun,[t0,tfinal],x0) %⼆三阶RK 法%%ode 系列数值求解形如 / = ( , )的微分⽅程组, 并绘图。
第四讲Matlab求解微分方程（组）理论介绍：Matlab求解微分方程（组）命令求解实例：Matlab求解微分方程（组）实例实际应用问题通过数学建模所归纳得到的方程，绝大多数都是微分方程，真正能得到代数方程的机会很少.另一方面，能够求解的微分方程也是十分有限的, 特别是高阶方程和偏微分方程（组）.这就要求我们必须研究微分方程（组）的解法：解析解法和数值解法.一.相关函数、命令及简介1.在Matlab中，用大写字母D表示导数，Dy表示y关于自变量的一阶导数, D2y表示y关于自变量的二阶导数，依此类推.函数dsolve用来解决常微分方程（组）的求解问题，调用格式为：X=dsolve（＜eqnl，,，eqn2函数dsolve用来解符号常微分方程、方程组，如果没有初始条件，则求出通解，如果有初始条件，则求出特解.注意，系统缺省的自变量为t2.函数dsolve求解的是常微分方程的精确解法，也称为常微分方程的符号解. 但是，有大量的常微分方程虽然从理论上讲，其解是存在的，但我们却无法求出其解析解，此时，我们需要寻求方程的数值解，在求常微分方程数值解方面，MATLAB具有丰富的函数，我们将其统称为solver,其一般格式为：［T,Y］=solver（odefun,tspan,yO）说明：（1 ）solver 为命令ode45、ode23、odel 13、odel5s、ode23s、ode23t、ode23tb、odel5i 之一.（2）odefun是显示微分方程），=f （t,y）在积分区间tspan =［心心］上从心到“用初始条件儿求解.（3）如果要获得微分方程问题在其他指定时间点bG©…心上的解，则令（span = 『“，•••『/■］（要单调的）.（4）因为没有一种算法可以有效的解决所有的ODE问题，为此，Matlab提供T多种求解器solver,对于不同的ODE问题，采用不同的solver.程（组）的初值问题的解的Matlab常用程序，其中：ode23采用龙格-库塔2阶算法，用3阶公式作误差估计来调节步长，具有低等的精度.。
matlab求解微分方程sir模型摘要：一、微分方程SIR模型的简介1.SIR模型的基本概念2.SIR模型的数学表达式二、MATLAB求解SIR模型的方法1.MATLAB的ODE45函数2.MATLAB的ode23tb函数三、实例演示求解SIR模型1.初始条件的设定2.参数设置3.求解过程4.结果分析正文：一、微分方程SIR模型的简介SIR模型是描述传染病在人群中传播的一个数学模型，它由三个方程式组成，分别是S(t)表示易感者人数，I(t)表示感染者人数，R(t)表示康复者人数。
SIR模型的基本概念是易感者与感染者之间的接触会导致疾病传播，感染者康复后会获得免疫力。
SIR模型的数学表达式如下：dS/dt = -βSIdI/dt = βSI - γIdR/dt = γI其中，β表示每个感染者在单位时间内接触并传染给易感者的概率，γ表示感染者康复或死亡的概率。
二、MATLAB求解SIR模型的方法MATLAB提供了求解常微分方程的ODE45和ode23tb函数。
其中，ODE45是一种具有四阶和五阶龙格库塔法的数值积分方法，适用于求解非刚性问题；ode23tb是一种具有二阶和三阶龙格库塔法的数值积分方法，适用于求解刚性问题。
三、实例演示求解SIR模型下面我们通过一个实例来演示如何用MATLAB求解SIR模型。
1.初始条件的设定我们设定初始条件为：S(0)=1000, I(0)=1, R(0)=0。
2.参数设置我们设定参数为：β=0.4, γ=0.3。
3.求解过程我们使用MATLAB的ODE45函数来求解SIR模型，输入参数为：fun，即SIR模型的数学表达式，以及初始条件和参数。
[~, params] = ode45(@(t, S) [-0.4*S*I; 0.4*S*I - 0.3*I; 0.3*I], [0, 20], [1000, 1, 0]);4.结果分析我们可以通过params返回的结果来分析SIR模型的解。
