def test_create_lscp_realisation():
"""Test that create method is working at least."""
N = 256
c = 1.4
Fs = 8.0
m = 12.0
d = 0.4
x, H, Rx = utils.create_lscp_realisation(N, c, Fs, m, d)
assert len(x) == N
assert_allclose(H.shape, (N, N))
assert_allclose(Rx.shape, (N, N))
def test_create_lscp_realisation():
"""Test that create method is working at least."""
N = 256
c = 1.4
Fs = 8.0
m = 12.0
d = 0.4
x, H, Rx = utils.create_lscp_realisation(N, c, Fs, m, d)
assert len(x) == N
assert_allclose(H.shape, (N, N))
assert_allclose(Rx.shape, (N, N))
def generate_lsp_processes_figures():
"""Generates Figure 1 and Figure 2 from the paper."""
N = 256
Fs = 7.0
t_vector = np.arange(-(N/2), N/2, 1) / Fs
X, Y = np.meshgrid(t_vector, t_vector)
random_seed = 1
x_1, _, Rx_1 = utils.create_lsp_realisation(N, 1.1, Fs, random_seed)
x_2, _, Rx_2 = utils.create_lsp_realisation(N, 3.0, Fs, random_seed)
x_3, _, Rx_3 = utils.create_lscp_realisation(N, 1.1, Fs, 2.0, -2.0, random_seed)
x_4, _, Rx_4 = utils.create_mlsp_realisation(N, (1.1, 20.0), Fs, random_seed)
Rxs = [
('$R_x^{LSP}, c=1.1$', Rx_1),
('$R_x^{LSP}, c=3$', Rx_2),
('$R_x^{LSCP}, c=1.1, m=2, d=-2$', Rx_3),
('$R_x^{MLSP}, c_1=1.1, c_2=20$', Rx_4)
]
for k in six.moves.range(1, 5):
ax = plt.subplot(2, 2, k)
ax.contour(X, Y, np.real(Rxs[k - 1][1]))
ax.axis('square')
ax.set_title(Rxs[k - 1][0])
ax.set_xlim([-5.0, 5.0])
ax.set_ylim([-5.0, 5.0])
ax.set_xlabel('t')
ax.set_ylabel('s')
plt.show()
xs = [
('$x^{LSP}, c=1.1$', x_1),
('$x^{LSP}, c=3$', x_2),
('$x^{LSCP}, c=1.1, m=2, d=-2$', x_3),
('$x^{MLSP}, c_1=1.1, c_2=20$', x_4)
]
for k in six.moves.range(1, 5):
ax = plt.subplot(2, 2, k)
ax.plot(t_vector, np.real(xs[k - 1][1]))
ax.set_title(xs[k - 1][0])
ax.set_xlim([-10.0, 10.0])
ax.set_ylim([-1.5, 1.5])
ax.set_xlabel('t')
plt.show()
def generate_lsp_processes_figures():
"""Generates Figure 1 and Figure 2 from the paper."""
N = 256
Fs = 7.0
t_vector = np.arange(-(N / 2), N / 2, 1) / Fs
X, Y = np.meshgrid(t_vector, t_vector)
random_seed = 1
x_1, _, Rx_1 = utils.create_lsp_realisation(N, 1.1, Fs, random_seed)
x_2, _, Rx_2 = utils.create_lsp_realisation(N, 3.0, Fs, random_seed)
x_3, _, Rx_3 = utils.create_lscp_realisation(N, 1.1, Fs, 2.0, -2.0,
random_seed)
x_4, _, Rx_4 = utils.create_mlsp_realisation(N, (1.1, 20.0), Fs,
random_seed)
Rxs = [('$R_x^{LSP}, c=1.1$', Rx_1), ('$R_x^{LSP}, c=3$', Rx_2),
('$R_x^{LSCP}, c=1.1, m=2, d=-2$', Rx_3),
('$R_x^{MLSP}, c_1=1.1, c_2=20$', Rx_4)]
for k in six.moves.range(1, 5):
ax = plt.subplot(2, 2, k)
ax.contour(X, Y, np.real(Rxs[k - 1][1]))
ax.axis('square')
ax.set_title(Rxs[k - 1][0])
ax.set_xlim([-5.0, 5.0])
ax.set_ylim([-5.0, 5.0])
ax.set_xlabel('t')
ax.set_ylabel('s')
plt.show()
xs = [('$x^{LSP}, c=1.1$', x_1), ('$x^{LSP}, c=3$', x_2),
('$x^{LSCP}, c=1.1, m=2, d=-2$', x_3),
('$x^{MLSP}, c_1=1.1, c_2=20$', x_4)]
for k in six.moves.range(1, 5):
ax = plt.subplot(2, 2, k)
ax.plot(t_vector, np.real(xs[k - 1][1]))
ax.set_title(xs[k - 1][0])
ax.set_xlim([-10.0, 10.0])
ax.set_ylim([-1.5, 1.5])
ax.set_xlabel('t')
plt.show()
