MATLAB 数学应用 微分方程 时滞微分方程 具有常时滞的DDE「建议收藏」

MATLAB 数学应用 微分方程 时滞微分方程 具有常时滞的DDE「建议收藏」本文讲述了如何使用dde23对具有常时滞的DDE(时滞微分方程)方程组求解。方程组为:y1′(t)=y1(t−1)y’_1(t)=y_1(t−1)y1′​(t)=y1​(t−1)y2′(t)=y1(t−1)+y2(t−0.2)y’_2(t)=y_1(t-1)+y_2(t-0.2)y2′​(t)=y1​(t−1)+y2​(t−0.2)y3′(t)=y2(t)y’_3(t)=y_2(t)y3′​(t)=y2​(t).t≤0的历史解函数是常量y1(t)=y2(t)=y3(t)=1y_1(t)=y

大家好,又见面了,我是你们的朋友全栈君。如果您正在找激活码,请点击查看最新教程,关注关注公众号 “全栈程序员社区” 获取激活教程,可能之前旧版本教程已经失效.最新Idea2022.1教程亲测有效,一键激活。

Jetbrains全系列IDE稳定放心使用

本文讲述了如何使用 dde23 对具有常时滞的DDE(时滞微分方程)方程组求解。

方程组为:

y 1 ′ ( t ) = y 1 ( t − 1 ) y’_1(t)=y_1(t−1) y1(t)=y1(t1)

y 2 ′ ( t ) = y 1 ( t − 1 ) + y 2 ( t − 0.2 ) y’_2(t)=y_1(t-1)+y_2(t-0.2) y2(t)=y1(t1)+y2(t0.2)

y 3 ′ ( t ) = y 2 ( t ) y’_3(t)=y_2(t) y3(t)=y2(t).

t≤0 的历史解函数是常量 y 1 ( t ) = y 2 ( t ) = y 3 ( t ) = 1 y_1(t)=y_2(t)=y_3(t)=1 y1(t)=y2(t)=y3(t)=1

方程中的时滞仅存在于 y 项中,并且时滞本身是常量,因此各方程构成常时滞方程组。

要在 MATLAB 中求解此方程组,您需要先编写方程组、时滞和历史解的代码,然后再调用时滞微分方程求解器 dde23,该求解器适用于具有常时滞的方程组。您可以将所需的函数作为局部函数包含在文件末尾,或者将它们作为单独的命名文件保存在 MATLAB 路径上的目录中。

编写时滞代码

首先,创建一个向量来定义方程组中的时滞。此方程组有两种不同时滞:

  • 在第一个分量 y 1 ( t − 1 ) y_1(t−1) y1(t1) 中时滞为 1。
  • 在第二个分量 y 2 ( t − 0.2 ) y_2(t−0.2) y2(t0.2) 中时滞为 0.2。

dde23 接受时滞的向量参数,其中每个元素是一个分量的常时滞。

lags = [1 0.2];

编写方程代码

现在,创建一个函数来编写方程的代码。此函数应具有签名 dydt = ddefun(t,y,Z),其中:

t 是时间(自变量)。

y 是解(因变量)。

Z(:,j) 用于逼近时滞 y ( t − τ j ) y(t−τ_j) y(tτj),其中常时滞 τ j τ_j τj 由 lags(j) 给定。

求解器会自动将这些输入传递给该函数,但是变量名称决定如何编写方程代码。在这种情况下:

  • Z ( : , 1 )   →     y 1 ( t − 1 ) Z(:,1) →  y_1(t−1) Z(:,1)y1(t1)
  • Z ( : , 2 )   →     y 2 ( t − 0.2 ) Z(:,2) →  y_2(t−0.2) Z(:,2)y2(t0.2)
function dydt = ddefun(t,y,Z)
  ylag1 = Z(:,1);
  ylag2 = Z(:,2);

  dydt = [ylag1(1); 
          ylag1(1)+ylag2(2); 
          y(2)];
end

注意:所有函数都作为局部函数包含在示例的末尾。

编写历史解代码

接下来,创建一个函数来定义历史解。历史解是时间 t ≤ t 0 t≤t_0 tt0 的解。

function s = history(t)
  s = ones(3,1);
end

求解方程

最后,定义积分区间 [ t 0     t f ] [t_0  t_f] [t0tf] 并使用 dde23 求解器对 DDE 求解。

tspan = [0 5];
sol = dde23(@ddefun, lags, @history, tspan);

对解进行绘图

