clear all clf set(0,'defaultaxesfontsize',20); n1=-100:100; n1=n1+0.000000001; NL=12; wc=pi/2%+(pi/8); % wc=pi; n1=n1+0.5; % h1=sin(wc*(n1))./(pi*(n1)); w1=pi/4; w2=3*pi/4; w3=pi/2; %w3=5*pi/8; %w1=pi/8; w2=3*pi/8; w3=pi/2; % h1=sin(w1*(n1))./(pi*(n1))+ sin(w2*(n1))./(pi*(n1)); %h1=sin(w3*(n1))./(pi*(n1)); mi=0; %h1=8*((sin(w1*(n1+mi*NL))./(pi*(n1+mi*NL))).^2).*cos(w3*(n1+mi*NL)); %h1=3*16*j*(sin(w1*(n1+mi*NL))./(pi*(n1+mi*NL))).*(sin(w2*(n1+mi*NL))./(pi*(n1+mi*NL))).*sin(w3*(n1+mi*NL)); %h1=4*(sin(w1*(n1+mi*NL))./(pi*(n1+mi*NL))).*(sin(w2*(n1+mi*NL))./(pi*(n1+mi*NL))); h0=h1; for mi=1:150, % h1=h1+4*(sin(w1*(n1+mi*NL))./(pi*(n1+mi*NL))).*(sin(w2*(n1+mi*NL))./(pi*(n1+mi*NL))); % h1=h1+4*(sin(w1*(n1-mi*NL))./(pi*(n1-mi*NL))).*(sin(w2*(n1-mi*NL))./(pi*(n1-mi*NL))); % h1=h1+sin(wc*(n1+mi*NL))./(pi*(n1+mi*NL)) + sin(wc*(n1-mi*NL))./(pi*(n1-mi*NL)); % h1=h1+8*((sin(w1*(n1+mi*NL))./(pi*(n1+mi*NL))).^2).*cos(w3*(n1+mi*NL)); % h1=h1+8*((sin(w1*(n1-mi*NL))./(pi*(n1-mi*NL))).^2).*cos(w3*(n1-mi*NL)); % h1=h1+sin(w3*(n1+mi*NL))./(pi*(n1+mi*NL)) + sin(w3*(n1-mi*NL))./(pi*(n1-mi*NL)); h1=h1+sin(w1*(n1+mi*NL))./(pi*(n1+mi*NL)) + sin(w2*(n1+mi*NL))./(pi*(n1+mi*NL)); h1=h1+sin(w1*(n1-mi*NL))./(pi*(n1-mi*NL)) + sin(w2*(n1-mi*NL))./(pi*(n1-mi*NL)); %h1=h1+3*16*j*(sin(w1*(n1+mi*NL))./(pi*(n1+mi*NL))).*(sin(w2*(n1+mi*NL))./(pi*(n1+mi*NL))).*sin(w3*(n1+mi*NL)); %h1=h1+3*16*j*(sin(w1*(n1-mi*NL))./(pi*(n1-mi*NL))).*(sin(w2*(n1-mi*NL))./(pi*(n1-mi*NL))).*sin(w3*(n1-mi*NL)); end X8=fft(h1(101:101+NL-1),NL) % wincr=2*pi/NL; wsd=0:wincr:(2*pi-wincr); %Hs=abs(fft(h0,Nfft)); % Nfft=1032; Nd=Nfft/NL; wfincr=2*pi/Nfft; wdf=0:wfincr:(2*pi-wfincr); Hf=exp(j*wdf*100).*(fft(h0,Nfft)); ws=linspace(0,2*pi,Nfft); plot(ws,Hf,'Linewidth',3); hold on plot(wsd,Hf(1:Nd:end),'rx','MarkerSize',18,'Linewidth',3) %stem(wsd,Hf(1:Nd:end),'r','Linewidth',3) hold off % winc=2*pi/NL; % ideal=exp(j*0.5*winc*[0:NL-1]) % svd([X8.' ideal.']) % % wc1=5*pi/8; wc2=3*pi/8; % h1=(8/3)*(sin(wc1*(n1))./(pi*(n1))).*(sin(wc2*(n1))./(pi*(n1))); % for mi=1:150, % h1=h1+(8/3)*(sin(wc1*(n1+mi*NL))./(pi*(n1+mi*NL))).*(sin(wc2*(n1+mi*NL))./(pi*(n1+mi*NL))); % h1=h1+(8/3)*(sin(wc1*(n1-mi*NL))./(pi*(n1-mi*NL))).*(sin(wc2*(n1-mi*NL))./(pi*(n1-mi*NL))); % end % X8=fft(h1(101:101+NL),NL) % winc=2*pi/NL; % ideal=exp(j*0.5*winc*[0:NL-1]) % svd([X8.' ideal.'])