MATLAB学习笔记(一)求解三阶微分方程

一、求解三阶微分方程

        对于多变量三阶微分方程求解问题,这里介绍一种求解方法。

例题如下:

        eq?%5Cddot%7Bx%7D_%7B1%7D+%5Cdot%7Bx%7D_%7B1%7D+x_%7B1%7D-%5Cdot%7Bx%7D_%7B2%7D-x_%7B2%7D%3D0

eq?%5Cdddot%7Bx%7D_%7B2%7D+%5Cddot%7Bx%7D_%7B2%7D+%5Cdot%7Bx%7D_%7B2%7D+x_%7B2%7D-%5Cddot%7Bx_%7B1%7D%7D-%5Cdot%7Bx_%7B1%7D%7D-x_%7B1%7D%3D0

        对于以上方程,给定边界条件,eq?x_%7B2%7D%5Cmid%20_%7Bt%3D0%7D%3D0eq?%5Cdot%7Bx%7D_%7B2%7D%5Cmid%20_%7Bt%3D0%7D%3D0eq?%5Cddot%7Bx%7D_%7B2%7D%5Cmid%20_%7Bt%3D0%7D%3D0eq?%7Bx%7D_%7B1%7D%5Cmid%20_%7Bt%3D0%7D%3D0eq?%5Cdot%7Bx%7D_%7B1%7D%5Cmid%20_%7Bt%3D0%7D%3D0。求解eq?x_%7B1%7Deq?x_%7B2%7D的表达式。

二、解题步骤

(1)建立辅助变量h(t表示自变量),令eq?h%5Cleft%20%28%201%20%5Cright%20%29%3Dx_%7B1%7Deq?h%5Cleft%20%28%202%20%5Cright%20%29%3D%5Cdot%7Bx%7D_%7B1%7Deq?h%5Cleft%20%28%203%20%5Cright%20%29%3D%7Bx%7D_%7B2%7Deq?h%5Cleft%20%28%204%20%5Cright%20%29%3D%5Cdot%7Bx%7D_%7B2%7Deq?h%5Cleft%20%28%205%20%5Cright%20%29%3D%5Cddot%7Bx%7D_%7B2%7D

(2)根据等式阶数确定向量个数,将每阶的求导依次写出来,由(1)可知,eq?%5Cfrac%7B%5Cmathrm%7Bd%7D%20h%281%29%7D%7B%5Cmathrm%7Bd%7D%20t%7D%3Dh%5Cleft%20%28%202%20%5Cright%20%29eq?%5Cfrac%7B%5Cmathrm%7Bd%7D%20h%283%29%7D%7B%5Cmathrm%7Bd%7D%20t%7D%3Dh%5Cleft%20%28%204%20%5Cright%20%29eq?%5Cfrac%7B%5Cmathrm%7Bd%7D%20h%284%29%7D%7B%5Cmathrm%7Bd%7D%20t%7D%3Dh%5Cleft%20%28%205%20%5Cright%20%29

(3)在matlab中建立函数,将每阶的求导表示出来;

(4)对x1和x2进行求解,画图如下,

%一阶导数
figure(1)
%二阶导数
figure(2)

此方法可以有效得到解的图。

55e6f3853ec344f5a9d431b966b9f06d.png

x1

51733842b4624986ba475514b6f5f4cc.png

x2

 

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.hqwc.cn/news/298565.html

如若内容造成侵权/违法违规/事实不符,请联系编程知识网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

Python遥感影像深度学习指南(2)-在 PyTorch 中创建自定义数据集和加载器

在上一篇 文章中,我们Fast.ai 在卫星图像中检测云轮廓,检测物体轮廓被称为语义分割。虽然我们用几行代码就能达到 96% 的准确率,但该模型无法考虑数据集中提供的所有输入通道(红、绿、蓝和近红外)。问题在于,深度学习框架(如 Keras、Fast.ai 甚至 PyTorch)中的大多数语…

复分析——第1章——复分析准备知识(E.M. Stein R. Shakarchi)

