-
Notifications
You must be signed in to change notification settings - Fork 0
Expand file tree
/
Copy pathrun.py
More file actions
80 lines (58 loc) · 1.78 KB
/
Copy pathrun.py
File metadata and controls
80 lines (58 loc) · 1.78 KB
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
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
import numpy
import sympy
from sympy.abc import s
from sympy import exp, pi
from iplot import butterworth, reserve
import signal_examples
from signal_examples import time, noise, slider, tf, add, mul, delta
import signals
import wrapper
from wrapper import plot, complex_plot, ortskurve, bode, wok, impuls, sprung, harmonisch, response
wrapper.register(s)
signal_examples.register(s)
g = 1/(s**2-4)
gr = 1 + .2 * s
go = 100 * gr * g
gg = 96/100 * go / (1 + go)
response(tf(g, slider().p()).p())
# pt1 = 1 / (1 + s)
# pt2_frei = 1 / (s**2 + 1)
# pt2_leicht = 1 / ((s+.1)**2 + 1)
# pt2_kritisch = 1 / (s+1)**2
# pt2_stark = 1 / (s+10) / (s+.1)
# pt2_instabil = 1 / (s**2 - 1)
# bw = butterworth(1, s, wg=1)
# f = (s-3) / (s+.01)**2 / ((s+2)**2+1) # 1/((s+.05)**2+1)/((s+.5)**2+1)
# h = s/(s-2)/(s-8)
# f = 5 / (1 + 1.414 * s + s**2) / s
# r = .1 * (1 + .536 * s)
# g = 1 / (s-1) / (s+2)
# f2 = 1 / ((s+1)**2+1)
# f_aufg6_8 = 1/(s+8)/(s+2)
# f_aufg6_9 = 1/s/(s+1)
# f_aufg6_10 = 1/(s**2+10*s+11)
# f_aufg6_11 = s/(s**2+110*s+1000)
# f_aufg6_12 = 1e7*(s+10)**2/s**2/(s+100)/(s+1000)
# f_aufg6_13 = 4e4*(s+2)/(s+20)/(s+200)**2
# f_aufg6_15 = 1e4 / (s+10)/(s+20)/(s+30) * exp(-.1 * s)
# f_aufg6_15t = 1e4 / (s+10)/(s+20)/(s+30)
# f_test = 1/((s+1)**2+4) * exp(-.1*s)
# f_aufg6_17 = 7 /(s**2+2*s+2)/(s-3)
# f_aufg6_18 = 1 /(s+1)/s**3
# f_aufg7_2 = s/(s-2)/(s-8)
# f_aufg7_3 = (s+4)/(s-3)/(s**2+4*s+53)
# f_aufg7_4 = 1/(s+2)**2/(s+5)/(s+9)
# f_aufg7_5 = 1/(s-2)/(s**2+8*s+25)
# probe_4 = 1/(s+1)**2
# gz = .5 * (z+1) / z
# g = 1/(s**2-4)
# gr = 1 + .2 * s
# go = 100 * gr * g
# gg = 96/100 * go / (1 + go)
# gg = gg.expand().factor().nsimplify().expand().simplify()
# wok(gr * g)
# complex_plot(gg)
# print(gg)
# complex_plot(gg)
# usercontrol(gg, dt=.003, shown=1000)
# timeresponse(lambda t: 2*(t%2-1), gg)