WICHTIG: Der Betrieb von goMatlab.de wird privat finanziert fortgesetzt. - Mehr Infos...

Mein MATLAB Forum - goMatlab.de

Mein MATLAB Forum

  
Gast > Registrieren       Autologin?   

Partner:




Forum
      Option
[Erweitert]
  • Diese Seite per Mail weiterempfehlen
     


Gehe zu:  
Neues Thema eröffnen Neue Antwort erstellen

FFT-Daten ausgeben
 

01detlef
Forum-Anfänger

Forum-Anfänger

Beiträge: 47
Anmeldedatum: 07.05.10
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 16.08.2010, 21:06     Titel: FFT-Daten ausgeben
  Antworten mit Zitat      
Hallo,

ich berechne eine Frequenzanalyse mit dem folgenden Code
Code:

...
Y = fft(vL,NFFT)/L
f = Fs/2*linspace(0,1,NFFT/2+1);

subplot(3,1,3);
stem(f,2*abs(Y(1:NFFT/2+1)))
 


Nun möchte ich zusätzlich die Frequenzen und die dazugehörigen Amplituden in einer Tabelle ausgeben. Wie ist das möglich? Die Amplituden sind sehr klein, wie kann man diese dennoch anzeigen?

detlef
Private Nachricht senden Benutzer-Profile anzeigen


Fourier
Forum-Guru

Forum-Guru

Beiträge: 329
Anmeldedatum: 05.04.10
Wohnort: Karlsruhe
Version: 2009b
     Beitrag Verfasst am: 16.08.2010, 21:43     Titel:
  Antworten mit Zitat      
Hallo,

du kannst die Amplituden logaritmisch anzeigen lasse, ist auch üblich bei einer FFT-Analyse.
Code:

20*log(abs(Y))  % Betragsspektrum
 


gruss
_________________

Ein schlechter General ist besser als zwei gute.
Private Nachricht senden Benutzer-Profile anzeigen
 
01detlef
Themenstarter

 Forum-Anfänger

Forum-Anfänger

Beiträge: 47
Anmeldedatum: 07.05.10
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 17.08.2010, 08:05     Titel:
  Antworten mit Zitat      
Hallo,

ja das mache ich auch in einem anderen Bild! Aber wie kann man die konkreten Werte in einer Text-Datei ausgeben lassen (Frequenz /Amplitude)

detlef
Private Nachricht senden Benutzer-Profile anzeigen
 
Fourier
Forum-Guru

Forum-Guru

Beiträge: 329
Anmeldedatum: 05.04.10
Wohnort: Karlsruhe
Version: 2009b
     Beitrag Verfasst am: 17.08.2010, 08:36     Titel:
  Antworten mit Zitat      
Hi,

unter Annahme dass Y auch zeilenvektor ist:
Code:

dlmwrite(file,[f' Y'],'Delimiter',' ','Newline','pc')
 


gruss
_________________

Ein schlechter General ist besser als zwei gute.
Private Nachricht senden Benutzer-Profile anzeigen
 
01detlef
Themenstarter

 Forum-Anfänger

Forum-Anfänger

Beiträge: 47
Anmeldedatum: 07.05.10
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 17.08.2010, 09:54     Titel:
  Antworten mit Zitat      
Hallo,

also ich habe das nun so gemacht:
Code:
if (MOD (Y,2100)== 0)
fid2 = fopen('Ausgabewert.txt','wt');
%fprintf(fid2,'Hallo\n\n');

fprintf(fid2,'%f \n',Y);
fclose(fid2);
end;


habe aber zwei Probleme. Die Amplituden sind ganz kleine Werte und es werden nur Nullen ausgegaben, wie kann ich das ändern? Und wie kann ich die Bedingung einbauen, dass nur die Werte herausgeschrieben werden, die durch 2100 teilbar ohne Rest sind?

detlef
Private Nachricht senden Benutzer-Profile anzeigen
 
01detlef
Themenstarter

 Forum-Anfänger

Forum-Anfänger

Beiträge: 47
Anmeldedatum: 07.05.10
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 17.08.2010, 10:33     Titel:
  Antworten mit Zitat      
also
Code:

Y = fft(vL,NFFT)/L;
f = Fs/2*linspace(0,1,NFFT/2+1)

%FFT-Werte in txt-Datei schreiben

fid2 = fopen('Ausgabewert.txt','wt');
fprintf(fid2,'%d     %f \n',f,Y);
fclose(fid2);


