Quanten.de Diskussionsforum  

Zur?ck   Quanten.de Diskussionsforum > Plauderecke

Hinweise

Plauderecke Alles, was garantiert nichts mit Physik zu tun hat. Seid nett zueinander!

Antwort
 
Themen-Optionen Ansicht
  #21  
Alt 27.02.12, 19:37
Benutzerbild von richy
richy richy ist offline
Singularität
 
Registriert seit: 01.05.2007
Ort: karlsruhe
Beitr?ge: 4.170
Standard AW: Math DZGL Katalog

*************************
Schwingungs (harmonische) DZGL
*************************
x[1]:=cos(dx); dx ist die Schrittweite
y[1]:=sin(dx);
x[n+1]:=x[n]*c1 - y[n]*c2; # -> Cosinuswerte
y[n+1]:=y[n]*c1 + x[n]*c2; # -> Sinuswerte

Herleitung :

Eine exakte ! numerische Simulation der harmonischen Schwingung ergbt sich aus den Additionstheoremen.

Additionstheoreme
**************




Fuer eine Startwert x und eine Schrittweite dx erhaelt man aus der ersten Gleichung:

cos(x+dx)=cos(x)*cos(dx)-sin(dx)*sin(x)
Anm: Ziel ist es eine Substitution X(n)=cos(x+n*dx) durchzufuehren.

Die Werte fuer die Sinusfunktion koennte man ueber die Beziehung sin(x)^2+cos(x)^2=1 ausdruecken. Allerdings weist eine solche lteration aufgrund des Wurzelvorzeichens lediglich enen Gueltigkeitsbereich im Intervall 0-180 Grad auf :



Eleganter,allgemeiner und effizienter ist es die Sinusfunktion ebenfalls iterativ aus dem Additionstheorem zu berechnen :

cos(x+dx)=cos(x)*cos(dx)-sin(dx)*sin(x)
sin(x+dx)=sin(x)*cos(dx)+sin(dx)*cos(x)

dx stellt in einer iterativen Betrachtung eine Konstante dar =>

cos(x+dx)=c1*cos(x)-c2*sin(x)
sin(x+dx)=c1*sin(x)+c2*cos(x)

die Gleichungen sind exakt fuer :
c1=cos(dx)
c2=sin(dx)

Erstetz man C1 ud C2 durch Taylorapproximationen erhaelt man 1.Ordnung eine Variante der Euler Chromer Integration. (Nur eine formale Variante !)
c1=1
c2=0

Die "Merman Naeherung" entspricht eine Taylorapproximation 2.Ordung die sich auch graphisch herleiten laesst.
c1:=1-dx^2/2;
c2:=sqrt(dx^2-dx^4/4);

(Notwendig sind solche Approximationen nicht.)

Mittels Substitution erhalten wir ein DZL System fuer X(n),Y(n)
(Im folgenden in Kleinbuchstaben geschrieben)

Zitat:
Zitat von Quellcode
restart;

N:=360; # Anzahl Abtastwerte
P:=3; # Anzahl Perioden
M:=P*N;
dx:=evalf(2*Pi/N);

c1:=1-dx^2/2; #(Beispiel Menrman Naeherung)
c2:=sqrt(dx^2-dx^4/4); #(Beispiel Menrman Naeherung)

x[1]:=c1;
y[1]:=c2;

for n from 1 to M do
x[n+1]:=x[n]*c1 - y[n]*c2; # -> Cosinuswerte
y[n+1]:=y[n]*c1 + x[n]*c2; # -> Sinuswerte
od:

# DRUCKROUTINEN
druck1:=seq([i,y[i]],i=1..M):
druck2:=seq([i,sin(i*2*Pi/360)],i=1..M):
druck3:=seq([i,x[i]],i=1..M):
druck4:=seq([i,cos(i*2*Pi/360)],i=1..M):
plot([[druck1],[druck2],[druck3],[druck4]]);
Der Cordic Algo baut auf die vorgestellte DZGL auf.
Die Taylor (graphischen) Naeherungen sind nur von theoretischen Interesse, da die Methode der Additionstheoreme exakt ist. Diese stellt keine integrale Appproximation dar ! Bereits eine Taylor, Merman Approximation 2.ter Ordnung ist sehr genau.



Da die Approximation nach den Additionstheoremen exakt ist, stellt die Loesung der DZGL die harmonische Schwingung dar. Die Taylor Naeherungen ergeben dagegen sehr komplexe Ausdruecke fuer die Loesung.


(Kommentare bitte im Diskussionsthread :-)

Ge?ndert von richy (27.02.12 um 20:06 Uhr)
Mit Zitat antworten
  #22  
Alt 20.03.12, 10:39
Benutzerbild von richy
richy richy ist offline
Singularität
 
Registriert seit: 01.05.2007
Ort: karlsruhe
Beitr?ge: 4.170
Standard AW: Math DZGL Katalog

Weitere Spezielle Lukasfolgen (in P Fibonacci Form)
**************************************
Fuer P=2 wird die P Fibonacci Form auch als Pell Folge bezeichnet :

Gegeben sei folgende Differenzengleichung :

y(0)=1, y(1)=1
y(k+1)=p*y(k)+y(k-1)


Der Wert y(k+1)/y(k) konvergiert gegen den Wert folgender Iteration :

z(k+1)=p+1/z(k), z(0)=1

und damit gegen
z=1/2*(p+-Wurzel (p^2+4))

Fuer p=2 somit gegen 1+Wurzel(2)

y[2] := 3
y[3] := 7
y[4] := 17
y[5] := 41
y[6] := 99
y[7] := 239
y[8] := 577
y[9] := 1393
y[10] := 3363 ...

Weitere Attraktoren der P-Lukasfolge P=1,2,3,...

Mit Zitat antworten
Antwort

Lesezeichen

Themen-Optionen
Ansicht

Forumregeln
Es ist Ihnen nicht erlaubt, neue Themen zu verfassen.
Es ist Ihnen nicht erlaubt, auf Beitr?ge zu antworten.
Es ist Ihnen nicht erlaubt, Anh?nge hochzuladen.
Es ist Ihnen nicht erlaubt, Ihre Beitr?ge zu bearbeiten.

BB-Code ist an.
Smileys sind an.
[IMG] Code ist an.
HTML-Code ist aus.

Gehe zu


Alle Zeitangaben in WEZ +1. Es ist jetzt 17:19 Uhr.


Powered by vBulletin® Version 3.8.8 (Deutsch)
Copyright ©2000 - 2024, vBulletin Solutions, Inc.
ScienceUp - Dr. Günter Sturm