Fourierkoeffizienten berechnen

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

Partner:

Forum

Option

[Erweitert]

• Diese Seite per Mail weiterempfehlen

Gehe zu:

chrisi99

Forum-Anfänger


	Beiträge: 31

	Anmeldedatum: 09.04.08

	Wohnort: ---

	Version: ---

Verfasst am: 15.05.2008, 21:26 Titel: Fourierkoeffizienten berechnen

Hallo Leute!

ich würde gerne die Fourierkoeffizienten einer (ein)gegebenen Funktion berechnen.Leider funktioniert es noch nicht so richtig. Vielleicht kann mir jemand helfen?

Code:

[a,b]=fouriercoeff(f,n)

f_a = @(x) f(x).*cos(k*x);
f_b=@(x) f(x).*sin(k*x);

a=zeros(1,n);
b=zeros(1,n);

for k=1:n
a(k)=1/pi*quadl(f_a,-pi,pi);
b(k)=1/pi*quadl(f_b,-pi,pi);
end

Funktion ohne Link?

f soll als function handle per @ übergeben werden.

lg
Chris

nschlange

Ehrenmitglied


	Beiträge: 1.318

	Anmeldedatum: 06.09.07

	Wohnort: NRW

	Version: R2007b

Verfasst am: 16.05.2008, 21:31 Titel:

Hi,

ich würde das so machen:

Code:

function fouriermain
clear all;
clc;

f1=@(t) cos(t);
f2=@(t) sin(t);
f3=@(t) sin(t).*cos(t);

n=5;

[a0 a b]=fouriercoeff(f3,n);

disp('a_0=')
disp(a0)
disp('a_n')
disp(a)
disp('b_n')
disp(b)
end

function [a0 a b] = fouriercoeff(f,n,T)
if nargin==2
T=2*pi; % normalerweise 2Pi-periodisch
end

w=2*pi/T;

a0=2/T*quad(f,0,T);
ak=@(k) (2/T*quad(@(t) (f(t).*cos(k*w*t)),0,T));
bk=@(k) (2/T*quad(@(t) (f(t).*sin(k*w*t)),0,T));

a=zeros(1,n);
b=a;
for k=1:n
a(k)=ak(k);
b(k)=bk(k);
end

end

Funktion ohne Link?

Hoffentlich stimmt's auch Wink

Des Weiteren bin ich wie immer an Verbesserungsvorschlägen interessiert.
Und wenn noch jemand eine schöne Ausgabe dazu bastelt haben wir wieder was für die Skript-Ecke Cool

fouriermain.m

Beschreibung:

Download

Dateiname:

fouriermain.m

Dateigröße:

536 Bytes

Heruntergeladen:

970 mal

_________________

Viele Grüße
nschlange

"Chuck Norris ejakuliert fluessigen Stahl!"

chrisi99

Themenstarter

Forum-Anfänger


	Beiträge: 31

	Anmeldedatum: 09.04.08

	Wohnort: ---

	Version: ---

Verfasst am: 17.05.2008, 17:27 Titel:

ok, habe es jetzt so versucht wie du oben, leider noch immer fehlerhaft Sad

vielleicht findet jemand mein Problem....

Code:

function [a,b]=fouriercoeff(f,n)

a=zeros(1,n);
b=zeros(1,n);

a0=1/pi*quadl(f,-pi,pi);

for k=1:n
ak=@(k) (1/pi*quadl(@(x) (f(x) .*cos(k*x)),-pi,pi));
bk=@(k) (1/pi*quadl(@(x) (f(x).*sin(k*x)), -pi,pi));
a(k)=ak(k);
b(k)=bk(k);

end

Funktion ohne Link?

der Code funktioniert zwar, spuckt aber einfach nicht die richtige Koeffizienten aus (und damit meine ich nicht Abweichungen in der 10. Kommastelle, sondern einfach falsche Werte...)

lg

nschlange

Ehrenmitglied


	Beiträge: 1.318

	Anmeldedatum: 06.09.07

	Wohnort: NRW

	Version: R2007b

Verfasst am: 17.05.2008, 18:01 Titel:

Welche Funktion schmeißt Du denn darein?
Ist die auch 2pi-periodisch?

Hiermit klappt es:

Code:

f=@(x) 2+3*sin(2*x)-5*cos(3*x)
[a, b]=fouriercoeff(f,6)

Funktion ohne Link?

Außer daß Du a0 nicht zurückgibst.
_________________

Viele Grüße
nschlange

"Chuck Norris ejakuliert fluessigen Stahl!"

chrisi99

Themenstarter

Forum-Anfänger


	Beiträge: 31

	Anmeldedatum: 09.04.08

	Wohnort: ---

	Version: ---

Verfasst am: 17.05.2008, 18:08 Titel:

ja, sind 2pi-periodisch.

dann füge ich a0 und b0 in die Vektoren dazu....

Code:

a=[a0,a];
b=[0,b];

Funktion ohne Link?

komisch, bei mir sagt er immer, dass es fehlerhaft ist ....

ärgerlich!

nschlange

Ehrenmitglied


	Beiträge: 1.318

	Anmeldedatum: 06.09.07

	Wohnort: NRW

	Version: R2007b

Verfasst am: 17.05.2008, 18:12 Titel:

Was ist fehlerhaft? Fehlermeldung?
Welche Funktion packst Du denn konkret rein?
Welche Koeffizienten kommen konkret raus?
Welche sollten rauskommen?
_________________

Viele Grüße
nschlange

"Chuck Norris ejakuliert fluessigen Stahl!"


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.