matlab求解延迟微分方程_状态依赖时滞微分方程的动力学研究

matlab求解延迟微分方程_状态依赖时滞微分方程的动力学研究固定时滞的微分方程:满足下面的形式,也就是微分方程右边包含时滞部分,且时滞为常数。使用dde23函数求解:问题:       (1)微分方程定义:多了一个时滞部分创建myddefun.m文件,文件里的内容如下:functiondy=myddefun(t,y,Z)dy=[Z(1,1);Z(1,1)+Z(2,2);

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

Jetbrains全系列IDE稳定放心使用

固定时滞的微分方程:满足下面的形式,也就是微分方程右边包含时滞部分,且时滞为常数。

matlab求解延迟微分方程_状态依赖时滞微分方程的动力学研究

使用dde23函数求解:

问题:

matlab求解延迟微分方程_状态依赖时滞微分方程的动力学研究            matlab求解延迟微分方程_状态依赖时滞微分方程的动力学研究

(1)微分方程定义:多了一个时滞部分

创建myddefun.m文件,文件里的内容如下:

function dy = myddefun(t,y,Z)
dy=[
    Z(1,1);
    Z(1,1)+Z(2,2);
    y(2) ];

说明:其中Z表示时滞部分:即y(t-T),T是时滞。

Z(i,j)表示y(i)(t-T),即y(i)的时滞形式;j表示T选取第j个时滞值

(2)外部调用方程,输入参数求解

lags=[1,0.2];
history=[1;1;1];
tspan=[0,5];
sol = dde23(@myddefun,lags,history,tspan)

其中lags表示时滞数组,如上面输入了2个时滞,Z(1,1)表示选用第一个时滞T=1,Z(2,2)表示选取第二个时滞T=0.2;

history表示变量的初值,上面代码的意思是t=0时,y=[1,1,1];

tspan表示解的范围,即t的范围,上面表示求t在[0,5]范围内y的解。

(3)绘制解的图

plot(sol.x,sol.y) 

结果:

matlab求解延迟微分方程_状态依赖时滞微分方程的动力学研究


版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 举报,一经查实,本站将立刻删除。

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

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

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

(0)
blank

相关推荐

发表回复

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

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