THESIS-FINAL-CODE-2.rar_full_thesis simulation

%cos(2*pi/3*n) omegap=pi/4; Apmax=3; omegar=2/3*pi; Armin=50; T=5*pi*10^(-6); wp=[omegap/pi]; wr=[omegar/pi]; [N,wc]=buttord(wp,wr,Apmax,Armin,'s'); [b,a]=butter(N,wc,'low'); popsize=100; chromlength=10; pc=0.5; pm=0,01; pop=initpop(popsize,chromlength); for i=1:100 [objvalue]=calobjvalue(pop); fitvalue=calfitvalue(objvalue); [newpop]=selection(pop,fitvalue); [newpop]=crs(pop,pc); [newpop]=mutation(pop,pm); [bestindividual,bestfit]=best(pop,fitvalue); y(i)=max(bestfit); n(i)=i; pop5=bestindividual; x(i)=decodechrom(pop5,1,chromlength)*8/1023-4; pop=newpop; end figure lowpas sysl=tf(b,a) w=linspace(0,2*pi,500); h=freqz(b,a,w); length(h) mag=fftshift(20*log10(abs(h))); n=0:23; x=1+cos(pi/4*n)+cos(2*pi/3*n); X=fft(x,500); X1=fftshift(X); omega=2*pi/500*([0:499]); Y=fftshift(X.*h); figure resp figure firsttaska figure title(''); ylabel('');xlabel(''); plot(omega,abs(X1));axis([0,2*pi,0,35]); figure title(''); ylabel('Gain(dB)');xlabel('Frequency(rad)'); plot(w,mag);axis([0,2*pi,-100,5]); figure plot(omega,abs(Y));axis([0,2*pi,0,25]);