1import numpy as np
2from pylab import figure, clf, plot, xlabel, ylabel, xlim, ylim, title, grid, axes, show,semilogx, semilogy
3import scipy.fftpack
4import os.path
5import struct
6import argparse
7import sys
8from inspect import getsourcefile
9
10current_path = os.path.abspath(getsourcefile(lambda:0))
11current_dir = os.path.dirname(current_path)
12parent_dir = current_dir[:current_dir.rfind(os.path.sep)]
13
14sys.path.insert(0, parent_dir)
15
16import PatternGeneration.DebugTools as d
17
18f = "f16"
19
20inputPath = os.path.join(parent_dir,"Patterns","DSP","Filtering","BIQUAD","BIQUAD%s" % f.upper(),"BiquadInput1_%s.txt" % f )
21refPath = os.path.join(parent_dir,"Patterns","DSP","Filtering","BIQUAD","BIQUAD%s" % f.upper(),"BiquadOutput1_%s.txt" % f)
22outputPath= os.path.join(parent_dir,"Output","DSP","Filtering","BIQUAD","BIQUAD%s" % f.upper(),"Output_1.txt")
23
24
25
26
27inSig = d.readF16Pattern(inputPath)
28
29refSig = d.readF16Pattern(refPath)
30
31sig = d.readF16Output(outputPath)
32
33
34figure()
35plot(inSig)
36figure()
37plot(refSig)
38figure()
39plot(sig)
40
41#print(d.SNR(refSig,sig))
42
43#figure()
44#plot(np.unwrap(np.angle(refSig)))
45#figure()
46#plot(np.unwrap(np.angle(sig)))
47#figure()
48#plot(np.unwrap(np.angle(sig)) - np.unwrap(np.angle(refSig)))
49show()#