function[objvalue]=GAaim（popsizerhoapharhoa0pha0nfreq）
for?j1=1:popsize
????beta=0;
????for?j2=1:nfreq
????????beta=beta+（（rhoa（j2）-rhoa0（j1j2））/rhoa（j2））^2+（（pha（j2）-pha0（j1j2））/pha（j2））^2;
????end
????objvalue（j11）=sqrt（beta）;
end