毕业设计（论文）-基于小波变换的图像去噪方法研究报告附MATLAB程序.docx

基于小波变换的图像去噪方法的研究摘要图像在生成和传输的过程中会受到各种噪声的干扰,对信息的处理、传输和存储造成极大的影响。寻求一种既能有效地减小噪声,又能很好地保留图像边缘信息的方法,是人们一直追求的目标。小波分析是局部化时频分析,它用时域和频域联合表示信号的特征,是分析非平稳信号的有力工具。它通过伸缩、平移等运算功能对信号进行多尺度细化分析,能有效地从信号中提取信息。本文对基于小波变换的图像去噪方法进行了研究分析, 详细介绍了其去噪原理和算法,分析了去噪过程中参数的选取问题,并给出了一些选取依据;详细介绍了小波系数相关性去噪方法的原理和算法。最后对均值滤波、中值滤波和维纳滤波方法在高斯噪声下进行了分析比较,并给出了仿真实验结果，结果证明小波去噪十分有效，其结果好于其它3种滤波。关键词：小波变换，图像去噪，阈值，阈值函数1引言数字图像在我们日常生活中起着非常重要的作用，它与我们的日常生活息息相关，例如在卫星、电视、核磁共振、计算机视觉、地球信息系统以及天文学中应用非常广泛。但是一般情况下采集到的数字图像是含有噪声的。噪声1可以理解为“妨碍人们感觉器官对所接收的信源信息理解的因素”。图像在生成和传输的过程中灰受到各种噪声的干扰，对信号的处理、传输和存储造成极大的影响。数字图像之所以含有噪声这是因为在图像的采集、获取、编码和传输的过程中，所有的图像均不同程度地被可见或不可见的噪声“污染”。对于这种“污染”，如果信噪比(snr)低于一定水平，就会影响图像场景内容的表示，直接导致图像质量的下降。除了视觉质量上下降外，噪声还可能掩盖一些重要的图像细节，使图像的熵增大，从而对于图像数据的有效压缩起到了一定的妨碍作用。对于图像在采集、获取过程造成的“污染”，我们虽然尽量提高硬件设备以获取质量更高的图像，但图像传感器的截止频率总是有一定的，受硬件水平和价格的限制，且图像在编码和传输过程中造成的“污染”，必需采取有效的降噪技术才能提高图像的质量。对图像进行去噪最初主要是在空域内进行的，图像空域去噪方法很多，主要是通过各种滤波器对图像进行去噪。例如均值滤波器、顺序统计滤波器、维纳滤波器等。为了进一步提高去噪的效果，在变换域中进行降噪处理成为有效的方法，图像变换域去噪就是对图像进行某一种变换，然后将图像从时域变换到变换域中，再对变换域中的图像变换系数按照某种方法进行处理，最后再对处理后的系数按照某种方法进行反变换，这样就实现了将图像去除图像噪声的目的。将图像从时域转换到变换域的变换方法很多，例如傅立叶变换、小波变换等等。小波变换是在短时傅立叶变换的基础上发展起来的一种新型的变换方法。小波变换具有多分辨率分析的特点，在时域、频域都具有较强的表征信号局部特征的能力，因此基于小波分析的图像去噪技术已成为图像去噪的一个重要方法。本文采用小波阈值去噪的方法，从去噪的效果上比较了多种去噪方法的优劣，实验证明小波去噪在图像噪声处理中起到很好的效果。2小波变换概述2.1 小波变化去噪技术研究现状上个世纪八十年代mallet 提出了 mra(multi_resolution analysis)，并首先把小波理论运用于信号和图像的分解与重构，利用小波变换模极大值原理进行信号的奇异性检测，提出了交替投影算法用于信号重构，为小波变换用于图像处理奠定了基础1。后来，人们根据信号与噪声在小波变换下模极大值在各尺度上的不同传播特性，提出了基于模极大值去噪的基本思想。1992年，donoho和johnstone提出了“小波收缩”，它较传统的去噪方法效率更高。“小波收缩”被 donoho和johnstone证明是在极小化极大风险中最优的去噪方法，但在这种方法中最重要的就是确定阈值。1995年，stanford大学的学者d.l.donoho和i.m.johnstone提出了通过对小波系数进行非线性阈值处理来降低信号中的噪声2。从这之后的小波去噪方法也就转移到从阈值函数的选择或最优小波基的选择出发来提高去噪的效果。影响比较大的方法有以下这么几种：eero p.semoncelli和edward h.adelson提出的基于最大后验概率的贝叶斯估计准则确定小波阈值的方法3；elwood t.olsen等在处理断层图像时提出了三种基于小波相位的去噪方法：边缘跟踪法、局部相位方差阈值法以及尺度相位变动阈值法；学者kozaitis结合小波变换和高阶统计量的特点提出了基于高阶统计量的小波阈值去噪方法4；g.p.nason等利用原图像和小波变换域中图像的相关性用gcv(general cross- validation)法对图像进行去噪；hang.x和woolsey等人提出结合维纳滤波器和小波阈值的方法对信号进行去噪处理5，vasily strela等人将一类新的特性良好的小波(约束对)应用于图像去噪的方法6；同时，在19世纪60年代发展的隐马尔科夫模型(hidden markov model)，是通过对小波系数建立模型以得到不同的系数处理方法；后又有人提出了双变量模型方法7，它是利用观察相邻尺度间父系数与子系数的统计联合分布来选择一种与之匹配的二维概率密度函数。这些方法均取得了较好的效果，对小波去噪的理论和应用奠定了一定的基础。小波具有低墒性、多分辨率、去相关性、选基灵活性等特点。2.2 连续小波变换1. 小波基函数所谓小波(wavelet)，即存在于一个较小区域的波。小波函数的数学定义是：设为一平方可积函数，即，若其傅立叶变换满足：（2.1）时，则称为一个基本小波或小波母函数，并称上式是小波函数的可容许条件。根据小波函数的定义，小波函数一般在时域具有紧支集或近似紧支集，即函数的非零值定义域具有有限的范围，这即所谓“小”的特点;另一方面，根据可容许性条件可知，即直流分量为零，因此小波又具有正负交替的波动性。将小波母函数进行伸缩和平移，设其伸缩因子(亦称尺度因子)为，平移因子为，并记平移伸缩后的函数为，则:（2.2）并称为参数和小波基函数。由于和均取连续变换的值，因此又称为连续小波基函数，它们是由同一母函数经伸缩和平移后得到的一组函数系列。2. 连续小波变换将空间的任意函数在小波基下进行展开，称其为函数的连续小波变换cwt，变换式为：（2.3）当小波的容许性条件成立时，其逆变换为：（2.4）其中为的容许性条件。另外，在小波变换过程中必须保持能量成比例，即:（2.5）由cwt的定义可知，小波变换和傅立叶变换一样，也是一种积分变换，其中为小波变换系数。可见小波变换对函数在小波基上的展开具有多分辨率的特性，这种特性正是通过缩放因子和平移因子来得到的。一个一维函数的连续小波变换是一双变量的函数，变量比多一个，因此称连续小波变换是超完备的，因为它要求的存储量和它代表的信息量都显著增加了。对于变量超过一个的函数来说，这个变换的维数也将增加。若是一个二维函数，则它的连续小波变换是：（2.6）其中，表示在两个维度上的平移，二维连续小波逆变换为：（2.7）同样的方法可以推广到两个或两个以上的变量函数上。2.3 离散小波变换计算机中的图像信息是以离散信号形式存放的，所以需要将连续小波变换离散化。而最基本的离散化方法就是二进制离散，一般将这种经过离散化的小波及其变换叫做二进小波和二进变换。需要注意的是这里的离散化都是针对连续的尺度因子和连续平移因子的，而不是针对时间的。这儿限制尺度因子总是正数。1. 尺度与位移的离散化对连续小波基函数尺度因子和平移因子进行离散化可以得到离散小波变换，从而减少小波变换系数的冗余度。在离散化时通常对尺度因子和平移因子按幂级数进行离散化，即取（为整数，但一般都假定），得到离散小波函数为：（2.8）其对应系数为：（2.9）2. 二进制小波变换二进小波变换是一种特殊的离散小波变换，特别地令参数，则有。该二进尺度分解的原理在二十世纪三十年代由 littlewood 和 paley 在数学上进行了研究证明。离散小波变换为：（2.10）离散二进小波变换为：（2.11）二维离散小波变换：我们考虑二维尺度函数是可分离的情况，也就是：（2.12）设是与对应的一维小波函数，则有：（2.13）（2.14）（2.15）以上三式就建立了二维小波变换的基础。3. 图像的小波变换及mallat算法图像是二维信号，二维多分辨率分析与一维类似，但空间变成，一维中引入的尺度函数变为。设是的一个多分辨率分析，则可以证明，张量空间：构成的一个多分辨率分析，并且二维多分辨率分析的二维尺度函数为 (2.16)式中：是尺度函数(一维)。式(2.28)说明了二维尺度函数的可分离性。对于每一个,函数系构成的规范正交基，这里下标j，n，m的含义是：（2.17）我们将称为的可分离多分辨率分析。因、都是低通的尺度函数，所以平滑的低通空间。如果是一维多分辨率分析的正交小波基，则二维多分辨率分析的三个小波函数为：（2.18）对于每一个，它们的整数平移系为：（2.19）注意这里的上标只是索引而不是指数。它们构成了的规范正交基。因此以上的三个正交基中都至少包含一个带通的或，所以它们都是带通的。也就是说这三部分反应的都是细节信息。具体来说，函数系，是的正交归一基，其中均为整数，=1,2,3分别对应于水平、垂直和对角三个方向。对于任一二维图像信号，在分辨率下有： （2.20）上式表明，在分辨率上将图像分解成、和四个子图，其中代表原图像在分辨率上的近似（即图像的低频部分，不妨用ll来表示），则代表这种近似的误差（即图像的高频部分或“细节”部分）；对应于垂直方向的高频成分，即水平的边缘（细节）信息（不妨用lh表示）；对应于水平方向的高频成分，即垂直的边缘（细节）信息（不妨用hl表示）；则对应于对角方向的高频成分（不妨用hh表示）。图2.1形象地表示了二维图像的多分辨率小波分解。图中符号的上标表示图像的小波分解层数，图中示意了图像的2级小波分解。可以看到，在每一分解层上，图像均被分解为ll, lh, hl和hh的四个频带；下一层的分解仅对低频分量ll进行分解。f (x , y)ll11ll2lh2hl2hh2lh1hl1hh1(a) 图像的小波树形分解ll2hl2 lh2hh2hl1lh1hh1(b) 图像小波分解的塔型结构图2.1二维图像的小波分解按照mallat的快速算法，图像的小波分解算法如图2.2所示：ghghghx与滤波器x卷积h低通滤波器g高通滤波器从两列中取一列从两行中取一行图2.2图像的小波分解算法图2.2示意了图像的一步小波分解过程，可以看到：二维图像的小波分解可以对图像依次按行、按列与一维的低通(h)和高通(g)滤波器作卷积来实现，在卷积之后进行相应的降2采样。图像小波分解的重构算法如图2.3所示。hgxghhg与滤波器x卷积h低通滤波器g高通滤波器在相邻两列间插入一列零在相邻两行间插入一行零图2.3图像的小波重构算法图2.3示意了小波分解图像的一步重构过程，小波分解图像的重构是先对列或行进行升2采样(在相邻列或行间插入一零列或零行)，然后再按行、按列与一维的低通或高通滤波器进行卷积，这样递推下去便可重构原图像。二维图像的这种行、列可分离性简化了图像的小波变换。3小波阈值去噪3.1阈值函数的选取在阈值去噪中，阈值函数体现了对超过和低于阈值的小波系数模的不同处理策略以及不同估计方法。常用的阈值函数有硬阈值函数和软阈值函数两种8-9。硬阈值策略保留大于闽值的小波系数，而把小于阈值的小波系数都设定为零。软阈值策略把小于闽值的小波系数置零，把大于阈值的小波系数的绝对值减去阈值以除去噪声的影响。设a为元素小波系数，t为阈值，两种阈值函数的表达式分别如下：（1） 硬阈值函数：（3.1）（2） 软阈值函数：（3.2）-ttat-ta（a）硬阈值（b）软阈值硬阈值法得到的小波系数的连续性较差，重构信号可能出现突变或振荡现象；如图3.2（a）所示。软阈值法的到的小波系数的连续性好，但当小波系数较大时，得到的处理后的小波系数和实际的小波系数有一定的偏差，会导致重构结果的误差。如图3.2（b）所示。图3.2 两种阈值方法3.2阈值的选取小波阈值去噪方法除了阈值函数的选取，另一个关键因素是对阈值的具体估计。如果阈值较小，去噪后的图像信号与输入比较接近，但是残留了较多噪声。若阈值较大，则得到较多为零的小波系数，对于软阈值策略重建图像变得模糊。在小波域阈值去噪中，阈值的选取直接影响滤波效果。本文用了2种的阈值估计方法。（1） visushrink阈值也称通用阈值，是donoho和johnstone提出的，针对多维独立正态变量联合分布，在维数趋于无穷是得出的结论，是基于最小最大估计得出的最优阈值10。阈值t的选择满足（3.3）其中，n为信号长度，为噪声标准差。的估计公式如下：（3.4）（2） surehrink阈值也称stein无偏风险阈值，是一种基于stein的无偏似然估计的自适应阈值选择s53，它是针对软阈值函数得出的结论。sureshrink阈值的具体计算过程如下：求取信号的长度n；将某一层的小波系数的绝对值由小到大排列，得打一个新的向量，x=x1,x2,x3,xn，其中x1x2x320 ?nj=3, coif2小波3层分解噪声图j=4,coif2小波4层分解噪声图提取1,2,3层小波系数提取1,23,4层小波系数对横、竖、斜三个尺度进行软阈值去噪小波重构输出去噪图像本文采用了小波基为coif2的小波分解法，以及以上2中阈值去噪法，过程如下：图3.4本文小波去噪流程图4实验仿真及结论本文为说明小波去噪的有效性和优越性，对含有高斯白噪声的lenna图像进行消噪处理，其中噪声平均为0，方差为0.01，0.02，0.03，0.04，0.05，0.06。为了客观评价噪声图像的噪声污染度，以及各种去噪方法的效果，本文引用了峰值信噪比psnr：（4.1）式中i(i,j)和i(i,j)分别是有噪图像和原始理想图像在点（i,j）处的灰度值，m、n分别表示图像的行数和列数。由于在各方差噪声下，噪声图像的峰值信噪比均小于20，所以小波尺度取j=4。本文对方差为0.01噪声图像进行小波系数分解，得到下图：图4.1小波系数分解图本文分别用小波去噪、均值去噪、中值去噪和维纳去噪对8位的255灰度值的lena图像进行去噪处理，结果图如下所示：图4.2去噪效果图此外，本文也对多个方差的噪声进行去噪处理，并得出psnr值进行对比，如下表所示。表4.1 各种去噪方法的psnr值比较实验序号噪声方差噪声图像psnr均值滤波psnr中值滤波psnr维纳滤波psnr小波去噪psnr10.0119.534726.303826.918228.369528.535620.0216.924225.617825.482426.418927.080030.0315.432324.992824.523725.155326.131340.0414.348324.453023.624624.231425.411650.0513.597723.989023.041723.493824.829660.0612.971823.564022.437322.934824.3197通过以上结果分析可知，小波阈值去噪的效果较均值滤波、中值滤波和维纳滤波效果都好。参考文献1 杨福生.小波变换的工程分析与应用.北京:科学出版社,1999：150-176.2 阮秋琦.数字图像处理学m.第二版.北京.电子工业出版社.2007.1-23 d.l.donoho, i.m.johnstone. wavelet shrinkage: asymptopia j. r. stat. soc. b.1995, 57：301-369.4 d.l.donoho, i.m.johnstone. ideal time-frequency denoising. standford university,technical report, dept. of statistics,1994：3397-3415.5 张晔，黄秀明.小波变换及在图像处理中的小波特性分析j.中国图像图形学报，1997,2(7):480-484.6 elwood t.olsen, biquan lin. a wavelet phase filter for emission tomography. spie,1996，2491：829-839.7 kozaitis, basuhail. adaptive wavelet threshold selection using higher order statistics for signal denoising. spie,1998，3391：68-74.8 tu dan and shen jianjun.the design of wavelet domain wiener filter and its application in inage denoising.systems engineering and electronics.2011.23.4-79 shubhankar r,bani k.a bayesian transformation model for wavelet shrinkage. ieee trans. on ip.2003.12(12).1512-152010 donoho d l,johnstone i.adapting to unknown smoothness via wavelet shinkage. journal of amercan stat.assoc.1995.90.1201-122511 张磊，潘泉，张洪才，戴冠中，小波域滤波阈值参数的选取，电子学报， 2001，3，29(3).12 xu chen, zhao ruizhen, gan xiaobing. applicationalgorithm of wavelet analysism. beijing：science press，2004.13 chao rui, zhang ke, li yanjun. a wavelet transformbased image fusion methodj.acta electronica sinica，2004，32(5)：750-753.14 gai li-ping,wang gui-lian. technology of filtering in medical image processingj.chinesemedical equipmentjournal,2007,(06):50-51.15 gao qing-wei,li bin. an image de-noising method based on stationary wavelet transformj.journal of computer research and development,2002,(12):1689-1693.tsp题目程序说明1、主函数：run_tsp_ga.m（运行20次）2、次函数：tsp_ga.m备注：（1）、tsp_ga.m可以直接运行，其结果是一次遗传算法计算产生的解。（9198的注释符“%”去掉可以，显示一个最优解并画图）（2）、种群大小在该函数中修改。3、适应度计算函数：calculate_fitness.m4、选择函数：selection.m5、交叉函数：gene_cross.m6、变异函数：bianyi.m7、路程计算函数：lenghtcal.m关于数据存储的说明：数据没有存到外存，直接存在matlab的workspace里。20次计算产生20个历史最优值用nowminlenghtcun存储，其维数是20*100020次计算产生20个每代最优值用lenghtcun存储，其维数是20*100020次计算产生20个每代最优解路径用shortestpath存储，其维数是20*15程序如下主函数：run_tsp_ga.mticclcclear allshortestlen=0; %用来存储每次的最优值shortestpath=1 2 3 4 5 6 7 8 9 10 11 12 13 14 15; %用来存储每次的最优解lenghtcun=zeros(1,1000); %用来存储每次的1000代的每代最好值nowminlenghtcun=zeros(1,1000); for i=1:20i opt_path,nowminlenght,lenght=tsp_ga();%ga运行20次，计算性能if i=1shortestlen=min(nowminlenght);shortestpath=opt_path;lenghtcun=lenght;nowminlenghtcun=nowminlenght;elseshortestlen=shortestlen min(nowminlenght);shortestpath=shortestpath;opt_path;lenghtcun=lenghtcun;lenght;nowminlenghtcun=nowminlenghtcun;nowminlenght;endendaver_shortestlen=mean(shortestlen); %平均性能maxlen=max(shortestlen); %最差性能minlen x=min(shortestlen); %最好性能fangcha=var(shortestlen); %方差fprintf(20次随机，平均最短距离为%3.3fn,aver_shortestlen);fprintf(20次随机，最好性能为%dn,minlen);fprintf(20次随机，最差性能为%dn,maxlen);fprintf(20次随机，方差为%.3fn,fangcha);tocplot(1:1000,lenghtcun(x,:),b)hold onplot(1:1000,nowminlenghtcun(x,:),r)legend(每代最优值,历史最优值)xlabel(代数)ylabel(路程)title(无灾变情况下种群大小=40的最优解)子函数：tsp_ga.m（该子函数是单次ga算法，可直接当主函数运行）% *程序说明*% 作者：陈福沨% 日期：2014年1月9日% 父代种群：individium% 子代：nextgenerat% 下一代种群：nextgenera（竞争选择之后）% 当前最优值：nowminlenght% 每代最优值： lenght% 可更改种群大小 n% 该子程序为一次计算，可直接运行，将15,9198的注释符“%”去掉即显示一个最有解并画图%*function opt_path,nowminlenght,lenght=tsp_ga()clcclear all%tic% 城市间距离邻接矩阵d=0 66 72 16 74 25 60 24 78 44 63 17 67 70 -999; 66 0 7 49 8 67 10 42 13 24 6 49 6 5 22; 72 7 0 56 2 74 14 49 6 30 10 55 7 4 28; 16 49 56 0 58 27 44 9 62 28 47 2 51 54 38; 74 8 2 58 0 75 16 50 4 32 11 57 8 5 30; 25 67 74 27 75 0 64 36 80 52 67 29 71 71 49; 60 10 14 44 16 64 0 36 19 16 4 43 8 14 26; 24 42 49 9 50 36 36 0 55 20 39 8 43 47 35; 78 13 6 62 4 80 19 55 0 36 15 61 12 9 33; 44 24 30 28 32 52 16 20 36 0 20 27 24 29 27; 63 6 10 47 11 67 4 39 15 20 0 46 4 10 25; 17 49 55 2 57 29 43 8 61 27 46 0 50 53 39; 67 6 7 51 8 71 8 43 12 24 4 50 0 8 28; 70 5 4 54 5 71 14 47 9 29 10 53 8 0 24;-999 22 28 38 30 49 26 35 33 27 25 39 28 24 0;% 重要参数定义t=1000; %迭代代数1000代n=40; %种群个体数score=zeros(1,n); %用来存子代分数lenght=zeros(1,n); %用来存子代路程fitness=zeros(1,n); %用来存子代适应度lenght=0; %来存储每代种群中最好的个体nowminlenght=0; %存储当前最好值opt_path=1 2 3 4 5 6 7 8 9 10 11 12 13 14 15; %存储最优解zaibian=0; %灾变标记% 个体的表示和初始化种群individium=struct(path,fitness,lenght,pc,); %把个体用结构体表示：路径，适应度，路程，生存概率for n=1:n %对115随机排序，产生初始种群individium(n).path=randperm(15);end% 主函数nextgenerat=gene_cross(individium,n); %交叉：产生下一代nextgenera为下一代个体individiumt,lenght,score,fitness=calculate_fitness(nextgenerat,n,d); %计算子代的适应度和生存概率nextgenera=selection(individiumt,n,fitness); %选择：用生存概率选择产生下一代种群individiumt,lenght,score,fitness=calculate_fitness(nextgenera,n,d); %计算种群的适应度score=score;lenght=lenght;fitness=fitness;individium=nextgenera; %下一代种群lenght=min(lenght); nowminlenght=lenght;for t=2:t % 遗传算法：交叉变异选择 t %显示代数 nextgenerat=gene_cross(individium,n); %交叉：产生下一代nextgenera为下一代个体 individiumt,lenght,score,fitness=calculate_fitness(nextgenerat,n,d); %计算子代的适应度和生存概率 nextgenerat=bianyi(individiumt,n,1); %变异：变异概率分别为0.06和0.003，p=1是一个系数，用来方便调整变异率 individiumt,lenght,score,fitness=calculate_fitness(nextgenerat,n,d); %计算种群的适应度和生存概率 nextgenera=selection(individiumt,n,fitness); %选择：用生存概率选择产生下一代种群 individiumt,lenght,score,fitness=calculate_fitness(nextgenerat,n,d); score=score;score;lenght=lenght;lenght; fitness=fitness;fitness; individium=nextgenera; %*以下存储最优值和最优解*nowminlenght=nowminlenght min(lenght); %存储当前最优值 minlenght x=min(lenght); if minlenght=nowminlenght(t) zaibian=zaibian+1;elsezaibian=0;endifzaibian=200 for n=1:n %灾变：重新随机初始化individium(n).path=randperm(15);endzaibian=0;endend%opt_path%plot(1:t,lenght,b)%hold on%plot(1:t,nowminlenght,r)%legend(每代最优值,历史最优值)%xlabel(代数)%ylabel(路程)%title(种群大小=10的最优解)子函数：calculate_fitness.m（适应度函数）function individium,lenght,score,fitness=calculate_fitness(individium,n,d) %计算种群的每个个体适应度和生存概率score=zeros(1,n); %处理前的分数scoret=zeros(1,n); %处理后的分数lenght=zeros(1,n); %存储路程fitness=zeros(1,n); %存储适应度for n=1:n score(n)=lenghtcal(individium(n).path,d);individium(n).lenght=score(n);lenght(n)=score(n);end% 对路程值进行取反，并归一化为scorescore_max=max(score);score_min=min(score);ifscore_max=score_min score=(score_max-score)/(score_max-score_min); %由于该实例的存在负值的路程，且值越小适应度越大，归一化处理elseifscore_max0score=score./999;elsescore=-score./999;endendscore=sum(score);average=score/n;% 对score进行处理成scoret,舍弃得分小于平均值的个体for n=1:nif score(n)0 &seletlunpan(i) &selet=lunpan(i+1) nextgenera(n)=individiumt(i+1);endendendendnextgenera(1)=individiumt(x); %将精英保存下去子函数：gene_cross.m（交叉）function nextgenera=gene_cross(individium,n) %交叉算子nextgenera=struct(path,fitness,lenght,pc,);for n=2:2:n site=