()
| 1448 | """Cumulative integration with proper corrections.""" |
| 1449 | |
| 1450 | def plot_lundeby(): |
| 1451 | c0, c1, interIdx, BGL = lundebyParams |
| 1452 | fig = plt.figure(figsize=(10, 5)) |
| 1453 | ax = fig.add_axes([0.08, 0.15, 0.75, 0.8], polar=False, |
| 1454 | projection='rectilinear', xscale='linear') |
| 1455 | line = c1*inputSignal.timeVector + c0 |
| 1456 | ax.plot(inputSignal.timeVector, 10*np.log10(inputSignal.timeSignal**2), label='IR') |
| 1457 | ax.axhline(y=10*np.log10(BGL), color='#1f77b4', label='BG Noise', c='red') |
| 1458 | ax.plot(inputSignal.timeVector, line, label='Late slope', c='black') |
| 1459 | ax.axvline(x=interIdx/inputSignal.samplingRate, label='Truncation point', c='green') |
| 1460 | ax.grid() |
| 1461 | ax.set_xlabel('Time [s]') |
| 1462 | ax.set_ylabel('Amplitude [dBFS]') |
| 1463 | plt.title('{0:.0f} [Hz]'.format(band)) |
| 1464 | ax.legend(loc='best', shadow=True, fontsize='x-large') |
| 1465 | |
| 1466 | # timeSignal = inputSignal.timeSignal[:] |
| 1467 | # Substituted by SignalObj.crop in analyse function |
no test coverage detected