%《MATLAB及其在理工課程中的應用指南》第九章例9－1－7程序
% 方波經過低通濾波器后恢復
% 西安電子科技大學出版社出版  陳懷琛編著 1999年10月初版，2004年11月第二版
%
%連續信號的付利葉分析
%設已知信號的時域波形f(t),則其付利葉變換可表為
clear,tf=10; N=256;
t = linspace(0,tf,N); 	% 給出時間分割
w1 = linspace(eps,20,N);dw =20/(N-1);
f =[ones(1,N/2),zeros(1,N/2)]; 	    % 給出信號（此處是方波）
F1 = f*exp(-j*t'*w1)*tf/(N-1);		% 求付利葉變換
w = [-fliplr(w1),w1(2:N)]; 		    % 補上負頻率
F = [fliplr(F1),F1(2:N)]; 		    % 補上負頻率區的頻譜
w2 = w(N/2:3*N/2);				    % 取出中段頻率
F2 = F(N/2:3*N/2);			 	    % 取出中段頻譜
subplot(1,2,1),plot(w,abs(F),'linewidth',1.5),grid
set(gcf,'color','w')                % 置圖形背景色為白色
f1=F2*exp(j*w2'*t)/pi*dw; 			% 對中段頻譜求付利葉逆變換
subplot(1,2,2),plot(t,f,t,f1,'linewidth',1.5),grid