test.m

来自「量子遗传算法MATLAB源代码」· M 代码 · 共 19 行

clf(figure(1));clf(figure(2));
clear
for i=1:10
   [EMinB(i,:),HInd(i,:),HTar,w]=main();
end
[E,place]=min(EMinB(:,200));
clf(figure(1));clf(figure(2));
figure(1)
plot(w/pi,abs(HInd(place,:)),'k');grid on;hold on 
axis([0,1,0,1.2]);
plot(w/pi,abs(HTar));grid on
figure(2)
plot(w/pi,20*log10(abs(HInd(place,:))),'r');grid on;hold on 
axis([0,1,-80,5]);
%画等阶数MATLAB的巴特沃兹数字低通滤波器
figure(1)
[num,den]=butter(2*3,wp/pi);
Hb=freqz(num,den,w);
plot(w/pi,abs(Hb),'g');