matlab求解时滞微分方程「建议收藏」

matlab求解时滞微分方程「建议收藏」matlab求解时滞微分方程,dde23调用格式:sol=dde23(ddefun,lags,history,tspan);–ddefun函数句柄,求解微分方程y’=f(t,y(t),y(t-τ1),…,y(t-τk))必须写成下面形式:dydt=ddefun(t,y,Z);其中t对应当前时间t,y为列向量,近似于y(…

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

Jetbrains全系列IDE稳定放心使用

matlab求解时滞微分方程,dde23调用格式:

sol = dde23(ddefun,lags,history,tspan);

–ddefun函数句柄,求解微分方程y=f(t,y(t),y(t-τ1),…,y(t-τk))

        必须写成下面形式:

        dydt =ddefun(t,y,Z);

        其中t对应当前时间t,y为列向量,近似于y(t);Z(:,j)近似于y(t-τj)

–lags为延迟时间,为正常数。

            例:方程中包含y1(t-0.2)和y2(t-1),则可以表示为lags=[0.2,1]

–history t≤t0 时的状态变量的值

–tspan 积分区间 t0 = tspan(1),tf =tspan(end)。

 

看下面例子:

      假定系统状态方程为dXdt =- Ax(t) – Bx(t-0.23)+Bx(t-0.56); A= [3,-1,-1,-1;-1,2,0,-1;-1,0,2,-1;-1,-1,-1,3];B=A;

 

程序如下:

     (1) 编写延迟函数

function dx = ddefun(t,y,Z)
A = [3,-1,-1,-1;-1,2,0,-1;-1,0,2,-1;-1,-1,-1,3]; tau1= Z(:,1); tau2= Z(:,2); dx=-A*y-A*tau1-A*tau2
;

 

 (1) 编写主调函数

 

 

tau = [0.23,0.56];
y0 = [1,7,3,0];
sol = dde23(@ddefun, tau, y0, [0, 50]);

% plot the system states

plot(sol.x,sol.y);

*注意:该函数返回的sol中结构体sol.x和sol.y均为按行排列,与ode45等不同
如只显示一组数据plot(sol.x,sol.y(1,:));

 

matlab求解时滞微分方程「建议收藏」

 

  

转载于:https://www.cnblogs.com/rainrainbow/p/3958132.html

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

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

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

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

(0)
blank

相关推荐

  • apache中文乱码_文件名称乱码怎么解决

    apache中文乱码_文件名称乱码怎么解决RestSharp是一个第三方开源的Http模拟请求辅助类,其底层实现基于System.Net.HttpWebRequest,且不依赖于任何第三方控件。其github地址为:https://github.com/restsharp/RestSharp,start数可以说明该类库的知名度,当然侧面也可以证明它的确是一个比较好用的HTTP请求辅助类。一般情况下,RestSharp都工作的很好,只是当…

  • Shell循环读取文件

    Shell循环读取文件注:部分概念介绍来源于网络方法1:while循环中执行效率最高,最常用的方法。functionwhile_read_LINE(){whilereadLINEdoecho$LINEdone<$FILENAME}注释:这种方式在结束的时候需要执行文件,就好像是执行完的时候再把文件读进去一样。方法2:重定向法;管道法:cat$FILENAME|whilereadLINEfunctionWhile_read_LINE(){cat$FILENAME|..

  • Java学习之Spring框架入门篇

    Java学习之Spring框架入门篇0x00前言最近在学SSM的三大框架,上篇文章,已经更新了Mybatis的相关内容,那么这篇文章就来写写Spring的入门。0x01Spring概述S

    2021年12月12日
  • php实现微信扫码登陆,PHP实现微信开放平台扫码登录源码下载[通俗易懂]

    php实现微信扫码登陆,PHP实现微信开放平台扫码登录源码下载[通俗易懂]1、首先到微信开放平台申请https://open.weixin.qq.com/获取到appid和APPSECRET,前台显示页面如下varobj=newWxLogin({id:”login_container”,appid:”wxed782be999f86e0e”,scope:”snsapi_login”,redirect_uri:encodeURIComponent(“htt…

  • atof函数_log函数怎么比较大小

    atof函数_log函数怎么比较大小atof函数原型doubleatof(constchar*str);作用:将字符串转换为双精度浮点数(double).头文件:#include&lt;stdlib.h&gt;返回值: 返回转换后的浮点数,如果字符串str不能被转换为double,那么返回0.0函数说明:atof()会扫描茶树str字符串,跳过前面的空格字符,直到遇到数字或者正…

    2022年10月22日
  • UAVStack之文件数据归集「建议收藏」

    UAVStack之文件数据归集「建议收藏」UAVStack之文件数据归集

发表回复

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

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