仿真求助

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


本帖完毕
回帖
  • 等待你,消灭零评论
科研通是完全免费的文献互助平台,具备全网最快的应助速度,最高的求助完成率。 对每一个文献求助,科研通都将尽心尽力,给求助人一个满意的交代。
实时播报
刚刚
Jeffrey完成签到 ,获得积分10
刚刚
小蘑菇应助科研通管家采纳,获得10
刚刚
在水一方应助科研通管家采纳,获得10
刚刚
刚刚
1秒前
1秒前
无风发布了新的文献求助10
2秒前
4秒前
刻苦小凝发布了新的文献求助10
5秒前
坚果燕麦发布了新的文献求助10
5秒前
晴枫3648完成签到,获得积分10
5秒前
糖丸完成签到,获得积分10
8秒前
FashionBoy应助无风采纳,获得10
9秒前
科研通AI6.1应助csj采纳,获得10
10秒前
在水一方应助七页禾采纳,获得10
11秒前
科学实验站完成签到,获得积分10
11秒前
12秒前
12秒前
YY完成签到 ,获得积分10
13秒前
13秒前
phenory发布了新的文献求助10
14秒前
汉堡包应助余铸海采纳,获得10
14秒前
初景发布了新的文献求助10
14秒前
球球尧伞耳完成签到,获得积分10
15秒前
16秒前
krabs完成签到,获得积分10
16秒前
lemon完成签到 ,获得积分10
18秒前
涛1118完成签到,获得积分10
18秒前
lmhytr完成签到,获得积分10
18秒前
yuki完成签到,获得积分10
18秒前
xny发布了新的文献求助10
19秒前
19秒前
21秒前
芝士雪豹发布了新的文献求助30
21秒前
一颗葡萄完成签到 ,获得积分10
21秒前
英姑应助念l采纳,获得10
22秒前
linn发布了新的文献求助10
22秒前
23秒前
汉堡包应助流星雨采纳,获得10
24秒前
热门帖子
关注 科研通微信公众号,转发送积分 6403900
求助须知:如何正确求助?哪些是违规求助? 8222932
关于积分的说明 17427862
捐赠科研通 5456380
什么是DOI,文献DOI怎么找? 2883487
邀请新用户注册赠送积分活动 1859773
关于科研通互助平台的介绍 1701151
最新评论
学习啦加油 3小时前
学习了 同学们 5小时前
我是第一次用,好像是违规了,现在发布不了文献,能帮忙处理一下吗,以后会尽量多的提供信息,拜托了 6小时前
学习喽 7小时前
学习了 11小时前
学习啦 11小时前
评论会有积分吗 13小时前
学习学习 13小时前