[转载]matlab中求时滞常微分方程组
matlab中求时滞常微分方程组有两个函数，一个是dde23函数，求解时滞函数为常数的情形；另一个是ddesd函数，求解时滞函数为一般函数的情形。下边主要讲dde23函数，ddesd可以参考help文档。
以控制系统中的时滞常微分方程组为例：
代码如下：
function y=controsys()
miu=0.1;beta=0.2;omiga=1;xita=0.3;tao=5;%参数的赋值
sc=[0.6,0.2,0.2];%设定初始值
for i=5:5:100
sol=dde23(@eq,[1 1],@(t) eqhist(t,sc),[i-5 i]);
%eq是求解的方程组;eqhist是定义的时滞函数的值[i-tao i];
1]是S,I的时滞参数，[i-5 i]是求解区间
sc=sol.y(:,length(sol.y));
sc(3,:)=sc(1,:)*xita+sc(3,:);sc(1,:)=sc(1,:)*(1-xita);sc=sc';
%求解新的初始值
plot(sol.x,sol.y')
on%画出函数图象
&&& plot([i
i],[sol.y(:,length(sol.y)) sc'])
on%连接跳跃点
end%求解跳跃点
function s=eqhist(t,sc)
function dy=eq(t,y,Z)
miu=0.1;beta=0.2;omiga=1;xita=0.3;tao=5;
dy=[0.6 0.2 0.2];
Slag=Z(1);
Ilag=Z(2);
dy=[miu-beta*y(1).*y(2)-miu*y(1);
beta*y(1)*y(2)-beta*exp(-miu*omiga)*Slag*Ilag-miu*y(2);
beta*exp(-miu*omiga)*Slag*Ilag-miu*y(3)];
以上网友发言只代表其个人观点，不代表新浪网的观点或立场。典型混沌系统的Matlab仿真实现
The Simulation Reality of Typical Chaotic Systems with Matlab
简要介绍了三种典型的混沌系统:Logistic映射,Lorenz系统和Henon映射,利用Matlab软件,使用编程方法进行了仿真.
Wang Gaiyun
Ma Shujing
桂林电子科技大学计算机与控制学院,541004
年，卷(期)
机标分类号
广西区教育厅立项项目，桂林电子科技大学学科软环境项目
本文读者也读过
加载中，请稍候
．客服电话
&&8:00-11:30,13:00-17:00(工作日)查看: 60|回复: 0|关注: 0
求教大神，lorenz方程参数如果是跟x有关的变量怎么办
<h1 style="color:# 麦片财富积分
新手, 积分 5, 距离下一级还需 45 积分
发表于 7&天前
]{WB)Z]M]QFQ@11V5]Y]UZF.png (14.69 KB, 下载次数: 0)
7&天前 上传
上面这个是lorenz方程。下面这个是跟lorenz方程类似的一个方程中我取的bet,gam,sig值，但是这样编写程序是不对的，求大神帮助！！！万分感谢！
$P[LDZE_)Z0OBCUG~YZ(7.png (7.22 KB, 下载次数: 0)
7&天前 上传
站长推荐 /3
使用MATLAB和Simulink进行软件无线电设计
Powered by&#xe621; 上传我的文档
&#xe602; 下载
&#xe60c; 收藏
该文档贡献者很忙，什么也没留下。
&#xe602; 下载此文档
正在努力加载中...
基于Pyragas时滞反馈的统一Lorenz系统的混沌控制研究_[定稿]
下载积分：1300
内容提示：基于Pyragas时滞反馈的统一Lorenz系统的混沌控制研究_[定稿]
文档格式：PDF|
浏览次数：0|
上传日期： 05:57:02|
文档星级：&#xe60b;&#xe612;&#xe612;&#xe612;&#xe612;
全文阅读已结束，如果下载本文需要使用
&#xe71b; 1300 积分
&#xe602;下载此文档
该用户还上传了这些文档
基于Pyragas时滞反馈的统一Lorenz系统的混沌控制研究_
官方公共微信}