fg_02_02.m

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

% fg_02_02        Continuous and discrete least squares.
N=1000;
k=[0:N-1];
T=35;
ks=[0:T:N-T];
f=5*[0:124,125:-1:-124,-125:124,125:-1:-124,-125:-1]/125;
fs=f(ks+1);
fhat=40*sin(2*pi*k/500)/(pi^2);
%
sp_fig(1);
set(gcf,'position',[80,40,700,450]);
axis([0 1000 -5 5]);
set(gca,'fontsize',18);
set(gca,'position',[.13 .15 .9-.13 .9-.15]);
%
line(ks,fs,'color','k','linestyle','none','marker','o');
line(k,fhat,'color','k'); grid;
xlabel('Time (t)');
ylabel('f_n=f(nT) and l.s. approx.');