fg_02_03.m

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

% fg_02_03       Commands executed in (2.17).
f=[0.00 1.40 2.80 4.20 4.40 3.00 1.60 0.20 -1.20 -2.60 -4.00 -4.60];
m=[1:3];
n=[0:11];
T=35;
G=sin(2*pi*(n')*(2*m-1)*T/500);

figure(1); clf;
% Text box (2.17)
c=(G'*G)\(f*G)';
n=[0:12*35-1]';
arg=2*pi*n/500;
fhat=[sin(arg) sin(3*arg) sin(5*arg)]*c;
plot(35*[0:11],f,'ok',n,fhat); grid
line(n,fhat,'color','k');
% Do not include the next statement.
set(gca,'fontname','times','fontsize',18);
%
xlabel('t'); ylabel('signal, f(t)');