第一章 复分析准备知识 (Preliminaries to Complex Analysis) The sweeping development of mathematics during the last two centuries is due in large part to the introduction of complex numbers; paradoxically, this is based on the seemingly absurd no…

【DataGrip使用小技巧】2

【技巧1】要将活动编辑器与剪贴板内容进行比较,请右键单击编辑器中的任意位置,然后从上下文菜单中选择compare with Clipboard. 如图所示 剪贴板和编辑器的比较一目了然 【技巧2】若要创建表的备份副本,请将表拖到数据库树视图中同一数据源的…

2023.12.24 关于 Redis 中 String 类型内部编码 及 应用场景

目录 String 类型内部编码 3 种内部编码方式 String 类型应用场景 Cache 缓存 键名命名规则 计数(Counter) 共享会话(Session ) 手机验证码 总结 String 类型内部编码 3 种内部编码方式 int:用来表示 64 位 —…

VUE项目中使用jessibuca实现直播流播放

静态文件引入&#xff1a; 将静态文件存放到public文件夹下 在index.html中引入 具体实现&#xff1a; <template><div class"container" ref"container" style"width:100%;height:200px;"> </div> </template><s…

zabbix-邮件告警与钉钉机器人告警

zabbix-邮件告警 在部署邮件告警前需要先将zabbix部署成功 邮件告警需要得到邮箱的授权码&#xff0c;这里将以qq邮箱为例。 1.获取邮箱的授权码 2.编辑Email的相关信息 3.测试Email是否配置成功 测试成功。 4.为用户添加邮箱报警服务 5.创建主机群组、主机、监视项、触…

Maven高级篇

Maven依赖管理原则; 可选依赖&#xff1a;隐藏当前项目中的指定的包&#xff0c;如此&#xff0c;别的包引用当前包时&#xff0c;当前包中的指定包就被隐藏了&#xff0c;在别的包中无法看到隐藏的包 排除依赖&#xff1a;指定排除引用包中所包含的依赖&#xff0c;防止与当…

【数据结构和算法】删掉一个元素以后全为 1 的最长子数组

其他系列文章导航 Java基础合集数据结构与算法合集 设计模式合集 多线程合集 分布式合集 ES合集 文章目录 其他系列文章导航 文章目录 前言 一、题目描述 二、题解 2.1 方法一&#xff1a;滑动窗口 2.2 滑动窗口解题模板 三、代码 3.1 方法一&#xff1a;滑动窗口 四…

【每日一题】美丽塔 II

Tag 【单调栈】【数组】【2023-12-21】 题目来源 2866. 美丽塔 II 题目解读 题目意思相对明确&#xff0c;所谓的美丽塔数组就是山状数组&#xff0c;即有一个高度为 maxHeight[i] 的山峰&#xff0c;山峰两侧的高度要小于 maxHeight[i] 并且小于各自的允许高度。需要找出满…

台达A2-M伺服

地址: P3.00 从站地址0x01~0x7F【1~127】 P3.01 通讯速度UZYX【0403】 X:0【4800】1【9600】2【19200】3【38400】4【57600】5【115200】Z:0【125Kbit/s】1【250】2【500】3【750】4【1Mbit/s】 P3.02 通讯格式6【8N2】7【8E1】8【8O1】 P3.03 1通讯错误刹停…

每日一题——LeetCode141.环形链表

个人主页&#xff1a;白日依山璟 专栏&#xff1a;Java|数据结构与算法|每日一题 文章目录 1. 题目描述示例1&#xff1a;示例2&#xff1a;示例3&#xff1a;提示: 2. 思路3. 代码 1. 题目描述 给你一个链表的头节点 head &#xff0c;判断链表中是否有环。 如果链表中有某…

synchronized详解

synchronized详解 基本使用源码解析常见面试题好书推荐 基本使用 Java中的synchronized关键字用于在多线程环境下确保数据同步。它可以用来修饰方法和代码块 当一个线程访问一个对象的synchronized方法或代码块时&#xff0c;其他线程将无法访问该对象的其他synchronized方法…