| 查看: 1825 | 回复: 4 | ||
星志.Ding木虫 (小有名气)
|
[求助]
（十万火急）切换退化时滞系统（带约束条件的微分方程）的matlab作图程序
|
|
系统 E * (dx/dt) = Ai*x(t) + Bi*x(t-h) (i=1,2,即两个子系统) 其中矩阵E=[1 0;0 0] Ai，Bi 是任意实矩阵，比如取A1=[-8 2;1 -4] ，B1=[2 0;3 -1]，A2=[-6 -1;0 -5] ，B2=[-4 -2;1 -2] h=0.01是时滞. 切换时间序列任意给，比如：[0,50]运行第一个系统，[51,150]运行第二个系统。。。。即两个系统交替运行就行 方程也可以写成分量的形式如下（实际上就是一个带约束条件的微分方程作图程序）： 第一个系统： dx1/dt = -8*x1(t) + 2*x2(t) + 2*x1(t-h) 0 = x1(t) - 4*x2(t) + 3*x1(t-h) - x2(t-h) 第二个系统： dx1/dt = -6*x1(t) - x2(t) - 4*x1(t-h) - 2*x2(t-h) 0 = - 5*x2(t) + x1(t-h) - 2*x2(t-h) 其中x1表示x的第一个分量，x2表示x的第二个分量 跪求，谢谢了，实在不能做切换系统的，能够做出一个系统的运行图形性可以！感谢各位好心人！！ |
回复此楼» 猜你喜欢
科研孤儿太难了
已经有9人回复
-
时间戳又变了
已经有15人回复
-
时间戳变了，能看出什么问题？
已经有23人回复
-
只有每年这种时候来逛逛小木虫
已经有22人回复
-
今天系统多次维护，明天很可能放榜！
已经有10人回复
-
放榜前的不淡定
已经有12人回复
-
言之凿凿，明天官网会通知放榜？
已经有16人回复
-
重要消息，中午系统在维护
已经有12人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有20人回复
-
今天基金会出结果吗？20260819
已经有15人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 关于Matlab作图的问题
已经有14人回复
- 【【吐血求助】】Matlab微分方程作图-报错报错
已经有4人回复
- 科研作图时候matlab中使用latex
已经有3人回复
soliton923: 违规存档, 别灌水啊~~学科版没有沙发 2012-11-27 19:41:50
 |
2楼2012-11-27 16:42:24
aaron1988
木虫 (正式写手)
- 应助: 15 (小学生)
- 金币: 2261.1
- 红花: 3
- 帖子: 938
- 在线: 135.2小时
- 虫号: 413419
- 注册: 2007-06-26
- 专业: 化工系统工程
【答案】应助回帖
★ ★ ★
星志.Ding(sweety代发): 金币+3, 鼓励应助 2013-07-01 13:21:56
星志.Ding(sweety代发): 金币+3, 鼓励应助 2013-07-01 13:21:56
|
1. 这种情况劝你用Simulink来做。 有现成的杂化系统模块（但不是matlab原有的），你去找pws-piecewise system toolbox或者hybrid system。 2. 如果你硬要写m-file的话，有几种方法： i. 离散化用idpoly或者相关的命令。你会得到y（z)=G(z)u(z)这样的形式的传递方程。然后你定义u就是时间。要响应的时间为1，不要就为0. ii. 把两个系统分别写成 function， 然后用ode来解，具体点。就是每到下个时间点就解一次ODE 我感觉还缺少不少条件，所以不能给出更具体的方案了 |
3楼2013-06-30 00:23:55
pippi6
铁杆木虫 (著名写手)
工程和科学数值计算咨询
- 数学EPI: 6
- 应助: 413 (硕士)
- 贵宾: 0.002
- 金币: 7116.5
- 散金: 15
- 红花: 63
- 帖子: 1639
- 在线: 798.9小时
- 虫号: 2469437
- 注册: 2013-05-14
- 专业: 计算数学与科学工程计算
4楼2013-07-03 09:44:20
youandiandhe
铁杆木虫 (知名作家)
- 应助: 25 (小学生)
- 金币: 13706.4
- 散金: 5605
- 红花: 29
- 帖子: 8071
- 在线: 2670.7小时
- 虫号: 3483828
- 注册: 2014-10-18
- 专业: 控制理论与方法
5楼2017-05-21 14:00:54
40