混沌保密通信-C语言编程最佳范例

《混沌保密通信-C语言编程最佳范例》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、#include#include#include#defineMM33554432#defineA1103515245#defineC12345#defineRA1.0intx1=21,x2=1,next=11;intshengchengsuijishu1(){intr;x2=(A*x1+C)%MM;r=(A*x1)%MM;x1=x2;return(r);}doubleshengchengsuijishu(){next=shengchengsuijishu1();next=next*A+C;next=(in

2、t)(next/256)%65536;return(fabs((double)next/65536));}voidlorenz(doublesig,doubleb,doubler,double*x){inti;doubley[3];for(i=0;i<3;i++)y[i]=*(x+i);*x=-sig*(y[0]-y[1]);*(x+1)=r*y[0]-y[1]-y[0]*y[2];*(x+2)=y[0]*y[1]-b*y[2];}voidrun_kut(doublesig,doubleb,doubler,double*x,doublet){inti;double

3、y[4][3];for(i=0;i<3;i++)y[0][i]=*(x+i);lorenz(sig,b,r,y[0]);for(i=0;i<3;i++)y[1][i]=*(x+i)+t*y[0][i]/2;lorenz(sig,b,r,y[1]);for(i=0;i<3;i++)y[2][i]=*(x+i)+t*y[1][i]/2;lorenz(sig,b,r,y[2]);for(i=0;i<3;i++)y[3][i]=*(x+i)+t*y[2][i];lorenz(sig,b,r,y[3]);for(i=0;i<3;i++)*(x+i)+=t*(y[0][i]+

4、2*y[1][i]+2*y[2][i]+y[3][i])/6;}doubleffunction(intn,double*param,double*mx1,doublet){inti;doublerr=0.0,sig,b,r,x[3],sx1[600];sig=10;b=8/3.0;r=*(param);x[0]=*mx1;sx1[0]=*mx1;x[1]=*(param+1);x[2]=*(param+2);for(i=1;i

5、x1+i)-sx1[i]);returnrr;}voidant_rea(double*x,double*x1){//*x=(30+0.0)*((*x)+50)/100.0+0.0;//30//*(x+1)=(10+0.0)*(*(x+1)+50)/100.0+0.0;//10*(x+1)=(10+10.0)*(*(x+1)+50)/100.0-10.0+(*(x1+1));*(x+2)=(10.0+10.0)*(*(x+2)+50)/100.0-10.0+(*(x1+2));*(x)=(50+0.0)*(*(x)+50)/100.0-0.0+(*(x1));}vo

6、idrea_ant(double*x,double*x1){//*x=100*((*x)-0.0)/(30+0.0)-50.0;//*(x+1)=100*(*(x+1)-0.0)/(10+0.0)-50.0;*(x+1)=100*(*(x+1)+10.0-(*(x1+1)))/(10+10.0)-50.0;*(x+2)=100*(*(x+2)+10.0-(*(x1+2)))/(10+10.0)-50.0;*(x)=100*(*(x)+0.0-(*(x1)))/(50+0.0)-50.0;}voidant_opt(intNN,double*mx1,doublet,d

7、ouble*x1){constintM=1000,N=3,K=10,S=1000;inti,j,k,js,qd,jishu=0;doublea=200,b=0.5,f,tm=0.0,midle,range=50,fi=0.075,r[M],v[M],x[M][N],p[M][N],p1[M][N],y[M],fitness[M],ex[N],para[3];//charname[50];for(i=0;i<3;i++)para[i]=*(x1+i);para[0]=0.0;for(i=0;i

8、0;;qd