仿真求助

搜索技巧 2317 0
jiao
jiao 47 来自贵州 发表于:2025-03-19 19:55:29

 

有没有大哥大姐帮忙看一下这个matlab代码为什么是错的,其中k不进行循环,但是打印出来它是迭代的,这是为什么呢?求求了

Xi = [1.5, 0;0 ,1.2];

mu=0.1;    %时间步长

T=100; % 最大时间步

N=100;%仿真步数

t=0:mu:(N-1)*mu;%时间向量 

                          % 更新状态

                          fprintf('k = %d\n', k); % 打印 k 的值

for period=1:1

  for k = 6:T-2

        f1=@(k)0.1*tanh(k);% 定义激活函数

        f2=@(k)0.2*tanh(k);% 定义激活函数

        u1 =0;                      %定义输入

        u2 =[0.2*cos(k),0.1*sin(k)];%定义输入

%%%%%%%%%%%%%%%%

    c1 =@(k) 1.90 + 0.80*tanh(k);

    c2 =@(k) 1.85+0.45*tanh(2*k);

    d1 =@(k) 1 + tanh(k);

    d2 =@(k) 1.10 + 0.42*tanh(2*k);

    a11=@(k) 0.15 + 0.08 * sin(k);

    a12=@(k) 0.10 + 0.17 * sin(2*k);

    a21=@(k) 0.16 + 0.10 * sin(3*k);

    a22=@(k) 0.35 + 0.01 * sin(4*k);

    b11=@(k) 0.17 + 0.02 * cos(k);

    b12=@(k) 0.10 + 0.20 * cos(2*k);

    b21=@(k) 0.15 + 0.50 * cos(3*k);

    b22=@(k) 0.06 + 0.13 * cos(4*k);

%%%%%%%%%%% 

    Cshang=[c1(k),0;0,c2(k)];

    Dshang=[d1(k),0;0,d2(k)];

    Ashang=[a11(k),a12(k);a21(k),a22(k)];

    Bshang=[b11(k),b12(k);b21(k),b22(k)];

%%%%%%%%%%%

    Cxing=Cshang-2*eye(2);

    Dxing=Dshang+eye(2)-Cshang;

    Axing=Ashang;

    Bxing=Bshang;

            x=zeros(2,N);

            h1=1;    %定义时滞h(k)

            h2=3;

            % h= randi([1, min(5,k-1)]);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

            % x1_delay = k-n1;  % 确保索引至少为 1

            % x2_delay = k-n2;

   x(1,k+2)= -Cxing(1,1)*x(1,k+1)-Dxing(1,1)*x(1,k) ...

            +Axing(1,1)*f1(x(1,k))+Axing(1,2)*f2(x(2,k)) ...

            +Bxing(1,1)*f1(x(1,k-h1))+Bxing(1,2)*f2(x(2,k-h2))+u1;

 

   x(2,k+2)= -Cxing(2,2)*x(2,k+1)-Dxing(2,2)*x(2,k) ...

            +Axing(2,1)*f1(x(1,k))+Axing(2,2)*f2(x(2,k)) ...

            +Bxing(2,1)*f1(x(1,k-h1))+Bxing(2,2)*f2(x(2,k-h2))+u1;

 

     fprintf('k = %d\n', k); % 打印 k 的值

  end 

  %subplot(2, 1, 1);

  plot(t, x(1, 1:N), 'c', 'LineWidth', 2);hold on;

  plot(t, x(2, 1:N), 'r', 'LineWidth', 2);grid on;

  xlabel('time t');

  ylabel('x(t)');

  title('x_1(t) and x_2(t) state of u1');

end

  

       fprintf('k = %d\n', k); % 打印 k 的值


本帖完毕
回帖
  • 等待你,消灭零评论
科研通是完全免费的文献互助平台,具备全网最快的应助速度,最高的求助完成率。 对每一个文献求助,科研通都将尽心尽力,给求助人一个满意的交代。
实时播报
zzmyyds发布了新的文献求助10
刚刚
1秒前
suxuan完成签到,获得积分10
1秒前
慕青应助Spring采纳,获得10
1秒前
LYJ发布了新的文献求助10
2秒前
Rong0618发布了新的文献求助10
2秒前
广发牛勿完成签到,获得积分10
2秒前
2秒前
CC完成签到,获得积分20
2秒前
巩志成完成签到,获得积分10
2秒前
思源应助Bminor采纳,获得10
2秒前
3秒前
YC完成签到 ,获得积分10
3秒前
科研通AI6.4应助哈哈采纳,获得15
4秒前
震动的医学生完成签到,获得积分20
5秒前
蛋黄啵啵完成签到 ,获得积分10
5秒前
认真自行车完成签到 ,获得积分10
5秒前
虎桔发布了新的文献求助10
6秒前
学术垃圾完成签到,获得积分10
6秒前
亚历山大关注了科研通微信公众号
6秒前
biubiubiu发布了新的文献求助10
6秒前
丘比特应助妮妮采纳,获得30
7秒前
zl_tf完成签到,获得积分10
7秒前
希望天下0贩的0应助114514采纳,获得30
8秒前
8秒前
8秒前
达克赛德发布了新的文献求助10
8秒前
小蘑菇应助咔嚓采纳,获得10
9秒前
无限彤完成签到,获得积分20
9秒前
9秒前
Unfair发布了新的文献求助10
9秒前
上官若男应助细腻的易真采纳,获得10
10秒前
JamesPei应助fxs采纳,获得10
10秒前
NexusExplorer应助科研通管家采纳,获得10
10秒前
xxxy应助科研通管家采纳,获得10
10秒前
完美世界应助科研通管家采纳,获得10
10秒前
深情安青应助科研通管家采纳,获得10
10秒前
顾矜应助科研通管家采纳,获得10
10秒前
酷波er应助frl采纳,获得10
10秒前
情怀应助科研通管家采纳,获得10
11秒前
热门帖子
关注 科研通微信公众号,转发送积分 6070005
求助须知:如何正确求助?哪些是违规求助? 7901866
关于积分的说明 16335449
捐赠科研通 5210951
什么是DOI,文献DOI怎么找? 2787112
邀请新用户注册赠送积分活动 1769943
关于科研通互助平台的介绍 1648020
最新评论
非常便捷实用的文献互助平台,解决了很多论文下载不了的情况,希望平台越来越好 1小时前
搜一下5730图书馆,注册个会员,从知网入口进去下载 2小时前
已捐赠,祝越来越好! 5小时前
这个网站真的好方便,有求必应,速度极快,文献下载再也不愁了,推荐大家收藏使用。 18小时前
非常棒,帮大忙了 18小时前
办法不多,,, 22小时前
非常棒 1天前
00哈哈哈哈好坏 1天前