%《MATLAB在電子信息課程中的應用》第七章例7.9程序q709
% 離散時間傅立葉變換
% 電子工業出版社出版  陳懷琛 吳大正 高西全合著 2001年10月

disp('八點時間信號的離散時間傅立葉變換')
x0=sin(2*pi* [1:8]/8) *5;			 % x0是8點行向量
dt=2*pi/8;
w=linspace(-2*pi,2*pi,1000)/dt;	 % w是1000點行向量
X0 = dtft(x0,w)*dt;				% 求得頻率響應X0
subplot(3,1,1),plot(w,abs(X0)),grid,shg		% 畫圖
disp('重復N次的八點時間信號的離散時間傅立葉變換')
N=input('N= ');						% 用鍵盤輸入延拓周期數
x1=reshape(x0'*ones(1,N),1,N*length(x0));	% 延拓后的時域信號x1
X1 = dtft(x1,w)*dt;				% 求x1的頻率響應X1
subplot(3,1,2),plot(w,abs(X1)),grid,shg		% 畫圖
%disp('重復100次的八點時間信號的離散時間傅立葉變換')
%x1=reshape(x0'*ones(1,100),1,100*length(x0));	% 延拓后的時域信號x1
%X1 = dtft(x1,w)/100*dt;				% 求x1的頻率響應X1
%subplot(4,1,3),plot(w,abs(X1)),grid,shg		% 畫圖
disp('重復無窮次的八點信號的離散時間傅立葉變換-傅立葉級數')
pause,X2=fft(x0*dt);							% 離散傅立葉變換
w1=2*pi* [0:length(x0)-1]/length(x0);			 % 離散頻點向量
subplot(3,1,3),stem([-w1,w1],[abs(X2),abs(X2)]),grid,shg
axis([min(w),max(w),0,max(abs(X2))]),grid,shg
%gtext('八點時間信號的離散時間傅立葉變換')
%gtext('重復N次的八點時間信號的離散時間傅立葉變換')
%gtext('重復100次的八點時間信號的離散時間傅立葉變換')
%gtext('重復無窮次的八點信號的離散傅立葉變換')