www.mikrocontroller.net

Forum: Digitale Signalverarbeitung / DSP FFT eines codierten Signals


Autor: Dennis Klein (kaleef)
Datum:

Bewertung
0 lesenswert
nicht lesenswert
Hallo zusammen!

Ich habe folgendes Problem und ich hoffe jemand von euch könnte mir 
dabei helfen:

Gegeben ist ein binärer Code (Abfolge von +1 und -1)der mit einem 
Trägersignal gefaltet wird. Dieses Signal wiederum, wird mit seinem 
Spiegelbild gefaltet um eine Delta-Ähnliche funktion raus zu bekommen.
Schaut euch mal mein Matlab-code an, dann wird es euch hoffentlich 
deutlicher:


code = [1 1 1 1 1 -1 -1 1 1 -1 1 -1 1];

Fs = 50000;   %Abtastfreqeunz
T = 1/Fs;
L = 1000;      %Länge
t = (0:L-1)*T;
S = sin(2*pi*t*50); %Trägersignal


%Faltung des Codes mit dem Tägersignal
K = [];
for i = 1:size(code,2)
   K = [K, conv(code(i),S)];
end

D = conv(K,fliplr(K)); %erneute Faltung


figure;subplot(4,1,1);stem(code);
       subplot(4,1,2);plot(S);
       subplot(4,1,3);plot(K)
       subplot(4,1,4);plot(D)
%----------------------------------------------------------------------

Mein problem: ich versuch das FFT-Spektrum meines Ergebnisses (also D) 
darstellen zu lassen, es funktioniert aber irgenwie nicht. Vielleicht 
könnte es an der Abtastfrequenz liegen?! Hat vielleicht jemand eine 
Idee?


Das FFT hab ich wie folgt berechnet:

Y= fft(D,size(D,2));
Py=Y.*conj(Y)/size(D,2);
f=Fs/size(D,2)*(1:size(D,2)/2);
figure;
plot(f,Py(1:size(D,2)/2));
%-------------------------------


habe es aber auch so probiert:


NFFT = 2^nextpow2(size(D,2));
Y = fft(D,NFFT)/size(D,2);
f = Fs/2*linspace(0,1,NFFT/2);

figure;plot(f,2*abs(Y(1:NFFT/2))


Ich hoffe ich konnte mich klar ausdrücken...danke im voraus

Antwort schreiben

Die Angabe einer E-Mail-Adresse ist freiwillig. Wenn Sie automatisch per E-Mail über Antworten auf Ihren Beitrag informiert werden möchten, melden Sie sich bitte an.

Wichtige Regeln - erst lesen, dann posten!

  • Groß- und Kleinschreibung verwenden
  • Längeren Sourcecode nicht im Text einfügen, sondern als Dateianhang

Formatierung (mehr Informationen...)

  • [c]C-Code[/c]
  • [avrasm]AVR-Assembler-Code[/avrasm]
  • [code]Code in anderen Sprachen, ASCII-Zeichnungen[/code]
  • [math]Formel in LaTeX-Syntax[/math]
  • [[Titel]] - Link zu Artikel
  • Verweis auf anderen Beitrag einfügen: Rechtsklick auf Beitragstitel,
    "Adresse kopieren", und in den Text einfügen




Bild automatisch verkleinern, falls nötig
Bitte das JPG-Format nur für Fotos und Scans verwenden!
Zeichnungen und Screenshots im PNG- oder
GIF-Format hochladen. Siehe Bildformate.
Hinweis: der ursprüngliche Beitrag ist mehr als 6 Monate alt.
Bitte hier nur auf die ursprüngliche Frage antworten,
für neue Fragen einen neuen Beitrag erstellen.

Mit dem Abschicken bestätigst du, die Nutzungsbedingungen anzuerkennen.