freqcalcu.m

量子遗传算法MATLAB源代码

% 计算个体频域
% FreqCalcu
function [HInd,A]=FreqCalcu(n,k,m,w,HTar,individualR)
% n量子位长度
% k系统单元数
% m频域抽样点数
% w数字频率范围
% Compute the frequency samples of the DTFT
for i=1:k
    sos(i,:) = [1 individualR(1,i) individualR(2,i) 1 individualR(3,i) individualR(4,i)];
end
[num,den] = sos2tf(sos);
G = freqz(num,den,w);
%计算A
A=sum(abs(G).*abs(HTar))/sum(abs(G).*abs(G));
%个体频域
HInd=A*abs(G);