s3fig4a.m

来自「Industrial Mathematics」· M 代码 · 共 18 行

% Routine 3.1      Bode plots
clear all;
fS= 1000;
T = 1/fS;
f = 0:fS/2;
o = 2*pi*f;
w = exp(-i*o*T);
N = (1 + w)/2;
D = 1
H = N./D;
r = abs(H);
phi = angle(H);
plot(f,r,'k-')
%plot(f,20*log10(r));
% plot(f,360*phi/(2*pi))