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