|
23 | 23 | cnse = np.convolve(nse, r)*dt
|
24 | 24 | cnse = cnse[:len(t)]
|
25 | 25 |
|
26 |
| -s = 0.1*np.sin(2*np.pi*t) + cnse # the signal |
| 26 | +s = 0.1*np.sin(4*np.pi*t) + cnse # the signal |
27 | 27 |
|
28 |
| -fig, axes = plt.subplots(nrows=3, ncols=2) |
| 28 | +fig, axes = plt.subplots(nrows=3, ncols=2, figsize=(7, 7)) |
29 | 29 |
|
30 | 30 | # plot time signal:
|
31 |
| -axes[0, 0].plot(t, s) |
| 31 | +axes[0, 0].set_title("Signal") |
| 32 | +axes[0, 0].plot(t, s, color='C0') |
32 | 33 | axes[0, 0].set_xlabel("Time")
|
33 | 34 | axes[0, 0].set_ylabel("Amplitude")
|
34 | 35 |
|
35 | 36 | # plot different spectrum types:
|
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) |
| 37 | +axes[1, 0].set_title("Magnitude Spectrum") |
| 38 | +axes[1, 0].magnitude_spectrum(s, Fs=Fs, color='C1') |
| 39 | + |
| 40 | +axes[1, 1].set_title("Log. Magnitude Spectrum") |
| 41 | +axes[1, 1].magnitude_spectrum(s, Fs=Fs, scale='dB', color='C1') |
| 42 | + |
| 43 | +axes[2, 0].set_title("Phase Spectrum ") |
| 44 | +axes[2, 0].phase_spectrum(s, Fs=Fs, color='C2') |
| 45 | + |
| 46 | +axes[2, 1].set_title("Angle Spectrum") |
| 47 | +axes[2, 1].angle_spectrum(s, Fs=Fs, color='C2') |
| 48 | + |
40 | 49 | axes[0, 1].remove() # don't display empty ax
|
41 | 50 |
|
42 | 51 | fig.tight_layout()
|
|
0 commit comments