CoCalc -- 435.sagews

All published worksheets from http://sagenb.org

Project: sagenb.org published worksheets

Views: ¹⁶⁸⁷³³
Image: ubuntu2004

var('t f')

a = 1500*t+1
b = 5*sin(1400*t)-cos(14000*t+pi/6)
c = sin(40000*t)+2.5

sample_rate = 10000
samples = 20
bits = 8
total_time = (samples / sample_rate)
resol = 256
range = 5

reso_range = range/resol

lookup = {}

def adc(x,f):
    x = floor(x * sample_rate)/sample_rate
    if x in lookup:
        return lookup[x]
    else:
        n = max( min( floor( f(x)/reso_range ) ,255),0)*reso_range
        lookup[x] = n
        return n

def adc2(x,f):
    x = floor(x * sample_rate)/sample_rate
    if x in lookup:
        return lookup[x]
    else:
        n = max( min( floor( f(x)/reso_range ) ,255),0)
        lookup[x] = n
        return n

def getTable(f):
    s = "<table border=1><thead><tr><th>Sample</th><th>Time (sec.)</th><th>v(t) (volts)</th><th>ADC value</th><th>Voltage read by ADC</th></tr></thead><tbody>"
    for i in xrange(0,20):
        t = i / 10000
        v = f(t)
        d = adc2(t,f)
        e = d * reso_range
        s = s + "<tr><td>%d</td><td>%f</td><td>%f</td><td>%d</td><td>%f</td></tr>" % (i,t,v,d,e)
    s = s + "</tbody></table>"
    return s

lookup = {}
show(plot(a(t),(0,total_time))+plot(lambda x: adc(x,a),(0,total_time),color="orange"),ymin=0,ymax=5,figsize=(8,10),dpi=75)

html(getTable(a))

Sample	Time (sec.)	v(t) (volts)	ADC value	Voltage read by ADC
0	0.000000	1.000000	51	0.996094
1	0.000100	1.150000	58	1.132812
2	0.000200	1.300000	66	1.289062
3	0.000300	1.450000	74	1.445312
4	0.000400	1.600000	81	1.582031
5	0.000500	1.750000	89	1.738281
6	0.000600	1.900000	97	1.894531
7	0.000700	2.050000	104	2.031250
8	0.000800	2.200000	112	2.187500
9	0.000900	2.350000	120	2.343750
10	0.001000	2.500000	128	2.500000
11	0.001100	2.650000	135	2.636719
12	0.001200	2.800000	143	2.792969
13	0.001300	2.950000	151	2.949219
14	0.001400	3.100000	158	3.085938
15	0.001500	3.250000	166	3.242188
16	0.001600	3.400000	174	3.398438
17	0.001700	3.550000	181	3.535156
18	0.001800	3.700000	189	3.691406
19	0.001900	3.850000	197	3.847656

lookup = {}
show(plot(c(t),(0,total_time))+plot(lambda x: adc(x,c),(0,total_time),color="orange"),ymin=0,ymax=5,figsize=(8,10),dpi=75)

lookup = {}
html(getTable(b))

Sample	Time (sec.)	v(t) (volts)	ADC value	Voltage read by ADC
0	0.000000	-0.866025	0	0.000000
1	0.000100	1.043245	53	1.035156
2	0.000200	2.365261	121	2.363281
3	0.000300	2.027593	103	2.011719
4	0.000400	1.668638	85	1.660156
5	0.000500	2.896683	148	2.890625
6	0.000600	4.600232	235	4.589844
7	0.000700	4.775020	244	4.765625
8	0.000800	3.835106	196	3.828125
9	0.000900	3.911727	200	3.906250
10	0.001000	5.304134	255	4.980469
11	0.001100	5.974470	255	4.980469
12	0.001200	4.925392	252	4.921875
13	0.001300	3.853468	197	3.847656
14	0.001400	4.333640	221	4.316406
15	0.001500	5.208722	255	4.980469
16	0.001600	4.517448	231	4.511719
17	0.001700	2.760256	141	2.753906
18	0.001800	2.081190	106	2.070312
19	0.001900	2.723773	139	2.714844

lookup = {}
show(plot(b(t),(0,total_time))+plot(lambda x: adc(x,b),(0,total_time),color="orange"),figsize=(8,10),dpi=75)

lookup = {}
html(getTable(c))

Sample	Time (sec.)	v(t) (volts)	ADC value	Voltage read by ADC
0	0.000000	2.500000	128	2.500000
1	0.000100	1.743198	89	1.738281
2	0.000200	3.489358	178	3.476562
3	0.000300	1.963427	100	1.953125
4	0.000400	2.212097	113	2.207031
5	0.000500	3.412945	174	3.398438
6	0.000600	1.594422	81	1.582031
7	0.000700	2.770906	141	2.753906
8	0.000800	3.051427	156	3.046875
9	0.000900	1.508221	77	1.503906
10	0.001000	3.245113	166	3.242188
11	0.001100	2.517702	128	2.500000
12	0.001200	1.731745	88	1.718750
13	0.001300	3.486628	178	3.476562
14	0.001400	1.978449	101	1.972656
15	0.001500	2.195189	112	2.187500
16	0.001600	3.420026	175	3.417969
17	0.001700	1.602072	82	1.601562
18	0.001800	2.753823	140	2.734375
19	0.001900	3.066108	156	3.046875