Mercurial > pygpibtoolkit / file comparison
comparison: pygpibtoolkit/HP3562A/trace_decoder.py
pygpibtoolkit/HP3562A/trace_decoder.py
- changeset 46
- 235686915f92
- parent 40
- 1bbea188a7e5
- child 53
- 8e32c806fcdd
equal
deleted
inserted
replaced
| 45:cafc62915d39 | 46:235686915f92 |
|---|---|
| 112 print e | 112 print e |
| 113 sys.exit(1) | 113 sys.exit(1) |
| 114 | 114 |
| 115 if options.displayheader: | 115 if options.displayheader: |
| 116 print format_header(header, HEADER, 100) | 116 print format_header(header, HEADER, 100) |
| 117 f0 = header['Start freq value'] | |
| 118 dx = header['Delta X-axis'] | |
| 119 n = header['Number of elements'] | |
| 120 x = numpy.linspace(f0, f0+dx*n, len(data)) | |
| 121 y = data.copy() | |
| 117 if options.plot: | 122 if options.plot: |
| 118 f0 = header['Start freq value'] | |
| 119 dx = header['Delta X-axis'] | |
| 120 n = header['Number of elements'] | |
| 121 x = numpy.linspace(f0, f0+dx*n, len(data)) | |
| 122 y = data.copy() | |
| 123 | |
| 124 import pylab | 123 import pylab |
| 125 if options.ymode != "lin": | 124 if options.ymode != "lin": |
| 126 minv = min(y[y>0]) | 125 minv = min(y[y>0]) |
| 127 y[y==0] = minv | 126 y[y==0] = minv |
| 128 y = numpy.log10(y) | 127 y = numpy.log10(y) |
| 130 y = y*10 | 129 y = y*10 |
| 131 pylab.ylabel('db') | 130 pylab.ylabel('db') |
| 132 pylab.grid() | 131 pylab.grid() |
| 133 pylab.plot(x, y) | 132 pylab.plot(x, y) |
| 134 pylab.xlabel('frequency') | 133 pylab.xlabel('frequency') |
| 135 pylab.show() | 134 pylab.show() |
| 136 | 135 return header, x, y |
| 137 | 136 |
| 138 if __name__ == "__main__": | 137 if __name__ == "__main__": |
| 139 main() | 138 h, x, y = main() |
