清晨好,您是今天最早来到科研通的研友!由于当前在线用户较少,发布求助请尽量完整地填写文献信息,科研通机器人24小时在线,伴您科研之路漫漫前行!

仿真求助

搜索技巧 2286 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
2秒前
牛仔完成签到 ,获得积分10
6秒前
chichenglin完成签到 ,获得积分0
6秒前
传奇3应助ybwei2008_163采纳,获得10
8秒前
FashionBoy应助ybwei2008_163采纳,获得10
8秒前
Lucas应助科研通管家采纳,获得10
9秒前
活力的珊完成签到 ,获得积分10
15秒前
YZY完成签到 ,获得积分10
21秒前
27秒前
量子星尘发布了新的文献求助10
28秒前
盘尼西林发布了新的文献求助10
30秒前
周大福完成签到 ,获得积分10
30秒前
qianci2009完成签到,获得积分0
36秒前
Gary完成签到 ,获得积分10
41秒前
Enyiqi001完成签到 ,获得积分10
48秒前
Fe_Al_Po完成签到,获得积分0
52秒前
52秒前
melody完成签到 ,获得积分10
57秒前
57秒前
貔貅完成签到 ,获得积分10
1分钟前
mictime完成签到,获得积分10
1分钟前
研友_Ze2oV8完成签到 ,获得积分10
1分钟前
Orange应助淡水美人鱼采纳,获得10
1分钟前
1分钟前
研友_LN25rL完成签到,获得积分10
1分钟前
junjie完成签到 ,获得积分10
1分钟前
lzy完成签到,获得积分10
1分钟前
1分钟前
Lucas应助盘尼西林采纳,获得10
1分钟前
1分钟前
ybwei2008_163发布了新的文献求助10
1分钟前
1分钟前
yyyyy发布了新的文献求助30
1分钟前
话说dota完成签到 ,获得积分10
1分钟前
zhangxiaoqing完成签到,获得积分10
1分钟前
GankhuyagJavzan完成签到,获得积分10
1分钟前
cadcae完成签到,获得积分10
2分钟前
李健应助科研通管家采纳,获得10
2分钟前
科研通AI6.4应助yyyyy采纳,获得20
2分钟前
隐形荟完成签到 ,获得积分10
2分钟前
热门帖子
关注 科研通微信公众号,转发送积分 6066511
求助须知:如何正确求助?哪些是违规求助? 7898785
关于积分的说明 16322787
捐赠科研通 5208390
什么是DOI,文献DOI怎么找? 2786268
邀请新用户注册赠送积分活动 1769013
关于科研通互助平台的介绍 1647813
最新评论
这个网站真的好方便,有求必应,速度极快,文献下载再也不愁了,推荐大家收藏使用。 8小时前
非常棒,帮大忙了 9小时前
办法不多,,, 13小时前
非常棒 14小时前
00哈哈哈哈好坏 15小时前
非常好的平台 15小时前
非常好 20小时前
养龙虾的人太多了,把阿里云的资源挤兑的厉害,我们也被影响了 22小时前