Map 映射

clear all;

format rat

clc;

x=0.5;

x1=[x];

tx1=[0];

mu0=1.5;

for i=1:200

if x<0.5

tx=mu0*x;

else

tx=(1-x)*mu0;

end

x1=[x1 x tx ];

tx1=[tx1 tx tx];

x=tx;

xout(i)=x;

end

plot(x1,tx1)

axis([0 1 0 1])

hold on

line([0 1],[0 1])

line([0 0.5],[0 0.75])

line([1 0.5],[0 0.75])

xlabel(x)

ylabel(T_x)