Springen naar inhoud

FouriercoŽfficiŽnt bepalen


  • Log in om te kunnen reageren

#1

Andy

    Andy


  • >250 berichten
  • 294 berichten
  • Ervaren gebruiker

Geplaatst op 15 mei 2005 - 15:38

Hellow

kheb een dik probleem!v

voor een toepassing (golfvergelijking-achtige toestand) heb ik nood aan een fourierontbinding van

f(x)=sin(1/2*x)^3

echter, tis niet zomaar die functie, de functie gegeven ligt aan de basis...
de functie moet periode 2*Pi hebben en oneven zijn... (alleja, -f(-x) =f(x) , das voorwaarde, nogal problemen met even oneven toestanden, wissel ze soms wel eens om)

dus

ff(x) = herhaling van het stukje sin(1/2*x)^3 tussen -Pi en Pi

tot daar: geen problemen

alleja, via berekeningen (die normaal moeten kloppen) kom ik uit:

f(x) = D+(sum(C*a*sin(lambda*x), n = 1 .. infinity))

dus D moet nul zijn aangezien f(x) oneven is

bijgevolg blijven C*a nog over als fouriercoeffn

ik reken vervolgens uit:

C*a := 2*int(f(x)*sin(lambda*x), x = 0 .. L)/L

met als gevolg

Ca := 1/2*sin(Pi*n)*(4*n^3-4*n+6*n^2-3)/(n*(n^3+2*n^2-n-2)*Pi)

sin(Pi*n) = 0 want n is geheel getal!

dus al mijn fouriercoeff zijn nul

maw mijn functie is nul!
terwijl duidelijk gegeven is dat die verschilt van nul

Waar zit mijn fout???

dank
mvg
Andy

Dit forum kan gratis blijven vanwege banners als deze. Door te registeren zal de onderstaande banner overigens verdwijnen.

#2

Bart

    Bart


  • >5k berichten
  • 7224 berichten
  • VIP

Geplaatst op 15 mei 2005 - 17:12

De functie is 4 ;) periodiek en niet 2 :shock:

Maar of daar het probleem zit. Ik kwam er namelijk ook niet uit.
If I have seen further it is by standing on the shoulders of giants.-- Isaac Newton

#3

Andy

    Andy


  • >250 berichten
  • 294 berichten
  • Ervaren gebruiker

Geplaatst op 15 mei 2005 - 18:14

weljah, daar ben ik mij bewust van, maar ik definieer een nieuwe functie die 2Pi continu is ipv 4Pi door gewoonweg stuk te nemen (nl stuk rond -Pi..Pi) en dat naast elkaar te plakken. ik construeer zo wel discontinuiteiten, maar we hadden gezien dat dit geen probleem geeft, je krijgt dan het Gibbs effect, maar het zou moeten lukken... ik weet niet of je maple hebt, maar de functie die ik wil plotten is maw
f:=(x)-> piecewise(sin(1/2*x)^3,x<Pi,sin(1/2*(x-3*Pi))^3,x<3*Pi) enz...

#4

Andy

    Andy


  • >250 berichten
  • 294 berichten
  • Ervaren gebruiker

Geplaatst op 16 mei 2005 - 15:37

nuja, had blijkbaar wel wat fouten gemaakt... geen antwoord meer nodig:
ge kunt het ontbinden in 2 termen met sinus... kweetniet juist meer welke twee maar tkomt erop neer dat het in orde is
toch bedankt voor reactie!

mvg
Andy





0 gebruiker(s) lezen dit onderwerp

0 leden, 0 bezoekers, 0 anonieme gebruikers

Ook adverteren op onze website? Lees hier meer!

Gesponsorde vacatures

Vacatures