arg(Xb)で求まります。計測したサンプリングデータについて離散的に上記演算を行う場合は、1周期ではなく、10周期以上フーリエ変換したほうが良いでしょう。信号周期がサンプリング周期の整数倍でないと、切取窓による誤差が出ますから、窓関数を考慮するかフーリエ変換長を長くして、切取誤差の影響を小さくします。信号b(t)の位相が、時々刻々と変化する場合は、フーリエ変換による方法は適していません。このような場合は、ヒルベルト変換フィルタを用いて瞬時位相を求めるほうが良いでしょう。この方法に興味があるのでしたら、補足で要求してください。信号の角周波数は既知なので、これをω0とします。