Vergleich Signal und Schätzung

[mazze32 (Gast)]

Hier mein Lösungsvorschlag

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-');

und die entsprechende figure

http://img207.imageshack.us/img207/7711/grafgg5.jpg

ist alles beachtet worden? und richtig?

danke für die hilfe


mfg

[Waschtl (Gast)]

Ich hab nicht dein ganzne Code durchgeschaut od. die Augabe aber eins 
solltest du wissen: DFT ist ne Ortogonale Transformation. D.h. alle 
basis Funktionen sind senkrecht aufeinander.
==> inv(F) == F'
mit nem konstanten Faktor noch korrigieren.