simulink建模之求解微分方程

simulink建模之求解微分方程

0.前言

simulink建模的基础其实就是利用加减乘除等数学基本四则运算,有时加上积分和微分对数学公式进行表示,因此理解这些数学公式背后的物理意义,才是重点。所以各行各业都有门槛,但是通用的数学知识却是基础中的基础,重点中的重点。那么如果已经知道物理含义的前提下,结合simulink图形化的公式这个工具表示出来,可以方便快速定位问题,解决问题,快速进行分析,这项能力就显得很重要。后续还可以自动生成代码,结合各种编程语言,动态链接库等知识点,完成工程师的需求。本博客重点所记录为利用微分方程式转化为simulink模型(即分析问题抽象为数学模型的能力)以及给出simulink模型快速分析出数学公式(即看懂模型优化数学公式的能力),此两种能力相辅相成,若同时修炼内力大增,堪比九阳神功,可与算法和数据结构相比光辉,但学习难度应该相对较少(个人片面观点)。

1.题目如下

如前言所述,建模的基础能力,就是对数学公式的转化能力(个人狭隘理解)。从解微分方程开始入手,因为微分方程涉及高等数学知识,基本上包含很多重要算法的雏形。废话到此结束,请看下文分解。
在这里插入图片描述

2.分析题目

第一种看法:如下图所示,如果采取微分的方法,按步骤一步步的去微分,发现无法继续下去
在这里插入图片描述
网上查阅资料对比可知:就算成功变好了模型,也会出现代数环的问题(PS:就算连续系统加了memory模块一样不满足情况,所以请问找到初始条件到哪里设置??)。如下图
在这里插入图片描述

第二种看法:建立x关于x的一阶导数和二阶导数的表示关系,但是此题是关于x二次方程也不适用,simulink不存在双向箭头,如下图
在这里插入图片描述
第三种看法:建立X"关于x和X’的关系,如下图:
在这里插入图片描述
生成的效果图如下:
在这里插入图片描述

3.总结此题

(人长的丑,字也写的差…)具体如下图所示:
在这里插入图片描述

4.举一反三

重头戏来了,解微分方程的套路如下:
第一步:建立简便易求解的y = f(y’,y"…)的此类表达形式
第二步:尽量多使用积分模块,少使用微分模块,关键是建立等价联系,方程组的左边等于右边的关系!
如下题,屡试不爽:
在这里插入图片描述

5.后续补充用传递函数的方式求解微分方程

待续。。。。。

  • 61
    点赞
  • 279
    收藏
    觉得还不错? 一键收藏
  • 23
    评论
要在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.001和0.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 ]

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 23
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值