Vergleich Signal und Schätzung

T=1;   % % Zeitkonstante T=1

N=128; %  Stützstellenzahl

fp=32; % Größe des spektralen Fensters

f0=fp/N; % --> Erfassen des Spektrums an N Stützstellen

         %     im Raster f0=fp/N

t0=1/fp; 

tp=1/f0;

t=0:t0:tp-t0;

f=-fp/2:f0:fp/2-f0

U=1./(1+j*2*pi*f*T);

D=f0.*ifftshift(U);

F=dftmtx(N);

d=N.*inv(F)*conj(D');

stem(t,d,'r'); % Darstellung des geschätzten Zeitsignals

xlabel('t in sec'); % Achsenbeschriftung

ylabel ('d(n*t_0); u_p(t)'); % Achsenbeschriftung

t=linspace (0,4.*T,500); % lieanar gestufter Argumentenvektor

U0=1;

ut=U0*exp(-t/T);

ut(1)=U0/2;

hold on;

plot(t,ut,'b-');