funchirp1.m

d大学gatbs工具箱所有算例下载及调试说明；数学建模讲义；数学建模先进算法讲义

function [q,qq]=funchirp1(fmin,fmax,smax,fs,y,t)
% fmin:搜寻最小偏置频率；fmax:搜寻最大偏置；smax:搜寻最大频偏；fs：采样频率；y：信号；t:时间；perc:保留
% 与最大值相差比例以内数据；No：确定鞍点搜寻点数；
% q 维度定义：1分量序号，2幅值，3频偏，4斜率，5能量，6幅角；
if fmax>fs/2
    fmax=fs/2;
end
ly=length(y);
NIND=40;               %个体数目(Number of individuals)
MAXGEN=50;             %最大遗传代数(Maximum number of generations)
NVAR=2;                %变量数目
PRECI=25;              %变量的二进制位数(Precision of variables)
GGAP=0.9;   
%建立区域描述器(Build field descriptor)
FieldD=[rep([PRECI],[1,NVAR]);rep([fmin -smax;fmax smax],[1,NVAR]);rep([1;0;1;1],[1,NVAR])];
gen=0;                                                 
trace=zeros(MAXGEN, 2);                                %遗传算法性能跟踪初始值
x=bs2rv(Chrom, FieldD);                                %初始种群十进制转换
ObjV=Shubert(x(:,1),x(:,2));                           %计算初始种群的目标函数值
while gen<MAXGEN
    FitnV=ranking(ObjV);                               %分配适应度值(Assign fitness values)
    SelCh=select('sus',Chrom,FitnV,GGAP);              %选择
    SelCh=recombin('xovsp',SelCh,0.7);                 %重组
    SelCh=mut(SelCh);                                  %变异
    x=bs2rv(SelCh,FieldD);                             %子代十进制转换
    ObjVSel=Shubert1(x(:,1),x(:,2));
    [Chrom ObjV]=reins(Chrom,SelCh,1,1,ObjV,ObjVSel);  %重插入
    gen=gen+1;
    [Y, I]=min(ObjV);
    Y,bs2rv(Chrom(I,:),FieldD)                         %输出每一次的最优解及其对应的自变量值
    trace(gen,1)=min(ObjV);                            %遗传算法性能跟踪
    trace(gen,2)=sum(ObjV)/length(ObjV);
    if(gen==50)                                        %迭代数为50时画出目标函数值分布图
        figure(2);
        plot(ObjV);hold on;
        plot(ObjV,'b*');grid;
    end
end
f=Y;
a1=x(1,1);
b1=x(1,2);
q(1,1)=1;
q(1,2)=abs(f);
q(1,3)=a1;
q(1,4)=b1;
h=y-abs(f)*cos(2*pi*((t-min(t))*a1+b1*(t-min(t)).^2)-angle(f));
q(1,5)=(y*y'-b*b');
q(1,6)=-angle(f);
end