FFT Matlab

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

Partner:

Forum

Option

[Erweitert]

• Diese Seite per Mail weiterempfehlen

Gehe zu:

ThKo26

Forum-Century


	Beiträge: 184

	Anmeldedatum: 21.09.17

	Wohnort: ---

	Version: 2015b

Verfasst am: 17.08.2019, 20:00 Titel: FFT Matlab

Hallo zusammen,

es gibt schon mehrere Themen mit FFT, ich habe auch schon viele verschiedene Sachen versucht, komme aber "glaube ich" nicht auf die richtige Lösung.

Zu meinem Code:

Code:

%% Fourier Zerlegung eines periodischen Rechtecksignals
clc
clear
close all
%% Werte eingeben
a = 1; % Impulshöhe
xend = 10; % max. x-Wert
deltax = 0.01; % x-Inkrement sampling period
fend = 9; % max. Oberwelle (Skalar oder Vektor möglich)

%% Rechnungen
dim = 1+(xend / deltax); % Anzahl an Zwischenwerten
x = 0:deltax:xend; % Time Vector
f1 = figure('Name', 'Rechteck'); % Graphikbestimmung
y = zeros(1,dim); % Damit Startwert: Überall y = 0;
igv = 1;
ym = zeros(round(fend/2),length(x));
s1 = subplot(1,3,1);
for i = 1:2:fend
y = y + (sin(i *x)/i);
ym(igv,:) = (sin(i *x)/i);
p11 = plot(x,ym(igv,:));
hold all
igv = igv+1;
end
grid on
xlabel('Time (sec)');
ylabel('y = f(x)');

%% Graphikausgabe
s2 = subplot(1,3,2);
y = ((4 * a)/pi) .* y; % Normierung
p = plot(x,y); % graphische Darstellung der Ergebnisse
hold all;grid on;
xlabel('Time (sec)');
ylabel('y = f(x)');

f2 = figure(2);
set(gcf,'color',[1 1 1],'position',[10 10 770 520])
% hf = stem3(zeros(1,31),fspec(1:31), ampspec(1:31),'fill');
grid on
hold all
igj = 1;
hs = zeros(length(1:2:fend),1);
for igh = 1:2:fend
hs(igj) = plot3(x(1:n), igh.*ones(1,n), ym(igj,1:n),'linewidth',2);
igj = igj+1;
end
hs1 = plot3(x(1:n), zeros(1,n), y(1:n),'-r','linewidth',3);

xlabel('Time [s]','fontsize',18)
ylabel('Frequency [Hz]','fontsize',18)
zlabel('Amplitude [-]','fontsize',18)

l2 = plot3([0,xend],[0,0],[0,0],'Color','black','linewidth',2);
l3 = plot3([0,0],[0,xend],[0,0],'Color','black','linewidth',2);
view([-45 25]);

Funktion ohne Link?

Ich habe eine Fourierreihe abgebildet. nun möchte die fft auf y anwenden

Hierzu folgender Code:

Code:

%% FFT
N = 1024; % gewünschte FFT-Länge (N=2^x, sonst wird der DFT-Algorithmus verwendet!)
yfft = fft(y,N); % FFT
Fs = 1/deltax; % sampling frequency Abtastung

% FREQUENZ
fspec1 = linspace(0,Fs,N);
fspec2 = linspace(0,Fs,N/2+1)
% AMPLITUDE
ampspec = 2.*abs(yfft)/N; % Amplitude
ampspec2 = ampspec(1:N/2+1)

Funktion ohne Link?

Wenn ich nun das Frequenzband zur Amplitue plotte, sollte sich nach meinem Verständnis die Überhöhungen bei den Frequenzen 1,3,5,7,9 sein. Dies ist in der figure f2 zu erkennen. oder ist dies nur der Fall, wenn ich alle getrennt von einander betrachte, als die ym step by step analysiere.

Wo liegt mein Denkfehler. Gerne ausführliche Informationen.

Vielen Dank.

Grüße
Thomas

ThKo26

Themenstarter

Forum-Century


	Beiträge: 184

	Anmeldedatum: 21.09.17

	Wohnort: ---

	Version: 2015b

Verfasst am: 18.08.2019, 19:59 Titel:

Hallo zusammen,

ich habe einige Antworten gefunden.

Die Frequenzen sind nicht bei 1,3,5,..., da in diesem Fall der Sinus mit 2*pi mulitpliziert werden muss.

Grüße
Thomas


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 |

RSS

Hosted by:

Copyright © 2007 - 2024 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.