%項目二 連續系統的頻域分析
%線性系統為;H(jw)=1/(jw+1);輸入信號為周期矩形波
%畫雙邊幅度頻譜
tau_T = 1/4;   %占空比
n = -20:20;
t = -3:0.01:3;
A = 2;
T = 2;
f_t = A*0.5*(square(pi*t,25)+1);  %占空比為25%的方波
F_n = tau_T*sinc(tau_T*n*pi).*exp(-j*tau_T*n.*pi);
Fn_max = max(abs(F_n));
Fn_min = min(abs(F_n));

%畫輸入信號f（t）
subplot(2,1,1),plot(t,f_t);
axis([-3 3 -1 3]);xlabel('t');ylabel('f(t)');
title('周期矩形波');

%畫輸入信號f（t）的幅度頻譜|F(jw)|
subplot(2,1,2),plot(n,abs(F_n)); 
axis([-16 16 Fn_min-1 Fn_max+1]);
xlabel('w');ylabel('|F(jw)|');
title('周期矩形波雙邊帶頻譜|F(jw)|');
%畫系統函數|H（jw）|
w = -20:0.01:20;
H_w = 1./(j*w+1); %注意此處要用點除 為數組的除法
figure(2);
plot(w,abs(H_w),'color',g);
xlabel('w');ylabel('|H(jw)|');
title('系統函數幅度頻譜|H(jw)|');

%畫輸出信號y（t）和其幅度譜|Y(jw)|