当前位置：文档之家 › 微分方程求特解的公式

微分方程求特解的公式

  • 格式：docx
  • 大小：36.89 KB
  • 文档页数：2

微分方程求特解的公式

微分方程是数学中的重要概念，广泛应用于自然科学、工程技术和经济学等领域。求解微分方程的特解是解决实际问题的关键步骤之一。本文将介绍微分方程求特解的公式。

一、一阶线性常微分方程的特解公式

对于一阶线性常微分方程形如：dy/dx + P(x)y = Q(x)，其中P(x)和Q(x)是已知函数，则可以得到特解公式为：

y = e^(-∫P(x)dx) * [∫Q(x)e^(∫P(x)dx)dx + C]，

其中C为任意常数。

二、二阶常系数齐次线性微分方程的特解公式

对于二阶常系数齐次线性微分方程形如：ay'' + by' + cy = 0，其中a、b、c是已知常数，则可以得到特解公式为：

1. 当方程的特征方程有两个不相等的实根r1和r2时，特解公式为：

y = C1e^(r1x) + C2e^(r2x)，其中C1和C2为任意常数。

2. 当方程的特征方程有两个相等的实根r1=r2=r时，特解公式为：

y = C1e^(rx) + C2xe^(rx)，其中C1和C2为任意常数。

3. 当方程的特征方程有两个共轭复根α±βi时，特解公式为：

y = e^(αx)(C1cos(βx) + C2sin(βx))，其中C1和C2为任意常数。

三、二阶非齐次线性微分方程的特解公式 对于二阶非齐次线性微分方程形如：ay'' + by' + cy = f(x)，其中a、b、c是已知常数，f(x)是已知函数，则可以得到特解公式为：

1. 根据待定系数法，特解形式可以根据f(x)的类型选择。

* 当f(x)是常数时，特解形式为y = k，其中k是常数。

* 当f(x)为多项式时，特解形式为y = P(x)，其中P(x)是与f(x)次数相同的多项式。

* 当f(x)为三角函数时，特解形式为y = A sin(mx) + B cos(mx)，其中A和B是待定常数，m是f(x)的角频率。

* 当f(x)为指数函数时，特解形式为y = Ae^(kx)，其中A和k是待定常数。

* 当f(x)为幂函数时，特解形式为y = Ax^n，其中A和n是待定常数。

2. 将特解形式代入非齐次微分方程并解方程组，可以得到特解公式中的待定常数。

本文简单介绍了微分方程求特解的公式，包括一阶线性常微分方程的特解公式和二阶常系数齐次非齐次微分方程的特解公式。通过应用这些公式，可以解决许多实际问题，并在科学研究与工程技术中发挥重要作用。了解这些公式对于深入理解微分方程和应用数学方法具有重要意义。

合集下载

常微分方程通解、特解、所有解的区别与联系

常微分方程通解、特解、所有解的区别与联系

第３Ｏ卷第２期　

２０１４年４月　大　学　数　学　

ＣＯＬＬＥＧＥ　ＭＡＴＨＥＭＡＴＩＣＳ　Ｖ０１．３Ｏ，№．２　Ａｐｒ．２０１４　

常微分方程通解、特解、所有解的区别与联系　

刘雄伟，　王　晓　

（国防科技大学理学院，湖南长沙４１００７３）　

［摘　要］通过具体实例分析、讨论了高等数学中常微分方程的通解、特解和微分方程的所有解之间的区　别与联系，并对高等数学教材中二阶线性微分方程的降阶法与二阶常系数非齐次线性微分方程特解求解过　程中的作法进行了说明．　［关键词］高等数学；常微分方程；通解；定解条件；特解　［中图分类号］０１７５．１　［文献标识码］ｃ　［文章编号］１６７２—１４５４（２０１４）０２—００８８—０３　

在国内的高等数学教材与常微分方程教材中，对于　阶常微分方程　

Ｆ（ｚ，　，五ｄｙ，器，…，　）一０，　㈩　

关于通解、特解通常定义如下ｎ　］：　

如果包含有　个相互独立的任意常数Ｃ　，Ｃ　一，Ｃ　的关系式　

（　，Ｙ，Ｃ１，Ｃ２，…，Ｃ　）一０　（２）　

确定的函数　一　（　，Ｃ　，Ｃ　，…，Ｃ　）是（１）的解，则称（２）为（１）的通解．通过初值条件　

ｙ（ｘｏ）一Ｙ１，Ｙ　（　ｏ）＝＝＝Ｙ２，…，Ｙ‘　（ｚｏ）一Ｙ　（３）　确定通解中的任意常数后，得到的解　（　，　）一０称为（１）的特解．　

对于非数学专业的学生来说，在高等数学课程的常微分方程学习过程中，通解、特解、所有解的定义　

及相互关系，二阶齐次线性微分方程降阶法和二阶常系数非齐次线性微分方程所设特解的处理方式，一　直是学生在学习过程中容易产生困惑的内容．为了更好的理解和求解常微分方程，关于常微分方程的通　

解、特解和微分方程的所有解的基本概念和求解过程中的一些处理方式有几点需要特别注意．　（ｉ）通解并不一定包含微分方程的所有解．　例如［１］，Ｙ—ｓｉｎ（ｘ＋Ｃ）是微分方程　

丁ｄｙ＝　（４）　丁一Ｖ　１一　

的通解，但Ｙ一土１也是（４）的解，显然不包含在通解中，即通解不包含方程的所有解．因此，在做练习时，　

实验二MATLAB 求微分方程的解

