仿真求助

搜索技巧 117 0
jiao
jiao 1 来自贵州 发表于: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
1秒前
liam发布了新的文献求助10
1秒前
1秒前
李健应助yw采纳,获得10
2秒前
结实的以冬完成签到,获得积分10
3秒前
ff发布了新的文献求助20
3秒前
科研通AI5应助丫丫采纳,获得10
3秒前
归尘发布了新的文献求助10
3秒前
科研通AI5应助迅速的网络采纳,获得10
3秒前
4秒前
zzz发布了新的文献求助10
6秒前
好运来完成签到,获得积分10
6秒前
爆米花应助黄黄大可爱采纳,获得10
6秒前
传奇3应助玩命的猫咪采纳,获得10
7秒前
上官若男应助zjy03259采纳,获得10
8秒前
明亮若枫完成签到,获得积分10
8秒前
zyue发布了新的文献求助10
9秒前
小蘑菇应助YuLu采纳,获得10
9秒前
安息香发布了新的文献求助20
9秒前
10秒前
bkagyin应助背后的傥采纳,获得10
11秒前
华仔应助chiahaokuo采纳,获得10
12秒前
猪猪hero应助科研通管家采纳,获得10
13秒前
尤瑞米应助科研通管家采纳,获得10
13秒前
13秒前
充电宝应助科研通管家采纳,获得10
13秒前
深情安青应助科研通管家采纳,获得10
13秒前
猪猪hero应助科研通管家采纳,获得10
13秒前
13秒前
FashionBoy应助科研通管家采纳,获得10
13秒前
猪猪hero应助科研通管家采纳,获得10
13秒前
猪猪hero应助科研通管家采纳,获得10
13秒前
赘婿应助拾意采纳,获得10
13秒前
科研通AI5应助科研通管家采纳,获得30
13秒前
迟大猫应助科研通管家采纳,获得10
13秒前
13秒前
13秒前
14秒前
天真书南完成签到,获得积分10
15秒前
科目三应助桃核采纳,获得10
16秒前
热门帖子
关注 科研通微信公众号,转发送积分 3590184
求助须知:如何正确求助?哪些是违规求助? 3158575
关于积分的说明 9520512
捐赠科研通 2861581
什么是DOI,文献DOI怎么找? 1572636
邀请新用户注册赠送积分活动 738020
科研通“疑难数据库(出版商)”最低求助积分说明 722636
最新评论
真好用,效率太高啦 11小时前
是什么差不多啊,是峰的形状还是峰面积,有截图看看 13小时前
也可以是一个浓度的菌,最后的cDNA梯度稀释,然后qPCR,用ct值和od生成标曲,就是简单的线性拟合。qPCR引物就选择细菌的管家基因,也就是普通做qPCR时的内参。 13小时前
非常好用感谢 1天前
好的谢谢 1天前
好滴 谢谢站长回复 1天前
正常文献求助,在最下面选择“仅需补充材料”,得到的便是了 1天前
把没贴壁的洗下去留贴壁的继续长不就好了么,这个能说是长满了么?这个是刚复苏的细胞把,应该状态不好,或者再养养看会不会贴壁的多一些,想换液可以离心一下把细胞还留着吧 1天前