二维光子晶体禁带的遗传优化算法MATLAB源码

function [Zp,Yp,Xp,LC1,LC2]=PBGGA(M,N,Pm)
%--------------------------------------------------------------------------
% PBGGA.m
% Photonic Band Gap Genetic Algorithm
% 二维光子晶体禁带的遗传优化算法
% GreenSim团队原创作品，转载请注明
% 欢迎访问GreenSim――算法仿真团队→http://blog.sina.com.cn/greensim
%--------------------------------------------------------------------------
% 函数功能
% 在平面波展开法的基础上，使用遗传算法寻找具有最大光子禁带的二维光子晶体结构
%--------------------------------------------------------------------------
% 主要参考文献
% 龚春娟,胡雄伟.快速遗传算法优化计算二维光子晶体[J].
% 半导体学报.2006,27(6):1098-1102
%--------------------------------------------------------------------------
% 输入参数列表
% M 遗传进化迭代次数
% N 种群规模(取偶数)
% Pm 变异概率
%--------------------------------------------------------------------------
% Zp 最大禁带
% Yp 最大禁带的区间上下限
% Xp 光子晶体结构决策变量（10×10的01矩阵）
% LC1 收敛曲线1，各代最优个体适应值的记录
% LC2 收敛曲线2，各代群体平均适应值的记录
%--------------------------------------------------------------------------
%第一步：变量初始化
epsa=11.56;%介质介电常数
epsb=1;%空气介电常数
Xp=zeros(10,10);%最优结构