实验目标:

大家都是初学者嘛,既然会点进这种标题进来的话。但是我们这里也不准备解释什么是常微分方程了,直接进入本文实验的终极目标:

告诉你

(y+x)dy+(x−y)dx=0

(y+x)dy+(x-y)dx=0
然后问这个微分方程组的解是什么?

有没有一脸懵13的感觉ahhh,这tm哪里是变量分离方程了?emmm它确实不是,只不过我们这里也不准备花太多时间在变量分离方程上面(哎,要想解好变量分离方程建议大家去查查数学分析里面的求积分的方法更靠谱,这篇文章更多是在讲如何用变量变换把较一般的方程变成变量分离方程。)

先解释一下变量分离方程吧

emmm大家都知道万丈高楼平地起咯,虽然我觉得没有必要解释但是还是照顾一下初学者吧。变量分离方程大概就是长这个样子的:

dydx=f(x)φ(y)

\frac{dy}{dx}=f(x)\varphi(y)

有什么特点呢?特点就是这个方程特别容易解呀~很容易就可以将这个方程化为:

1φ(y)dy=f(x)dx

\frac{1}{\varphi(y)}dy=f(x)dx
这样的形式,然后我们再两边同时积分:

∫1φ(y)dy=∫f(x)dx

\int \frac{1}{\varphi(y)}dy=\int f(x)dx

接下来就是在考察大家的微积分知识了。

ok我们接下来进入正题,说一说一些常见的需要我们用些小手段进行才能进行变量分离的微分方程。

形如dydx=g(yx)\frac{dy}{dx}=g(\frac{y}{x})的形式转化为微分方程

其实这种方程叫齐次微分方程,我们用什么手段对其进行变换呢?让我思考一下。。。

先想好一个变换函数(这名字我自己想的),比如说

u=yx

u=\frac{y}{x}
然后我们在把这个变换函数给变换一个形式:

y=ux

y=ux
然后我们再对两边同时求导

dydx=u′⋅x+u

\frac{dy}{dx}=u'·x+u
看!我们又得到了一个dydx\frac{dy}{dx}的表达式,我们现在把它代入我们的标题中,就有:

u′⋅x+u=g(u)

u'·x+u=g(u)
这里已经得到一个变量分离函数啦~不造大家有米有看出来呀?没有看出来的话,我再化简一点点就有:

dudx=g(u)−ux

\frac{du}{dx}=\frac{g(u)-u}{x}
这不就是了吗?大家应该还记得怎么解吧···不记得的话自己看上面吧。

