srenblad
/
artiq
forked from M-Labs/artiq
1
0
Fork 0
Code Pull Requests Activity

test/fir: look at overshoot behavior

This commit is contained in:
Robert Jördens 2017-06-12 20:04:58 +02:00
parent 6ac9d0c41e
commit 39a1dcbb3d
1 changed files with 18 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

24
artiq/gateware/test/dsp/fir.py
View File

@ -22,17 +22,23 @@ class Transfer(Module):
yield yield
yi[:] = (yield from [(yield o) for o in self.dut.o]) yi[:] = (yield from [(yield o) for o in self.dut.o])
def run(self, samples, amplitude=1.): def run(self, samples, amplitude=1., seed=None):
if seed is not None:
np.random.seed(seed)
w = 2**(self.dut.width - 1) - 1 w = 2**(self.dut.width - 1) - 1
x = np.round(np.random.uniform( x = np.round(np.random.uniform(
-amplitude*w, amplitude*w, samples)) -amplitude*w, amplitude*w, samples))
y = np.empty_like(x) y = self.run_data(x)
run_simulation(self, [self.drive(x), self.record(y)],
vcd_name="fir.vcd")
x /= w x /= w
y /= w y /= w
return x, y return x, y
def run_data(self, x):
y = np.empty_like(x)
run_simulation(self, [self.drive(x), self.record(y)],
vcd_name="fir.vcd")
return y
def analyze(self, x, y): def analyze(self, x, y):
fig, ax = plt.subplots(3) fig, ax = plt.subplots(3)
ax[0].plot(x, "c-.", label="input") ax[0].plot(x, "c-.", label="input")
@ -81,7 +87,7 @@ class UpTransfer(Transfer):
def _main(): def _main():
if True: if True:
coeff = fir.halfgen4_cascade(8, width=.4, order=8) coeff = fir.halfgen4_cascade(2, width=.4, order=8)
dut = fir.ParallelHBFUpsampler(coeff, width=16) dut = fir.ParallelHBFUpsampler(coeff, width=16)
# print(verilog.convert(dut, ios=set([dut.i] + dut.o))) # print(verilog.convert(dut, ios=set([dut.i] + dut.o)))
tb = UpTransfer(dut) tb = UpTransfer(dut)
@ -91,7 +97,13 @@ def _main():
# print(verilog.convert(dut, ios=set(dut.i + dut.o))) # print(verilog.convert(dut, ios=set(dut.i + dut.o)))
tb = Transfer(dut) tb = Transfer(dut)
x, y = tb.run(samples=1 << 10, amplitude=.5) if True:
x, y = tb.run(samples=1 << 10, amplitude=.5, seed=0x1234567)
else:
x = np.zeros(100)
x[:50] = 1 << 8
x[50:] = 1 << 13
y = tb.run_data(x)
tb.analyze(x, y) tb.analyze(x, y)
plt.show() plt.show()