hdlcjd.asv

混沌matlab源代码 实际例子,仅供参考

clear
n=432;
t=2;
for m=1:15
  h=m;   
x=[11.89 	5.86 	3.29 	6.14 	6.64 	64.02 	40.18 	98.83 	107.05 	140.35 	57.54 	28.66... 
20.62 	6.14 	3.96 	2.26 	3.32 	113.01 	161.78 	771.38 	557.28 	407.12 	60.91 	32.41...  
30.53 	13.92 	9.29 	16.87 	15.78 	22.52 	98.57 	214.00 	73.61 	127.76 	37.58 	33.75... 
15.64 	11.79 	4.93 	16.41 	17.92 	37.84 	154.81 	433.90 	482.11 	199.27 	62.99 	32.68...  
24.88 	8.11 	5.14 	5.42 	5.46 	33.44 	142.76 	281.23 	181.96 	122.94 	29.03 	11.78...  
7.15 	5.31 	3.59 	5.70 	15.53 	83.46 	221.24 	210.52 	48.47 	34.82 	37.58 	49.01...  
16.31 	4.85 	3.37 	2.49 	5.06 	44.84 	99.90 	358.91 	80.61 	117.05 	55.21 	25.82...  
19.28 	6.35 	6.32 	6.56 	6.03 	62.47 	121.87 	138.47 	116.64 	32.41 	31.36 	14.52...  
6.54 	3.86 	3.16 	2.93 	10.31 	64.28 	122.94 	184.27 	233.54 	60.00 	192.07 	45.26...  
13.02 	8.53 	9.19 	9.82 	44.46 	223.95 	380.33 	621.39 	285.12 	84.64 	30.59 	12.56...  
13.02 	8.01 	4.96 	4.59 	9.78 	111.20 	77.67 	171.69 	97.98 	34.02 	41.73 	9.29...  
8.09 	6.07 	5.30 	3.97 	5.62 	29.55 	47.94 	109.55 	197.77 	34.82 	17.65 	9.72...  
6.54 	4.46 	2.87 	4.20 	1.23 	7.08 	40.71 	136.33 	114.05 	45.53 	14.39 	6.67...  
4.82 	4.12 	1.82 	3.27 	6.70 	277.34 	130.17 	228.74 	177.81 	86.24 	15.47 	9.78...  
6.62 	5.16 	3.46 	3.03 	2.14 	41.21 	64.01 	116.51 	239.76 	152.13 	20.30 	11.57...  
8.49 	4.43 	2.17 	2.67 	5.54 	13.12 	21.56 	72.85 	34.47 	73.39 	8.50 	5.14...  
3.48 	2.41 	1.23 	1.11 	3.00 	48.99 	126.42 	129.37 	103.68 	15.64 	6.97 	5.22...  
4.37 	3.97 	2.09 	3.37 	1.31 	25.48 	51.96 	80.89 	54.17 	85.17 	11.61 	6.86...  
10.45 	3.73 	6.59 	2.62 	1.90 	17.76 	26.17 	145.97 	147.74 	110.35 	60.13 	16.98...  
18.94 	2.93 	1.98 	4.38 	7.45 	54.43 	149.72 	66.69 	58.84 	18.96 	6.66 	4.50...  
3.62 	2.41 	2.06 	3.76 	34.55 	146.19 	235.97 	190.70 	246.50 	79.01 	22.97 	9.19...  
11.46 	5.16 	3.56 	4.07 	3.83 	7.62 	310.69 	366.94 	316.22 	428.54 	58.84 	20.68...  
24.67 	12.52 	5.52 	6.40 	3.46 	32.92 	60.80 	119.99 	311.04 	103.65 	12.39 	16.55...  
12.37 	4.33 	3.40 	4.30 	3.67 	37.07 	30.80 	79.28 	103.94 	23.19 	11.51 	4.61...  
2.41 	1.24 	1.77 	2.80 	0.70 	7.54 	17.84 	71.51 	41.21 	47.14 	7.18 	5.52...  
3.08 	0.88 	1.55 	1.19 	27.86 	206.84 	380.33 	229.54 	118.45 	194.99 	22.37 	20.89...  
12.94 	3.27 	1.47 	3.76 	4.42 	61.17 	81.96 	220.43 	329.18 	133.65 	91.76 	23.06...  
15.83 	9.59 	8.04 	3.37 	19.50 	22.01 	68.83 	52.76 	54.43 	85.71 	27.73 	13.26...  
10.79 	6.27 	3.80 	2.62 	10.02 	13.27 	35.62 	126.42 	191.03 	51.43 	11.28 	6.45...  
3.91 	4.38 	3.54 	2.46 	7.63 	68.69 	210.52 	130.97 	106.27 	14.22 	8.04 	8.54...  
6.24 	1.74 	1.04 	1.66 	3.64 	78.80 	149.45 	257.13 	241.06 	177.31 	105.75 	29.73...  
13.50 	7.44 	6.59 	3.11 	15.40 	74.39 	164.99 	256.59 	98.76 	25.10 	34.47 	14.44...  
8.04 	5.05 	2.60 	2.70 	2.46 	51.06 	256.59 	182.67 	277.34 	218.56 	24.88 	13.07...  
9.29 	5.60 	5.73 	4.46 	2.38 	316.22 	624.07 	332.12 	154.48 	45.80 	41.47 	18.64...  
23.06 	8.89 	3.43 	1.89 	15.29 	41.21 	167.67 	514.25 	425.09 	62.14 	123.64 	37.23...  
21.37 	23.87 	22.44 	4.22 	32.94 	99.01 	283.91 	261.41 	247.54 	129.10 	63.24 	21.35];	

for i=1:(n-(m-1)*t)
    for j=1:((m-1)*t+1)
        y(i,j)=x(i+j-1);
    end
end
r(n-(m-1)*t,n-(m-1)*t)=0;
for i=1:(n-(m-1)*t)
    for j=1:(n-(m-1)*t)
        for k=1:(m-1)*t+1
            r(i,j)=r(i,j)+(y(i,k)-y(j,k))^2;
        end
    end
end
R=sqrt(r);
dmax=0;
for i=1:(n-(m-1)*t)
    for j=1:(n-(m-1)*t)
        if R(i,j)>dmax 
         dmax=R(i,j);
        end
    end
end
d=dmax
for i=1:15
    r0(h,i)=i*dmax/15;
end
  
for k=1:15
    b=0;
    for i=1:(n-(m-1)*t)
        a=0;
        for j=1:(n-(m-1)*t)
            a=r0(h,k)-R(i,j);
            if a>0 b=b+1;
            end
        end
    end
    C(h,k)=b/((n-(m-1)*t)^2);
end
end
c=C;
w=r0;