[Playlisten] [Impressum und Datenschutzerklärung]

17.3 FFT in MATLAB(R), Window (Fensterfunktion), Hann


CC-BY-NC-SA 3.0

Tempo:

Anklickbares Transkript:

beiHerrn Foyer fing es damit an das er sich überlegt hat das wenn er eineperiodische Funktion nimmtman Herausforderungenperiodische Funktion zu zeichnenKommadass man dieService erkenntdas soll eine Periode sein das man die zerlegenkann in sieben gasförmige Schwingungenderselben Frequenzoder eines ganzzahligenvielfachendieser Frequenzvielleicht mag hierirgendwasvielleichtmarkiert diese Sinus Schwingungdrin seinmit der Phasenlageals Grundfälledann die Nuss große Schwingung gegen die gemischt passende Phasenlagedoppelte Fremdsprachefrequenzund so weiter und sofort alles zusammengesetztgibt dann zum Schluss netterweise wieder die Originalschwingungdas ließ sichin Sinus und Kosinus etwas billiger hinschreiben als mit den wie hoch dieFunktionendie schreib noch mal hin zur Wiederholungalter der andere ist meinemein Signal sollten sagen ich mein Funktion sei mein Signal ausgemessenesgern von der Zeit abhängtmit der Periode Tnungroß T mit Periode groß T setzt sich an alsReihe eine unendlich lange Summe das ist die Fourierreihewas wird auf summiert diese netten Funktionenjedoch zwo Pidie Sohn zu vielfacherFrequenzähmsagt das Vielfache die Frequenz gegenüber der Grundfälle erhöht ist und es auch hier noch mal die Zeit durch die Periodenlängedamit auch hinhaut??Zeit durch Periodenlängeist in einer Periodenlängeeinmal durch jede von null eins wächst um einsähm malmacht N Durchgänge pro Periodenlängezwei Pi den FallumdrehungenN Malpunkt S I sorgt dafür dass es wirklich Sinus und Kosinus werden mit der Euler Formelso konnte hinten zustandedavor stehen die Foyerkoeffizientenjetzt die komplexen ProjektkoeffizientenC genanntdie anderen die ursprünglichen Hardbasefür großen Sinus auch schon gesehendas ist der Ansatz und der funktioniert tatsächlichmit ein paar Körnchen Salz wenn diese Fonds Ursprünge hat aberansonsten hat sie wirklich heute anders verglichenanich kann meine Funktionausschlimmstenfallsunendlich viel von diesen Funktionen zusammenbauenetwas überraschend ist vielleicht das mir positive und negative Frequenzen hat der was drin ?? E hoch zwei Pi maldreizehnC durch die ich aber auch in gleichminus dreizehn dabei das Brauchtum wirklich Sinus und Kosinus wieder bilden zu können steht ja Kosinus Bussi mal Sinus mit dem Minus steht ein großesMinus ihm als Sinusdurch positive und negative Frequenzen schaff ich eswirklich ?? wieder großes Sinus zu kriegenkam und dieseFoyerkoeffizientenrichte man dannbillig mit einem Integral das war das integralüber eine Periode zum Beispiel von null bis T könnte genauso nehmen von minusgroß T halbes Plus die halbe oder irgendwas anders Hauptsache Periode langähmeins durch die Periodenlängedass das was wird werden Mittelwertwie bei ?? gibt's quer beim Effektivwertwurzelmittelwertquadratesteht auch einst durch die Länge über die ich dann integriereähm was wird integriert Skalarproduktdiese Funktion hier das war so komisch das die dann komplex korrigiert werden musste in diesemeigenwilligen Skalarproduktdes minus her durch die mal meine Originalfunktionvon DDTsind viele Tees glaube ich ähmdas hier mal T UnterstrichTT Unterstrich Daten unterstrichendazu machendas ist meine Integrationsvariablefunktioniert das bei der Foyerreihemeiner Originalfunktionwird zusammengesetztaus wie hoch die mal SoundsofunktionendieAmplituden und die Phasen stecken hier diesen Fourierkoeffizientendrines ist der spannende Vorgang hier aus der Originalfunktiondie Foyerkoeffizientenzu bestimmenden Beistrich welche Frequenzen drin sindmit welcher Rasen lange drin sind das Diktieren den verborgenund das war der Weg dahin müsste diese integral jetzt ausrechnen könnte man für jedes enden sagenwas das CN istnutzen Abschnitt drei in der Warenwelthabe ich keine Funktionenanbeliebigeher fein aufgelöstenZeitabständenich habe nicht so eine Funktion wie sie hier steht eine Kurve die ich durchzeichnen kann sowie das Signal in der Warenwelt nicht ausin der heutigen Warenwelt sieht das Signal so aus das es gestempelt ist für die Nummer fünfzehnich habe einenMesswert zum Zeitpunktvon minus null Sekundenich habe einen Messwert vom ersten Zeitpunkt eine Millisekunde an Messwert zwei Millisekunden drei Millisekunden vier Millisekunden fünf hundert und sechs Millionen sieben ?? und so weiterso sieht das Signal in der Warenwelt aus Punkt sie haben keine durchgezogeneKurve diese jetzt beliebig unter dem Mikroskop vergrößern könnensondernist denn bitte Wertebedeuteten einen Schritt weiter dass mandannauch nicht gleichmäßig abtastet sondern der abtastet da genau abtastet wussten bisschen schwieriger wirdKomma dass man ?? so weit treibendie typischedas typische Signal sieht so was siein irgend einem Bildungssystem kriegenoder bei der Messwerterfassungsowie das typische Signal aus ja meine Folge an Messwerten die allein festen Zeitabständengenommen sindkeine durchgezogeneKurveund damitist dieses integral hierMakulaturdieses integral zwangsläufigzu einer Summe werden man sich das überlegtin sie dir nehmendannan dir zwangsläufig eine Summe stehenich hab auch nicht mehr unendlich viele Frequenzen drin die möglich sind weil sie diese Auflösungder begrenzten Sicht kann ich so Schwingung der drin habenundich kann in dieser Wellenform nichtso eine Schwingung drin habenkönnte mit diesen Punkten gar nicht messendann die Folge ist wenn man das macht für echte Signaledas wird eine Summeund diese Summe hier wirdkeine rein mehr werden sondern ganz normale handelsüblicheSummeüber endlich wiederdiesen Dekaden auch noch eine Summe über endlich vieleam?? guckt sich also periodischeSignaleangemesseneperiodische Signale das solche Fälle noch einmalins periodisch wird für das heißt mir geht ?? wieder von vorne losKlammer zu machenanmodischhier wird es natürlichmit dem Ende hier weitergehendas interessiert mich gemessenen Signale zwar periodischegemessenen Signaledie Sache wird deutlich einfacher ich füge sie nicht die Formen vor so richtig spannendan dieser Stellehiervon die?? Foyerreihezu bestimmen?? das wird wesentlich einfacherPunkt das zeige ich einfach mal mitmatlap wir bauen also irgend ein Signalich nenne das malkreativ Fkames hatte genau das ?? nur aus Nullen bestehteineZeileeine Zeile tausend vierundzwanzigSpalten tausend zwanzig nur hintereinanderdas soll eine Periode meines Signals sein das endlich sowieso langweilig setzt irgendwas rein sagen wir eine Stelle dreizehn steht nach einzelnen Und-Zeichen werdende Stellezeigte sich auch einst Zensus etwas spannender das Signal des Coursera anguckendas ?? natürlich jetzt nur eine Periode was soll ich unendlich viele Perioden angeblich gebe eine Periode von diesem Signal an der Strecke zwischender weiteren zwei Stellen ist es einsansonsten ist es nurwieder sie typischerweisegezeichnetwird ist vielleichtetwasirreführendder sieht das ja so schönmit geraden Linienverbindet das irdische mit geraden Linienamsollte man sich nicht so vorstellenals Fieberkurveund sollte sich das eher vorstellen oder was das aus demRahmen soll sich das eher so vorstellendieeinzelnen Nadelnzu bestimmten Zeiten habe ich Ergebnisseund ansonsten nichts?? und hier bei der Nummer dreizehndas Ergebnis da oben und dann wieder null null null null nullweiteren ??ähm wieder ganz viele Nullenund den hier bei der zweiundvierzigdas Ergebnis oben so würde ich mir das vollständig durchgezogeneLinienandas trifft es nicht so ganz ich kenne ja nicht die Werte dazwischen ich kenne nicht den Wert an der Stelleeinundvierzig Komma sieben neun acht oder was auch immer ich weiß bestelle zwei vierzig ist einzelne Stelle einundvierzigstesnull und dazwischen habe ich kein Messwertalso lieber so vorstellennicht mit durchgezogeneLinienmit dem Signal gehe ich nunin diesen ganzen Formalismus reinführedie Foyeranalyseaus frage mich okay welche Schwingungen sind denn jetzt in diesem Signal mit seinen beiden Klicks periodisch wiederholttrenne ich sonst auch noch mal zeigen wie das dann mit der periodischen Wiederholung aussiehthabendas Karree auch locker hinkriegen ich baue eine Matrixeine Matrix die einfach mehrereWiederholungen davon hatalso mein Signal ist eine Matrix mit einerWahl ?? mit einer Zeile und tausend vierundzwanzig Spaltendieser sich einfachsiebenmal hintereinanderdas ?? war die Periode siehtin den sie dann siebenmal in diePeriode so ist das gemeint das ist das Signal was analysiert werden sollBeistrich unsinnigsiebenmal das ganze anzugeben die eine Periode bestimmter Signale komplettnurjedwede Funktion heißt RFTFestfoyerzwarjeder noch amerikanischer sprechen das für Transferdann effektivvon meinem Fnun?? anders herauskommtergebenPunkt man muss das ?? vom Typ herauskommtsehen das davonvon der Sorte Herrdes selber dieselbe Anzahl wieder rauskommt ich geh mit tausend in zwanzig Werten reinigt jedoch tausend zwanzig werden rausdas hier sind jetzt die frühen Koeffizientenwas in dem Gel steht FFCist das was passiertwenn mandiese Formel rechnet aber mit solchen?? gesendeten Signalen das hier das integral zu Summe wird diese Summe wird endlicheigentlich heißt das erst mal deinen DFTdie diskrete Fourier TransformationDFTdiskrete Fouriertransformationwarenaber man rechnetdasselbe Ergebnisals FFC-Phasegeht uns vor und alle Leute nennen es dann nur noch FF die Avus korrekterweise DFT heißen also wasähmdas Resultat ist die DFT und FFC ist eigentlich eine Art dieses Resultatauszurechnendie meisten Leute nennen es dann aber noch RFT Fast Fourier Transformeigentlich ein Algorithmusum solche Signaleschnell in ihreSinus vermengen Bestandteile zu zerlegen dass es FF Themensind der Stelle eingebaut ?? ich kriege die Koeffizientenwieder rausnundiese Koeffizientenguckenwas auchähm wenn man sie Koeffizienten angucken gegenüber kleineren Strecken die sind natürlichkomplette Zahlenverstehen oder komplette Zahlen was kein großes Wunder ist eine Minute nachdenkenhier steht jaeine üble komplexe Zahl mal meine Funktion wenn das wenig komplex wärees komischer sich die mit tausend vier zwanzig Fällen Zahlen rein und kommen mit tausend vier zwanzigkomplexen Zahlen wieder rauskönnen zumindest den Betrag Schallplatten?? PlotAppsvon Gsieht das dann ausdas wäre wenn sie wollendas Spektrumwas ist an Frequenzen drinnen unabhängig von der Phaseguck ich mir einfach an was in Frequenzen drin ist bestimmte sind drinnen bestimmte sind anscheinend nicht Rinder sind wieder hier dieses schon deutlich trennen die da wir deutlich nichtund so weiterdas sind die enthaltenenFrequenzenhaben ich nehme gleich mal etwas vernünftigesSignal hier mit den beiden Peaks das ist nicht so richtig lustigich nehme mal ein Signalan einen Sinus für einige Fälle mal guckenKomma daoh je zweimal Piauf jeden Fallich hätte gerne eine Sinus schwammige Fälle die?? zweimal dreizehn mal durch es auf den tausend vierundzwanzigalso Frequenz dreizehnmalhier müssen jetzt meine Zeit Wertereihendurch die tausend vierundzwanzigund die Zeitwertesindvon null bis Uppsala von null bis tausend dreiundzwanzigin Einzelschrittendas sollte jetzteine Sinuswelleseinmit tausend vier in zwanzig SamplesSelektor hier tausend vier und zwei ??Tempels eine Sinuswelle der dreizehn fachenFrequenzder GrundwellenlängenbestimmtPunkt an Bordbleibtdas dreizehnZiffer bis plausibelokay alsoeine Sinus verliert und damit gehe ich nun eine weitere sinnvolle ?? vor das wäre das gemessene Signalwird dadurch periodischer weiter links und rechts bis das männliche dieses Signal immer wieder immer wiederallem was seine Sinuswelleder dreizehnfache Frequenz bleibtuns das mit derTransformationam FFCvomelfterdaich Platte dasErgebnisdas Ergebniszu sehenaucheigenwillig ?? zwei Klicks in meinem SpektrumArm der eine peakder eine peak ist bei dreizehnwie sich das gehört ich geh mit Frequenz dreizehn reinen erwarte ich auch gefälligst das mir die Projekt Information acht Ober dreizehntes was drinnen diese Koeffizienten sind bei dreizehn nennenswert und ansonsten ?? sie nicht dazu nennenswertund sie kriegenungewolltund unvermeidbarauch noch eine Spiegelung hier am oberen Endesozusagen minus dreizehndas ist tausend vierundzwanzigminus dreizehn vom oberen Endedreizehn abwärtsandas sind die negativen Frequenzen hier aufgetaucht sind das ehrlicher vorgeführt und hiermit Sinus und Kosinus zu bildenbrauche ich die negativen Frequenzenwenn ich den Kosinus der ?? den Sinus der dreizehn fachen Frequenz haben will Komma hier C dreizehn und zehn minus dreizehnBeistrich hier nicht nurnach Sinus und Kosinus rausund hierbei der fast für jede Information schlecht sich das so nieder dass wir hierdie Frequenz minus dreizehn beitausend vier zwanzig minus dreizehn sind das ist der zweite peakhatte zwei Ticks die mir dann sagen was in Frequenzen drinnen istnundas heißt mit diesem Verfahren wieder so aus konnte man noch letzte Frequenzanalysemachen Sie Schweißen einfach nur Signal das sie analysieren wollen?? Schwingung ?? Signal was sie analysieren wollen in sonderFFC und dann ablesenlassen Frequenzen drinnen wardie starkannicht ganz es gibt ein Problem das ich zumindest mal andeuten wollte für den Lückentextsechzehnwird es danach ?? wenn sich hier eine krumme Frequenz habenimmer nicht dreizehn sondern nicht nehmendreizehn Kommakeine Ahnungirgendwas krummes vernichten krumme Frequenz habenach denselben ProzessFourier Transformationund Plattenden sie dass wir hier einenziemlich dicken Fuß kriechtder Sitz nicht hübschauf einer Stellemir das ?? ?? der Wagen damit dass er zwischen dreizehn und vierzehn antik sitztund sonst nichtses wird fürchterlich auswaren wenn ich die Frequenz war nicht ganz so hoch wähle dann wird's noch etwas deutlicher ich geh mal hier was mache Komma nicht dreizehn ich nehme mal dreidas drei Komma drei vier facheselber Protzpersönlich ihrerTransformationlautennurdas ist eher ungeschicktdass das hier so weit wirddas wohl Frequenzen drinnen scheinbar sind hohe Frequenzen drinnen dieser gar nicht haben will irgendwer dichtet mir hier hohe Frequenzen dareinkönnen zwar die Wellenformanguckenum zu sehenwas da passiertjedochden zurückwarum das Ärger geben könnte Platten wir mal die OriginalquelleFlotte Sieber wird jetzt durchgezogen nicht als einzelne Punkte sondern durchgezogen auch wenn ich eben gesagt habe das erkannt und IdeendieseWelle wird jetztFoyer analysiertnimmt dieses Ding finster ist hängt es unendlich oft aneinandervorstellen was dann passiertirgendeine Idee was wird passieren wenn ich das unendlich oft aneinanderhier ist die Periode nicht ganz durch es geht hier sofort wieder bei Null weiter einen Seitensprungdrindas nervt sobald diese Frequenz hier nichtganzzahliges Vielfaches ist original Frequenzeinen Sprung?? underscore Zimmermann guckenerreichte sieben ?? miteinander daich hänge mal die Originalwelle paar Mal hintereinander das man es sehen kannwie das denn aussieht was er jetztauch für das bisschen viel des?? soll dannwenn ich hier mal reinsuchewennes mir gelingt so hier sehen Sie was passiert das ist das Ende der ersten Periodehier kommt die nächste Periodean dieser Sprung hier ist keine gute Idee da stecken schön hohe Frequenzen drindie abrupt das ?? Übergang ist hoffentlichnicht ganz unanschaulich abrupten Übergang zum Aerofrequenzenhat ?? in seiner Wellehier kommen extremere besser weil ich da gar nicht drin haben will?? das ist der Ärgerund die Lösung ist folgendedie Lösung nennt sich Fenster Punkt das ist für die Nummer sechzehn für den Lückentext sechzehnanich nehmemeineQuelleund schneide da ein Fenster Rauswindungenund das blende ich ein und auseinerFensterfunktionPunkt das gelingt mir jetztdas ist die Idee mal eine FensterfunktionSchneider aus meiner Welle ein Stück raus zum Beispiel tausend vier zwanzig Samples das ist üblicheGröße fürFFC-Steinezum Beispiel tausend vier zwanzig raus aberdiesen Sprung hier möchte ich vermeidenund um den zu vermeidenblende ich ein und blende wieder ausmit einerFensterfunktionvon der Gipsmindestens Dutzenddie billigste Fensterfusion die man haben kann ist das man so Kosinus für micheinbindet und wieder ausblendetsanft reinund sanft wieder Punkt das nennt sich Hanandist eingebaut wenn wir Hantausend vier zwanzig sagenKomma dass sofortdas Essen wurde zum weich ein Ausblendenbeziehe einfach diese Kurve mit einem Signalund habe dann ebenweiche Einblendung von Signal und weiche Ausblendungund damit gehe ich dann ist die Transformationalso wir bilden erstmaligjedes Mal F einsähm eine Funktion Fmaldie FensterfunktionBindungsempfang schonandere sind Bleckmannoder Hemminghabenkommen jeweils die besten Eigenschaften wärennicht so spannendwichtig ?? Idee zu Filderstadtein Problem bestehtokay das ist jetzt meine Sinuswellemit der Sonne zu vielfachen Frequenzund jetztFensterein und ausgeblendetdurch MultiplikationKlammer zusage ich mir so vergessen das es erst eine Matrix istbesser dasunsnochmals das aneinsso sieht das dann aus können nur mit Mühe hier die Sinus Schwingung erkennen siehst ausgeblendetam Anfang angewendet am Anfang ausgeblendet am Ende wird habe ich kein Problem wenn das unendlich oft miteinander kommtdie Gitter schön sauber mit null raus und der Konzern sauber wieder mit null rein keine Sprünge mehr wenn ich das mit aufeinander hängedamit keinen Ärger mit der Fourier Transformationdie ja so tut als ob das unendlich oft ineinanderalso die für jede Summation voneinsBegehung und kannnach?? gibt es sie angucken sie sehen das dashoffentlich bedeutend schmaler geworden ist als das es eben waramweil mandendiesen Randeffektunterdrücktperfekt ganz lange drüber theoretisierendas geht nicht perfektan Komma kann diesen Randeffektschon unterdrückenzumindestkomplett aberweitgehend unterdrückt mit solchen Fensterfunktion?? ich sollte auch mal die Rückstandsformationüberhaupt zeigenalso aus meiner Originalfunktionwird jetztaus mein Original tausend vier zwanzig Daten so was mein Original tausend vier zwanzig Datenwerden jetzt tausend vierundzwanzigkomplexe Foyerkoeffizientennatürlich geht das zurück auchdie FFCinvers Fourier Transform in das Fass Fourier Transform die FFCaufdieG war das Ergebnis und dann müsste jetzt wieder das Original rauskommenFFCmit der FFC hinmit der IFCzurück das Jahr ist im Endeffekt diese Foyerreihezu bildenhaben dieser Teil hier ist die Fourier Transformationvon der Originalfunktionzu den Koeffizientendieser Teil ist die inverse Fouriertransformationvon den KoeffizientenzurOriginalfunktionsie beides ineinander machen soll das Original wieder rauskommendrei H anguckenanSiedler meine ich ganz unplausibel ausnur das von F einsvon meinem F eins billig die FouriertransformationGGwird sie wieder zurück von dem G die inverse Fourier TranslationRHwenn ich den Plattesehen sie in der Tatdas scheint das wieder zu seiner zuvor reingefüttert haben?? hier findet die Fourierreihe auf diese Weise mit dem FFC und demdie FFC aber in der Praxis geht's eben nochekligeRandnotizenhiermit den Fensterfunktion