|
25 | 25 |
|
26 | 26 | s = 0.1*np.sin(2*np.pi*t) + cnse # the signal
|
27 | 27 |
|
28 |
| -fig, axx = plt.subplots(3, 2) |
| 28 | +fig, axes = plt.subplots(nrows=3, ncols=2) |
29 | 29 |
|
30 | 30 | # plot time signal:
|
31 |
| -axx[0, 0].plot(t, s) |
32 |
| -axx[0, 0].set_xlabel("Time $t$") |
33 |
| -axx[0, 0].set_ylabel("Signal $s(t)$") |
| 31 | +axes[0, 0].plot(t, s) |
| 32 | +axes[0, 0].set_xlabel("Time") |
| 33 | +axes[0, 0].set_ylabel("Amplitude") |
34 | 34 |
|
35 | 35 | # plot different spectrum types:
|
36 |
| -axx[1, 0].magnitude_spectrum(s, Fs=Fs) |
37 |
| -axx[2, 0].phase_spectrum(s, Fs=Fs) |
38 |
| -axx[0, 1].remove() # don't display empty ax |
39 |
| -axx[1, 1].magnitude_spectrum(s, Fs=Fs, scale='dB') |
40 |
| -axx[2, 1].angle_spectrum(s, Fs=Fs) |
| 36 | +axes[1, 0].magnitude_spectrum(s, Fs=Fs) |
| 37 | +axes[1, 1].magnitude_spectrum(s, Fs=Fs, scale='dB') |
| 38 | +axes[2, 0].phase_spectrum(s, Fs=Fs) |
| 39 | +axes[2, 1].angle_spectrum(s, Fs=Fs) |
| 40 | +axes[0, 1].remove() # don't display empty ax |
41 | 41 |
|
42 | 42 | fig.tight_layout()
|
43 | 43 | plt.show()
|
0 commit comments