最近同学毕设需要求解循坏摆的微分方程,我在帮忙过程中学习了一下常微分方程的解析解和数值解的求法,在此分享。

以下讲解遵循Matlab官方文档提供的方程和写法。

(强烈建议大家有问题多看官方文档,非常有用)

介绍一下核心函数

常用形式:S = dsolve(eqn,cond)

功能介绍:解微分方程eqn,其中eqn是一个符号方程,cond是初值条件(边界条件)。使用diff和==表示微分方程。例如,diff(y,x)==y表示方程dy/dx=y。通过将eqn指定为这些方程的向量来求解微分方程组。

1. 一阶常微分方程求解(无初值)

方程:dy/dt=ay

代码:

syms y(t) a; %syms 定义方程用的变量和未知常量(可有可无)。y(t)表示y是关于t的函数,求解t。

eqn=diff(y,t)==a*y; %diff(y,t)即dy/dt,默认为1阶。

S=dsolve(eqn) %dsolve()为求解函数,将方程放入即可

说明:首先定义方程所用变量和未知常量,然后列出微分方程,最后dsolve()求解。

结果:

S =

C2*exp(a*t)

分析:C2为未定的常数项,需要提供一个初值,才能确定。

2. 二阶常微分方程求解(无初值)

方程:dy²/dt²=ay

代码:

syms y(t) a; %syms 定义方程用的变量和未

Matlab求时变微分方程组解,Matlab求常微分方程组的解析解相关推荐

  1. matlab常系数线性矩阵微分方程组,基于Matlab常系数线性微分方程组的求解

    ·基础数学· 基于 Matlab 常系数线性微分方程组的求解* 严水仙 (赣南师范大学 数学与计算机科学学院,江西 赣州 341000) 摘 要: 在常微分方程课程教学中,常系数线性微分方程组可以通过 ...

  2. matlab解二阶微分方程组ode,MATLAB解含参数方程、矩阵方程、二阶微分方程组

    1.如下公式: 其中x.y.m.n为参数,a.b为未知数,利用MATLAB求解方程: syms x y m n a b; [a,b]=solve('x=m*cos(a)+n*cos(a+b)','y= ...

  3. 龙格库塔法解微分方程组的matlab程序,MATLAB实例源码教程:龙格库塔法求解微分方程组源代码实例.doc...

    MATLAB实例源码教程:龙格库塔法求解微分方程组源代码实例.doc MATLAB实例源码教程龙格库塔法求解微分方程组源代码实例题目用经典 Runge-Kutta方法求下列一阶微分方程组的近似解y1 ...

  4. matlab求非线性函数的解,MATLAB求解非线性方程(转)

    1.解方程 最近有多人问如何用matlab解方程组的问题,其实在matlab中解方程组还是很方便的,例如,对于代数方程组Ax=b(A为系数矩阵,非奇异)的求解,MATLAB中有两种方法: (1)x=i ...

  5. matlab 一元线性微分,怎样用matlab求一元线性函数极值

    clear; %%%%% 清变量 clc; %%%%% 清屏幕 syms x; %%%% 定变量 f = sqrt(5^2+(8-x)^2)*28.8+sqrt(15^2+(5+x)^2)*7.2; ...

  6. matlab离散点数字微分,matlab数值积分与微分详解.ppt

    * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 格式: [I,n]= ...

  7. matlab 二阶非线性微分方程组,二阶非线性常微分方程的打靶法matlab实现.doc

    二阶非线性常微分方程的打靶法matlab实现.doc 二阶非线性常微分方程的打靶法1.问题:试用打靶法求二阶非线性常微分方程亮点边值的数值解:要求用Matlab编程计算,请给出一些例子,验证你的算法与 ...

  8. matlab求解一阶微分方程组

    在matlab里运用ode计算一阶微分方程组 可以根据需求选择合适的ode 举个例 程序为:function dy = odefun(t,y) dy = zeros(3,1); % 一个列向量 dy( ...

  9. Matlab求解常微分方程组

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

最新文章

  1. iOSunicode转中文
  2. 错误:Subquery returns more than 1 row 表示子查询返回了多行数据
  3. 读《我编程,我快乐,程序员的职业规划之道》有感
  4. Oracle数据库案例整理-Oracle系统执行时故障-Shared Pool内存不足导致数据库响应缓慢...
  5. python源码编译 mingw_在windows上用gcc(mingw32)从命令行编译Cython扩展
  6. 静态html js文件上传,js实现动态添加上传文件页面
  7. 【数据库系统原理】数据库设计
  8. python time、datetime模块学习使用
  9. Web Dev领域:2017精彩事件和 2018预测
  10. java基本类型运算溢出_JAVACard 基本数据类型的运算及溢出问题
  11. Balrum 汉化指引
  12. 学习笔记(5):2020华为HCIA/HCNA/数通/路由交换/实验/视频/教程/持续更新赠题库-HCIA数通IP地址编址及报头详解
  13. 使用BarTender连接Excel打印标签
  14. Unity进阶之ET网络游戏开发框架 01-下载、运行
  15. SpringBoot - Unable to start ServletWebServerApplicationContext due to missing ServletWeb
  16. 金蝶云·星空python插件示例代码
  17. supervisor备忘
  18. Material Design学习
  19. ORACLE EBS常用表及查询语句(二)
  20. SAP商超订单统一管理系统

热门文章

  1. N、Z、Q、R 分别代表什么
  2. nil Foundation blueprint模块代码解析
  3. Zephyr DMA
  4. 解决模拟器不能上网问题
  5. 七个习惯之六:统合综效
  6. Java黑皮书课后题第5章:5.3(将千克转换成磅)编写程序,显示下面的表格(注意:1kg为2.2磅)
  7. 好用速度又快,多御安全浏览器下载安装教程
  8. jQuery与JS实现AJAX实例!(附AJAX教程、路线图)
  9. 初中生怎样学习编程语言?
  10. 【Hyper-V】Windows10家庭版安装Hyper-V虚拟机