实验二MATLAB 求微分方程的解

实验二 微分方程求解

一、问题背景与实验目的

实际应用问题通过数学建模所归纳而得到的方程，绝大多数都是微分方程，真正能得到代数方程的机会很少．另一方面，能够求解的微分方程也是十分有限的，特别是高阶方程和偏微分方程（组）．这就要求我们必须研究微分方程（组）的解法，既要研究微分方程（组）的解析解法（精确解），更要研究微分方程（组）的数值解法（近似解）．

对微分方程（组）的解析解法(精确解)，Matlab 有专门的函数可以用，本实验将作一定的介绍．

本实验将主要研究微分方程(组)的数值解法（近似解），重点介绍 Euler 折线法．

二、相关函数（命令）及简介

1．dsolve('equ1','equ2',„)：Matlab 求微分方程的解析解．equ1、equ2、„为方程（或条件）．写方程（或条件）时用 Dy 表示y 关于自变量的一阶导数，用用 D2y 表示 y 关于自变量的二阶导数，依此类推．

2．simplify(s)：对表达式 s 使用 maple 的化简规则进行化简．

例如：

syms x

simplify(sin(x)^2 + cos(x)^2)

ans=1

3．[r,how]=simple(s)：由于 Matlab 提供了多种化简规则，simple 命令就是对表达式 s 用各种规则进行化简，然后用 r 返回最简形式，how 返回形成这种形式所用的规则．

例如：

syms x

[r,how]=simple(cos(x)^2-sin(x)^2)

r = cos(2*x)

how = combine

4．[T,Y] = solver(odefun,tspan,y0) 求微分方程的数值解．

说明：

(1) 其中的 solver为命令 ode45、ode23、ode113、ode15s、ode23s、ode23t、ode23tb 之一．

(2) odefun 是显式常微分方程：00)(),(ytyytfdtdy

通过MATLAB求二阶全微分方程解析解

通过MATLAB求二阶全微分方程解析解

1. 关于二阶全微分方程a𝒅𝒅𝒅𝒅𝒅+𝒅𝒅𝒅𝒅𝒅+𝒅𝒅=𝒅𝒅𝒅⁡(𝒅𝒅) ，不同的a,b,c取值会求出不同的解析解，解析解又是由齐次解和特解组成。其中，齐次解由特点方程决定，而特解的决定因素那么比较复杂。

2. 关于二阶全微分方程的分析，咱们大致分为三种情形：

b^2-4ac>0(两个不同的实根)

b^2-4ac=0（两个相同的重根）

b^2-4ac<0（两个不同的复数根）

对三种情形进行MATLAB编程，分析齐次解和特解后，再改变W的值，观看解析解的转变

3. b^2-4ac>0的情形

STEP1：求解析解

s1=dsolve('D2y+3*Dy+2*y=0','y(0)=2,Dy(0)=0','t');

s2=dsolve('D2y+3*Dy+2*y=sin(t)','y(0)=2,Dy(0)=0','t');

s3=dsolve('D2y+3*Dy+2*y=sin(2*t)','y(0)=2,Dy(0)=0','t');

s4=dsolve('D2y+3*Dy+2*y=sin(5*t)','y(0)=2,Dy(0)=0','t');

s5=dsolve('D2y+3*Dy+2*y=sin(13*t)','y(0)=2,Dy(0)=0','t');

s6=dsolve('D2y+3*Dy+2*y=sin(25*t)','y(0)=2,Dy(0)=0','t');

STEP2：绘制图形

（1）求w=1情形下的通解和齐次解

t=1::10;

s1=4*exp(-t)-2*exp(-2*t) %general solution

s2=-3/10*cos(t)+1/10*sin(t)-11/5*exp(-2*t)+9/2*exp(-t) %special

solution

subplot(2,1,1);

plot(t,s2);

xlabel('t')

微分方程数值解 打靶法求边值问题

微分方程数值解 打靶法求边值问题

微分方程数值解 实验七

学号 班级 姓名 指导教师

实验题目 打靶法求边值问题 评 分

实验要求：

（1） 掌握打靶法求边值问题

问题： 用打靶法求如下边值问题

12132124-2,,(1)0,(2)0,[1,2]yyytyeyyt

2、实验内容：（求解过程，包括简单的分析，求解程序代码，所得结果）

程序代码：

（1）先猜测初值，调用T7函数龙格库塔算法计算出该初值下的边界值，与题中的边界值相减，计算差值，既F(s)=yb-y(b).

function y=f(a)

%利用龙格库塔算法计算F(s)=yb-y(b)

qj=[1,2];

cz=[0,a];

[t,w]=T7(qj,cz,0.01);

y=0-w(end,2);

end

（2）在ydot函数中输入题目中的微分方程，用龙格库塔算法计算预设初值下的边界值。

function [t,w]=T7(qj,cz,h)

%求解实验7问题

%qj为所求区间

%cz=[y1,y2]为y1,y2的初值

%h为步长

%eg：[t,w]=T7([qj,cz,0.01);

n=round((qj(2)-qj(1))/h);

t=zeros(n,1);

w=zeros(n,length(cz));

t(1)=qj(1);

w(1,:)=cz;

for i=1:n t(i+1)=t(i)+h;

w(i+1,:)= RungeKutta(t(i),w(i,:),h);

end

end

function y = RungeKutta(t,w,h)

%龙格—库塔法求解微分方程

%h为步长

k1= ydot(t,w);

k2= ydot(t+h/2,w+h/2*k1);

k3= ydot(t+h/2,w+h/2*k2);