线性代数笔记24——微分方程和exp(At)

  原文:https://mp.weixin.qq.com/s/COpYKxQDMhqJRuMK2raMKQ

  微分方程指含有未知函数及其导数的关系式,解微分方程就是找出未知函数。未知函数是一元函数的,叫常微分方程;未知函数是多元函数的,叫做偏微分方程。常微分方程有时也简称方程。微分方程是一门复杂的学科,对于常微分方程来说,可以使用特征值和特征向量的知识求解。

  

  相关前置知识:

    微分方程:单变量微积分11——常微分方程和分离变量

    泰勒公式:单变量微积分30——幂级数和泰勒级数 

    泰勒公式在0点展开的原因:多项式函数能够拟合非线性问题原理

    求逆矩阵:线性代数笔记8——求解逆矩阵

    求行列式:线性代数20——行列式和代数余子式

    特征值和特征向量:线性代数22——特征值和特征向量

      矩阵对角化:线性代数笔记23——矩阵的对角化和方幂

 

常微分方程的一般解法

  根据概念构造一个常微分方程:

  目标是求得原函数u=u(t)的具体形式。通过积分求解:

  这就是最终答案的通解,C是任意常数。实际上这种解法就是利用了不定积分的知识:

  

  如果du/dt=u,可以使用分离变量法的求解方式:

  也就是说,当函数的导数是函数本身的时候,这个函数就是型如Aet的函数,由于A=eC是任意常数,所以经常用C代替A,写成u=Cet的形式。

 

  同理,对于du/dt=λu,微分方程的解是u(t)=Ceλt。当t=0时:

  由于C是任意常数,因此可以取C=u(0),得到u(t)= u(0)eλt,这样做可以去掉常数C。在实际问题中,u可以表示关于时间t的函数,对于时间来说,通常可以把t=0看作初始条件。

常微分方程与矩阵

  现在将常微分方程扩展为常微分方程组,u1=u1(t),u2=u2(t),初

要在MATLAB中求解二元偏微分方程,可以按照以下步骤进行操作。 1. 打开MATLAB,在命令窗口中输入 "pdepe" 函数,并按回车键。这个函数用于求解偏微分方程。 2. 根据提示,输入偏微分方程的系数。系数可以根据具体的方程进行设置。 3. 输入偏微分方程的边界条件初始条件。这些条件需要根据具体的问题进行设置。 4. 输入求解的时间范围空间范围。 5. 根据提示,进行可选的设置,如相对容差绝对容差。 6. 根据具体的方程类型,选择使用有限元法或有限差分法进行求解。 7. 运行程序,MATLAB会自动计算并给出偏微分方程的解。 参考文献: 点击第7个图标(显示PDE字样),按提示输入偏微分方程的系数即可。在这里笔者求解波动方程:∂2u∂2t=∇u. 本课件内容首先介绍了MATLAB进行数学建模的方法,给出了优化求解方程组求解的示例,阐述了数学建模的思想;然后介绍了MATLAB在信号处理方面的应用,演示了音频图像的读取、分析处理过程;最后讲解了使用MATLAB进行...。 第四行第五行表示相对容差绝对容差,笔者查看了Matlab帮助中心,大概了解到这两个参数似乎与浮点数0的截断精度有关,太小的话会延长计算时间,如果你想了解更多,笔者把链接提供上来Absolute tolerance - MATLAB & Simulink - MathWorks 中国,假如我们对计算精度没有要求的话,使用默认值就可以了。这里笔者为了演示使用了0.0010.0001。如果想跟着一起做,那么笔者把方程的代码也放上来:第一个是atan(cos(pi/2*x)),第二个是3*sin(pi*x).*exp(cos(pi*y))。<span class="em">1</span><span class="em">2</span><span class="em">3</span> #### 引用[.reference_title] - *1* *3* [Matlab偏微分方程快速上手:使用pde有限元工具箱求解二维偏微分方程](https://blog.csdn.net/weixin_47006934/article/details/113524513)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v93^chatsearchT3_1"}}] [.reference_item style="max-width: 50%"] - *2* [Matlab基础应用学习笔记.md](https://download.csdn.net/download/weixin_52057528/88284511)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v93^chatsearchT3_1"}}] [.reference_item style="max-width: 50%"] [ .reference_list ]
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

我是8位的

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值