rationallogist1.m

这么代码是用Matlab编写的

clear all;
key0=0.1;
key1=4;
M=256;
N=256;
for(i=1:M)
    key0=key1*key0*(1-key0);
    a(i)=key0;
end
key3=0.2;
key2=4;
M=256;
N=256;
for(i=1:M)
    key3=key2*key3*(1-key3);
    b(i)=key3;
end
key4=0.3;
key2=4;
M=256;
N=256;
for(i=1:M)
    key4=key2*key4*(1-key4);
    c(i)=key4;
end
t=0.4;
w0=0.2;
w1=0.5;
w2=0.3;
w=(1-t)^2*w0+2*t*(1-t)*w1+t^2*w2;
for(i=1:M)
    P(i)=(1-t)^2*a(i)*w0+2*t*(1-t)*b(i)*w1+t^2*c(i)*w2;
    d(i)=P(i)/w;
end
% figure(1)
% plot(a);
% figure(2)
% plot(b);
% figur
plot(d);