【调剂】北京石油化工学院2024年16个专业接受调剂
查看: 2415  |  回复: 6
【悬赏金币】回答本帖问题,作者ygj2015将赠送您 15 个金币

ygj2015

新虫 (小有名气)

[求助] 离散时滞系统的matlab实现已有1人参与

离散时滞系统
x(k+1)=Ax(k)+Bx(k-d(k))

matlab实现:
A1=[-0.4 -0.3;0.1 0.15];
A2=[0.1 0.25;-0.9 -0.1];
N=150;
x=zeros(2,N+1);
x(:,1)=[-0.1;0.1];
for k=1:N
        d=[k-(0.1-0.1*sin(k))]
        x(:,k+1)=A1*x(:,k)+A2*x(:,k-d);
    end
figure
    tout=0:N;
    x1=x(1, ;
   x2=x(2, ;
    plot(tout,x1,tout,x2)  
     xlabel('t');
    ylabel('x(t)');
    legend('x_1(t)','x_2(t)','location','northeast')
运行报错如下:
位置 2 的索引无效。数组索引必须为正整数或逻辑值。

出错 delayUntitled2 (line 44)
    x(:,k+1)=A2*x(:,k)+Ad2*x(:,k-d)

该如何修改此程序,望高人指点。
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hrbnu0828

木虫 (正式写手)

k-d一定是个大于0的整数?

发自小木虫IOS客户端
人不能因为担心不能够永远而拒绝开始
2楼2018-10-15 01:11:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

googon

木虫 (著名写手)


【答案】应助回帖

感谢参与,应助指数 +1
你的这个问题一看就清楚明白,d不可能都是正整数,因为你用了sin,建议你用四舍五入法处理此问题,还有d不能大于k,索引值不能为负,需要更精细的解答,请Q搜matlab编程爱好者

发自小木虫Android客户端
3楼2018-10-16 07:09:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ygj2015

新虫 (小有名气)

引用回帖:
3楼: Originally posted by googon at 2018-10-16 07:09:36
你的这个问题一看就清楚明白,d不可能都是正整数,因为你用了sin,建议你用四舍五入法处理此问题,还有d不能大于k,索引值不能为负,需要更精细的解答,请Q搜matlab编程爱好者
...

感谢回复,明白了,谢谢啊
4楼2018-10-16 08:28:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ygj2015

新虫 (小有名气)

引用回帖:
3楼: Originally posted by googon at 2018-10-16 07:09:36
你的这个问题一看就清楚明白,d不可能都是正整数,因为你用了sin,建议你用四舍五入法处理此问题,还有d不能大于k,索引值不能为负,需要更精细的解答,请Q搜matlab编程爱好者
...

按照四舍五入处理后,运行通了,但是图中的初值和理论对不上,不知道什么原因?
5楼2018-10-16 23:14:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

googon

木虫 (著名写手)


【答案】应助回帖

引用回帖:
5楼: Originally posted by ygj2015 at 2018-10-16 23:14:18
按照四舍五入处理后,运行通了,但是图中的初值和理论对不上,不知道什么原因?...

加群详谈,把相关文献发给我

发自小木虫Android客户端
6楼2018-10-17 06:51:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ygj2015

新虫 (小有名气)

引用回帖:
6楼: Originally posted by googon at 2018-10-17 06:51:24
加群详谈,把相关文献发给我
...

已经运行通过了,谢谢你。
7楼2018-10-18 10:05:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ygj2015 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] with editor 两个月了,什么原因? +6 yiersan9 2024-04-24 12/600 2024-04-26 02:18 by yiersan9
[论文投稿] Chemical Engineering Journal投稿3周了,一直显示With editor状态。这是送审了吗? 10+4 yifeng11 2024-04-20 11/550 2024-04-25 22:19 by f798710538
[找工作] 普通院校药学硕士,做合成的,感觉找不到工作 +8 pom戴墨镜 2024-04-24 15/750 2024-04-25 20:50 by pom戴墨镜
[论文投稿] Nature一直在编辑手里,考虑好几天了,是悬了吗 +12 彩虹初见 2024-04-24 12/600 2024-04-25 19:21 by 雪径踏青
[教师之家] 刚博士毕业,入职普通一本高校 +18 保持肌肤 2024-04-20 19/950 2024-04-25 18:30 by lyfbangong
[论文投稿] 催稿后秒拒 +10 lizhengke06 2024-04-19 12/600 2024-04-25 17:50 by 追分筝
[考博] 取博导收留 5+3 zzb777888 2024-04-20 8/400 2024-04-25 17:40 by 那年的旧城
[基金申请] 基金开始函评了吗? +11 wych1103 2024-04-25 11/550 2024-04-25 16:54 by 未岸亡灵
[教师之家] 博士论文被抄袭 +23 和尚敲小木鱼 2024-04-22 39/1950 2024-04-25 15:48 by 和尚敲小木鱼
[硕博家园] 聊天 +11 暮色恋伊人 2024-04-22 12/600 2024-04-25 13:53 by UCTS
[考博] 真的好想读博! +14 wangzhe_bs 2024-04-22 17/850 2024-04-25 11:36 by 庭前花未开
[博后之家] 南京大学-广州大学联合招聘博士后 欢迎广大优秀人才!!! +4 黑魔变身啾 2024-04-20 12/600 2024-04-25 11:18 by dodonaomi
[电化学] 耗材发问 +4 Happy C 2024-04-22 4/200 2024-04-25 11:03 by 普通小虫
[考研] 0854-0855调剂 +7 shangannum1 2024-04-21 10/500 2024-04-25 01:31 by 啊廖sh
[考研] 博士生招生!! +5 2458827884 2024-04-19 5/250 2024-04-24 22:28 by 不懂的风清
[基金申请] 化学结清有情发出来了 +3 starboy7286 2024-04-20 3/150 2024-04-24 15:52 by 嘿,黑贝
[教师之家] 河海大学常州校区新进的讲师待遇如何?考核要求如何? +3 德国小球迷 2024-04-19 18/900 2024-04-22 19:23 by ZHONGWU_U
[考博] 青岛科技大学顾玉超教授课题组招收博士生 +3 yanqingwang 2024-04-19 3/150 2024-04-22 15:21 by 理科生yes
[论文投稿] 编辑是选国外的好还是国内的好。 +8 lizhengke06 2024-04-20 8/400 2024-04-22 08:58 by cuiyunjian
[硕博家园] 材料工程,二本硕士,值得读吗 +3 简风23 2024-04-20 6/300 2024-04-21 12:19 by 简风23
信息提示
请填处理意见