Wieso klappt es so nicht, die Daten von f in der ersten Spalte darzustellen und die von Y in der zweiten?

detlef
Private Nachricht senden Benutzer-Profile anzeigen
 
Fourier
Forum-Guru

Forum-Guru

Beiträge: 329
Anmeldedatum: 05.04.10
Wohnort: Karlsruhe
Version: 2009b
     Beitrag Verfasst am: 17.08.2010, 11:04     Titel:
  Antworten mit Zitat      
Code:

fprintf(fid2,'%d     %f \n',[f;Y]);

 


wobei Y auch reihenvektor ist, ansonsten Y'
_________________

Ein schlechter General ist besser als zwei gute.
Private Nachricht senden Benutzer-Profile anzeigen
 
01detlef
Themenstarter

 Forum-Anfänger

Forum-Anfänger

Beiträge: 47
Anmeldedatum: 07.05.10
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 17.08.2010, 11:08     Titel:
  Antworten mit Zitat      
Es kommt dann immer die Fehlermeldung
??? Error using ==> vertcat
CAT arguments dimensions are not consistent.

Error in ==> FFT_res1_Variante2 at 44
fprintf(fid2,'%f %f \n',[f;Y']);

detlef
Private Nachricht senden Benutzer-Profile anzeigen
 
01detlef
Themenstarter

 Forum-Anfänger

Forum-Anfänger

Beiträge: 47
Anmeldedatum: 07.05.10
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 17.08.2010, 11:14     Titel:
  Antworten mit Zitat      
Das Problem ist, dass nicht spaltenweise ausgegeben wird, sondern erst f und dann danach Y' !?
Private Nachricht senden Benutzer-Profile anzeigen
 
01detlef
Themenstarter

 Forum-Anfänger

Forum-Anfänger

Beiträge: 47
Anmeldedatum: 07.05.10
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 17.08.2010, 11:22     Titel:
  Antworten mit Zitat      
Okay,

ich habs jetzt hinbekommen!

Wie kann man nun noch angeben, dass die Ausgabe als Potenzschreibweise angegeben wird, also X * 10^Y ? Meine Zahlen sind alle sehr klein und
0.000000000045 möchte ich lieber als Potenz schreiben!

danke detlef
Private Nachricht senden Benutzer-Profile anzeigen
 
01detlef
Themenstarter

 Forum-Anfänger

Forum-Anfänger

Beiträge: 47
Anmeldedatum: 07.05.10
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 17.08.2010, 18:25     Titel:
  Antworten mit Zitat      
Hallo,

so also nochmal konkret: Es läuft alles wunderbar, aber die Ausgabe nicht mir zu unübersichtlich und deshalb möchte ich, dass die Werte als Potenz gespeichert werden. Also nicht 0.001 sondern 1*10^-3 soll in die txt-Datei geschrieben werden.

Wie macht man das?

detlef
Private Nachricht senden Benutzer-Profile anzeigen
 
denny
Supporter

Supporter


Beiträge: 3.853
Anmeldedatum: 14.02.08
Wohnort: Ulm
Version: R2012b
     Beitrag Verfasst am: 17.08.2010, 18:38     Titel:
  Antworten mit Zitat      
Hallo

so:
Code:

fprintf('%.2e\n',0.001)
 
Private Nachricht senden Benutzer-Profile anzeigen
 
Neues Thema eröffnen Neue Antwort erstellen



Einstellungen und Berechtigungen
Beiträge der letzten Zeit anzeigen:

 Du kannst Beiträge in dieses Forum schreiben.
Du kannst auf Beiträge in diesem Forum antworten.
Du kannst deine Beiträge in diesem Forum nicht bearbeiten.
Du kannst deine Beiträge in diesem Forum nicht löschen.
Du kannst an Umfragen in diesem Forum nicht mitmachen.
Du kannst Dateien in diesem Forum posten
Du kannst Dateien in diesem Forum herunterladen
.





 Impressum  | Nutzungsbedingungen  | Datenschutz | FAQ | goMatlab RSS Button RSS

Hosted by:


Copyright © 2007 - 2025 goMatlab.de | Dies ist keine offizielle Website der Firma The Mathworks

MATLAB, Simulink, Stateflow, Handle Graphics, Real-Time Workshop, SimBiology, SimHydraulics, SimEvents, and xPC TargetBox are registered trademarks and The MathWorks, the L-shaped membrane logo, and Embedded MATLAB are trademarks of The MathWorks, Inc.