iGEM竞赛2016年南开大学团队matlab微分方程代码

这是一份整理于2012年的南开大学数据库原理100题，希望能帮助到学习数据库原理的同学们。

借助 Matlab 工具，探索求解微分方程的方法。本教程涵盖解析解和数值解的求解技巧，并提供实例和实验作业，加深理解。

提供微分方程解代码

偏微分方程的数值解，有时候真的挺头疼的。是 2D 扩散问题或者泊松方程，自己从头撸一套还挺费劲的。这份matlab 微分方程代码就蛮实用的，直接用Jacobi 方法来二维扩散（热传导）和泊松问题，结构也清晰，代码也不复杂。适合你快速上手试试数值方法的思路。代码参考了 PowerPoint 的演示 PDF，还比较详细。虽然有两个辅助文件Utilities.cu和Utilities.cuh没放上来，但不影响你理解主要逻辑。Jacobi 迭代部分挺经典的，搞懂这块再拓展 Gauss-Seidel、SOR 啥的都容易。建议你顺手看看这几个相关资源：Matlab 微分方程求解、泊松方程的数值解法，有些思

阐述了Matlab在解决微分方程及数学建模中的应用实例。

微分方程的解法一直是建模里绕不开的话题，MATLAB的工具箱是真的挺给力，适合新手入门。数学实验里的第四个任务就是搞定微分方程的求解，用MATLAB来做还挺省事的，不光能数值解，连符号解也能整。像ode45这种函数，用起来挺顺手的。只要定义好微分方程、初始条件和时间范围，一行代码就能跑出结果。如果你习惯看代码示例，可以看看这个基本示例，讲得还蛮清楚的，连图都画了。要是你对建模比赛感兴趣，国赛微分方程类获奖论文也可以瞄一眼，看看人家是怎么建模和解题的。实在搞不懂符号解和数值解区别？别急，这篇符号解法文章可以帮你理清思路。如果你经常写代码，建议写个通用模板，比如：function dydt = m

MATLAB 使用 Runge-Kutta-Fehlberg 方法解 ODE 问题，以有限个点进行计算，点间距由解本身决定。 可使用 ode23 求解 2-3 阶常微分方程组，使用 ode45 使用 4-5 阶 Runge-Kutta-Fehlberg 方法。 例如，在命令行中使用 ode45 函数代替 solver，其中 x' 是 x 的微分，而非 x 的转置。

matlab 的微分方程解法资源挺丰富的，尤其是对常微分方程的数值方法比较全，适合平时搞建模、做控制系统仿真的同学参考。文章不只是讲原理，还配了 MATLAB 实现，代码也挺清晰。比如欧拉法、Adams 方法这些常见套路，基本都能找到，而且用的语言你一看就懂，不绕弯子。如果你是新手，建议先从欧拉法的那篇开始，思路简单，代码也好上手。

微分方程的解法在 MATLAB 里算是标配了，适合搞科研或者工程建模的同学。用起来最顺手的还是 ode45，稳定还好上手。你只要把微分方程写进一个函数文件，比如dxdt.m，再在主脚本里调一下就行了，响应也快，代码也简单。嗯，还有一个完整的例子打包好了，叫dxdt_solve，里面包括怎么写方程、怎么调ode45，还有怎么画图结果，比较适合刚上手或者想回顾用法的你。微分方程像dy/dt = f(t, y)这种的常见，不管是做电路仿真、生物种群建模，还是控制系统设计，都离不开这类模型。你就把它写进dxdt.m，函数结构也挺直观的：function dydt = dxdt(t, y) dydt