私は MEG 信号を扱っており、空間漏れの問題に対処するために、時系列の直交化を実行したいと考えていました。理解を深めるために、正弦波で作業を開始しました。Python でコードを書き、直交化するために [spectral python][1] モジュールの関数を使用しました。次のコードを実行して信号をプロットすると、添付の図のような結果が得られます。直交信号も元の信号のように見えるはずです。大幅に変更する必要はありません。信号を直交化するより良い方法があれば教えてください。[ここに画像の説明を入力][2] これが私のコードです:
import matplotlib.pyplot as plt
import spectral
from scipy import signal
from scipy.signal import hilbert
#Computing Sine wave#
t = np.linspace(-0.02, 0.05, 1000)
sig = 325 * np.sin(2*np.pi*50*t)
sig1 = 200 * np.sin(2*np.pi*50*t)
signal = np.zeros((2, 1000))
signal[0] = sig
signal[1] = sig1
#Orthogonalisation
ortho = spectral.orthogonalize(signal)
#Plotting
plt.figure()
plt.plot(signal[0], 'Blue')
plt.plot(signal[1], 'Green')
plt.title('Signal')
plt.figure()
plt.plot(ortho[0], 'Red')
plt.plot(ortho[1], 'Black')
plt.title('Ortho')```
[1]: https://www.spectralpython.net/class_func_ref.html#orthogonalize
[2]: https://i.stack.imgur.com/GKVZ9.png