微分方程几种求解方法
第五章 控制系统仿真
§5.2 微分方程求解方法
以一个自由振动系统实例为例进行讨论。 如下图1所示弹簧-阻尼系统，参数如下： M=5 kg, b=1 N.s/m, k=2 N/m, F=1N
F
图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(*4.0)2(*2.02.0)2()2()2()1(x x x x x x x && 下面就可以进行程序的编制。
%写出函数文件myfun1.m
function xdot=myfun1(t,x)
xdot=[x(2);0.2-0.2*x(2)-0.4*x(1)];
% 主程序wffc1.m
t=[0 30];
x0=[0;0];
[tt,yy]=ode45(@myfun1,t,x0);
plot(tt,yy(:,1),':b',tt,yy(:,2),'-r') hold on
plot(tt,0.2-0.2*yy(:,2)-0.4*yy(:,1),'-k') legend('位移','速度',’加速度’)
title('微分方程的解 x(t)') xlabel(‘时间t ’)
ylabel(‘输出’)
二、方法2：
F kx x b x m =++?
??
251)()()(2++==s s s F s X s G
%用传递函数编程求解ksys1.m num=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)'; figure(6)
bode(num,den,w)
grid on
[m p]=bode(num,den,w); [gm,pm,wpc,wgc]=margin(sys) figure(7)
margin(sys)
figure(8)
nyquist(sys)
figure(9)
nichols(sys)
方法3：Simulink 中传递函数模块求解
方法4：Simulink 中积分模块求解 F kx x b x m =++???
125=++???x x x
x x x 4.02.02.0--=?
?? x''x'x u(t)x_t
To Workspace Scope 1s Int21s Int1
0.2
Gs
0.4
G20.2
G1
各种类型的微分方程及其相应解法教程文件
各种类型的微分方程及其相应解法 专业班级：交土01班 姓名：高云 学号：1201110102 微分方程的类型有很多种，解题时先判断微分方程是哪种类型，可以帮助我们更快解题，所以我们有必要归纳整理一下各类型（主要是一阶和二阶）的微分方程及其相应解法。 一、一阶微分方程的解法 1.可分离变量的方程 dx x f dy y g )()(=,或)()(y g x f dx dy = 其特点是可以把变量x 和y 只分别在等式的两边，解法关键是把变量分离后两边积分。 例1.求微分方程ydy dx y xydy dx +=+2的通解. 解 先合并dx 及dy 的各项,得dx y dy x y )1()1(2-=- 设,01,012≠-≠-x y 分离变量得 dx x dy y y 1112-=- 两端积分??-=-dx x dy y y 1112得 ||ln |1|ln |1|ln 2 112C x y +-=- 于是 2212)1(1-±=-x C y 记,21C C ±=则得到题设方程的通解 .)1(122-=-x C y 2.齐次方程 （1）)(x y f dx dy = （2） )(c by ax f dx dy ++=（a ，b 均不等于0） 例2求解微分方程.2222xy y dy y xy x dx -=+- 解 原方程变形为=+--=2222y xy x xy y dx dy ,1222?? ? ??+--??? ??x y x y x y x y 令,x y u =则,dx du x u dx dy +=方程化为,1222u u u u dx du x u +--=+ 分离变量得?? ????-+--??? ??--112212121u u u u ,x dx du = 两边积分得 ,ln ln ln 2 1)2ln(23)1ln(C x u u u +=----
081数值计算方法—常微分方程(组)
科学计算—理论、方法 及其基于MATLAB 的程序实现与分析 微分方程（组）数值解法 §1 常微分方程初值问题的数值解法 微分方程（组）是科学研究和工程应用中最常用的数学模型之一。如揭示质点运动规律的Newton 第二定律： ()()()?????'='==0 00022x t x x t x t F dt x d m （１） 和刻画回路电流或电压变化规律的基尔霍夫回路定律等，但是，只有一些简单的和特殊的常微分方程及常微分方程组，可以求得用公式给出的所谓“解析解”或“公式解”，如一阶线性微分方程的初值问题： () ()0 0y y t f ay dt dy =+= （２） 的解为： ()()()τττd f e y e t y t t a at ?-+=00 （３） 但是，绝大多数在实际中遇到的常微分方程和常微分方程组得不到“解析解”，因此，基于如下的事实：
１、绝大多数的常微分方程和常微分方程组得不到（有限形式的）解析解； ２、实际应用中往往只需要知道常微分方程（组）的解在（人们所关心的）某些点处的函数值（可以是满足一定精度要求的近似值）； 如果只需要常微分方程（组）的解在某些点处的函数值，则没有必要非得通过求得公式解，然后再计算出函数值不可，事实上，我们可以采用下面将介绍的常微分方程（组）的初值问题的数值解法，就可以达到这一目的。 一般的一阶常微分方程（组）的初值问题是指如下的一阶常微分方程（组）的定解问题： ()()0 00,y t y t t t y t F dt dy f =≤≤= （７） 其中 ()()()()???? ?? ? ??=t y t y t y t y n 21 （８） ()()()()???? ?? ? ??=y t f y t f y t f y t F n ,,,,21 （９） 常微分方程（组）的初值问题通常是对一动态过程（动态系统、动力系统）演化规律的描述，求解常微分方程（组）的初值问题就是要了解和掌握动态过程演化规律。 §１.1 常微分方程（组）的Cauch 问题数值解法概论
各类微分方程的解法大全
各类微分方程的解法 1.可分离变量的微分方程解法 一般形式:g(y)dy=f(x)dx 直接解得∫g(y)dy=∫f(x)dx 设g(y)及f(x)的原函数依次为G(y)及F(x),则G(y)=F(x)+C为微分方程的隐式通解 2.齐次方程解法 一般形式:dy/dx=φ(y/x) 令u=y/x则y=xu,dy/dx=u+xdu/dx,所以u+xdu/dx=φ(u),即du/［φ(u)-u］=dx/x 两端积分,得∫du/［φ(u)-u］=∫dx/x 最后用y/x代替u,便得所给齐次方程的通解 3.一阶线性微分方程解法 一般形式:dy/dx+P(x)y=Q(x) 先令Q(x)=0则dy/dx+P(x)y=0解得y=Ce－ ∫P(x)dx,再令y=u e－∫P(x)dx代入原方程解得u=∫Q(x) e∫P(x)dx dx+C,所以y=e－∫P(x)dx［∫Q(x)e∫P(x)dx dx+C］ 即y=Ce－∫P(x)dx +e－ ∫P(x)dx∫Q(x)e∫P(x)dx dx为一阶线性微分方程的通解 4.可降阶的高阶微分方程解法 ①y(n)=f(x)型的微分方程 y(n)=f(x) y(n-1)= ∫f(x)dx+C1 y(n-2)= ∫［∫f(x)dx+C1］dx+C2 依次类推,接连积分n次,便得方程y(n)=f(x)的含有n个任意常数的通解②y”=f(x,y’) 型的微分方程 令y’=p则y”=p’,所以p’=f(x,p),再求解得p=φ(x,C1) 即dy/dx=φ(x,C1),所以y=∫φ(x,C1)dx+C2 ③y”=f(y,y’) 型的微分方程
令y’=p则y”=pdp/dy,所以pdp/dy=f(y,p),再求解得p=φ(y,C1) 即dy/dx=φ(y,C1),即dy/φ(y,C1)=dx,所以∫dy/φ(y,C1)=x+C2 5.二阶常系数齐次线性微分方程解法 一般形式:y”+py’+qy=0,特征方程r2+pr+q=0 6.二阶常系数非齐次线性微分方程解法 一般形式: y”+py’+qy=f(x) 先求y”+py’+qy=0的通解y0(x),再求y”+py’+qy=f(x)的一个特解y*(x) 则y(x)=y0(x)+y*(x)即为微分方程y”+py’+qy=f(x)的通解 求y”+py’+qy=f(x)特解的方法: ①f(x)=P m(x)eλx型 令y*=x k Q m(x)eλx［k按λ不是特征方程的根,是特征方程的单根或特征方程的重根依次取0,1或2］再代入原方程,确定Q m(x)的m+1个系数 ②f(x)=eλx［Pｌ(x)cosωx+P n(x)sinωx］型 令y*=x k eλx［Q m(x)cosωx+R m(x)sinωx］［m=max﹛ｌ,n﹜,k按λ+iω不是特征方程的根或是特征方程的单根依次取0或1］再代入原方程,分别确定Q m(x)和R m(x)的m+1个系数
微分方程几种求解方法
第五章 控制系统仿真 §5.2 微分方程求解方法 以一个自由振动系统实例为例进行讨论。 如下图1所示弹簧-阻尼系统，参数如下： M=5 kg, b=1 N.s/m, k=2 N/m, F=1N F 图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(*4.0)2(*2.02.0)2()2()2()1(x x x x x x x && 下面就可以进行程序的编制。 %写出函数文件myfun1.m function xdot=myfun1(t,x) xdot=[x(2);0.2-0.2*x(2)-0.4*x(1)]; % 主程序wffc1.m t=[0 30]; x0=[0;0]; [tt,yy]=ode45(@myfun1,t,x0); plot(tt,yy(:,1),':b',tt,yy(:,2),'-r') hold on plot(tt,0.2-0.2*yy(:,2)-0.4*yy(:,1),'-k') legend('位移','速度',’加速度’)
各类微分方程的解法大全
创作编号：BG7531400019813488897SX 创作者：别如克* 各类微分方程的解法 1.可分离变量的微分方程解法 一般形式:g(y)dy=f(x)dx 直接解得∫g(y)dy=∫f(x)dx 设g(y)及f(x)的原函数依次为G(y)及F(x),则G(y)=F(x)+C为微分方程的隐 式通解 2.齐次方程解法 一般形式:dy/dx=φ(y/x) 令u=y/x则y=xu,dy/dx=u+xdu/dx,所以u+xdu/dx=φ(u),即du/［φ(u)-u］ =dx/x两端积分,得∫du/［φ(u)-u］=∫dx/x 最后用y/x代替u,便得所给齐次方程的通解 3.一阶线性微分方程解法 一般形式:dy/dx+P(x)y=Q(x) 先令Q(x)=0则dy/dx+P(x)y=0解得y=Ce－ ∫P(x)dx,再令y=u e－∫P(x)dx代入原方程 解得u=∫Q(x) e∫P(x)dx dx+C,所以y=e－∫P(x)dx［∫Q(x)e∫P(x)dx dx+C］ 即y=Ce－∫P(x)dx +e－ ∫P(x)dx∫Q(x)e∫P(x)dx dx为一阶线性微分方程的通解 4.可降阶的高阶微分方程解法 ①y(n)=f(x)型的微分方程 y(n)=f(x) y(n-1)= ∫f(x)dx+C1
y(n-2)= ∫［∫f(x)dx+C1］dx+C2 依次类推,接连积分n次,便得方程y(n)=f(x)的含有n个任意常数的通解②y”=f(x,y’) 型的微分方程 令y’=p则y”=p’,所以p’=f(x,p),再求解得p=φ(x,C1) 即dy/dx=φ(x,C1),所以y=∫φ(x,C1)dx+C2 ③y”=f(y,y’) 型的微分方程 令y’=p则y”=pdp/dy,所以pdp/dy=f(y,p),再求解得p=φ(y,C1) 即dy/dx=φ(y,C1),即dy/φ(y,C1)=dx,所以∫dy/φ(y,C1)=x+C2 5.二阶常系数齐次线性微分方程解法 一般形式:y”+py’+qy=0,特征方程r2+pr+q=0 6.二阶常系数非齐次线性微分方程解法 一般形式: y”+py’+qy=f(x) 先求y”+py’+qy=0的通解y0(x),再求y”+py’+qy=f(x)的一个特解y*(x) 则y(x)=y0(x)+y*(x)即为微分方程y”+py’+qy=f(x)的通解 求y”+py’+qy=f(x)特解的方法: ①f(x)=P m(x)eλx型 令y*=x k Q m(x)eλx［k按λ不是特征方程的根,是特征方程的单根或特征方程的重根依次取0,1或2］再代入原方程,确定Q m(x)的m+1个系数 ②f(x)=eλx［Pｌ(x)cosωx+P n(x)sinωx］型
各种类型的微分方程及其相应解法
各种类型的微分方程及其相应解法 专业班级：交土01班 姓名：高云 学号：1201110102 微分方程的类型有很多种，解题时先判断微分方程是哪种类型，可以帮助我们更快解题，所以我们有必要归纳整理一下各类型（主要是一阶和二阶）的微分方程及其相应解法。 一、一阶微分方程的解法 1.可分离变量的方程 dx x f dy y g )()(=,或)()(y g x f dx dy = 其特点是可以把变量x 和y 只分别在等式的两边，解法关键是把变量分离后两边积分。 例1.求微分方程ydy dx y xydy dx +=+2的通解. 解 先合并dx 及dy 的各项,得dx y dy x y )1()1(2-=- 设,01,012≠-≠-x y 分离变量得 dx x dy y y 1112-=- 两端积分??-=-dx x dy y y 1112得 ||ln |1|ln |1|ln 2 112C x y +-=- 于是 2212)1(1-±=-x C y 记,21C C ±=则得到题设方程的通解 .)1(122-=-x C y 2.齐次方程 （1）)(x y f dx dy = （2） )(c by ax f dx dy ++=（a ，b 均不等于0） 例2求解微分方程.2222xy y dy y xy x dx -=+- 解 原方程变形为=+--=2222y xy x xy y dx dy ,1222?? ? ??+--??? ??x y x y x y x y 令,x y u =则,dx du x u dx dy +=方程化为,1222u u u u dx du x u +--=+ 分离变量得?? ????-+--??? ??--112212121u u u u ,x dx du = 两边积分得 ,ln ln ln 2 1)2ln(23)1ln(C x u u u +=----
计算方法_微分方程数值解
120 第6章 常微分方程初值问题数值解法 6.1 问题的描述和基本概念 1、常微分方程初值问题 ● 一般形式 0(,)()y f x y y a y '=??=? 式中(,)f x y 已知，0()y a y =称为初值条件. ● 初值问题的数值方法和数值解 求函数()y y x =在若干离散点k x 上的近似值 (0,1,)k y k = 的方法称为初值问题的数值方法，而 称(0,1,)k y k = 为初值问题的数值解.
121 2. 建立数值解法的思想与方法 用离散化方法将初值问题化为差分方程, 然后再求解. 设节点为 011n n a x x x x +=<<<<< 距离1k k k h x x +=-称为步长. 求数值解一般是从0y 开使逐次顺序求出12,,y y . 初值问题的解法有单步法和多步法两种： ● 单步法：计算1k y +时只用到k y 一个值； ● 多步法：计算1k y +时要用1,,,k k k l y y y -- 多个值。 数值解法还有显格式和隐格式之分。
122 微分方程离散化方法主要有 数值微分法，数值积分法和Taylor 展开法 1) 数值微分法 由'()(,())k k k y x f x y x =，用数值微分的2点前差公式代替'()k y x ，得近似离散化方程 记1k k h x x +=-，做k k ，“”，得差分方程 1(,)k k k k y y f x y h +-= 即 1(,)k k k k y y hf x y +=+ （Euler 公式） 由初值条件0()y y a =及Euler 公式可求出数值解 12,,,,n y y y .Euler 公式是显式单步法.
常微分方程解题方法总结.doc
常微分方程解题方法总结 来源：文都教育 复习过半， 课本上的知识点相信大部分考生已经学习过一遍 . 接下来， 如何将零散的知 识点有机地结合起来， 而不容易遗忘是大多数考生面临的问题 . 为了加强记忆， 使知识自成 体系，建议将知识点进行分类系统总结 . 著名数学家华罗庚的读书方法值得借鉴， 他强调读 书要“由薄到厚、由厚到薄”，对同学们的复习尤为重要 . 以常微分方程为例， 本部分内容涉及可分离变量、 一阶齐次、 一阶非齐次、 全微分方程、 高阶线性微分方程等内容， 在看完这部分内容会发现要掌握的解题方法太多， 遇到具体的题 目不知该如何下手， 这种情况往往是因为没有很好地总结和归纳解题方法 . 下面以表格的形 式将常微分方程中的解题方法加以总结，一目了然，便于记忆和查询 . 常微分方程 通解公式或解法 ( 名称、形式 ) 当 g( y) 0 时，得到 dy f (x)dx ， g( y) 可分离变量的方程 dy f ( x) g( y) 两边积分即可得到结果； dx 当 g( 0 ) 0 时，则 y( x) 0 也是方程的 解 . 解法：令 u y xdu udx ，代入 ，则 dy 齐次微分方程 dy g( y ) x dx x u g (u) 化为可分离变量方程 得到 x du dx 一 阶 线 性 微 分 方 程 P ( x)dx P ( x) dx dy Q(x) y ( e Q( x)dx C )e P( x) y dx
伯努利方程 解法：令 u y1 n，有 du (1 n) y n dy ， dy P( x) y Q( x) y n（n≠0,1）代入得到du (1 n) P(x)u (1 n)Q(x) dx dx 求解特征方程：2 pq 三种情况： 二阶常系数齐次线性微分方程 y p x y q x y0 二阶常系数非齐次线性微分方程 y p x y q x y f ( x) (1)两个不等实根：1, 2 通解： y c1 e 1x c2 e 2x (2) 两个相等实根：1 2 通解： y c1 c2 x e x (3) 一对共轭复根：i , 通解： y e x c1 cos x c2 sin x 通解为 y p x y q x y 0 的通解与 y p x y q x y f ( x) 的特解之和. 常见的 f (x) 有两种情况： x （ 1）f ( x)e P m ( x) 若不是特征方程的根，令特解 y Q m ( x)e x；若是特征方程的单根，令特 解 y xQ m ( x)e x；若是特征方程的重根， 令特解 y*x2Q m (x)e x； （2）f (x) e x[ P m ( x) cos x p n ( x)sin x]
郑州大学研究生课程数值分析复习---第八章 常微分方程数值解法
郑州大学研究生课程（2012-2013学年第一学期）数值分析 Numerical Analysis 习题课 第八章常微分方程数值解法
待求解的问题：一阶常微分方程的初值问题/* Initial-Value Problem */: ?????=∈=0 )(] ,[),(y a y b a x y x f dx dy 解的存在唯一性（“常微分方程”理论）：只要f (x , y ) 在[a , b ] ×R 1 上连续，且关于y 满足Lipschitz 条件，即存在与x , y 无关的常数L 使 对任意定义在[a , b ] 上的y 1(x ) 和y 2(x ) 都成立，则上述IVP 存在唯一解。 1212|(,)(,)||| f x y f x y L y y ?≤?一、要点回顾
§8.2 欧拉(Euler)法 通常取(常数),则Euler 法的计算格式 h h x x i i i ==?+1?? ?=+=+) (),(001x y y y x hf y y i i i i i =0,1,…,n ( 8.2 )
§8.2 欧拉(Euler)法(1) 用差商近似导数 )) (,()()()()(1n n n n n n x y x hf x y x y h x y x y +=′+≈+?? ?=+=+) (),(01a y y y x hf y y n n n n 差分方程初值问题向前Euler 方法h x y x y x y n n n ) ()()(1?≈ ′+)) (,() ()(1n n n n x y x f h x y x y ≈?+))(,()(n n n x y x f x y =′
数值分析_第五章_常微分方程数值解法
图５畅２ 令珔h ＝h λ，则y n ＋１＝１＋珔 h ＋１２珔h ２ ＋１６珔h ３＋１２４ 珔 h ４y n ．由此可知，绝对稳定性区域在珔h ＝h λ复平面上满足 ｜１＋珔 h ＋１２珔h ２＋１６珔h ３＋１２４珔h ４ ｜≤１的区域，也就是由曲线 １＋珔h ＋ １２珔h ２＋１６珔h ３＋１２４ 珔h ４＝e i θ 所围成的区域．如图５畅２所示． 例22 用Euler 法求解 y ′＝－５y ＋x ，y （x ０）＝y ０， x ０≤x ≤X ． 从绝对稳定性考虑，对步长h 有何限制？ 解 对于模型方程y ′＝λy （λ＜０为实数）这里λ＝抄f 抄y ＝－５．由 ｜１＋h λ｜＝｜１－５h ｜＜１ 得到对h 的限制为：０＜h ＜０畅４． 四、习题 １畅取步长h ＝０畅２，用Euler 法解初值问题 y ′＝－y －x y ２ ， y （０）＝ １． （０≤x ≤０畅６）， ２畅用梯形公式解初值问题 y ′＝８－３y ， （１≤x ≤２），
取步长h＝０畅２，小数点后至少保留５位． ３畅用改进的Euler公式计算初值问题 y′＝１x y－１x y２， y（１）＝０畅５， １＜x＜１畅５， 取步长h＝０畅１，并与精确解y（x）＝ x １＋x比较． ４畅写出用梯形格式的迭代算法求解初值问题 y′＋y＝０， y（０）＝１ 的计算公式，取步长h＝０畅１，并求y（０畅２）的近似值，要求迭代误差不超过１０－５． ５畅写出用四阶经典Runge唱Kutta法求解初值问题 y′＝８－３y， y（０）＝２ 的计算公式，取步长h＝０畅２，并计算y（０畅４）的近似值，小数点后至少保留４位． ６畅证明公式 y n＋１＝y n＋h９（２K１＋３K２＋４K３）． K１＝f（x n，y n）， K２＝f x n＋h２，y n＋h２K１， K３＝f x n＋３４h，y n＋３４h K２， 至少是三阶方法． ７畅试构造形如 y n＋１＝α（y n＋y n－１）＋h（β０f n＋β１f n－１）
各种类型的微分方程及其相应解法教学文案
各种类型的微分方程及其相应解法
各种类型的微分方程及其相应解法 专业班级：交土01班 姓名：高云 学号：1201110102 微分方程的类型有很多种，解题时先判断微分方程是哪种类型，可以帮助我们更快解题，所以我们有必要归纳整理一下各类型（主要是一阶和二阶）的微分方程及其相应解法。 一、一阶微分方程的解法 1.可分离变量的方程 dx x f dy y g )()(=,或)()(y g x f dx dy = 其特点是可以把变量x 和y 只分别在等式的两边，解法关键是把变量分离后两边积分。 例1.求微分方程ydy dx y xydy dx +=+2的通解. 解 先合并dx 及dy 的各项,得dx y dy x y )1()1(2-=- 设,01,012≠-≠-x y 分离变量得 dx x dy y y 1112-=- 两端积分??-=-dx x dy y y 1112得 ||ln |1|ln |1|ln 2 112C x y +-=- 于是 2212)1(1-±=-x C y 记,21C C ±=则得到题设方程的通解 .)1(122-=-x C y 2.齐次方程 （1） )(x y f dx dy = （2） )(c by ax f dx dy ++=（a ，b 均不等于0） 例2求解微分方程.2222xy y dy y xy x dx -=+- 解 原方程变形为=+--=2222y xy x xy y dx dy ,1222?? ? ??+--??? ??x y x y x y x y
令,x y u =则,dx du x u dx dy +=方程化为,1222u u u u dx du x u +--=+ 分离变量得??? ???-+--??? ??--1122 121 21 u u u u ,x dx du = 两边积分得 ,ln ln ln 21 )2ln(23 )1ln(C x u u u +=---- 整理得 .)2(1 2/3Cx u u u =-- 所求微分方程的解为 .)2()(32x y Cy x y -=- 3.一阶线性微分方程 ?+??==+-])([),()()()(C dx e x Q e y x Q y x p dx dy dx x p dx x p 其通解为 例3. x y dx dy x sin 2=+， ππ1 )(=y ； 解 将方程改写为 x x y x dx dy sin 2=+， 这里x x p 2)(=，x x x q sin )(=，故由求解公式得 )sin (1 sin 22 2 ??+=??? ????+?=-xdx x C x dx e x x C e y dx x dx x 22sin cos x x x x x C +-=. 由初值条件ππ1 )(=y ，得0=C . 所以初值问题的解为 2cos sin x x x x y -= 例4.设非负函数()f x 具有一阶导数，且满足1 200()()()x f x f t dt t f t dt =+??，求 函数()f x ． 解：设120()A t f t dt =?，则0()()x f x f t dt A =+?，两边对x 求导，得 ()()()x f x f x f x Ce '=?=，由已知(0)()x f A C A f x Ae =?=?= 又 11222004 ()()1t A t f t dt t Ae dt A e ==?=+??，则 24 ()1x f x e e =+
解微分方程方法
MATLAB解微分方程(2011-07-15 17:35:25) 转载▼ 分类：matlab学习标签： 教育
先说明一下最常用的ode45调用方式，和相应的函数文件定义格式。 [t,x]=ode45(odefun,tspan,x0); 其中，Fun就是导函数，tspan为求解的时间区间（或时间序列，如果采用时间序列，则必须单调），x0为初值。 这时，函数文件可以采用如下方式定义 function dx=odefun(t,x) 对于上面的小例子，可以用如下的程序求解。
2.终值问题 tspan可以是递增序列，也可以为递减序列，若为递减则可求解终值问题。 [t,x]=ode45(@zhongzhiode,[3,0],[1;0;2]);plot(t,x) function dx=zhongzhiode(t,x) dx=[2*x(2)^2-2; -x(1)+2*x(2)*x(3)-1; -2*x(2)+2*x(3)^2-4]; 结果如下 3.odeset options = odeset('name1',value1,'name2',value2,...) [t,x]=solver(@fun,tspan,x0,options) 通过odeset设置options 第一，通过求解选项的设置可以改善求解精度，使得原本可能不收敛的问题收敛。options=odeset('RelTol',1e-10);
第二，求解形如M(t,x)x'=f(t,x)的方程。 例如，方程 x'=-0.2x+yz+0.3xy y'=2xy-5yz-2y^2 x+y+z-2=0 可以变形为 [1 0 0][x'] [-0.2x+yz+0.3xy] [0 1 0][y']=[2xy-5yz-2y^2 ] [0 0 1][z'] [x+y+z-2 ] 这样就可以用如下的代码求解该方程 function mydae M=[1 0 0;0 1 0;0 0 0]; options=odeset('Mass',M); x0=[1.6,0.3,0.1]; [t,x]=ode15s(@daedot,[0,1.5],x0,options);plot(t,x) function dx=daedot(t,x) dx=[ -0.2*x(1)+x(2)*x(3)+0.3*x(1)*x(2); 2*x(1)*x(2)-5*x(2)*x(3)-2*x(2)*x(2); x(1)+x(2)+x(3)-2]; 4.带附加参数的ode45
偏微分方程数值解法
《偏微分方程数值解法》 课程设计 题目: 六点对称差分格式解热传导方程的初边 值问题 姓名: 王晓霜 学院: 理学院 专业: 信息与计算科学 班级: ０91１01２ 学号: ０９１10１218 指导老师：翟方曼 ２０12年12月14
日 一、题目 用六点对称差分格式计算如下热传导方程的初边值问题 222122,01,01(,0),01 (0,),(1,),01x t t u u x t t x u x e x u t e u t e t +???=<<<≤?????=≤≤??==≤≤??? 已知其精确解为 2(,)x t u x t e += 二、理论 1.考虑的问题 考虑一维模型热传导方程 （1.１) )(22x f x u a t u +??=??,T t ≤<0 其中a 为常数。)(x f 是给定的连续函数。（1.１)的定解问题分两类: 第一,初值问题（Cauch ｙ 问题)：求足够光滑的函数()t x u ,,满足方程（１.1）和初始条件： (1.２) ()()x x u ?=0,， ∞<<∞-x 第二,初边值问题（也称混合问题）:求足够光滑的函数()t x u ,，满足方程(１.１)和初始条件： ()13.1 ()()x x u ?=0,， l x l <<- 及边值条件 ()23.1 ()()0,,0==t l u t u ， T t ≤≤0 假定()x f 和()x ?在相应的区域光滑,并且于()0,0，()0,l 两点满足相容条件,则上述问题有唯一的充分光滑的解。 现在考虑边值问题（１.1),（1.３)的差分逼近 取 N l h = 为空间步长,M T =τ为时间步长,其中N ,M 是自然数,
求解偏微分方程三种数值方法
数值模拟偏微分方程的三种方法介绍 （有限差分方法、有限元方法、有限体积方法） I．三者简介 有限差分方法（Finite Difference Methods）是数值模拟偏微分方程最早采用的方法，至今仍被广泛使用。该方法包括区域剖分和差商代替导数两个步骤。首先将求解区域划分为差分网格，用有限个网格节点代替连续的求解区域。其次，利用Taylor级数展开等方法将偏微分方程中的导数项在网格节点上用函数值的差商代替进行离散，从而建立以网格节点上的值为未知量的代数方程组。该方法是一种直接将微分问题变为代数问题的近似数值解法，数学概念直观，表达简单，是发展较早且十分成熟的数值方法。 差商代替导数后的格式称为有限差分格式，从格式的精度来考虑，有一阶格式、二阶格式和高阶格式。从差分的空间离散形式来考虑，有中心格式和迎风格式。对于瞬态方程，考虑时间方向的离散，有显格式、隐格式、交替显隐格式等。目前常见的差分格式，主要是以上几种格式的组合，不同的组合构成不同的差分格式。差分方法主要适用于结构网格，网格的大小一般根据问题模型和Courant 稳定条件来决定。 有限元方法（Finite Element Methods）的基础是虚位移原理和分片多项式插值。该方法的构造过程包括以下三个步骤。首先，利用虚位移原理得到偏微分方程的弱形式，将计算区域划分为有限个互不重叠的单元（三角形、四边形、四面体、六面体等），在每个单元上选择合适的节点作为求解函数的插值点，将偏微分方程中的变量改写成由各变量或其导数的节点值与所选用的分片插值基函数组成的线性表达式，得到微分方程的离散形式。利用插值函数的局部支集性质及数值积分可以得到未知量的代数方程组。 有限元方法有较完善的理论基础，具有求解区域灵活（复杂区域）、单元类型灵活（适于结构网格和非结构网格）、程序代码通用（数值模拟软件多数基于有限元方法）等特点。有限元方法最早应用于结构力学，随着计算机的发展已经渗透到计算物理、流体力学与电磁学等各个数值模拟领域。
(完整word版)偏微分方程数值解法答案
1. 课本2p 有证明 2. 课本812,p p 有说明 3. 课本1520,p p 有说明 4. Rit2法，设n u 是u 的n 维子空间，12,...n ???是n u 的一组基底，n u 中的任一元素n u 可 表为1n n i i i u c ?==∑ ，则,11 11()(,)(,)(,)(,)22j n n n n n n i j i j j i j j J u a u u f u a c c c f ???=== -=-∑∑是12,...n c c c 的二次函数，(,)(,)i j j i a a ????=，令 () 0n j J u c ?=?，从而得到12,...n c c c 满足1 (,)(,),1,2...n i j i j i a c f j n ?? ?===∑，通过解线性方程组，求的i c ，代入1 n n i i i u c ?==∑， 从而得到近似解n u 的过程称为Rit2法 简而言之，Rit2法：为得到偏微分方程的有穷维解，构造了一个近似解，1 n n i i i u c ? == ∑， 利用,11 11()(,)(,)(,)(,)22j n n n n n n i j i j j i j j J u a u u f u a c c c f ???===-=-∑∑确定i c ，求得近似解n u 的过程 Galerkin 法：为求得1 n n i i i u c ? == ∑形式的近似解，在系数i c 使n u 关于n V u ∈，满足(,)(,) n a u V f V =，对任 意 n V u ∈或（取 ,1j V j n ?=≤≤） 1 (,)(,),1,2...n i j i j i a c f j n ???===∑的情况下确定i c ，从而得到近似解1 n n i i i u c ?==∑的过程称 Galerkin 法为 Rit2-Galerkin 法方程： 1 (,)(,)n i j i j i a c f ???==∑ 5. 有限元法：将偏微分方程转化为变分形式，选定单元的形状，对求解域作剖分，进而构 造基函数或单元形状函数，形成有限元空间，将偏微分方程转化成了有限元方程，利用 有效的有限元方程的解法，给出偏微分方程近似解的过程称为有限元法。 6. 解：对求解区间进行网格剖分，节点01......i n a x x x x b =<<<<=得到相邻节点1,i i x x -
常微分方程数值解
第四章常微分方程数值解 [课时安排]6学时 [教学课型]理论课 [教学目的和要求] 了解常微分方程初值问题数值解法的一些基本概念，如单步法和多步法，显式和隐式，方法的阶数，整体截断误差和局部截断误差的区别和关系等；掌握一阶常微分方程初值问题的一些常用的数值计算方法，例如欧拉（Euler）方法、改进的欧拉方法、龙贝－库塔（Runge-Kutta）方法、阿达姆斯（Adams）方法等，要注意各方法的特点及有关的理论分析；掌握构造常微分方程数值解的数值积分的构造方法和泰勒展开的构造方法的基本思想，并能具体应用它们导出一些常用的数值计算公式及评估截断误差；熟练掌握龙格－库塔（Ｒ－Ｋ）方法的基本思想，公式的推导，Ｒ－Ｋ公式中系数的确定，特别是能应用“标准四阶Ｒ－Ｋ公式”解题；掌握数值方法的收敛性和稳定性的概念，并能确定给定方法的绝对稳定性区域。[教学重点与难点] 重点：欧拉方法，改进的欧拉方法，龙贝－库塔方法。 难点：R—K方法，预估-校正公式。 [教学内容与过程] 4.1 引言 本章讨论常微分方程初值问题 (4.1.1) 的数值解法，这也是科学与工程计算经常遇到的问题，由于只有很特殊的方程能用解析方法求解，而用计算机求解常微分方程的初值问题都要采用数值方法.通常我们假定(4.1.1)中 f(x,y)对y满足Lipschitz条件，即存在常数L＞0，使对，有 (4.1.2) 则初值问题(4.1.1)的解存在唯一. 假定(4.1.1)的精确解为，求它的数值解就是要在区间上的一组离散点 上求的近似.通常取 ，h称为步长，求(4.1.1)的数值解是按节点的顺序逐步 推进求得.首先，要对方程做离散逼近，求出数值解的公式，再研究公式的局部截
计算方法6_微分方程
习题6 6.1 试用三种方法导出线性二步方法 122+++=n n n hf y y 6.2 用Taylor 展开法求三步四阶方法类，并确定三步四阶显式方法. 6.3 形如 ∑=++=k i k n k j n j f h y 0βα 的k 阶方法称为Gear 方法，试确定一个三步Gear 方法，并给出其截断误差主项。 6.4 试用显式Euler 法及改进的Euler 法 )],(),([2 11n n n n n n n hf y t f y t f h y y +++=++ 6.5 给出线性多步法 ])13()3[(4 )1(212n n n n n f f h y y y +++=--++++αααα 为零稳定的条件，并证明该方法为零稳定时是二阶收敛的. 6.6 给出题(6.5)题中1=α时的公式的绝对稳定域. 6.7 指出Heun 方法 0 0 0 0 1/3 1/3 0 0 2/3 0 2/3 0 1/4 0 3/4 的相容阶，并给出由该方法以步长h 计算初值问题(6.45)的步骤. 6.8 试述刚性问题的基本特征，并给出s 级Runge-Kutta 方法为A -稳定的条件. 6.9 设有???=='00 )(),(y x y y x f y ，试构造形如 )()(11011--++++=n n n n n f f h y y y ββα 的二阶方法，并推导其局部截断误差首项。
6.10设有常微分方程初值问题???=='00 )(),(y x y y x f y 的单步法)],(2),([3 111+++++=n n n n n n y x f y x f h y y ，证明该方法是无条件稳定的。
微分方程的普通解法
微分方程的解法 1. 微分方程的基本概念 常微分方程, 微分方程的阶, 微分方程的解、通解, 初始条件和特解的概念。 2. 一阶微分方程 掌握变量可分离的微分方程及一阶线性微分方程的解法。 会解齐次方程和贝努利方程并从中领会变量代换求解微分方程的思想。 3. 可降阶的高阶方程 会 )()(x f y n =,),(y x f y '='',),(y y f y '=''的降阶解法。 4. 二阶线性微分方程 理解二阶线性微分方程解的结构。 掌握二阶常系数线性齐次微分方程的解法，了解高阶常系数线性齐次微分方程的解法。 会求非齐次项形如 x m e x P λ)(, )sin )(cos )((x x P x x P e n l x ω+ωλ 的二阶常系数非齐次线性微分方程的解法。
5.例题 例 验证函数 21 2 + =Cx y 是微分方程012=+-'y y x 的解。 解 将 212+ =Cx y 和Cx y 2='代入012=+-'y y x 的左边得 0112222=+--Cx Cx , 所以 21 2 + =Cx y 是方程012=+-'y y x 的解。 例 求微分方程 x y y 212-= '的通解。 解 这是可分离变量的微分方程, 分离变量得x dx y dy 212 =-,解 此方程如下: 1ln ln 2111ln 21C x y y +=+-11+-=?y y Cx . 即得通解为 )1(1+=-y Cx y . 例 求微分方程2 2 x xy y y -='的通解。 解 这是齐次方程,即12 -??? ??= x y x y dx dy ，令 x y u =u dx du x dx dy +=?
微分方程数值解法答案
包括基本概念，差分格式的构造、截断误差和稳定性，这些内容是贯穿整个教材的主线。解答问题关键在过程，能够显示出你已经掌握了书上的内容，知道了解题方法。这次考试题目的类型：20分的选择题，主要是基本概念的理解，后面有五个大题，包括差分格式的构造、截断误差和稳定性。 习题一 1． 略 2． y y x f -=),(，梯形公式：n n n n n n y h h y y y h y y )121(),(2111+-+=+- =+++，所以0122)1(01])121[()121()121(y h h y h h y h h y h h n h h n n n +--+--+-+=+-+==+-+= ，当0→h 时， x n e y -→。 同理可以证明预报-校正法收敛到微分方程的解. 3． 局部截断误差的推导同欧拉公式; 整体截断误差: ? ++++++-++≤1 ),())(,(11111n n x x n n n n n n n dx y x f x y x f R εε 11)(++-++≤n n n y x y Lh R ε，这里R R n ≤ 而111)(+++-=n n n y x y ε，所以 R Lh n n += -+εε1)1(，不妨设1 微分方程求解Newly compiled on November 23, 2020 求解微分方程 ：简单地说，就是去微分（去掉导数），将方程化成自变量与因变量关系的方程（没有导数）。 近来做毕业设计遇到微分方程问题，搞懂后，特发此文，来帮广大同学，网友。 1.最简单的例子： x dx dy 2= ——————》 C x y +=2 求微分方程 xy dx dy 2=的通解。 解 方程是可分离变量的，分离变量后得 两端积分 ： ,2??=xdx y dy 得： ,ln 12C x y += 从而 ： 2112 x C C x e e e y ±=±=+。 又因为 1C e ±仍是任意常数，可以记作C 2 x Ce y =。 非齐次线性方程 求方程25 )1(12'+=+-x x y y 的通解. 解：非齐次线性方程。 先求对应的齐次方程的通解。 01 2=+-x y dx dy ， 1 2+=x dx y dy ， 用常数变易法：把C 换成)(x u ，即令 2)1(+=x u y （1） 则有 )1(2)1('2+++=x u x u dx dy ， 代入原方程式中得 21 )1('+=x u ， 两端积分，得 C x u ++=23)1(3 2。 再代入（1）式即得所求方程通解 ])1(32[)1(23 2C x x y +++=。 法二： 假设待求的微分方程是: )()(x Q y x P dx dy =+ 我们可以直接应用下式 得到方程的通解，其中， 1 2)(+-=x x P ， 25 )1()(+=x x Q 代入积分同样可得方程通解 ])1(32[)1(23 2C x x y +++=， 2.微分方程的相关概念：(看完后你会懂得各类微分方程) 即得齐次方程通解。 ，代替分离变量，积分后将，，，则设的函数，解法：，即写成程可以写成 齐次方程：一阶微分方称为隐式通解。 得：的形式，解法：为：一阶微分方程可以化可分离变量的微分方程 或 一阶微分方程：u x y u u du x dx u dx du u dx du x u dx dy x y u x y y x y x f dx dy C x F y G dx x f dy y g dx x f dy y g dy y x Q dx y x P y x f y -=∴=++====+====+='??)()(),(),()()()()()()(0 ),(),(),(???一 阶线性微分方程： 全微分方程： 二阶微分方程： 二阶常系数齐次线性微分方程及其解法：微分方程求解