%《MATLAB在電子信息課程中的應用》第六章例6.11程序q611
% 方波的頻譜分析
% 電子工業出版社出版  陳懷琛 吳大正 高西全合著 2001年10月

% 連續信號的付利葉分析
% 演示(1).頻域樣本點數可取得與時域樣本點數不同；
% （2).若要求的頻譜太寬，而時域樣本點數又取得太少，會發生頻率泄漏
% 設已知信號的時域波形f(t),則其付利葉變換可表為

clear,tf=10; 
N = input('取時間分隔的點數N= ');
dt = 10/N;t = [1:N]*dt; 	% 給出時間分割
f =[ones(1,N/2),zeros(1,N/2)]; 	% 給出信號（此處是方波）
wf = input('需求的頻譜寬度wf= ');
Nf = input('需求的頻譜點數Nf= ');
w1 =linspace(0,wf,Nf);dw=wf/(Nf-1);
F1 = f*exp(-j*t'*w1)*dt;		% 求付利葉變換
w = [-fliplr(w1),w1(2:Nf)]; 		% 補上負頻率
F = [fliplr(F1),F1(2:Nf)]; 		% 補上負頻率區的頻譜
subplot(1,2,1),plot(t,f,'linewidth',1.5),grid
axis([0,10,0,1.1])
subplot(1,2,2),plot(w,abs(F),'linewidth',1.5),grid