project_2.m

来自「本文件夹共三个matlab文件和一个word文档」· M 代码 · 共 36 行

%项目二 连续系统的频域分析
%线性系统为;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)|