好的,有了这个变换—求导—重新代入的这几步小手段,我们可以轻而易举的搞定dydx=g(yx)\frac{dy}{dx}=g(\frac{y}{x})这种形式的微分方程,现在我们可以解决啥问题呢,比如dydx=yx+tanyx\frac{dy}{dx}=\frac{y}{x}+tan\frac{y}{x},比如xdydx+2xy−−√=y(x<0)x\frac{dy}{dx}+2\sqrt{xy}=y\quad (x
用以上所述的方法,这些基本上都只要做一些化简,或者根本就不用化简就可以搞定的微分方程。大家应该没有这需求,我也懒得在这里演示了。

接下来我们看看另外一种更为一般的且可以变量分离的微分方程:

dydx=a1x+b1y+c1a2x+b2y+c2

\frac{dy}{dx}=\frac{a_1x+b_1y+c_1}{a_2x+b_2y+c_2}
这样的形式。
显然我们之前所说的只不过是这个方程的一个特例而已。
那么这个样子的微分方程还可以用我们之前的套路吗?老实说不太行,我们来具体分析一下吧。
第一种情况:a1a2=b1b2=c1c2=k\frac{a_1}{a_2}=\frac{b_1}{b_2}=\frac{c_1}{c_2}=k的话那么有

dydx=k

\frac{dy}{dx}=k
接下来怎么做应该不用我说了吧
第二种情况:

a1a2=b1b2=k≠c1c2

\frac{a_1}{a_2}=\frac{b_1}{b_2}=k\neq \frac{c_1}{c_2}
这个时候怎么办呢?
构造变量分离方程u=a2x+b2yu=a_2x+b_2y
求导数有

dudx=a2+b2dydx

\frac{du}{dx}=a_2+b_2\frac{dy}{dx}
将原式代入进来:

dudx=a2+b2ku+c1u+c2

\frac{du}{dx}=a_2+b_2\frac{ku+c1}{u+c2}
这个东西是啥?它的右边只有uu而没有xx了啊,它就是一个变量分离方程。
好的,接下来看最后一种情况,这也是我们无法用以前的变换—求导—重代入套路实现的一种情况:

a1a2≠b1b2

\frac{a_1}{a_2}\neq \frac{b_1}{b_2}
那我们怎么做呢,我们要想办法消去常数项,用给x和y分别加上一个常数的方法x和y分别加上一个常数的方法,为什么要消去常数项呢?我们来看看消去之后的效果,方程就可以变成:

dYdX=a1X+b1Ya2X+b2Y

\frac{dY}{dX}=\frac{a_1X+b_1Y}{a_2X+b_2Y}
这个是啥呀,这个就是g(YX)g(\frac{Y}{X})的形式呀,分子分母同时除以X就可以了。
那么现在问题来了,我们怎么样才能得到这个形式呢?不妨假设X=x−α,Y=y−βX=x-\alpha,Y=y-\beta,则原式可以变成:

dYdX=a1x+b1y−a1α−b1β+c1a2x+b2y−a2α−b2β+c2

\frac{dY}{dX}=\frac{a_1x+b_1y-a_1\alpha-b_1\beta+c_1}{a_2x+b_2y-a_2\alpha-b_2\beta+c_2}
然后求解出使得这个式子的常数项0的α和β\alpha和\beta就可以了。
求解出来之后就按照我们刚才所讲的思路用dydx=g(yx)\frac{dy}{dx}=g(\frac{y}{x})的解题方法进行求解就可以了。

常微分实验(1.1) 变量分离方程与变量变换相关推荐

  1. 2.1 变量分离方程与变量变换

    2.1.1变量分离方程 形如 dydx=f(x)ϕ(y)\frac{dy}{dx}=f(x)\phi(y)dxdy​=f(x)ϕ(y) 的方程,称为变量分离方程,这里的f(x),ϕ(y)f(x),\p ...

  2. 【未完成】常微分实验3.3:解 连续的初值 可微性定理

    日常回顾 回顾我们之前所学,我们已经知道了如果对于一个待解的微分方程一个初值点,且在这个初值点的附近有利普希茨条件成立的话,那么在这个初值点的附近的内部(注意是附近的内部,附近的内部,初值点附近的内部 ...

  3. 考研:研究生考试(十五天学完)之《高等数学上/下册》研究生学霸重点知识点总结之考试内容各科占比及常考知识重点梳理(函数极限连续、一元/多元函数微分学/积分学、常微分函数、向量代数与空间几何、无穷级数)

    考研:研究生考试(十五天学完)之<高等数学上/下册>研究生学霸重点知识点总结之考试内容各科占比及常考知识重点梳理(函数极限连续.一元/多元函数微分学/积分学.常微分函数.向量代数与空间几何 ...

  4. Matlab求解常微分方程组

    求解这个常微分方程组. 初始条件为              其中ε取0.01,a是有上限的参数,求解方程的目的其实是找出a的临界值. syms y(t) for i = [0:0.5:1.5,1.7 ...

  5. 一阶微分方程的物理意义_MIT—微分方程笔记24 一阶常微分方程组

    18.03 微分方程 Differential Equations 第四单元 一阶常微分方程组 Unit 4 First-Order Systems 第24讲 一阶常微分方程组 第25讲 常系数齐次线 ...

  6. python求解微分方程组_python – SymPy / SciPy:求解具有不同变量的常微分方程组...

    我是SymPy和Python的新手,我目前正在使用Python 2.7和SymPy 0.7.5,其目标是: a)从文本文件中读取微分方程组 b)解决系统问题 我已经阅读了this question和t ...

  7. matlab中常微分方法,MATLAB解常微分方程组的解法(好东西要共享)

    1:问题 常微分方程的初值问题的标准数学表述为:y'=f(t,y),a<=t<=b,y(a)=y(0) :我们要求解的任何高阶常微分方程都可以用替换法化为上式所示的一阶形式,其中y为向量, ...

  8. python解常微分方程组

    美国气象学家洛伦兹(E.N.Lorenz)是混沌理论的奠基者之一.20世纪50年代末到60年代初,他的主要工作目标是从理论上进行长期天气预报研究.他在使用计算机模拟天气时意外发现,对于天气系统,哪怕初 ...

  9. matlab时域微分,实验一 MATLAB用于时域分析.doc

    实验一 MATLAB用于时域分析 实验一 MATLAB用于时域分析 一.实验目的 学会LTI系统模型的建立: 掌握相关函数的调用格式及作用: 掌握连续信号的基本运算. 二.实验设备 计算机一台,装有M ...

最新文章

  1. 《Redis设计与实现》简读
  2. 开发者需要了解的WebKit
  3. 最长不重复字符串python_Python简单实现查找一个字符串中最长不重复子串的方法...
  4. 浅谈.Net异步编程的前世今生----APM篇
  5. oracle clob截取_Oracle数据库设计规范建议
  6. python物理引擎模拟三体_一个物理引擎能不能模拟少量粒子之间的力?
  7. [原创]个人工具 - 对APK极限压缩并对齐的工具(58.ReExtremeZipAndAlignAPK)
  8. protobuf3 自定义option_Protobuf3语法详解
  9. 【java学习之路】(java SE篇)(练习)常用类和异常的简单练习
  10. python中stacked_栈式自动编码器(Stacked AutoEncoder)
  11. BP神经网络代码示例
  12. 只有PHP大牛才能读懂的内涵图
  13. 计算机软件企业申请商标,软件商标注册申请流程
  14. 用C语言读取txt文件并统计字符数
  15. 【历史上的今天】1 月 25 日:电子游戏起源;《吃豆人》作者出生;“蠕虫王”问世
  16. 微信小程序云开发项目——多肉植物销售小程序
  17. 软件设计领域没有银弹,但代码大师MaxKanat-Alexander的建议绝对能给你带来启发...
  18. UOS系统适配-常用开发工具安装
  19. html水平排列三个元素,且三者占宽度相同
  20. HSV色彩空间筛选 2021-10-06

热门文章

  1. 掌财社寒山:白马股和黑马股有什么区别?
  2. Excel如何比较两表格部分字段差异
  3. 微信小程序抓包——小白自用版
  4. 快速傅里叶变换-FFTW库的使用-参考和翻译官方文档
  5. 线程池面试题灵魂三问
  6. ValueError问题的解决
  7. Android 高效加载大图片
  8. 把字节数B转换为KB,MB,GB的方法
  9. html5 drm,DRM迈入HTML5倒计时:为什么DRM是邪恶的
  10. docker pull很慢解决办法