解结构体 sol 具有字段 sol.x 和 sol.y,这两个字段包含求解器在这些时间点所用的内部时间步和对应的解。(如果您需要在特定点的解,可以使用 deval 来计算在特定点的解。)

绘制三个解分量对时间的图。

plot(sol.x,sol.y,'-o')
xlabel('Time t');
ylabel('Solution y');
legend('y_1','y_2','y_3','Location','NorthWest');

在这里插入图片描述

局部函数

此处列出了 DDE 求解器 dde23 为计算解而调用的局部辅助函数。您也可以将这些函数作为它们自己的文件保存在 MATLAB 路径上的目录中。

function dydt = ddefun(t,y,Z) % equation being solved
  ylag1 = Z(:,1);
  ylag2 = Z(:,2);

  dydt = [ylag1(1); 
          ylag1(1)+ylag2(2); 
          y(2)];
end
%-------------------------------------------
function s = history(t) % history function for t <= 0
  s = ones(3,1);
end
%-------------------------------------------
版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 举报,一经查实,本站将立刻删除。

发布者:全栈程序员-用户IM,转载请注明出处:https://javaforall.cn/188171.html原文链接:https://javaforall.cn

【正版授权,激活自己账号】: Jetbrains全家桶Ide使用,1年售后保障,每天仅需1毛

【官方授权 正版激活】: 官方授权 正版激活 支持Jetbrains家族下所有IDE 使用个人JB账号...

(0)


相关推荐

  • devtools工具如何使用_devtool制作插件

    devtools工具如何使用_devtool制作插件7devtool快速参考目录7devtool快速参考7.1获得帮助7.2工作区层结构7.3向工作区层添加新配方7.4提取现有配方的来源7.5同步一个配方的提取源树7.6修改现有配方7.7编辑现有配方7.8更新配方7.9查看配方升级状态7.10升级配方7.11重置配方7.12建立你的配方7.13建立你的形象7.14在目标机器上部署你的软件7.15从目标机器上删除您的软件7.16在替代位置创建工作空间层7.17获取工作区中配方的状态

  • 详解 Java 内部类

    详解 Java 内部类前言内部类在Java里面算是非常常见的一个功能了,在日常开发中我们肯定多多少少都用过,这里总结一下关于Java中内部类的相关知识点和一些使用内部类时需要注意的点。从种类上说,内部类可以分为四类:普通内部类、静态内部类、匿名内部类、局部内部类。我们来一个个看:普通内部类这个是最常见的内部类之一了,其定义也很简单,在一个类里面作为类的一个字段直接定义就可以了,例:pub…

    2022年10月11日
  • form layui vue 和_layui是基于vue的吗?「建议收藏」

    form layui vue 和_layui是基于vue的吗?「建议收藏」layui不是基于vue的。layui是一款绝对开源和免费的前端UI框架,遵循原生“HTML/CSS/JS”的书写与组织形式;它虽然外在极简,但是内容丰富,里面包含众多组件从核心代码到API都非常适合界面的快速开发。layui不是基于vue的。layui(谐音:类UI)是一款采用自身模块规范编写的前端UI框架,遵循原生HTML/CSS/JS的书写与组织形式,门槛极低,拿来即用。lay…

  • java读写锁实现原理_java可重入锁原理

    java读写锁实现原理_java可重入锁原理一、synchronized和ReentrantLock的对比到现在,看到多线程中,锁定的方式有2种:synchronized和ReentrantLock。两种锁定方式各有优劣,下面简单对比一下:1、synchronized是关键字,就和if…else…一样,是语法层面的实现,因此synchronized获取锁以及释放锁都是Java虚拟机帮助用户完成的;ReentrantLock是类层面的…

  • java实现xml文件CRUD

    java实现xml文件CRUD

    2021年12月31日
  • [Python] ImportError: DLL load failed … 找不到指定的模块 此类问题解决方法

    [Python] ImportError: DLL load failed … 找不到指定的模块 此类问题解决方法最近升级Python项目,由Python2.7升级到Python3.8.3,项目使用了PySide2,对于较新的Python3.8.3,PySide2可能存在些许不兼容问题,环境配置完成后,出现一连串的ImportError:DLLloadfailed找不到指定模块对于很多Python开发者来说,这类问题最为头疼,不知道如何下手解决。我是在virtualenv虚拟环境下配置Python3.8.3的开发环境。OS:Windows7x64各种依赖安装完毕后

发表回复

您的电子邮箱地址不会被公开。

关注全栈程序员社区公众号