仿真求助

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


本帖完毕
回帖
  • 等待你,消灭零评论
科研通是完全免费的文献互助平台,具备全网最快的应助速度,最高的求助完成率。 对每一个文献求助,科研通都将尽心尽力,给求助人一个满意的交代。
实时播报
zhenhan发布了新的文献求助10
1秒前
1秒前
ana完成签到,获得积分10
2秒前
nav发布了新的文献求助10
2秒前
2秒前
2秒前
Nankdream完成签到,获得积分10
2秒前
king48615发布了新的文献求助20
2秒前
Nole应助csj采纳,获得30
3秒前
妮妮完成签到 ,获得积分10
4秒前
Khoilism完成签到,获得积分10
5秒前
111发布了新的文献求助10
5秒前
5秒前
Heron完成签到,获得积分10
6秒前
Rai完成签到,获得积分10
6秒前
冰冰发布了新的文献求助10
6秒前
Akim应助NJK采纳,获得10
6秒前
6秒前
6秒前
yueyue完成签到,获得积分10
7秒前
7秒前
7秒前
7秒前
dyq完成签到,获得积分20
7秒前
科研通AI6.2应助小超人采纳,获得10
8秒前
acale完成签到,获得积分10
8秒前
kwangil发布了新的文献求助10
11秒前
acale发布了新的文献求助10
12秒前
12秒前
JamesPei应助游阿游采纳,获得10
12秒前
矮小的行云完成签到,获得积分10
13秒前
13秒前
13秒前
CodeCraft应助Jackson_Cai采纳,获得10
14秒前
影子完成签到 ,获得积分10
14秒前
14秒前
研友_VZG7GZ应助SILICON采纳,获得10
15秒前
15秒前
爱笑的猪zz完成签到,获得积分10
15秒前
16秒前
热门帖子
关注 科研通微信公众号,转发送积分 7702561
求助须知:如何正确求助?哪些是违规求助? 9261044
关于积分的说明 20030071
捐赠科研通 7278224
什么是DOI,文献DOI怎么找? 3294245
关于科研通互助平台的介绍 2449697
邀请新用户注册赠送积分活动 2300900
最新评论
赞,稳定用了好几年了 16分钟前
非常好用 1小时前
非常好用 16小时前
你可以应助获得积分,很快的 18小时前
你好,我才看到你这个问题,解决了。我回忆一下,应该是要将沉淀进行真空的完全干燥,我是用干燥桶的负压进行干燥的,干燥一个小时左右到沉淀出现干涸的那种裂痕,然后加入到乙二醇中间去90℃搅拌就行了 21小时前
我的灵感是撰写一篇scoping review,尽管也是属于review类型,但是,是掌握了一种新的较为严谨的研究方法,让我觉得我从事的工作还是具有一定的学术门槛的。 21小时前
非常好用 21小时前
我也想问 23小时前