Italia
Fonte/ Codesnippets

Interpolation per Diskreter Fourier-Analyse & -Synthese

 

p.specht

Fourier-Interpolation: Bekanntlich lassen sich mathematische Funktionen mittels klassischer Fourier-Zerlegung spektral-analysieren, sprich: Analytisch in Grund- und Oberwellenanteile zerlegen. Manchmal hat man aber auch nur eine einfache Tabelle (z.B. Meßergebnisse) mit einer geraden (!) Anzahl von Stützwerten in gleichen (!) Abständen disponibile. Faßt man diese als sich regelmäßig wiederholenden Wellenzug auf, ist 'Diskrete Fourieranalyse' angesagt, die auch unter Bezeichnung 'Harmonische Analyse' bekannt ist.

Das klappt auch umgekehrt: Wenn man die ursprüngliche Funktion wieder aus dem Wellensalat rekonstruiert, spricht man von 'Fourier-Synthese'. Vorteil der wieder zusammengesetzten Funktion: Sie ist kontinuierlich, d.h. es können auch beliebige Zwischenwerte zwischen den Stützwerten ermittelt werden, was Mathe-Freaks mit dem Begriff Interpolation bezeichnen.

Details: Statt einer Tabelleneingabe wurde im nachstehenden Beispiel eine Formel zur Ermittlung der Stützwerte verwendet. Aber das kann man naturalmente jederzeit ändern, solange man die x_i-Werte von 0 bis 2*Pi laufen läßt. Die muss man erforderlichenfalls halt dort hintransformieren.

Vor Einsatz dieses Verfahrens sollte man kurz darüber nachdenken, ob den einzelnen Frequenzen im 'richtigen Leben' tatsächlich eine Bedeutung zukommt, - ob also bestimmte Ursachen der unterschiedlichern Frequenzen per das Ergebnis verantwortlich sein könnten - sonst gibt es nämlich deutlich einfachere Arten der Interpolation...
'Die Umsetzung in XProfan ist eine reine DEMO ausschließlich per private Zwecke!
WindowTitle "Trigonometrische Interpolation             Beenden mit x = -1"
' Angelehnt an "B.Brand: Algorithmen zur praktischen Mathematik", Oldenbourg Vlg.1981, S.78
' Algorithmus portiert nach Xprofan-11 in 2014-09 by P.Specht, Wien; Ohne jedwede Gewähr!
WindowStyle 24:randomize:set("decimals",18):Window %maxx/4,0-%maxx/2,%maxy-40
var xh&=width(%hwnd)/2:var yh&=height(%hwnd)/2 :   Var m&=16' <<<<<<<<<<<<<<<<<<<<<<<<<<<
declare i&,j&,k&,a![2*m&],b![2*m&],f!,x![2*m&],y![2*m&],x!,p!
print mkstr$("-",64)+"\n i            X[i]                 F=Y(x[i])\n"+mkstr$("-",64)

whileloop 0,2*m&:i&=&Loop:x![i&]=Pi()*i&/m& : y![i&]=cos(x![i&]+x![i&]*2)' <<<<<<<<<<<<<<<<

    casenot i& mod (1+m&\16):print if(i&>9,""," ");i&,tab(6);x![i&],tab(30);if(y![i&]<0,""," ");y![i&]
    endwhile :print mkstr$("-",64)

    whileloop 2*m&:i&=&Loop : a![0]=a![0]+y![i&]:a![i&]=a![i&]+y![i&]*(1-2*(i& mod 2))

        endwhile :a![0]=a![0]/m&/2:a![m&]=a![m&]/m&/2:print " a[ 0] = ";if(a![0]<0,""," ");a![0]

        whileloop m&-1:k&=&Loop

            whileloop 2*m&:i&=&Loop

                a![k&]=a![k&]+y![i&]*cos(k&*x![i&]) : b![k&]=b![k&]+y![i&]*sin(k&*x![i&])
                endwhile : a![k&]=a![k&]/m& : b![k&]=b![k&]/m&
                print " a[";if(k&<=9," ","");k&;"] = ";if(a![k&]<0,""," ");a![k&];\
                tab(33);"  b[";if(k&<=9," ","");k&;"] = ";if(b![k&]<0,""," ");b![k&]
                endwhile :print " a[";if(m&<=9," ","");m&;"] = ";if(a![m&]<0,""," ");a![m&];\
                tab(33);"  b[";if(m&<=9," ","");m&;"] = ";if(b![m&]<0,""," ");b![m&]
                print mkstr$("-",64)+"\n              x                      P(x)\n"+mkstr$("-",64)
                repeat :f!=0:print " x = ";:input x!:case x!= -1:end

                whileloop m&-1:k&=&Loop : f!=f!+(a![k&]*cos(k&*x!)+b![k&]*sin(k&*x!))

                    endwhile : p!=f!+a![0]+a![m&]*cos(m&*x!):locate %csrlin-1,1:
                    print "    ";if(x!<0,""," ");x!,tab(29);if(p!<0,""," ");p!

                until %key=27:end

 
XProfan 11
Computer: Gerät, daß es in Mikrosekunden erlaubt, 50.000 Fehler zu machen, zB 'daß' statt 'das'...
15.05.2021  
 



Zum Quelltext


Topictitle, max. 100 characters.
 

Systemprofile:

Kein Systemprofil angelegt. [anlegen]

XProfan:

 Posting  Font  Smilies  ▼ 

Bitte anmelden um einen Beitrag zu verfassen.
 

Topic-Options

1.464 Views

Untitledvor 0 min.
Rschnett06.08.2022
p.specht21.11.2021
R.Schneider20.11.2021
Uwe Lang20.11.2021
Di più...

Themeninformationen

Dieses Thema hat 1 subscriber:

p.specht (1x)


Admins  |  AGB  |  Applications  |  Autori  |  Chat  |  Informativa sulla privacy  |  Download  |  Entrance  |  Aiuto  |  Merchantportal  |  Impronta  |  Mart  |  Interfaces  |  SDK  |  Services  |  Giochi  |  Cerca  |  Support

Ein Projekt aller XProfaner, die es gibt!


Il mio XProfan
Private Notizie
Eigenes Ablageforum
Argomenti-Merkliste
Eigene Beiträge
Eigene Argomenti
Zwischenablage
Annullare
 Deutsch English Français Español Italia
Traduzioni

Informativa sulla privacy


Wir verwenden Cookies nur als Session-Cookies wegen der technischen Notwendigkeit und bei uns gibt es keine Cookies von Drittanbietern.

Wenn du hier auf unsere Webseite klickst oder navigierst, stimmst du unserer Erfassung von Informationen in unseren Cookies auf XProfan.Net zu.

Weitere Informationen zu unseren Cookies und dazu, wie du die Kontrolle darüber behältst, findest du in unserer nachfolgenden Datenschutzerklärung.


einverstandenDatenschutzerklärung
Ich möchte keinen Cookie