仿真求助

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


本帖完毕
回帖
  • 等待你,消灭零评论
更新
PDF的下载单位、IP信息已删除 (2025-6-4)

科研通是完全免费的文献互助平台,具备全网最快的应助速度,最高的求助完成率。 对每一个文献求助,科研通都将尽心尽力,给求助人一个满意的交代。
实时播报
圈圈黄发布了新的文献求助30
1秒前
渔夫完成签到,获得积分10
4秒前
天空之城完成签到,获得积分10
4秒前
5秒前
是我不得开心妍完成签到 ,获得积分10
5秒前
英姑应助77采纳,获得10
6秒前
shengwang完成签到,获得积分10
7秒前
的法国队完成签到,获得积分10
7秒前
安白发布了新的文献求助10
8秒前
9秒前
啊阿阿阿沐完成签到,获得积分10
10秒前
11秒前
7udo发布了新的文献求助10
12秒前
treasure完成签到 ,获得积分10
12秒前
迷人幻巧发布了新的文献求助10
14秒前
董帅发布了新的文献求助10
15秒前
不知月明是故乡完成签到 ,获得积分10
16秒前
一彤完成签到,获得积分10
17秒前
英俊的铭应助风清扬采纳,获得10
18秒前
18秒前
小蘑菇应助One采纳,获得10
18秒前
小丹小丹完成签到 ,获得积分10
18秒前
19秒前
7udo完成签到,获得积分10
21秒前
77发布了新的文献求助10
22秒前
Zhusy发布了新的文献求助30
22秒前
wills应助迷人幻巧采纳,获得10
22秒前
彭于晏应助迷人幻巧采纳,获得10
22秒前
朴素羊完成签到 ,获得积分10
23秒前
bkagyin应助楼梯口无头女孩采纳,获得10
24秒前
CipherSage应助杠赛来采纳,获得10
25秒前
26秒前
Zero完成签到,获得积分10
27秒前
hiter完成签到,获得积分10
29秒前
番茄炒蛋完成签到,获得积分10
29秒前
一只小学弱完成签到,获得积分10
30秒前
董帅完成签到,获得积分10
30秒前
风清扬发布了新的文献求助10
30秒前
kk完成签到 ,获得积分10
32秒前
未碎冰蓝完成签到,获得积分20
33秒前
热门帖子
关注 科研通微信公众号,转发送积分 5565699
求助须知:如何正确求助?哪些是违规求助? 4650686
关于积分的说明 14692512
捐赠科研通 4592693
什么是DOI,文献DOI怎么找? 2519716
邀请新用户注册赠送积分活动 1492102
关于科研通互助平台的介绍 1463316
最新评论
感谢!这是非常有用的改进 10小时前
感谢平台 12小时前
感谢平台 13小时前
还好吧 21小时前
请多一点包容,大家都是从不会过来的 22小时前
可能确实行不通,也可能没思考 1天前
感谢平台 1天前
我还以为我导师问的 吓死我了 1天前