examp35.m

来自「ADSP TOOLBOX: Version 2.0 and gui m-file」· M 代码 · 共 22 行

clc,echo on 
%EXAMPLE 35 
x='3*exp(-2*t)';              % Expression for signal x(t) 
v=fseries(x, 2);              % Fourier series results to k=32 are in v

t=0:0.02:2;                   % Time array for 1 period 
ak=v(:,2);bk=v(:,3);          % Trig coefficients from matrix v 
sm=ak(1);                     % Start reconstruction with dc value
echo off 
for k=1:30;                   % Loop for reconstruction
sm=sm+ak(k+1)*cos(k*pi*t)+bk(k+1)*sin(k*pi*t); 
if k==7,sm7=sm;end            % Save intermediate result if k=7 
end                           % End of loop. Prompt reappears
echo on
%for k=1:30;                  % Loop for reconstruction
%sm=sm+ak(k+1)*cos(k*pi*t)+bk(k+1)*sin(k*pi*t); 
%if k==7,sm7=sm;end           % Save intermediate result if k=7 
%end                          % End of loop. Prompt reappears
plot(t,eval(x),t,sm7,t,sm)    % Overplot results 
echo off %end of example