sorce codecleartau=1;t=-2:0.001:2;x(2001)=0;y=exp(-abs(t)/tau);b=-1:0.001:1;n=1;for a=-1:0.001:1for f=-10:0.1:10x(n)=x(n)+exp(-abs(f)/tau)*exp(-j*2*pi*f*a)*0.1;endx(n)=abs(x(n));n=n+1;endsubplot(2,1,1),plot(t,y),xlabel('t'),ylabel('x(t)'),title('tau=1')subplot(2,1,2),plot(b,x),xlabel('f'),ylabel('x(f)')tau 값을 변화 하면서 시뮬레이션 해보았다. 위에서 보는듯이tau 값이 커질수록 time domain에서는 그래프가 퍼지게 되는데 frequency도메인에서는 그래프가 점점 좁아지는 것을 확인할수 있었다.cleartau=1;t=-2:1:2;y=exp(-abs(t)/tau);p=-3:0.001:3;n=1;x(6001)=0;for a=-3:0.001:3for b=-6:1:6x(n)=x(n)+exp(-abs(b)/tau)*exp(-j*2*pi*a*b);endx(n)=abs(x(n));n=n+1;endsubplot(2,1,1);stem(t,y,'^'),ylabel('Xs(t)'),xlabel('t'),title('fs=1')subplot(2,1,2);plot(p,x),ylabel('Xs(f)');xlabel('f')fs가 커질수록 aliasing이 줄어드는 것을 확인할수 있다cleartau=1;t=-3:1:3;x=exp(-abs(t)/tau);p(7)=0;n=0;for k=1:1:7for n=1:1:7p(k)=p(k)+x(n)*exp(-j*2*pi/7*k*(n-1));endp(k)=abs(p(k));n=n+1;endk=1:1:7;subplot(2,1,1);stem(t,x,'^');xlabel('t');ylabel('Xsp(t)'),title('7point dft')subplot(2,1,2);stem(k,p);xlabel('k');ylabel('x(k)');point수가 늘어날수록 원 영상과 비슷한 모양의 그래프가 구현된다.