COMSOL Multiphysics多物理场仿真软件也提供了求救常微分方程(ODE)和偏微分方程(PDE)的接口,下面详细介绍一下。

(1)建立模型,选择模型向导–>零维–>数学–>全局常微分和微分代数方程(ge),选择研究,选择瞬态,点击完成
(2)在组件下面可以看到刚刚添加的全局常微分和微分代数方程(ge),在右边栏,全局方程那里输入需要求解的函数。
以上图电路实例来说,现有RLC串联电路,假设R、L、C的参数都为1,电容电压:uC(t)uC(t),电容电流:i(t)i(t),电感电压:uL(t)uL(t),电阻电压:uR(t)uR(t),输入电压uin=sin(2πft)uin=sin(2πft)则有下面的关系式:
uC(t)i(t)uL(t)uR(t)uin=uC(t)=CduC(t)dt=Ldi(t)dt=LCdu2C(t)dt2=Ri(t)=RCduC(t)dt=uC(t)+uL(t)+uR(t)
uC(t)=uC(t)i(t)=CduC(t)dtuL(t)=Ldi(t)dt=LCduC2(t)dt2uR(t)=Ri(t)=RCduC(t)dtuin=uC(t)+uL(t)+uR(t)

化简上式可得:
du2C(t)dt2+duC(t)dt+uC(t)=uin
duC2(t)dt2+duC(t)dt+uC(t)=uin
,换成常见的数学形式为
f′′+f′+f=sin(2πft)
f″+f′+f=sin(2πft)
要想在COMSOL求解该二阶微分方程,可以先用之前的方法在python里面求解,python代码如下
输出结果如下:
输入电源参数设置为2πf=1002πf=100,初值设定为y(0)=0y(0)=0,y′(0)=0y′(0)=0,由结果可知是一个衰减震荡的过程。

(3)步骤(2)得到的COMSOL设置结果如下:
更多资讯请关注公众号!

COMSOL求解常微分方程相关推荐

  1. matlab微分方程组边值,matlab求解常微分方程边值问题的方法

    matlab求解常微分方程边值问题的方法 Matlab 求解常微分方程边值问题的方法:bvp4c 函数 常微分方程的边值问题,即 boundary value problems ,简称 BVP 问题, ...

  2. MATLAB求解常微分方程

    MATLAB求解微分方程_Falcon的博客-CSDN博客_matlab微分方程求解 matlab求解常微分方程(组)---dsolve.ode系列函数详解(含例程)_假电工的真的博客-CSDN博客_ ...

  3. python解常微分方程_Python-sympy.dsolve求解常微分方程(组)

    这里分别介绍怎么利用sympy.dsolve求解常微分方程和常微分方程组. #首先利用sympy.dsolve求解单个的常微分方程: #代码 from sympy import Function, d ...

  4. matlab解微分方程组_MATLAB编程入门 求解常微分方程 通解 特解 数值解

    一.求微分方程的通解 clear all; f1=dsolve('Dy-y=sin(x)')     %默认情况下以t为自变量求解   求出的C2为积分常数 f2=dsolve('Dy-y=sin(x ...

  5. python解常微分方程龙格库_求解常微分方程组初值问题的龙格库塔法分析及其C代码...

    求解常微分方程组初值问题的 龙格库塔法分析及其 C 代码 1 .概 述 由高等数学的知识可知,一些特殊类型的常微分方程(组)能够求出给定初 始值的解析解, 而在科学与工程问题中遇到的常微分方程 (组) ...

  6. matlab 龙格-库塔 法求解常微分方程

    最近学习分室模型,里面碰到了用matlab 龙格-库塔 法求解常微分方程 研究了一阵子终于明白到底怎么实现了: 1. matlab 新建.m文件,编写龙格-库塔法求解函数 function [x,y] ...

  7. 利用sympy库求解常微分方程:dsolve()函数

    [小白从小学Python.C.Java] [计算机等级考试+500强双证书] [Python-数据分析] [sympy库的使用] 使用sympy库求解常微分方程 .dsolve()方法 选择题 下列说 ...

  8. 四阶龙格-库塔法求解常微分方程的初值问题

    算法原理和程序框图 龙格-库塔法是一种求其准确解y(x)在一系列点xi处y(xi)的近似值yi的方法,yi称为数值解.经典的四阶龙格库塔法方程如下: y'=f(t,y),y(t0)=y0输出按如下求解 ...

  9. Python求解常微分方程——sympy

    [常微分方程简介] 方程中未知量是函数而不是变量,且未知量涉及未知函数的导数的方程称为微分方程. 常微分方程(ordinary differential equation, ODE)是一类特殊情况,未 ...

  10. Python求解常微分方程

    sympy.numpy.scipy.matplotlib是强大的处理数学问题的库,可以执行积分.求解常微分方程.绘图等功能,其开源免费的优势可以与MATLAB媲美. 一阶常微分方程 from symp ...

最新文章

  1. 黑客带来了病毒攻击软件还是文化?
  2. 用安卓手机解锁目标检测模型YOLOv5,识别速度不过几十毫秒!
  3. springmvc如何访问静态文件,例如jpg,js,css
  4. 图像间距pitch是什么?(linesize、stride)(指图像中的一行图像数据所占的存储空间的长度)
  5. 2.1 《数据库系统概论》关系数据结构及形式化定义(关系、关系模式、关系数据库)
  6. 【贪心】逃跑(jzoj 1748)
  7. 深源恒际医疗票据OCR落地九省市 服务范围覆盖过半市场
  8. pythonrgb高精度浮点运算类型_python实现RGB字符串,按24位对齐后输出对应Integer行数字...
  9. 一条SQL的奇妙旅行
  10. python判断字符串回文_python实现--判断回文字符串、回文链表、回文数
  11. 关于原型污染漏洞的完整指南
  12. uploadify 3.1 /3.2 在Firefox,Chrome 浏览器下,进入页面请求两次问题解决办法。
  13. gflags的使用(转载)
  14. .NET中栈和堆的比较
  15. vue 基于网易云API的短信验证码登录(axios封装)
  16. GPS测量中的数据格式——RINEX格式说明
  17. Servlet 请求转发
  18. Kibana关联ES查询数据
  19. filebeat7.7.0相关详细配置预览- processors - add_fields
  20. 降噪蓝牙耳机对比测评:南卡和OPPO哪款降噪耳机好?

热门文章

  1. MaxProxy代理可以替代911s5代理吗?还有哪些软件可以使用
  2. Go语言开发第1课-环境搭建及简单程序入门
  3. leetcode + 腾讯面试
  4. 图文讲解flashfxp教程 flashfxp详细教程
  5. 火影手游服务器维护多久,《火影忍者》手游结束终极内测的停服公告
  6. 使命召唤手游显示服务器停服,使命召唤手游停服了吗 是手游还是端游
  7. C# HttpClientHelper
  8. [系统安全] 二十五.WannaCry勒索病毒分析 (1)Python复现永恒之蓝漏洞实现勒索加密
  9. “U盘写保护”什么解除
  10. 电影文件的合并与分割