fg_08_14.m

英文书《Digital Signal Processing with Examples in MATLAB》附带的MATLAB实例

% fg_08_14    Signals f and g for system identification demo.
K=64;
randn('seed',0);
f=randn(1,K);
c=exp(-.1);
d=filter([0 c*sin(.1*pi)], [1 -2*c*cos(.1*pi) c^2], f);

sp_fig(1,9,5.5);
subplot(2,1,1);
set(gca,'fontsize',16);
line(0:K-1,f,'color','r','marker','o'); grid on;
axis([0 K -3 3]);
xlabel('k'); ylabel('Input f_k');

subplot(2,1,2);
set(gca,'fontsize',16);
line(0:K-1,d,'color','b','marker','o'); grid on;
axis([0 K -3 3]);
xlabel('k'); ylabel('Plant output d_k');