% dm09301
% 周期方波信號頻譜分析
function  CTFS_SQ
%  繪制并觀察周期方波信號頻譜   
%  Nf：傅里葉級數展開的項數
%  an：各次諧波余弦項系數
display('Please input the value of Nf ');
Nf = input('Nf = ');
an = zeros(Nf+1,1);
an(1) = 0;
for i = 1:Nf 
    an(i+1) = 2/(i*pi)*sin(i*pi/2);                     %計算系數an
    cn(i+1) = abs(2/(i*pi)*sin(i*pi/2));                %計算幅度譜
end
t = -5:0.01:5;
x = square(pi*(t+0.5));                                 %用square函數求方波信號
subplot(211);
plot(t,x);                                              %繪制方波信號
axis([-5 5 -1.5 1.5]);
title('周期方波信號','Fontsize',8);
xlabel('t   (單位:s)', 'Fontsize',8);
subplot(212);
k = 0:Nf;
stem(k,cn);
hold on;
plot(k,cn);
title('幅度頻譜','Fontsize',8);
xlabel('諧波次數', 'Fontsize',8);
% End