MATLAB绘制Logistic混沌映射的分叉图

Logistic混沌映射函数:

x ( i + 1 ) = u ⋅ x ( i ) ⋅ [ 1 − x ( i ) ] , ( i = 0 , 1 , … , n ) x(i+1) = u · x(i) ·[1 - x(i)],(i=0,1,…,n) x(i+1)=ux(i)[1x(i)],i=0,1,,n)

其中 0 ≤ x ( i ) ≤ 1 , u ≤ 4 0 \leq x(i) \leq1,u \leq4 0x(i)1u4

clc,clear,close all;
x0 = 0.5;
u = 0:0.001:4;  %自定义步长改变迭代次数
xn = ft_logistic(x0,length(u)-1,u);

function [yn] = ft_logistic(x0,num_len,u)
    xn = x0;
    for j = 1:num_len
        u_i = u(j);
        for i = 1:j
            xn(i+1) = u_i*xn(i).*(1-xn(i));
        end
        plot(u_i,xn(j),'b.',"linewidth",1.2)
        hold on
        grid on
        yn(1,j) = xn(j);
    end
end

或将绘图函数移出循环速度更快:

clc,clear,close all;
x0 = 0.5;
u = 0:0.001:4;  %自定义步长改变迭代次数
xn = ft_logistic(x0,length(u)-1,u);
plot(u(2:end),xn,'b.',"linewidth",1.2)
grid on

function [yn] = ft_logistic(x0,num_len,u)
    xn = x0;
    for j = 1:num_len
        u_i = u(j);
        for i = 1:j
            xn(i+1) = u_i*xn(i).*(1-xn(i));
        end
        yn(1,j) = xn(j);
    end
end

在这里插入图片描述
在这里插入图片描述

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

Wayne_Fine

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值