仿真求助

搜索技巧 2294 0
jiao
jiao 来自贵州 发表于: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 的值


本帖完毕
回帖
  • 等待你,消灭零评论
科研通是完全免费的文献互助平台,具备全网最快的应助速度,最高的求助完成率。 对每一个文献求助,科研通都将尽心尽力,给求助人一个满意的交代。
实时播报
世佳何完成签到,获得积分10
刚刚
阿司匹林完成签到,获得积分10
1秒前
2秒前
kroll完成签到,获得积分10
2秒前
orixero应助tys0713104采纳,获得10
3秒前
4秒前
SDNUDRUG完成签到,获得积分10
5秒前
Rosaline完成签到 ,获得积分10
5秒前
蝎子莱莱启动完成签到,获得积分10
7秒前
可靠的南露完成签到,获得积分10
7秒前
sevenhill应助科研通管家采纳,获得10
8秒前
sevenhill应助科研通管家采纳,获得10
8秒前
科研通AI2S应助科研通管家采纳,获得10
8秒前
lyf完成签到 ,获得积分10
11秒前
量子星尘发布了新的文献求助10
11秒前
11秒前
悦耳含灵完成签到,获得积分10
11秒前
含蓄的孤丝完成签到 ,获得积分10
12秒前
瘦瘦的枫叶完成签到 ,获得积分10
13秒前
SCI硬通货完成签到 ,获得积分10
15秒前
LeOpard完成签到,获得积分10
17秒前
岁月旧曾谙完成签到,获得积分10
17秒前
19秒前
文献博士完成签到 ,获得积分10
20秒前
南攻完成签到,获得积分10
20秒前
O_O完成签到 ,获得积分10
21秒前
端庄代荷完成签到 ,获得积分10
22秒前
不能玩一下午吗应助kroll采纳,获得10
23秒前
麻雀完成签到 ,获得积分10
24秒前
无奈山雁完成签到 ,获得积分10
25秒前
無迪发布了新的文献求助10
25秒前
义气天真完成签到,获得积分10
27秒前
尊敬亦寒完成签到,获得积分10
28秒前
zcm1999完成签到,获得积分10
28秒前
WW完成签到 ,获得积分10
32秒前
hdhuang完成签到,获得积分10
33秒前
canghong完成签到,获得积分10
35秒前
cong完成签到 ,获得积分10
35秒前
忽晚完成签到 ,获得积分10
36秒前
量子星尘发布了新的文献求助10
38秒前
热门帖子
关注 科研通微信公众号,转发送积分 6066665
求助须知:如何正确求助?哪些是违规求助? 7898967
关于积分的说明 16322911
捐赠科研通 5208397
什么是DOI,文献DOI怎么找? 2786304
邀请新用户注册赠送积分活动 1769013
关于科研通互助平台的介绍 1647813
最新评论
这个网站真的好方便,有求必应,速度极快,文献下载再也不愁了,推荐大家收藏使用。 10小时前
非常棒,帮大忙了 10小时前
办法不多,,, 14小时前
非常棒 16小时前
00哈哈哈哈好坏 16小时前
非常好的平台 16小时前
非常好 22小时前
养龙虾的人太多了,把阿里云的资源挤兑的厉害,我们也被影响了 23小时前