仿真求助

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


本帖完毕
回帖
  • 等待你,消灭零评论
科研通是完全免费的文献互助平台,具备全网最快的应助速度,最高的求助完成率。 对每一个文献求助,科研通都将尽心尽力,给求助人一个满意的交代。
实时播报
可爱的函函应助ljy采纳,获得10
刚刚
司空元正完成签到,获得积分10
刚刚
FrozNineTivus完成签到,获得积分10
刚刚
学术渣渣发布了新的文献求助10
1秒前
钱途无量发布了新的文献求助20
1秒前
1秒前
糊涂的猕猴桃完成签到,获得积分20
2秒前
Ada1989完成签到,获得积分10
2秒前
penxyy应助blind采纳,获得30
2秒前
2秒前
2秒前
3秒前
eaton完成签到 ,获得积分10
3秒前
3秒前
4秒前
zqzqz完成签到,获得积分10
4秒前
找啊找完成签到,获得积分10
4秒前
4秒前
月亮发布了新的文献求助10
4秒前
传奇3应助清新的半烟采纳,获得10
5秒前
TQQT发布了新的文献求助10
5秒前
科研通AI2S应助代纤绮采纳,获得10
5秒前
ironsilica发布了新的文献求助10
5秒前
hhh完成签到 ,获得积分10
6秒前
LX发布了新的文献求助10
6秒前
fenghua完成签到,获得积分20
6秒前
zdw完成签到,获得积分10
6秒前
易方完成签到,获得积分10
6秒前
7秒前
7秒前
7秒前
8秒前
不是细菌完成签到,获得积分20
8秒前
能干冰菱完成签到,获得积分10
8秒前
可爱的函函应助gong采纳,获得10
8秒前
所所应助shuo采纳,获得10
8秒前
勇敢牛牛发布了新的文献求助10
8秒前
9秒前
希望天下0贩的0应助asd采纳,获得10
9秒前
9秒前
热门帖子
关注 科研通微信公众号,转发送积分 6067510
求助须知:如何正确求助?哪些是违规求助? 7899493
关于积分的说明 16326750
捐赠科研通 5209257
什么是DOI,文献DOI怎么找? 2786461
邀请新用户注册赠送积分活动 1769287
关于科研通互助平台的介绍 1647853
最新评论
这个网站真的好方便,有求必应,速度极快,文献下载再也不愁了,推荐大家收藏使用。 12小时前
非常棒,帮大忙了 12小时前
办法不多,,, 16小时前
非常棒 18小时前
00哈哈哈哈好坏 18小时前
非常好的平台 18小时前
非常好 1天前
养龙虾的人太多了,把阿里云的资源挤兑的厉害,我们也被影响了 1天前