blob: 3f1e815f4c719fced77b825093677cdeaafbe7a3 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
|
NOTE OSZILLATOR
L X.K=X.J+Y.J*DT
N X=2
L Y.K=(Y.J+DT*(F/M-X.J*OMEGANULLQUADRAT.J))/(1+GAMMA.J*DT)
N Y=3
C M=5
NOTE
NOTE linearer harmonischer Oszillator mit BETA=0 und F=0
NOTE
NOTE gedaempfter Oszillator mit BETA<>0 und F=0
NOTE
NOTE allgemeiner Oszillator mit BETA<>0 und F=f(TIME)
C BETA=0.5
A GAMMA.K=BETA/M
C F=0
C K=1
A OMEGANULLQUADRAT.K=K/M
NOTE hier heisst eine Konstante"K". DYNAMO verwechselt das nicht mit .K !!
C DT=0.01
NOTE DT WIRD EXTRA SO KLEIN GEWAEHLT; DAMIT DIE ANNAEHERUNG GUT IST
NOTE
NOTE DAS GEHT AUF KOSTEN DER RECHENZEITEN !!!
C LENGTH=68
C PLTPER=1
PLOT Y
|
