python_MM/curveplot1.py at master · Isaac2k1/python_MM · GitHub

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
#!/usr/bin/env python3
# runs in Python 3.7.4 on MacOS
# curveplot1.py
import matplotlib.pyplot as plt
import numpy as np
from math import *

"""
sON     =  "abs(   (0>sin(200*pi*Ch(""t"")+A )) AND   (0>sin(200*pi*Ch(""t"")+B )))"
sOFF    =  "abs(NOT(0>sin(200*pi*Ch(""t"")+A )) OR NOT(0>sin(200*pi*Ch(""t"")+B )))"
sSINrec =  "abs(sin(100*pi*Ch(""t"")))"
sDECH   = "pp1* (arctan(tan(100*pi*Ch(""t""))) + pp2)"

 sFormula1 = "C*(" & sON  & "*" & sSINrec & "+" & sOFF & " * " & sDECH & ") + O"

Sub setDCI(iCurve,sSymbols,sValues,sUnit,Pts,t0,dt)
'      sSymbols = array("Ampl" , "Start"     ,  "Stop"         , "pp1"     , "pp2"     , "sUnit", "O"   , "Pts", "t0", "dt")
       sSymbols = array("C"    , "A"         ,  "B"            , "pp1"     , "pp2"     , "sUnit", "O"   , "Pts", "t0", "dt")
select case iCurve
 Case 1 sValues = array("0"    , pi/3600*3575, pi/360000*56285 , -0.010378 , -94.79208 , "A"    ,  1000 , 10001,  0.0, 2E-5)
 Case 2 sValues = array("1000" , pi/3600*3575, pi/360000*56285 , -0.010378 , -94.79208 , "A"    ,  0    , 10001,  0.0, 2E-5)
"""
A = pi / 360000 * 416285  # 416285
B = pi / 360000 * 357500
C = 1000
pp1 = -0.010378
pp2 = -94.79208
O = 0

Fs = 1000000
f = 100
sample = 20000
t = np.arange(sample)

sON = (0 < np.sin(200 * pi * t / Fs + A)) * (0 > np.sin(200 * pi * t / Fs + B))
sOFF = 1 - sON
sSINrec = abs(np.sin(100 * pi * t / Fs))
sDECH = pp1 * (np.arctan(np.tan(100 * pi * t / Fs)) + pp2)

# sFormula1 = "C*(" & sON  & "*" & sSINrec & "+" & sOFF & " * " & sDECH & ") + O"
y1 = C * sSINrec + O
y2 = C * sDECH + O
# y3      = C * ( sSINrec * sON + sDECH * (1 - sON )) + O
y3 = C * (sON * (sSINrec - sDECH) + sDECH) + O
# y3      =

plt.plot(t, y1, color="b")
plt.plot(t, y2, color="g")
plt.plot(t, y3, color="r")
plt.xlabel("sample(n)")
plt.ylabel("voltage(V)")
plt.show()